东华大学MATLAB数学实验第二版答案.docx
东华大学MATLAB数学实验第二版答案数学实验答案Chapter1Page20,ex1(5)等于exp(1),exp(2);exp(3),exp(4)(7)3=1*3,8=2*4(8)a为各列最小值，b为最小值所在的行号(10)1>=4,false,2>=3,false,3>=2,ture,4>=1,ture(11)答案表明：编址第2元素知足不等式(30>=20)和编址第4元素知足不等式(40>=10)(12)答案表明：编址第2行第1列元素知足不等式(30>=20)和编址第2行第2列元素知足不等式(40>=10)Page20,ex2(1)a,b,c的值尽管都是1，但数据类型分别为数值，字符，逻辑，注意a与c相等，但他们不等于b(2)double(fun)输出的分别是字符a,b,s,(,x,)的ASCII码Page20,ex3>>r=2;p=;n=12;>>T=log(r)/n/log(1+*p)Page20,ex4>>x=-2:2;f=x.4-2.x;>>fmin,min_index=min(f)最小值最小值点编址>>x(min_index)ans=最小值点>>f1,x1_index=min(abs(f)求近似根-绝对值最小的点f1=x1_index=24>>x(x1_index)ans=>>x(x1_index)=;f=x.4-2.x;删去绝对值最小的点以求函数绝对值次小的点>>f2,x2_index=min(abs(f)求另一近似根-函数绝对值次小的点f2=x2_index=65>>x(x2_index)ans=>>z=magic(10)z=929918156774515840988071416735557644148188202254566370478587192136062697128869325296168755234172476839042492633652358289914830323966796139597293138457210129496783537444653111810077843643502759>>sum(z)>>sum(diag(z)>>z(:,2)/sqrt(3)>>z(8,:)=z(8,:)+z(3,:)Chapter2Page45ex1先在编辑器窗口写下列M函数，保存为functionxbar,s=ex2_1(x)n=length(x);xbar=sum(x)/n;s=sqrt(sum(x.2)-n*xbar2)/(n-1);例如>>x=81706551766690876177;>>xbar,s=ex2_1(x)Page45ex2s=log(1);n=0;whilesek=k+1;F(k)=F(k-1)+F(k-2);x=F(k)/F(k-1);enda,x,k计算至k=21可知足精度当前位置：文档视界东华大学MATLAB数学实验第二版答案东华大学MATLAB数学实验第二版答案x=2*sin(fai).*cos(theta);y=2*sin(fai).*sin(theta);z=2*cos(fai);surf(x,y,z)(7)x=linspace(0,pi,100);y1=sin(x);y2=sin(x).*sin(10*x);y3=-sin(x);plot(x,y1,x,y2,x,y3)page45,ex7x=:;y=*(x>+x.*(x=*(x1);p=p+b*exp(-y.2-6*x.2).*(x+y>-1).*(x+y>A=123;456;780;C=2-5-22;-5-24-56;-22-56-16;>>X=lyap(A,C)X=Chapter3Page65Ex1>>a=1,2,3;b=2,4,3;a./b,a.b,a/b,abans=ans=221ans=一元方程组x2,4,3=1,2,3的近似解ans=000000矩阵方程1,2,3x11,x12,x13;x21,x22,x23;x31,x32,x33=2,4,3的特解Page65Ex2(1)>>A=41-1;32-6;1-53;b=9;-2;1;>>rank(A),rank(A,b)A,b为增广矩阵ans=3ans=3可见方程组唯一解>>x=Abx=(2)>>A=4-33;32-6;1-53;b=-1;-2;1;>>rank(A),rank(A,b)ans=3ans=3可见方程组唯一解>>x=Abx= (3)>>A=41;32;1-5;b=1;1;1;>>rank(A),rank(A,b)ans=2ans=3可见方程组无解>>x=Abx=最小二乘近似解(4)>>a=2,1,-1,1;1,2,1,-1;1,1,2,1;b=123'%注意b的写法>>rank(a),rank(a,b)ans=3ans=3rank(a)=rank(a,b)>abans=110一个特解Page65Ex3>>a=2,1,-1,1;1,2,1,-1;1,1,2,1;b=1,2,3'>>x=null(a),x0=abx=x0=11通解kx+x0Page65Ex4>>x0='a=;>>x1=a*x,x2=a2*x,x10=a10*x>>x=x0;fori=1:1000,x=a*x;end,xx=>>x0='>>x=x0;fori=1:1000,x=a*x;end,xx=>>v,e=eig(a)v=e=>>v(:,1)./xans=成比例，讲明x是最大特征值对应的特征向量Page65Ex5用到公式>>B=6,2,1;,1,;3,;x=25520'>>C=B/diag(x)C=>>A=eye(3,3)-CA=>>D=171717'x=ADx=Page65Ex6(1)>>a=41-1;32-6;1-53;det(a),inv(a),v,d=eig(a)ans=-94ans=v=d=000000(2)>>a=11-1;02-1;-120;det(a),inv(a),v,d=eig(a)ans=1ans=v=+-+d=000+000-(3)>>A=5765;71087;68109;57910A=5765710876810957910>>det(A),inv(A),v,d=eig(A)ans=1ans=v=d=000000000000(4)(以n=5为例)方法一三个forn=5;fori=1:n,a(i,i)=5;endfori=1:(n-1),a(i,i+1)=6;endfori=1:(n-1),a(i+1,i)=1;enda方法二一个forn=5;a=zeros(n,n);a(1,1:2)=56;fori=2:(n-1),a(i,i-1,i,i+1)=156;enda(n,n-1n)=15;a方法三不用for
