南邮MATLAB数学实验答案(全).docx

《南邮MATLAB数学实验答案(全).docx》由会员分享，可在线阅读，更多相关《南邮MATLAB数学实验答案(全).docx（12页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、南邮MATLAB数学实验答案(全)第一次练习教学要求：熟练把握Matlab软件的基本命令和操作，会作二维、三维几何图形，能够用Matlab软件解决微积分、线性代数与解析几何中的计算问题。补充命令vpa(x,n)显示x的n位有效数字，教材102页fplot(f(x),a,b)函数作图命令，画出f(x)在区间a,b上的图形在下面的题目中m为你的学号的后3位1-9班或4位10班以上1.1计算30sinlimxmxmxx-与3sinlimxmxmxx-symsxlimit(902*x-sin(902*x)/x3)ans=366935404/3limit(902*x-sin(902*x)/x3,inf)

2、ans=01.2cos1000xmxye=，求ysymsxdiff(exp(x)*cos(902*x/1000),2)ans=(46599*cos(451*x)/500)*exp(x)/250000-(451*sin(451*x)/500)*exp(x)/2501.3计算221100xyedxdy+?dblquad(x,y)exp(x.2+y.2),0,1,0,1)ans=2.13941.4计算4224xdxmx+?symsxint(x4/(9022+4*x2)ans=(91733851*atan(x/451)/4-(203401*x)/4+x3/121.5(10)cos,xyemxy=求sy

3、msxdiff(exp(x)*cos(902*x),10)ans=-356485076957717053044344387763*cos(902*x)*exp(x)-3952323024277642494822005884*sin(902*x)*exp(x)1.60x=的泰勒展式(最高次幂为4).symsxtaylor(sqrt(902/1000+x),5,x)ans=-(9765625*451(1/2)*500(1/2)*x4)/82743933602+(15625*451(1/2)*500(1/2)*x3)/91733851-(125*451(1/2)*500(1/2)*x2)/40680

4、2+(451(1/2)*500(1/2)*x)/902+(451(1/2)*500(1/2)/5001.7Fibonacci数列nx的定义是121,1xx=12,(3,4,)nnnxxxn-=+=用循环语句编程给出该数列的前20项要求将结果用向量的形式给出。x=1,1;forn=3:20x(n)=x(n-1)+x(n-2);endxx=Columns1through1011235813213455Columns11through208914423337761098715972584418167651.8对矩阵211020411000Am?-?=?-?，求该矩阵的逆矩阵，特征值，特征向量，行列式

5、，计算6A，并求矩阵,PDD是对角矩阵，使得1APDP-=。A=-2,1,1;0,2,0;-4,1,902/1000;inv(A)ans=0.41070.0223-0.455400.500001.8215-0.4554-0.9107eig(A)ans=-0.5490+1.3764i-0.5490-1.3764i2.0000det(A)ans=4.3920 P,D=eig(A)P=%特征向量0.3245-0.3078i0.3245+0.3078i0.2425000.97010.89440.89440.0000D=-0.5490+1.3764i000-0.5490-1.3764i0002.0000

6、P*D6*inv(P)%A6的值ans=15.366112.1585+0.0000i-5.8531064.0000023.4124-5.8531+0.0000i-1.61961.9作出如下函数的图形注：先用M文件定义函数，再用fplot进行函数作图：1202()12(1)12xxfxxx?=?-1cossinxtytzt=?=?=?22cos2sinxtytzt=?=?=?t=-10:0.01:10;x1=cos(t);y1=sin(t);z1=t;plot3(x1,y1,z1);holdonx2=cos(2*t);y2=sin(2*t);z2=t;plot3(x2,y2,z2,m);grid

7、ontitle(第1.10题图)得下列图：第1.10题图1.11已知422134305,203153211ABm-?=-=-?-?，在MATLAB命令窗口中建立A、B矩阵并对其进行下面操作：(1)计算矩阵A的行列式的值det()A(2)分别计算下列各式：1122,*,.*,TABABABABABAA-解：A=4,-2,2;-3,0,5;1,5*902,3;B=1,3,4;-2,0,3;2,-1,1;det(A)ans=-1172882*A-Bans=7-70-407090215A*Bans=1210127-14-7-9013013537A.*Bans=4-6860152-45103A*inv(

8、B)ans=1.0e+003*-0.000000.00200.00000.00160.00011.0311-0.9016-1.4167inv(A)*Bans=0.34630.57670.53830.0005-0.0006-0.0005-0.19220.34600.9230A*Aans=2490124-7225569-135231352822561Aans=4-31-2045102531.12已知22()2()xfx-=分别在下列条件下画出)(xf的图形： (1)/600m=,分别为0,1,1-(在同一坐标系上作图);(2)0=,分别为1,2,4,/100m(在同一坐标系上作图).(1)x=-5:0.1:5;h=inline(1/sqrt(2*pi)/s*exp(-(x-mu).2/(2*s2);y1=h(0,902/600,x);y2=h(-1,902/600,x);y3=h(1,902/600,x);plot(x,y1,b,x,y2,m,x,y3,y)