数值分析matlab作业

数值分析matlab作业

ID:1665842

大小:297.00 KB

页数:20页

时间:2017-11-12

数值分析matlab作业_第1页
数值分析matlab作业_第2页
数值分析matlab作业_第3页
数值分析matlab作业_第4页
数值分析matlab作业_第5页
资源描述:

《数值分析matlab作业》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库

1、目录第一章2第1题2第3题2第4题:3第二章6第1题6第2题7第4题8第三章11第1题11第五章13第1题13第七章16第1题:16第八章18第1题1820第一章第1题程序:formatlongi=1;y1=sqrt(5);y2=sqrt(y1);whiley1-y2>10^-8y1=5.^(1./(2.^i));y2=5.^(1./(2.^(i+1)));i=i+1;endy1,y2,i运行结果:y1=1.000000011991247y2=1.000000005995624i=28第3题程序:x=-pi:0.01:pi;y=sin(x);y2=x;y5=x;y10=x;fori=1:2

2、y2=(-1)^i.*(x.^(2*i+1))./(factorial(2*i+1))+y2;endfori=1:5y5=(-1)^i.*(x.^(2*i+1))./(factorial(2*i+1))+y5;endfori=1:10y10=(-1)^i.*(x.^(2*i+1))./(factorial(2*i+1))+y10;end20plot(x,y,'-b',x,y2,'-g',x,y5,':r',x,y10,':k');legend('y','y2','y5','y10')运行结果:第4题:程序:①等分100份x=-10:0.2:10;y=-10:0.2:10;[x,y]=mes

3、hgrid(x,y);z=exp(-abs(x))+cos(x+y)+1./(x.^2+y.^2+1);z=exp(-abs(x))+cos(x+y)+1./(x.^2+y.^2+1);mesh(x,y,z)20运行结果:②等分200份x=-10:0.1:10;y=-10:0.1:10;[x,y]=meshgrid(x,y);z=exp(-abs(x))+cos(x+y)+1./(x.^2+y.^2+1);z=exp(-abs(x))+cos(x+y)+1./(x.^2+y.^2+1);mesh(x,y,z)运行结果:20③等分300份x=-10:0.05:10;y=-10:0.05:10

4、;[x,y]=meshgrid(x,y);z=exp(-abs(x))+cos(x+y)+1./(x.^2+y.^2+1);z=exp(-abs(x))+cos(x+y)+1./(x.^2+y.^2+1);mesh(x,y,z)运行结果:20第二章第1题程序:列主元素求解A=[11111;12345;1361015;14102035;15153570];[L,U,P]=lu(A)L=1.000000001.00001.00000001.00000.50001.0000001.00000.75000.75001.000001.00000.25000.7500-1.00001.0000U=1.

5、00001.00001.00001.00001.00002004.000014.000034.000069.000000-2.0000-8.0000-20.5000000-0.5000-2.37500000-0.2500P=1000000001001000001001000Inv(A)求逆:inv(A)ans=5.0000-10.000010.0000-5.00001.0000-10.000030.0000-35.000019.0000-4.000010.0000-35.000046.0000-27.00006.0000-5.000019.0000-27.000017.0000-4.000

6、01.0000-4.00006.0000-4.00001.0000左除命令AEE=[10000;01000;00100;00010;00001];AEans=5-1010-51-1030-3519-410-3546-276-519-2717-41-46-41第2题程序:20n=8;V=220;R=27;A=zeros(n,n);c=zeros(n,1);c(1)=V/R;A(1,1)=2;fori=1:n-1A(i,i+1)=-2;A(i+1,i)=-2/A(i,i);A(i+1,i+1)=5+2*A(i+1,i);end;fori=2:nc(i)=c(i)-A(i,i-1)*c(i-

7、1);end;x(n)=c(n)/A(n,n);fork=n-1:-1:1x(k)=(c(k)+2*x(k+1))/A(k,k);end;x运行结果:x=8.14784.07372.03651.01750.50730.25060.11940.0477第4题程序:(1)小题A=[10787;7565;86109;75910];det(A),cond(A),eig(A)运行结果:ans=1ans=2.9841e+003ans

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。