数学实验第二章答案[篇]

1、数学实验第二章答案3篇 以下是网友分享的关于数学实验第二章答案的资料3篇，希望对您有所帮助，就爱阅读感谢您的支持。数学实验第二章答案（一）Matlab数学实验（胡良剑）1.x=81,70,65,51,76,66,90,87,61,77;a=mean(x)b=std(x)2.s=0;forn=0:100s=s+log(1+n);ifs100break;endendn3.F=1,1;forn=3:100F(n)=F(n-1)+F(n-2);endFform=2:100a=F(m)/F(m-1);ifabs(a-(1+sqrt(5)/2)break;endendm4.clear;K=0;fori=1

2、:1000000K=K+sqrt(3)/(2);endKK=0;i=1;whileiK=K+sqrt(3)/(2);i=i+1;endKi=1:1000000;K=sum(sqrt(3)./(2.);K5.(1)t=*2324;T=*1*716;plot(t,T);(2)dlmwrite(Tem.txt,t;T,delimiter,b);(3)6.(1)subplot(1,2,1);fplot(x)x*sin(x-x-2),-2,2,r);x=-2:0.1:2;y=x.*sin(x.-x-2);subplot(1,2,2);plot(x,y,b:);(2)ezplot(x)/4+(y)/9-1

3、,-2,2,-3,3);(3)x,y=meshgrid(-3:0.1:3);z=x.+y.;mesh(x,y,z);(4)xa=-3:0.1:3;ya=-3:0.2:13;x,y=meshgrid(xa,ya);z=x.+3*x.+y.-2*x-2*y-2*x.*y+6;mesh(x,y,z);(5)t=0:pi/360:2*pi;x=sin(t);y=cos(t);z=cos(2*t);plot3(x,y,z);gridon;(6)aa=0:pi/180:2*pi;bb=0:pi/360:pi/2;a,b=meshgrid(aa,bb);x=2*sin(b).*cos(a);y=2*sin(

4、b).*sin(a);z=2*cos(b);mesh(x,y,z);(7)x=0:pi/180:pi;y1=sin(x);y2=sin(x).*sin(10*x);y3=-sin(x);plot(x,y1,b-,x,y2,r+-,x,y3,gx-);xlabel(x);ylabel(y);axis(0pi+1-1.51.5);text(pi/2,1.1,yl=sinx);text(pi/(3.5),0.5,y2=sinxsin10x);text(pi/2,-1.1,y3=-sinx);legend(y1,y2,y3);7.x=-3.3:0.1:3.3;y=-1.1*(x=-1.1&x1.1);

