实验内容作业

1、1程序：E=eye(3);R=rand(3,2);O=zeros(2,3);S=diag(1,2);A=E,R;O,S; .B=A*A,C=E R+R*S;O S*S,a=B=C结果：B = 1.0000 0 0 0.8205 1.0586 0 1.0000 0 1.7873 2.4395 0 0 1.0000 0.1158 0.0296 0 0 0 1.0000 0 0 0 0 0 4.0000 C = 1.0000 0 0 0.8205 1.0586 0 1.0000 0 1.7873 2.4395 0 0 1.0000 0.1158 0.0296 0 0 0 1.0000 0 0 0 0

2、 0 4.0000 a = 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 12.程序：A=magic(4),s1=sum(A),s2=sum(A'),.s3=sum(diag(A),s4=sum(diag(fliplr(A) 结果：A = 16 2 3 13 5 11 10 8 9 7 6 12 4 14 15 1 s1 = 34 34 34 34 s2 = 34 34 34 34 s3 = 34 s4 = 34 3.程序：x=-2:0.5:2;y1=-1./(1+x.2);y2=exp(-x.2/2);y3=sin(2*x);y4=s

3、qrt(4-x.2);a=y1+y2,b=y1.*y2,c=y3./y1,d=(5*y4-y1)./y2.2结果：a = -0.0647 0.0170 0.1065 0.0825 0 0.0825 0.1065 0.0170 -0.0647 b = -0.0271 -0.0999 -0.3033 -0.7060 -1.0000 -0.7060 -0.3033 -0.0999 -0.0271 c = -3.7840 0.4586 1.8186 1.0518 0 -1.0518 -1.8186 -0.4586 3.7840 d = 10.9196 65.6748 24.9002 13.4597 1

4、1.0000 13.4597 24.9002 65.6748 10.91966.程序1： x=1:0.01:3;y1=x.2;y2=x.3;y3=x.4;y4=x.5;plot(x,y1,x,y2,x,y3,x,y4)gtext('x2'),gtext('x3'),gtext('x4'),gtext('x5')结果1： 程序2：x=1:0.01:3;y1=x.2;y2=x.3;y3=x.4;y4=x.5;subplot(2,2,1),plot(x,y1),title('x2')subplot(2,2,2),plot

5、(x,y2),title('x3')subplot(2,2,3),plot(x,y3),title('x4')subplot(2,2,4),plot(x,y4),title('x5')gtext('x2'),gtext('x3'),gtext('x4'),gtext('x5')结果： 7程序1：x=1:0.01:3;y1=x.2;y2=8*log(x);plot(x,y1,x,y2),gtext('x2'),gtext('8Inx')结果1： 故方程x

6、2 = 8lnx的近似根为x1=1.1957,x2=2.9348 程序2：x=-4:0.01:4;y=4*sin(x)-x-2;z=0*x;plot(x,y,x,z),axis(-4 4 0 1),gtext('4sinxx2')结果2： 故方程4sinx x 2 = 0的近似根为：x1=1.8476,x2=0.7623,x3=-2.9117.8.程序：>> x=-10:0.01:10;y=x;X,Y=meshgrid(x,y);>> Z=X.2-Y.2;>> mesh(X,Y,Z)结果： 10.程序1：建立函数M文件如下：function

7、m=jiecheng(n)p=1;for i=1:n p=p*i;endm=p; 结果1：m=jiecheng(5) m = 120 程序2：建立函数M文件如下：function q=zuhe(n,m)p=1;if (n>=m && n>=0 & m>=0) m=min(n-m,m); for i=1:m p=p*(n-i+1)/i; endendq=p;结果2：q=zuhe(5,3) q = 10 程序3：建立函数M文件如下：function p=qiuhe(p1,p2)a=length(p1);b=length(p2);if a<b p1=zeros(1,b-a),p1;else