《【7A文】大学MATLAB数学实验第二版答案》由会员分享，可在线阅读，更多相关《【7A文】大学MATLAB数学实验第二版答案（20页珍藏版）》请在金锄头文库上搜索。
1、【MeiWei81-优质实用版文档】数学实验答案Chapter1Page20,eG1(5)等于eGp(1),eGp(2);eGp(3),eGp(4)(7)3=1G3,8=2G4(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,eG2(1)a,b,c的值尽管都是1，但数据类型分别为数值，字符，逻辑，注意a与c相等
2、，但他们不等于b(2)double(fun)输出的分别是字符a,b,s,(,G,)的ASCII码Page20,eG3r=2;p=0.5;n=12;T=log(r)/n/log(1+0.01Gp)Page20,eG4G=-2:0.05:2;f=G.4-2.G;fmin,min_indeG=min(f)最小值最小值点编址G(min_indeG)ans=0.6500最小值点f1,G1_indeG=min(abs(f)求近似根-绝对值最小的点f1=0.0328G1_indeG=24G(G1_indeG)ans=-0.8500G(G1_indeG)=;f=G.4-2.G;删去绝对值最小的点以求函数绝对值
3、次小的点f2,G2_indeG=min(abs(f)求另一近似根-函数绝对值次小的点f2=0.0630G2_indeG=65G(G2_indeG)ans=1.2500Page20,eG5z=magic(10)z=929918156774515840988071416735557644148188202254566370478587192136062697128869325296168755234172476839042492633652358289914830323966796139597293138457210129496783537444653111810077843643502759su
4、m(z)sum(diag(z)z(:,2)/sqrt(3)z(8,:)=z(8,:)+z(3,:)Chapter2Page45eG1先在编辑器窗口写下列M函数，保存为eg2_1.mfunctionGbar,s=eG2_1(G)n=length(G);Gbar=sum(G)/n;s=sqrt(sum(G.2)-nGGbar2)/(n-1);例如G=81706551766690876177;Gbar,s=eG2_1(G)Page45eG2s=log(1);n=0;whilesek=k+1;F(k)=F(k-1)+F(k-2);G=F(k)/F(k-1);enda,G,k计算至k=21可满足精度Pa
5、ge45eG4clear;tic;s=0;fori=1:1000000s=s+sqrt(3)/2i;ends,toctic;s=0;i=1;whilei1.1)+G.G(G=-1.1)-1.1G(G1);p=p+bGeGp(-y.2-6GG.2).G(G+y-1).G(G+y=1);p=p+aGeGp(-0.75Gy.2-3.75GG.2+1.5GG).G(G+yA=123;456;780;C=2-5-22;-5-24-56;-22-56-16;G=lyap(A,C)G=1.0000-1.0000-0.0000-1.00002.00001.0000-0.00001.00007.0000Chap
6、ter3Page65EG1a=1,2,3;b=2,4,3;a./b,a.b,a/b,abans=0.50000.50001.0000ans=221ans=0.6552一元方程组G2,4,3=1,2,3的近似解ans=0000000.66671.33331.0000矩阵方程1,2,3G11,G12,G13;G21,G22,G23;G31,G32,G33=2,4,3的特解Page65EG2(1)A=41-1;32-6;1-53;b=9;-2;1;rank(A),rank(A,b)A,b为增广矩阵ans=3ans=3可见方程组唯一解G=AbG=2.38301.48942.0213(2)A=4-33;
7、32-6;1-53;b=-1;-2;1;rank(A),rank(A,b)ans=3ans=3可见方程组唯一解G=AbG=-0.4706-0.29410(3)A=41;32;1-5;b=1;1;1;rank(A),rank(A,b)ans=2ans=3可见方程组无解G=AbG=0.3311-0.1219最小二乘近似解(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=1010一个特解Page65EG3a=2,1,-1,1;1,2,1,-1;1,1,2,1;b=
8、1,2,3;G=null(a),G0=abG=-0.62550.6255-0.20850.4170G0=1010通解kG+G0Page65EG4G0=0.20.8;a=0.990.05;0.010.95;G1=aGG,G2=a2GG,G10=a10GGG=G0;fori=1:1000,G=aGG;end,GG=0.83330.1667G0=0.80.2;G=G0;fori=1:1000,G=aGG;end,GG=0.83330.1667v,e=eig(a)v=0.9806-0.70710.19610.7071e=1.0000000.9400v(:,1)./Gans=1.17671.1767成比
9、例，说明G是最大特征值对应的特征向量Page65EG5用到公式(3.11)(3.12)B=6,2,1;2.25,1,0.2;3,0.2,1.8;G=25520;C=B/diag(G)C=0.24000.40000.05000.09000.20XX0.01000.12000.04000.0900A=eye(3,3)-CA=0.7600-0.4000-0.0500-0.09000.8000-0.0100-0.1200-0.04000.9100D=171717;G=ADG=37.569625.786224.7690Page65EG6(1)a=41-1;32-6;1-53;det(a),inv(a),
10、v,d=eig(a)ans=-94ans=0.2553-0.02130.04260.1596-0.1383-0.22340.1809-0.2234-0.0532v=0.0185-0.9009-0.3066-0.7693-0.1240-0.7248-0.6386-0.41580.6170d=-3.05270003.67600008.3766(2)a=11-1;02-1;-120;det(a),inv(a),v,d=eig(a)ans=1ans=2.0000-2.00001.00001.0000-1.00001.00002.0000-3.00002.0000v=-0.57730.5774+0.0000i0.5774-0.0000i-0.57730.57740.5774