5、plot(x,y);title(y函数);xlabel(X);ylabel(Y);8.9.函数文件P.mfunctionPxy=P(x,y)Pxy=0.5457*exp(-0.75*y.-3.75*x.-1.5*x).*(x+y1).+0.7575*exp(-y.-6*x.).*(x+y-1).+0.5457*exp(-0.75*y.-3.75*x.+1.5*x).*(x+y命令文件x,y=meshgrid(-3:0.1:3);s=size(x);Pxy=zeros(s(1),s(2);fori=1:s(1)forj=1:s(2)Pxy(i,j)=P(x(i,j),y(i,j);endendm

6、esh(x,y,Pxy);10.A=123;456;780;C=2-5-22;-5-24-56;-22-56-16;X=lyap(A,C);X数学实验第二章答案（二）实验二条件控制与循环编程作业要求：1. 将作业题目的所有程序编辑在一个m 文件中，文件名格式：sy2zq20114321(第一部分sy2是实验课题代号, 第二部分zq 是名字张强拼音缩写，第三部分20114321是学号)2. 作业第一行注释标明m 文件名，例:%sy2zq20114321，每题用注释标明题号。3. 交作业请发到老师指定邮箱：周二班(1,2号)(3,4号)(5,6号)(7,8号)(9,0号)周三班(1,2号)(3,4

7、号)(5,6号)(7,8号)(9,0号)4.5.6.7. 请写电邮主题为：数学实验作业二姓名学号 作业请在当天完成发送。 作业发送尽量不用附件，可将m 文件内容直接贴在信中，便于老师审阅。作业第一行注明m 文件名：%sy2zq20114321第一大题：建立数组，构造函数求对应的函数值.1.1 x 1=(3, 5, -1, 2, 8, 12), y 1=3x 2+e x -x +2x 2=(-2, -1, 0, 1, 2, 3, 4), y 2=sin(x +2) -1， x 3+61.2第二大题：运用条件控制语句编写程序, 实现根据输入的数值特性, 按分段函数计算出相应的函数值；2.1 构造由

8、键盘输入x 的值，由分段函数y1确定函数值2x -sin 4x x 0 y 1=x 3x 运行时键盘分别输入值x = 2 x = -12由程序得出相应y 1的函数值5x +3x 0x 4x =0键盘分别输入值x = 45,x = -32时y 2的值.第三大题：构造函数式M 文件f3.m, 函数表达式为ln(2x ) x 0，并求f 3(4), f 3(0), f (-2) y 3=32x -x +5x 0第四大题：运用条件控制语句、循环语句等基本编程语句编写程序, 实现有规律元素的矩阵生成; 实现各类特殊函数及数列的生成;4.1 由for 循环构造矩阵A45 21A 4= 1 1 121111

9、521112521112521 11252111254.2 用for 循环由数组t=3,2,5,4,-2, -3生成的范德蒙矩阵。（范德蒙矩阵的第i 行是数组的i-1次方）4.3 用while 循环构造求调和级数时的和S15与S20。 1前n 项和，项数由键盘输入。并求出n=15,n=20k k =1n4.4 用循环求费波那契数列的前40个数，以四个数为一行排成104的数阵F（费波那契数列的第一项、第二项是1，从第三项起各项是其前两项的和）4.5 键盘输入a 的值与项数n , 编程实现如下n 项和Sn ：Sn =a +aa + +aa an 个a当a = 2，n = 6时求S n（即：2+22

10、+222+2222+22222+222222）第五大题:运用多分枝控制语句，实现自动转换成绩制式功能. 用switch 语句编程实现输入百分制的成绩，输出90100为优秀，7089为良好，6069为及格，60以下为不及格的等级制成绩。第六大题：求200以内的素数.第七大题:用switch 语句编程实现运输费用的计算：运输公司对用户计算运费，距离s 越远，每公里运费越低。折扣标准如下：s250s500s1000s2000s3000s 15%折扣设每公里每吨货物基本运费为p, 货物重为w ，距离为s, 折扣为d, 则编程计算总费用f . 其f 的计算模型为：%sy2hwj20102739%1.1x

11、1=3 5 -1 2 8 12,y1=3*x1.+exp(x1)-x1+2%1.2x2=-2:4 y2=(sin(x2+2)-1)./(x2.+6) f =p w s (1-d )x=input(x=)if x=0y1=2*x-sin(4*x)elseif xy1=exp(x)+xendy1x=2*pix =6.2832 y1 =12.5664%2.2x=input(x=)if x0y2=5*x+3;elseif xy2=3*(x);else x=0y2=4;endy2%3.1function y3=f3(x)y3=log(2*x)else y3=2*(x)-x+5 end%4.1for i=

12、1:6for j=1:6if i=ja4(i,j)=5;elseif abs(i-j)=1a4(i,j)=2;else a4(i,j)=1;endendenda4 %4.2for i=1:6for j=1:6switch j1/x;x=x+1;endsum %4.4a44(1)=0;a44(2)=1;for k=3:1:40a44(k)=a44(k-1)+a44(k-2); k=k+1;endk=1;for i=1:4for j=1:10b44(i,j)=a44(k);k=k+1;endendb44 %4.5a=input(a=);n=input(n=);sum=0;for i=1:nj=0;

13、for x=0:i-1;j=j*10+a;endsum=sum+j;endsum %5x=input(x=);switch fix(x/10)case9,10f=优秀case 8,7f=良好case 6f=及格otherwisef=不及格end %6N=200;x=2:N;for u=2:sqrt(N)n=find(rem(x,u)=0 & x=u); x(n)=;end %7p=input(基本运费=);w=input(货物重量=); s=input(距离=);a=fix(s/250);switch a case 0f=p*w*s case 1f=p*w*s*0.98 case 2,3 f=

14、p*w*s*0.95 case4,5,6,7 f=p*w*s*0.92 case8,9,10,11 f=p*w*s*0.9 otherwisef=p*w*s*0.85 end 数学实验第二章答案（三）matlab数学实验 第一至第四章答案(胡良剑版)第一章%Page20,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

15、=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=0.5;n=12; T=log(r)/n/log(1+0.01*p) T =11.5813%Page20,ex4 x=-2:0.05:2;f=x.-2.; fmin,min_index=min(f) fmin =-1.3907 %

16、最小值 min_index =54 %最小值点编址 x(min_index)ans =0.6500 %最小值点 f1,x1_index=min(abs(f) %求近似根-绝对值最小的点 f1 =0.0328 x1_index = 24 x(x1_index) ans =-0.8500 x(x1_index)=;f=x.-2.; %删去绝对值最小的点以求函数绝对值次小的点 f2,x2_index=min(abs(f) %求另一近似根-函数绝对值次小的点 f2 =0.0630 x2_index = 65 x(x2_index) ans =1.2500%Page20,ex5 z=magic(10)z

17、 =92 99 1 98 80 7 4 81 88 85 87 19 86 93 25 17 24 76 23 5 82 79 6 13 10 12 94 11 18 100 sum(z) ans =505 505 505 sum(diag(z) ans = 505 z(:,2)/sqrt(3) ans =57.1577 46.1880 46.7654 50.2295 53.6936 13.8564 2.8868 3.4641 6.9282 10.3923 z(8,:)=z(8,:)+z(3,:) z =92 99 1 8 15 14 16 20 22 21 3 2 9 83 90 89 91

18、 95 97 96 78 77 84 505 505 8 15 67 74 73 55 54 56 60 62 61 68 42 49 48 30 29 31 35 37 36 43 505 505 67 74 51 58 57 64 63 70 69 71 75 52 26 33 32 39 38 45 44 46 50 27 505 505 51 58 40 41 47 28 34 65 66 72 53 59 505 40 98 80 7 14 16 73 55 57 64 41 4 85 86 1781 87 93 2488 20 22 54 56 63 70 47 19 21 3 6

19、0 62 69 71 28 25 2 9 61 68 75 52 34 76 83 90 42 49 26 33 6523 5 82 89 91 48 30 32 39 66 83 87 101 115 119 83 87 101 115 119 10 12 94 96 78 35 37 44 46 5311 18 100 77 84 36 43 50 27 59第二章%Page 40 ex1先在编辑器窗口写下列M函数，保存为eg2_1.m function xbar,s=ex2_1(x) n=length(x); xbar=sum(x)/n;s=sqrt(sum(x.)-n*xbar)/(n

20、-1); 例如x=81 70 65 51 76 66 90 87 61 77; xbar,s=ex2_1(x) xbar = 72.4000 s =12.1124%Page 40 ex2 s=log(1);n=0; while ss=s+log(1+n); end m=n计算结果m=37%Page 40 ex3 clear;F(1)=1;F(2)=1;k=2;x=0; e=1e-8; a=(1+sqrt(5)/2; while abs(x-a)ek=k+1;F(k)=F(k-1)+F(k-2); x=F(k)/F(k-1); enda,x,k计算至k=21可满足精度%Page 40 ex4 c

21、lear;tic;s=0; for i=1:1000000 s=s+sqrt(3)/2; end s,toc tic;s=0;i=1;while is=s+sqrt(3)/2;i=i+1; end s,toc tic;s=0;i=1:1000000;s=sqrt(3)*sum(1./2.); s,toc%Page 40 ex5 t=0:24;c=15 14 14 14 14 15 16 18 20 22 23 25 28 . 31 32 31 29 27 25 24 22 20 18 17 16; plot(t,c) %Page 40 ex6 %(1)x=-2:0.1:2;y=x.*sin(x

22、.-x-2);plot(x,y) y=inline(x*sin(x-x-2);fplot(y,-2 2)%(2)参数方法t=linspace(0,2*pi,100);x=2*cos(t);y=3*sin(t); plot(x,y) %(3)x=-3:0.1:3;y=x;x,y=meshgrid(x,y); z=x.+y.; surf(x,y,z)%(4)x=-3:0.1:3;y=-3:0.1:13; x,y=meshgrid(x,y);z=x.+3*x.+y.-2*x-2*y-2*x.*y+6; surf(x,y,z) %(5)t=0:0.01:2*pi;x=sin(t);y=cos(t);z

23、=cos(2*t); plot3(x,y,z) %(6)theta=linspace(0,2*pi,50);fai=linspace(0,pi/2,20); theta,fai=meshgrid(theta,fai); 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)%page41, ex7 x=-1.5:0.05:1

24、.5;y=1.1*(x1.1)+x.*(x=-1.1)-1.1*(x%page41,ex8分别使用which trapz, type trapz, dir C:MATLAB7toolboxmatlabdatafun%page41,ex9 clear;close; x=-2:0.1:2;y=x;x,y=meshgrid(x,y); a=0.5457;b=0.7575;p=a*exp(-0.75*y.-3.75*x.-1.5*x).*(x+y1); p=p+b*exp(-y.-6*x.).*(x+y-1).*(x+yp=p+a*exp(-0.75*y.-3.75*x.+1.5*x).*(x+y%p

25、age41, ex10 lookfor lyapunov help lyap A=1 2 3;4 5 6;7 8 0;C=2 -5 -22;-5 -24 -56;-22 -56 -16; X=lyap(A,C) X =1.0000 -1.0000 -0.0000 -1.0000 2.0000 1.0000 -0.0000 1.0000 7.0000第三章%Chapter 3 %Exercise 1 a=1,2,3;b=2,4,3;a./b,a.b,a/b,ab ans =0.5000 0.5000 1.0000 ans =2 2 1 ans =0.6552 %一元方程组x2,4,3=1,2,3

26、的近似解 ans =0 0 00 0 0 0.6667 1.3333 1.0000%矩阵方程1,2,3x11,x12,x13;x21,x22,x23;x31,x32,x33=2,4,3的特解%Exercise 2(1) A=4 1 -1;3 2 -6;1 -5 3;b=9;-2;1; rank(A), rank(A,b) %A,b为增广矩阵 ans = 3ans =3 %可见方程组唯一解 x=Ab x = 2.3830 1.4894 2.0213%Exercise 2(2) A=4 -3 3;3 2 -6;1 -5 3;b=-1;-2;1; rank(A), rank(A,b) ans = 3

27、ans =3 %可见方程组唯一解 x=Ab x =-0.4706 -0.2941 0 %Exercise 2(3) A=4 1;3 2;1 -5;b=1;1;1; rank(A), rank(A,b) ans = 2 ans =3 %可见方程组无解 x=Ab x =0.3311-0.1219 %最小二乘近似解 %Exercise 2(4) a=2,1,-1,1;1,2,1,-1;1,1,2,1;b=1 2 3;%注意b的写法 rank(a),rank(a,b) ans = 3ans =3 %rank(a)=rank(a,b) ab ans = 10 10 %一个特解%Exercise 3 a=

28、2,1,-1,1;1,2,1,-1;1,1,2,1;b=1,2,3; x=null(a),x0=ab x =-0.6255 0.6255 -0.2085 0.4170 x0 = 1 0 1 0%通解kx+x0%Exercise 4 x0=0.2 0.8;a=0.99 0.05;0.01 0.95; x1=a*x, x2=a*x, x10=a*x x=x0;for i=1:1000,x=a*x;end,xx =0.8333 0.1667 x0=0.8 0.2; x=x0;for i=1:1000,x=a*x;end,x x = 0.8333 0.1667 v,e=eig(a) v =0.9806

29、 -0.7071 0.1961 0.7071 e =1.0000 0 0 0.9400 v(:,1)./xans =1.17671.1767 %成比例，说明x是最大特征值对应的特征向量%Exercise 5%用到公式(3.11)(3.12) B=6,2,1;2.25,1,0.2;3,0.2,1.8;x=25 5 20; C=B/diag(x) C =0.2400 0.4000 0.0500 0.0900 0.2000 0.0100 0.1200 0.0400 0.0900 A=eye(3,3)-C A =0.7600 -0.4000 -0.0500 -0.0900 0.8000 -0.0100

30、 -0.1200 -0.0400 0.9100 D=17 17 17;x=AD x =37.5696 25.7862 24.7690%Exercise 6(1) a=4 1 -1;3 2 -6;1 -5 3;det(a),inv(a),v,d=eig(a) ans = -94 ans =0.2553 -0.0213 0.0426 0.1596 -0.1383 -0.2234 0.1809 -0.2234 -0.0532 v =0.0185 -0.9009 -0.3066 -0.7693 -0.1240 -0.7248 -0.6386 -0.4158 0.6170 d =-3.0527 0 0

31、0 3.6760 0 0 0 8.3766%Exercise 6(2) a=1 1 -1;0 2 -1;-1 2 0;det(a),inv(a),v,d=eig(a) ans = 1ans =2.0000 -2.0000 1.0000 1.0000 -1.0000 1.0000 2.0000 -3.0000 2.0000 v =-0.5773 0.5774 + 0.0000i 0.5774 - 0.0000i-0.5773 0.5774 0.5774 -0.5774 0.5773 - 0.0000i 0.5773 + 0.0000i d =1.0000 0 0 0 1.0000 + 0.000

32、0i 0 0 0 1.0000 - 0.0000i%Exercise 6(3) A=5 7 6 5;7 10 8 7;6 8 10 9;5 7 9 10 A =5 7 6 5 7 10 8 7 6 8 10 9 5 7 9 10 det(A),inv(A), v,d=eig(A) ans = 1 ans =68.0000 -41.0000 -17.0000 10.0000 -41.0000 25.0000 10.0000 -6.0000 -17.0000 10.0000 5.0000 -3.0000 10.0000 -6.0000 -3.0000 2.0000v =0.8304 0.0933

33、0.3963 0.3803 -0.5016 -0.3017 0.6149 0.5286 -0.2086 0.7603 -0.2716 0.5520 0.1237 -0.5676 -0.6254 0.5209 d =0.0102 0 0 0 0 0.8431 0 0 0 0 3.8581 0 0 0 0 30.2887%Exercise 6(4)(以n=5为例) %关键是矩阵的定义 %方法一（三个for） n=5;for i=1:n, a(i,i)=5;end for i=1:(n-1),a(i,i+1)=6;end for i=1:(n-1),a(i+1,i)=1;end a%方法二（一个fo

34、r） n=5;a=zeros(n,n);a(1,1:2)=5 6;for i=2:(n-1),a(i,i-1,i,i+1)=1 5 6;end a(n,n-1 n)=1 5; a%方法三（不用for） n=5;a=diag(5*ones(n,1);b=diag(6*ones(n-1,1); c=diag(ones(n-1,1);a=a+zeros(n-1,1),b;zeros(1,n)+zeros(1,n);c,zeros(n-1,1) %下列计算 det(a) ans = 665 inv(a) ans =0.3173 -0.5865 1.0286 -1.6241 1.9489 -0.0977

35、 0.4887 -0.8571 1.3534 -1.6241 0.0286 -0.1429 0.5429 -0.8571 1.0286 -0.0075 0.0376 -0.1429 0.4887 -0.5865 0.0015 -0.0075 0.0286 -0.0977 0.3173 v,d=eig(a) v =-0.7843 -0.7843 -0.9237 0.9860 -0.92370.5546 -0.5546 -0.3771 -0.0000 0.3771-0.2614 -0.2614 0.0000 -0.1643 0.00000.0924 -0.0924 0.0628 -0.0000 -

36、0.0628-0.0218 -0.0218 0.0257 0.0274 0.0257d =0.7574 0 0 0 00 9.2426 0 0 00 0 7.4495 0 00 0 0 5.0000 00 0 0 0 2.5505 %Exercise 7(1) a=4 1 -1;3 2 -6;1 -5 3;v,d=eig(a)v =0.0185 -0.9009 -0.3066-0.7693 -0.1240 -0.7248-0.6386 -0.4158 0.6170d =-3.0527 0 00 3.6760 00 0 8.3766 det(v)ans =-0.9255 %v行列式正常, 特征向

37、量线性相关，可对角化 inv(v)*a*v %验算ans =-3.0527 0.0000 -0.00000.0000 3.6760 -0.0000-0.0000 -0.0000 8.3766 v2,d2=jordan(a) %也可用jordanv2 =0.0798 0.0076 0.91270.1886 -0.3141 0.1256-0.1605 -0.2607 0.4213 %特征向量不同d2 =8.3766 0 00 -3.0527 - 0.0000i 00 0 3.6760 + 0.0000i v2a*v2ans =8.3766 0 0.00000.0000 -3.0527 0.0000

38、0.0000 0.0000 3.6760 v(:,1)./v2(:,2) %对应相同特征值的特征向量成比例ans =2.44912.44912.4491 %Exercise 7(2) a=1 1 -1;0 2 -1;-1 2 0;v,d=eig(a)v =-0.5773 0.5774 + 0.0000i 0.5774 - 0.0000i-0.5773 0.5774 0.5774-0.5774 0.5773 - 0.0000i 0.5773 + 0.0000id =1.0000 0 00 1.0000 + 0.0000i 00 0 1.0000 - 0.0000i det(v)ans =-5.0

39、566e-028 -5.1918e-017i %v的行列式接近0, 特征向量线性相关，不可对角化 v,d=jordan(a)v =1 0 11 0 01 -1 0d =1 1 00 1 10 0 1 %jordan标准形不是对角的，所以不可对角化 %Exercise 7(3) A=5 7 6 5;7 10 8 7;6 8 10 9;5 7 9 10A =5 7 6 57 10 8 76 8 10 95 7 9 10 v,d=eig(A)v =0.8304 0.0933 0.3963 0.3803-0.5016 -0.3017 0.6149 0.5286-0.2086 0.7603 -0.271

40、6 0.55200.1237 -0.5676 -0.6254 0.5209d =0.0102 0 0 00 0.8431 0 00 0 3.8581 00 0 0 30.2887 inv(v)*A*vans =0.0102 0.0000 -0.0000 0.00000.0000 0.8431 -0.0000 -0.0000-0.0000 0.0000 3.8581 -0.0000-0.0000 -0.0000 0 30.2887%本题用jordan不行, 原因未知 %Exercise 7(4)参考6(4)和7(1), 略 %Exercise 8 只有(3)对称, 且特征值全部大于零, 所以是正

41、定矩阵. %Exercise 9(1) a=4 -3 1 3;2 -1 3 5;1 -1 -1 -1;3 -2 3 4;7 -6 -7 0 rank(a)ans =3 rank(a(1:3,:)ans =2 rank(a(1 2 4,:) %1,2,4行为最大无关组ans =3 b=a(1 2 4,:);c=a(3 5,:); bc %线性表示的系数ans =0.5000 5.0000-0.5000 1.00000 -5.0000 %Exercise 10 a=1 -2 2;-2 -2 4;2 4 -2 v,d=eig(a)v =0.3333 0.9339 -0.12930.6667 -0.3

42、304 -0.6681-0.6667 0.1365 -0.7327d =-7.0000 0 00 2.0000 00 0 2.0000 v*v1.0000 0.0000 0.00000.0000 1.0000 00.0000 0 1.0000 %v确实是正交矩阵 %Exercise 11%设经过6个电阻的电流分别为i1, ., i6. 列方程组如下%20-2i1=a; 5-3i2=c; a-3i3=c; a-4i4=b; c-5i5=b; b-3i6=0;%i1=i3+i4;i5=i2+i3;i6=i4+i5;%计算如下 A=1 0 0 2 0 0 0 0 0;0 0 1 0 3 0 0 0

43、0;1 0 -1 0 0 -3 0 0 0; 1 -1 0 0 0 0 -4 0 0;0 -1 1 0 0 0 0 -5 0;0 1 0 0 0 0 0 0 -3; 0 0 0 1 0 -1 -1 0 0;0 0 0 0 -1 -1 0 1 0;0 0 0 0 0 0 -1 -1 1;b=20 5 0 0 0 0 0 0 0; Abans =13.34536.44018.54203.3274-1.18071.60111.72630.42042.1467 %Exercise 12 A=1 2 3;4 5 6;7 8 0; left=sum(eig(A), right=sum(trace(A)le

44、ft =6.0000right =6 left=prod(eig(A), right=det(A) %原题有错, (-1)应删去left =27.0000right =27 fA=(A-p(1)*eye(3,3)*(A-p(2)*eye(3,3)*(A-p(3)*eye(3,3)fA =1.0e-012 *0.0853 0.1421 0.02840.1421 0.1421 0-0.0568 -0.1137 0.1705 norm(fA) %f(A)范数接近02.9536e-013第四章%Exercise 1(1)roots(1 1 1) %Exercise 1(2)roots(3 0 -4 0

45、 2 -1) %Exercise 1(3)p=zeros(1,24);p(1 17 18 22)=5 -6 8 -5;roots(p) %Exercise 1(4)p1=2 3;p2=conv(p1, p1);p3=conv(p1, p2);p3(end)=p3(end)-4; %原p3最后一个分量-4roots(p3) %Exercise 2fun=inline(x*log(sqrt(x-1)+x)-sqrt(x-1)-0.5*x);fzero(fun,2) %Exercise 3fun=inline(x-2);fplot(fun,-2 2);grid on;fzero(fun,-1),fz

46、ero(fun,1),fminbnd(fun,0.5,1.5) %Exercise 4fun=inline(x*sin(1/x),x);fplot(fun, -0.1 0.1);x=zeros(1,10);for i=1:10, x(i)=fzero(fun,(i-0.5)*0.01);end;x=x,-x %Exercise 5fun=inline(9*x(1)+36*x(2)+4*x(3)-36;x(1)-2*x(2)-20*x(3);16*x(1)-x(1)-2*x(2)-16*x(3),x);a,b,c=fsolve(fun,0 0 0) %Exercise 6fun=(x)x(1)-

47、0.7*sin(x(1)-0.2*cos(x(2),x(2)-0.7*cos(x(1)+0.2*sin(x(2);a,b,c=fsolve(fun,0.5 0.5) %Exercise 7clear; close; t=0:pi/100:2*pi;x1=2+sqrt(5)*cos(t); y1=3-2*x1+sqrt(5)*sin(t);x2=3+sqrt(2)*cos(t); y2=6*sin(t);plot(x1,y1,x2,y2); grid on; %作图发现4个解的大致位置，然后分别求解y1=fsolve(x(1)-2)2+(x(2)-3+2*x(1)-5,2*(x(1)-3)+(x

48、(2)/3)-4,4,-4) %Exercise 8(1)clear;fun=inline(x.*sin(x.-x-2);fplot(fun,-2 2);grid on; %作图观察x(1)=-2;x(3)=fminbnd(fun,-1,-0.5);x(5)=fminbnd(fun,1,2);fun2=inline(-x.*sin(x.-x-2);x(2)=fminbnd(fun2,-2,-1);x(4)=fminbnd(fun2,-0.5,0.5);x(6)=2feval(fun,x)%答案: 以上x(1)(3)(5)是局部极小，x(2)(4)(6)是局部极大，从最后一句知道x(1)全局最小

49、， x(2)最大。 %Exercise 8(2)clear;fun=inline(3*x.-20*x.+10);fplot(fun,-3 3);grid on;%作图观察x(1)=-3;x(3)=fminsearch(fun,2.5);fun2=inline(-(3*x.-20*x.+10);x(2)=fminsearch(fun2,-2.5);x(4)=3;feval(fun,x) %Exercise 8(3)fun=inline(abs(x-x-x-2);fplot(fun,0 3);grid on;%作图观察fminbnd(fun,1.5,2.5)fun2=inline(-abs(x-x

50、-x-2);fminbnd(fun2,0.5,1.5) %Exercise 9close;x=-2:0.1:1;y=-7:0.1:1;x,y=meshgrid(x,y);z=y./9+3*x.*y+9*x.+y.+x.*y+9;mesh(x,y,z);grid on;%作图观察fun=inline(x(2)/9+3*x(1)*x(2)+9*x(1)+x(2)+x(1)*x(2)+9); x=fminsearch(fun,0 0)%求极小值fun2=inline(-(x(2)/9+3*x(1)*x(2)+9*x(1)+x(2)+x(1)*x(2)+9); x=fminsearch(fun2,0

51、-5)%求极大值 %Exercise 10clear;t=0:24;c=15 14 14 14 14 15 16 18 20 22 23 25 28 .31 32 31 29 27 25 24 22 20 18 17 16;p2=polyfit(t,c,2)p3=polyfit(t,c,3)fun=inline(a(1)*exp(a(2)*(t-14).),a,t);a=lsqcurvefit(fun,0 0,t,c)%初值可以试探f=feval(fun, a,t)norm(f-c)%拟合效果plot(t,c,t,f) %作图检验fun2=inline(b(1)*sin(pi/12*t+b(2)+20,b,t);%原题修改f(x)+20 b