资源描述：
《南邮 数学实验参考答案(选题版)》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、1.3、symsxy;>>a=int(int(exp(x^2+y^2),x,0,1),y,0,1)a=(pi*erfi(1)^2)/41.7、n=20;fori=1:(n-2)a(1)=1;a(2)=1;a(i+2)=a(i+1)+a(i);enda'ans=112358132134558914423337761098715972584418167651.8、>>A=[-2,1,1;0,2,0;-4,1,303/1000];>>inv(A)ans=0.08930.1027-0.294600.500001.1786-0.29
2、46-0.5893>>eig(A)ans=-0.8485+1.6353i-0.8485-1.6353i2.0000>>[p,D]=eig(A)p=0.2575-0.3657i0.2575+0.3657i0.2425000.97010.89440.89440.0000D=-0.8485+1.6353i000-0.8485-1.6353i0002.0000>>det(A)ans=6.7880>>A^6ans=45.01944.7452-6.3718064.0000025.4870-6.371830.3452>>A.^6ans=
3、1.0e+003*0.06400.00100.001000.064004.09600.00100.00001.9、M文件定义如下：functiony=f(x)ifx>=0&&x<=1/2y=2*x;elseifx>1/2&&x<=1y=2-2*x;endend命令窗口执行：fplot(@f,[0,1])1.10、t=-8:0.1:8;x=cos(t);y=sin(t);z=t;plot3(x,y,z,'r');holdonx1=2*cos(t);y1=2*sin(t);z1=t;plot3(x1,y1,z1)gridon1
4、.11、>>A=[4-22;-305;15*3033];>>B=[134;-20-3;2-1-1];>>det(A)ans=-39418>>2*A-Bans=7-70-4013030317>>A*Bans=1210207-14-17-30230-4544>>A.*Bans=4-6860-152-1515-3>>A*B^-1ans=-0.4211-1.47370.7368-1.0000-2.0000-3.0000637.7368716.5789398.2105>>A^-1*Bans=0.34670.57630.99950.0
5、015-0.0017-0.0013-0.19200.3458-0.0003>>A^2ans=2430224-775819-453845437586>>A'ans=4-31-2015152531.12、symsx;fplot('(1/(sqrt(2*pi)*514/600))*exp(-((x)^2)/2)',[-3,3],'r')holdonfplot('(1/(sqrt(2*pi)*514/600))*exp(-((x-1)^2)/2)',[-3,3],'b')holdonfplot('(1/(sqrt(2*pi)*51
6、4/600))*exp(-((x+1)^2)/2)',[-3,3],'g')holdofflegend('u为0','u为-1','u为1')symsx;fplot('(1/(sqrt(2*pi)*1))*exp(-((x)^2)/2)',[-3,3],'r')holdonfplot('(1/(sqrt(2*pi)*2))*exp(-((x)^2)/2)',[-3,3],'b')holdonfplot('(1/(sqrt(2*pi)*4))*exp(-((x)^2)/2)',[-3,3],'--')holdonfplot(
7、'(1/(sqrt(2*pi)*5.14))*exp(-((x)^2)/2)',[-3,3],'g')holdoff1.15、ezplot('exp(x)-3*303*x.^2',[-10,10]);gridonfsolve('exp(x)-3*303*x.^2',0)ans=-0.0326第二次练习：2.1、f=inline('(x+7/x)/2');symsx;x0=3;fori=1:1:15x0=f(x0);fprintf('%g,%g',i,x0);end结果如下：1,2.666672,2.645833,2.6
8、45754,2.645755,2.645756,2.645757,2.645758,2.645759,2.6457510,2.6457511,2.6457512,2.6457513,2.6457514,2.6457515,2.645752.2、同2.1的方法，把f=inline('(x+7/x)/2')
