MATLAB实验练习题(计算机)南邮MATLAB数学实验大作业答案

1、MATLA练习题要求：抄题、写由操作命令、运行结果，并根据要求，贴上运行图。1、求ex 3x2 =0的所有根。(先画图后求解)(要求贴图)>> solve('exp(x)-3*xA2',0) ans =-2*lambertw(-1/6*3A(1/2)-2*lambertw(-1,-1/6*3A(1/2)-2*lambertw(1/6*3A(1/2)2、求下列方程的根。5_1) x 5x 1 =0a=solve('xA5+5*x+1',0);a=vpa(a,6)1.10447+1.05983*i-1.00450+1.06095*i-.199936-1.

2、00450-1.06095*i1.10447-1.05983*i>> fzero('x*sin(x)-1/2', 3) ans =2.9726>> fzero('x*sin(x)-1/2',-3) ans =-2.9726>> fzero('x*sin(x)-1/2',0) ans =-0.74083) sin xcosx x2 =0 所有根>> fzero('sin(x)*cos(x)-xA2',0) ans =>> fzero('sin(x)*cos(x)-x

3、A2',0.6) ans =0.70223、求解下列各题:1)limx _0x -sin x3x>> sym x;>> limit(x-sin(x)/xA3) ans =1/6x(10)2) y =e cosx,求 y>> sym x;>> diff(exp(x)*cos(x),10) ans =(-32)*exp(x)*sin(x)1/223） e ex dx（精确到17位有效数字）> > sym x;>> vpa(int(exp(xA2),x,0,1/2),17)ans =0.54498710418362222

4、4)dxx225 4x2> > sym x;> > int(xA4/(25+xA2),x)ans =125*atan(x/5) - 25*x + *人3/35)求由参数方程Ix ="后F所确定的函数的一阶导数dy与二阶导数 吗I y =arctantdxdx> > sym t;> > x=log(sqrt(1+tA2);y=atan；> > diff(y,t)/diff(x,t)ans =1/t6)设函数y=f(x)由方程xy +3，=3所确定，求y'(x)。> > syms x y;f=x*y+exp(

5、y)-exp(1);> > -diff(f,x)/diff(f,y)ans =-y/(x + exp(y), x7)e- sin 2xdxo> > syms x;> > y=exp(-x)*sin(2*x);> > int(y,0,inf)ans =2/58） 将。1十x在x=0展开（最（Wj次帚为 8）>> syms xf=sqrt(1+x);taylor(f,0,9)ans =-(429*xA8)/32768 + (33*xA7)/2048 - (21*xA6)/1024 + (7*xA5)/256 - (5*xA4)/128 +

6、 xA3/16 -xA2/8 + x/2 + 1isin(3)9) y =e x 求 y (2)> > syms x y;> > y=exp(sin(1/x);> > dy=subs(diff(y,3),x,2)dy =-0.5826 X2 .,、.，I、,.，10)求变上限函数 Ja+tdt对变量x的导致。'x> > syms a t;> > diff(int(sqrt(a+t),t,x,xA2)Warning: Explicit integral could not be found.ans =2*x*(xA2 + a)A

