《南邮MATLAB数学实验答案全》由会员分享，可在线阅读，更多相关《南邮MATLAB数学实验答案全（23页珍藏版）》请在金锄头文库上搜索。
1、-第一次练习教学要求：熟练掌握Matlab软件的根本命令和操作，会作二维、三维几何图形，能够用Matlab软件解决微积分、线性代数与解析几何中的计算问题。补充命令vpa(*,n)显示*的n位有效数字，教材102页fplot(f(*),a,b)函数作图命令，画出f(*)在区间a,b上的图形在下面的题目中为你的*的后3位1-9班或4位10班以上1.1 计算与 syms *limit(902*-sin(902*)/*3)ans =366935404/3limit(902*-sin(902*)/*3,inf)ans =01.2，求syms *diff(e*p(*)*cos(902*/1000),2)a
2、ns = (46599*cos(451*)/500)*e*p(*)/250000 - (451*sin(451*)/500)*e*p(*)/2501.3 计算dblquad(*,y) e*p(*.2+y.2),0,1,0,1)ans =2.13941.4 计算syms *int(*4/(9022+4*2)ans =(91733851*atan(*/451)/4 - (203401*)/4 + *3/121.5 syms *diff(e*p(*)*cos(902*),10)ans =1.6 给出在的泰勒展式(最高次幂为4). syms *taylor(sqrt(902/1000+*),5,*)
3、ans = -(9765625*451(1/2)*500(1/2)*4)/82743933602 +(15625*451(1/2)*500(1/2)*3)/91733851 -(125*451(1/2)*500(1/2)*2)/406802 + (451(1/2)*500(1/2)*)/902 +(451(1/2)*500(1/2)/5001.7 Fibonacci数列的定义是用循环语句编程给出该数列的前20项要求将结果用向量的形式给出。*=1,1;for n=3:20 *(n)=*(n-1)+*(n-2);end*=Columns 1 through 10 1 1 2 3 5 8 13 21
4、 34 55Columns 11 through 20 89 144 233 377 610 987 1597 2584 4181 67651.8 对矩阵，求该矩阵的逆矩阵，特征值，特征向量，行列式，计算，并求矩阵是对角矩阵，使得。A=-2,1,1;0,2,0;-4,1,902/1000;inv(A)ans =0.4107 0.0223 -0.4554 0 0.5000 0 1.8215 -0.4554 -0.9107eig(A)ans =-0.5490 + 1.3764i -0.5490 - 1.3764i 2.0000 det(A)ans =4.3920P,D=eig(A)P = %特征向
5、量0.3245 - 0.3078i 0.3245 + 0.3078i 0.2425 0 0 0.9701 0.8944 0.8944 0.0000 D = -0.5490 + 1.3764i 0 0 0 -0.5490 - 1.3764i 0 0 0 2.0000 P*D6*inv(P) %A6的值ans =15.3661 12.1585 + 0.0000i -5.8531 0 64.0000 0 23.4124 -5.8531 + 0.0000i -1.6196 1.9 作出如下函数的图形注：先用M文件定义函数，再用fplot进展函数作图：m文件:function y=fenduan(*)i
6、f *=1/2y=2*else *=1y=2-2*endend执行函数：fplot(fenduan,0,1); grid on title(第1.9题图) 得以下图：1.10 在同一坐标系下作出下面两条空间曲线要求两条曲线用不同的颜色表示12t=-10:0.01:10;*1=cos(t);y1=sin(t);z1=t;plot3(*1,y1,z1);hold on*2=cos(2*t);y2=sin(2*t);z2=t;plot3(*2,y2,z2,m);grid ontitle(第1.10题图)得以下图：1.11 ，在MATLAB命令窗口中建立A、B矩阵并对其进展以下操作：(1) 计算矩阵A
7、的行列式的值(2) 分别计算以下各式：解：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 -7 0 -4 0 7 0 9021 5A*Bans =12 10 12 7 -14 -7 -9013 0 13537A.*Bans = 4 -6 8 6 0 15 2 -4510 3A*inv(B)ans =1.0e+003 * -0.0000 0 0.0020 0.0000 0.0016 0.0001 1.0311 -0.9016 -1.4167inv(A)*Bans = 0.3463 0
8、.5767 0.5383 0.0005 -0.0006 -0.0005 -0.1922 0.3460 0.9230A*Aans = 24 9012 4 -7 22556 9 -13523 13528 22561Aans =4 -3 1 -2 0 4510 2 5 3 1.12 分别在以下条件下画出的图形：(1),分别为(在同一坐标系上作图);(2),分别为(在同一坐标系上作图).(1)*=-5:0.1:5;h=inline(1/sqrt(2*pi)/s*e*p(-(*-mu).2/(2*s2);y1=h(0,902/600,*);y2=h(-1,902/600,*);y3=h(1,902/60
9、0,*);plot(*,y1,b,*,y2,m,*,y3,y)grid ontitle(第1.12题) (2)z1=h(0,1,*);z2=h(0,2,*);z3=h(0,4,*); z4=h(0,902/100,*);plot(*,z1,*,z2,y,*,z3,m,*,z4, g)grid ontitle(第1.12题)z1=h(0,1,*);z2=h(0,2,*);z3=h(0,4,*); z4=h(0,902/100,*);1.13 作出的函数图形。*=-10:0.2:10;y=*;* Y=meshgrid(*,y);Z=902*.2+Y.4;mesh(*,Y,Z);title(第1.1
10、3题)1.14对于方程，先画出左边的函数在适宜的区间上的图形，借助于软件中的方程求根的命令求出所有的实根,找出函数的单调区间,结合高等数学的知识说明函数为什么在这些区间上是单调的,以及该方程确实只有你求出的这些实根。最后写出你做此题的体会。解：作图程序：注：*围的选择是经过试探而得到的*=-1.7:0.02:1.7;y=*.5-902/200*-0.1;plot(*,y);grid on; title(第1.14题)由图形观察，在*=-1.5,*=0,*=1.5附近各有一个实根solve(*5-902/200*-0.1)ans = 1.45733649359338702809415339266
11、24*i + 0.00554243549622792973270286414996580.0055424354962279297327028641499658 - 1.4573364935933870280941533926624*i三个实根的近似值分别为：-1.4517,-0.0222,1.4628由图形可以看出，函数在区间单调上升，在区间单调下降，在区间单调上升。syms *diff(*5-902/200*-0.1,*)结果为5*4-4.51solve(5*4-902/200) ans = -(451(1/4)*500(3/4)/500 (451(1/4)*500(3/4)/500 -(451(1/4)*500(3/4)*i)/500 (451(1/4)*500(3/4)*i)/500vpa(ans)ans =得到两个实根：-0.9745与0.9745可以验证导函数在内为正，函数单调上升导函数在内为负，函数单调下降导函数在内为正，函数单调上升根据函数的单调性，最多有3个实根。1.15 求的所有根。先画图后