7、(1/2) - (a + xF(1/2)x-3vz-1,、4、求点(1,1,4)到直线L: =-= 的距离-102>> M0=1,1,4;M1=3,0,1;M0M1=M1-M0;v=-1,0,2;d=norm(cross(M0M1,v)/norm(v)d =1.095411x J5、已知f(x) =-=e 2犬，分别在下列条件下画出f(x)的图形：(要求贴图).2 二二(1)仃=1时，N=0,1 1 ,在同一坐标系里作图> > syms x;>> fplot('(1/sqrt(2*pi)*exp(-(x)A2)/2)',-3,3,'r

8、')> > hold on>> fplot('(1/sqrt(2*pi)*exp(-(x-1)A2)/2)',-3,3,'y')> > hold on>> fplot('(1/sqrt(2*pi)*exp(-(x+1)A2)/2)',-3,3,'g')> > hold offD4D35 03D25 02 015D.1门 JD133匡 0(2)卜=0时,仃=1,2,4,在同一坐标系里作图。>> syms x;fplot('(1/sqrt(2*pi

9、)*exp(-(x)A2)/2)',-3,3,T) hold onfplot('(1/(sqrt(2*pi)*2)*exp(-(x)A2)/(2*2A2)',-3,3,y) hold onfplot('(1/(sqrt(2*pi)*4)*exp(-(x)A2)/(2*4A2)',-3,3,'g') hold off(要求贴图)0 _t _200 _u _2x = u 由耐,y = u cosCt), z = t/46、画下列函数的图形：x =u sint(1) «y = u costt z =一 工4>> ezmes

10、h('u*sin(t)','u*cos(t)','t/4',0,20Q2)(2) z =sin(xy) 0 £x £3,0 My 三3>> x=0:0.1:3;y=x;X Y=meshgrid(x,y);Z=sin(X*Y);>> mesh(X,Y,Z)x =sint (3 cosu)(3) «y =cost (3+cosu)I|z =sin u0 _t _2二0_u _2二ezmesh('sin(t)*(3+cos(u)','cos(t)*(3+cos(u)'

11、,'sin(u)',0,2*pi,0,2*pi)x =p+c-Dsfu), y = ci&sfl) (3-hcos(u),: = wi"力口.5-7、已知A =4-3-201-24-3 ,在 MATLAB 命令窗口中建立 A、B矩阵并1对其进行以下操作：(1)计算矩阵A的行列式的值det(A)>> A=4,-2,2;-3,0,5;1,5,3;>> det(A) ans =12 TB, AB ,A B,A2, AT-158(2)分别计算下列各式：2A - B, A* B,A.*>> A=4,-22-3Q5;1,5,3;B=1,

12、3,4;-2Q-3;2,-1,1;>> 2*A-Bans =7-70-40130115>> A*Bans =1210247-14-7-30-8>> A.*Bans =4-6860-152-53>> A*inv(B)ans =-0.0000-0.00002.0000-2.7143-8.0000-8.14292.42863.00002.2857>> inv(A)*Bans =0.48730.41141.00000.3671-0.43040.0000-0.10760.24680.0000>> A*Aans =2424-7319-

13、81336>> Aans =4-31-205253>>8、 在MATLAB中分别利用矩阵的初等变换及函数rank、函数inv求下列矩阵的秩:门(1) A= 3-6 3 2一5 4 0 ,求 rank(A)=?-11 2 4，>> A=1,-6,3,2;3,-5,4,0;-1,-11,2,4;>> rank(A)ans =(2) B =I1 20 10 02 00 2>> B=3,5,0,1;1,2,0,0;1,0,2,0;1,2,0,2>> inv(B)ans =2.0000-1.0000-1.00000-4.00002.

14、50002.0000-0.5000-0.00000.00000.50000-1.00000.50000.50000.50009、在MATLAB中判断下列向量组是否线性相关，并找出向量组«1 =(1132)T,% =(11 13)T,ot3 =(5 289)T,o(4=(1317)T 中的一个最大线性无关组。>> a1=1 1 3 2'a2=-1 1 -1 3'a3=5 -2 8 9'a4=-1 3 1 7'A= a1, a2 ,a3 ,a4 ;R jb=rref(A)a1 =1132 a2 =-11-13 a3 =5-289 a4 =-13

15、171.0000001.09091.7879-0.060601.00000001.00000jb 二1>> A(:,jb)ans =113210、在-11-13MATLAB5-289中判断下列方程组解的情况，若有多个解，写出通解。X x1 f 4x3 -2x4 = 0 (1)- x1 X2 X3 +2人=0| 3x1 , X2 , 7 X3 - 2 X4 = 0x1 -3x2 -12x3 6x4 =0>> A=1,-1,4,2;1,-1,-1,2;3,1,7,-2;1,-3,-12,6;>> rank(A)ans =3>> rref(A)ans

16、=1000010-200100000>> A=1,-1,4,2;1,-1,-1,2;3,1,7,-2;1,-3,-12,6;>> format ratn=4;RA=rank(A)RA =3>> if(RA=n)fprintf('%方程只有零解)elseb=null(A,'r') endb =0201>> syms k X=k*bX =02*k0 k2x1 3x2 x3 =4x1 -2x2 4x3 = -5(2)3x1 8x2 -2x3 =134x1 -x2 9x3 = -6> > A=2 3 1;1 -2 4;

17、3 8 -2;4 -1 9;b=4 -5 13 -6'B=A b;> > n=3;> > RA=rank(A)RA =2>> RB=rank(B)RB =rref(B) ans =102-101-1200000000>> format ratif RA=RB&RA=n % 判断有唯一解X=Abelseif RA=RB&RA<n % 判断有无穷解X=Ab %求特解C=null(A,'r') %求AX=0的基础解系 else X='equition no solve' % 判断无解 end

18、Warning: Rank deficient, rank = 2, tol =8.9702e-015.X =03/2-1/2C =-211-2 1 1 、-._ 111、求矩阵A= 02 0的逆矩阵A及特征值和特征向量。T 1 3A=-2 1 1;0 2 0;-4 1 3;>> a1=inv(A)a1 =1/2-3/21/201/20-21/21>> P,R=eig(A)P =-985/1393-528/2177379/125700379/419-985/1393-2112/2177379/1257-100020002A的三个特征值是：r1=-1 , r2=2, r3

19、=2。三个特征值分别对应的特征向量是P1=1 0 1；p2=1 0 4;p3=1 3 12 22-2 ,12、化方阵A= 25M为对角阵。2 Y 5>> A=2 2 -2;2 5 -4;-2 -4 5;P,D=eig(A)P =- 0.29810.89440.3333- 0.5963-0.44720.6667- 0.74540-0.6667D =1.00000001.000000010.0000>> B=inv(P)*A*PB =1.0000-0.00000.00000.00001.00000.0000-0.0000程序说明：010.0000所求得的特征值矩阵 D即为矩

20、阵A对角化后的角矩阵， D和A相似。22213、求一个正交变换，将二次型= 5x1 +5x2 +3x3 -2x1x2 +6x1x3 -6x2x3 化为标准型。>> A=5 -1 3;-1 5 -3;3 -3 3;>> syms y1 y2 y3 y=y1;y2;y3;P,D=eig(A)881/2158985/1393-780/1351-881/2158985/1393780/1351-780/1351-881/1079>>x=P*y(6A(1/2)*y1)/6 + (2A(1/2)*y2)/2 - (30(1/2)*丫3)/3(2A(1/2)*y2)/2

21、- (6A(1/2)*y1)/6 + (3A(1/2)*y3)/3-(3A(1/2)*y3)/3 - (2人(1/2)*3人(1/2)*丫1)/3>> f=y1 y2 y3*D*y-y1A2/2251799813685248 + 4*y2A2 + 9*y3A27xn 1=(xn -)/214、6位有效设«xn，数列 xn是否收敛？若收敛，其值为多少？精确到X =3数字。f=inline('(x+7/x)/2');>> x0=3;>> fo门=1:20x0=f(x0);fprintf('%g,%gn',i,x0);en

22、d1,2.666672,2.645833,2.645754,2.645755,2.645756,2.645757,2.645758,2.645759,2.6457510,2.6457511,2.6457512,2.6457513,2.6457514,2.6457515,2.6457516,2.6457517,2.6457518,2.6457519,2.6457520,2.64575该数列收敛于三，它的值是11115、设xn =1+下+二+.+ = , xn是否收敛？右收敛， 其值为多少？精确到 17位有效2 P3pnp数字。(注：学号为单号的取p = 7 ,学号为双号的取p = 8 )>

23、>f=inline('1/(xA8)');x0=0;for i=1:20x0=(x0+f(i);fprintf('%g , %.16fn',i,x0);end1 , 1.00000000000000002 , 1.00390625000000003 , 1.00405866579027594,1.00407392457933845 ,1.00407648457933846,1.00407707995351927 , 1.00407725342004488 , 1.00407731302468969 , 1.004077336255262610 , 1.00

24、4077346255262611 , 1.004077350920336512 , 1.004077353246016813 , 1.004077354471911514,1.004077355149515015 , 1.004077355539699316,1.004077355772530017 , 1.004077355915883518 , 1.004077356006628119 , 1.004077356065508520 , 1.0040773561045711>>16、求二重极限limx 1y 0ln(x ey).x2 y2>> clear> &g

25、t; syms x y;> > f=(log(x+exp(y)/sqrt(xA2+yA2);> > 仅=limit(f,'x',1);>>仅丫=1而让(仅'丫',0)fxy =Iog(2)17、已知 ex-xyz=0,求名。 Fx>> clearsyms x y z;>> F=exp(x)-x*y*z;>> Fx= diff(F, 'x')Fx =exp(x) - y*z>> Fz= diff(F, 'z')Fz =-x*y>> G=-

26、Fx/FzG =(exp(x) - y*z)/(x*y)22218、已知函数 f(x, y,z)=x +2y +3z +xy+3x 3y 6z ,求梯度。一：>> clearsyms x y z;> > f=xA2+2*yA2+3*zA2+x*y+3*x-3*y-6*z;> > dxyz=jacobian(f)dxyz =2*x + y + 3, x + 4*y - 3, 6*z - 6二：> > clear> > syms x y z;> > f=xA2+2*yA2+3*zA2+x*y+3*x-3*y-6*z;>

27、> gr=jacobian(f)gr =2*x + y + 3, x + 4*y - 3, 6*z - 6119、计算积分I =口一(2xy)dxdy,其中D由直线y=x与y = x2围成。 d 2>> A=int(int (2-x-y),'y',xA2,x),'x',0,1)/2A =11/120220、计算曲线积分fzyds,其中曲线 C:x = cost,y =sint,z = t tw0,2n。c x yclearsyms x y z tx=cos(t);y=sin(t);z=t;dx=diff(x,t);dy=diff(y,t);dz

28、=diff(z,t);ds=sqrt(dxA2+dyA2+dzA2);f=zA2/(xA2+yA2);I=int(f*ds,t,0,2*pi)I =(8*2人(1*piA3)/321、计算曲面积分f f(x + y + z)dS ,其中 S: z = Ja2 -x2 - y2。S> > clear> > syms x y z a;> > z=sqrt(aA2-xA2-yA2);> > f=x+y+z;>> I=int(int(f,'y',0,sqrt(aA2-xA2),'x',0,a)I=1/2*aA3

29、+1/4*aA3*pi+1/3*aA2*(aA2)A(1/2)+1/3*(-1/2-1/4*pi)*aA32x6.3322、求解二阶微分方程：y -10y 9y -e2x, y(0) =-,y(0) = y o> > clear>> syms x y;>> d_equa='D2y-10*Dy+9*y=exp(2*x)'d_equa =D2y-10*Dy+9*y=exp(2*x)>> Condit= 'y(0)=6/7,Dy(0)=33/7'Condit =y(0)=6/7,Dy(0)=33/7 >> y

30、1=dsolve( d_equa , Condit , 'x')y1 =exp(9*x)/2 - exp(2*x)/7 + exp(x)/2 一二 1 ，一23、求数项级数I =£的和。nj(n 1)> > clear> > syms n;> > f=1/(n*(n+1);> > I=symsum(f,n,1,inf)I =11 一24、将函数f (x)=展开为(x -3)的帚级数。 x> > clear> > syms x;> > f=1/x;>> taylor(f,10

31、,x,3) ans =(x - 3)A2/27 - x/9 - (x - 3)A3/81 + (x - 3)A4/243 - (x - 3)A5/729 + (x - 3)A6/2187 - (x - 3)A7/6561 +(x - 3)人8/19683 - (x - 349/59049 + 2/325、能否找到一个分式线性函数axb ,使它产生的迭代序列收敛到给定的数？用这种办法cx d近似计算2o>> f=inline('(2+xA2)/(2*x)');x1=2;for i=1:20x1=f(x1);fprintf('%g,%gn',i,x1);

32、end;1,1.5 2,1.416673,1.414224,1.414215,1.414216,1.414217,1.414218,1.414219,1.4142110,1.4142111,1.4142112,1.4142113,1.4142114,1.4142115,1.4142116,1.4142117,1.4142118,1.4142119,1.4142120,1.414212x,0 <x < 126、函数f(x)=«2的迭代是否会广生混沌?12(1 -x), x _1 2>> x1=0:0.05:0.5;y1=2*x1;x2=0.5:0.05:1;y2

33、=2*(1-x2);figureplot(x1,y1,x2,y2)gtext('2*x,)gtext('2*(1-x)')27、函数f (x) =ax(1 x) (0EXE1)称为Logistic映射，试从 蜘蛛网”图观察它取初值为x0 =0.5产生的迭代序列的收敛性，将观察记录填人下表，作出图形。若出现循环，请指出它的周期。(要求贴图)0.90.B Q7 HE 0.5 Q4 0.3 Q2 0.1IJ0 Cl1 D.2 D3D.5D.6 Q7 £1B0.91f=inline('3.3*x*(1-x)');x=linspace(1,202,202

34、);y=linspace(1,202,202);x(1)=0.5;y(1)=0;x(2)=x(1);y(2)=x(1);3工值for i=1:100x(1+2*i)=x(2*i);x(2+2*i)=f(x(1+2*i);y(1+2*i)=x(2+2*i);y(2+2*i)=y(1+2*i);endplot(x,y,'r');hold on;syms x y;y=x;ezplot(x,0,1);ezplot(f(x),0,1);axis(0,1,0,3.3/4);T=0.35hold offhold onf=inline('3.5*x*(1-x)');x=lins

35、pace(1,202,202);y=linspace(1,202,202);x(1)=0.5;y(1)=0;x(2)=x(1);y(2)=x(1);fo门=1:100x(1+2*i)=x(2*i);x(2+2*i尸f(x(1+2*i);y(1+2*i)=x(2+2*i);y(2+2*i)=y(1+2*i);endplot(x,y,'r');hold on;syms x y;y=x;4T=0.4ezplot(x,0,1);ezplot(f(x),0,1);axis(0,1,0,3.5/4);hold offhold onf=inline('3.56*x*(1-x)'

36、;);x=linspace(1,202,202);y=linspace(1,202,202);x(1)=0.5;y(1)=0;x(2)=x(1);y(2)=x(1);fo门=1:100x(1+2*i)=x(2*i);x(2+2*i)=f(x(1+2*i);y(1+2*i)=x(2+2*i);y(2+2*i)=y(1+2*i);endplot(x,y,'r');hold on;syms x y;y=x;ezplot(x,0,1);DI (1201 Q5口日 DSezplot(f(x),0,1);axis(0,1,0,3.56/4);hold offhold onf=inline(

37、'3.568*x*(1-x)');x=linspace(1,202,202);y=linspace(1,202,202);x(1)=0.5;y(1)=0;x(2)=x(1);y(2)=x(1);fo门=1:100x(1+2*i)=x(2*i);x(2+2*i)=f(x(1+2*i);y(1+2*i)=x(2+2*i);y(2+2*i)=y(1+2*i);endplot(x,y,'r');hold on;syms x y;y=x;ezplot(x,0,1);ezplot(f(x),0,1);axis(0,1,0,3.568/4);hold onf=inline(&

38、#39;3.6*x*(1-x)');x=linspace(1,202,202);y=linspace(1,202,202);x(1)=0.5;y(1)=0;x(2)=x(1);y(2)=x(1);fo门=1:100x(1+2*i)=x(2*i);x(2+2*i)=f(x(1+2*i);y(1+2*i)=x(2+2*i);y(2+2*i)=y(1+2*i);endplot(x,y,'r');hold on;syms x y;y=x;ezplot(x,0,1);ezplot(f(x),0,1);axis(0,1,0,3.6/4);hold offhold onf=inlin

39、e('3.84*x*(1-x)');x=linspace(1,202,202);y=linspace(1,202,202);x(1)=0.5;y(1)=0;x(2)=x(1);y(2)=x(1);fo门=1:100x(1+2*i)=x(2*i);x(2+2*i)=f(x(1+2*i);y(1+2*i)=x(2+2*i);y(2+2*i)=y(1+2*i);endplot(x,y,'r');hold on;syms x y;y=x;ezplot(x,0,1);ezplot(f(x),0,1);axis(0,1,0,3.84/4); hold off表 Logist

40、ic迭代的收敛性ot3.33.53.563.5683.63.84序列收敛情况不收敛不收不收敛不收敛不收敛不收敛I 敛 I 28、由函数f (x, y) =y 一sgnxj|bx _c|与g(x,y) =a _x构成的二维迭代 Martin迭代。现观 察其当a =4.5 25 55 120, b =2, c = 400时取初值为(0,0)所得到的二维迭代散点图有什 么变化。(要求贴图)function Martin (a,b,c N)f=(x,y)(y-sign(x)*sqrt(abs(a*x-c);g=(x)(a-x);m=0;0;for n=1:Nm(:,n+1)=f(m(1,n),m(2,

41、n),g(m(1,n); endplot(m(1,:),m(2,:),'kx');axis equalMartin(4.52555120,2,-300,500)书上62页(Xi(0),x20)t =(1,2)T ,求出平面映射(x1(n),x2n)T的通项，并画出这些点的散点图。A=4,2;1,3;t=;for i=1:20x=2*rand(2,1)-1;t(length(t)+1,1:2)=x;for j=1:40x=A*x;t(length(t)+1,1:2)=x;endendplot(t(:,1),t(:,2),'*')grid('on')

42、30、对 A=4 1及随机给出的.(0)(0).T，叱(Xi)，x2),观察数列 jn).该数列有极限吗？31、若该地区的天气分为三种状态：晴、阴、雨。对应的转移矩阵为:明天晴2=阴雨今天晴阴雨一3/41/21/41/81/41/21/81/41/4-且p(0) =(0.5,0.25.0.25)T ,试根据这些数据来求出若干天之后的天气状态，并找出其特点(取4位有效数字)。>> A1=3/4,1/2,1/4;1/8,1/4,1/2;1/8,1/4,1/4;p=0.5;0.25;0.25;fo门=1:20p(:,i+1)=A1*p(:,i);endPP =Columns 1 thro

43、ugh 70.50000.56250.59380.60350.60690.60810.60850.25000.25000.22660.22070.21850.21780.21750.25000.18750.17970.17580.17460.17410.1740Columns 8 through 140.60860.60870.60870.60870.60870.60870.60870.21740.21740.21740.21740.21740.21740.21740.17390.17390.17390.17390.17390.17390.1739Columns 15 through 210.60870.60870.60870.60870.60870.60870.60870.21740.21740.21740.21740.21740.21740.21740.17390.17390.17390.17390.17390.17390.173932、对于上例中的 A,求出矩阵 A的特征值与特征向量，并将特征向量与上例中的结论作对比。>> A=3/4 1/2 1/4;1/8 1/4 1/2;1/8 1/4 1/4;>> P,R=eig(A)P =