matlab实验.doc

实验一 1、计算y=x3+(x-0.98)2/(x+1.35)3-5*(x+1/x)当x=2，x=4的值。2、计算cos（pi/3）-（9-2(1/2)）(1/3)。3、已知a=3，A=4,b=a2,B=b2-1,c=a+2*B+c,求C。4、创建一个3*3矩阵，然后用矩阵编辑器将其扩充为4*5矩阵。5、创建一个4*4矩阵魔方阵和相应的随机矩阵，将两个矩阵并接起来，然后提取任意两个列向量。7设A=1,4,8,13;-3,6,-5,-9;2,-7,-12,-8,B=5,4,3,-2;6,-2,3,-8;-1,3,-9,7,求C=A*B,D=A.*B。8、求（s2+2）（s+4）（s+1）/（s3+s+1）的“商”及“余”多项式。9、建立矩阵A,然后找出在10，20q区间的元素的位置。10、创建一个有7个元素的一维数组，并做如下处理：1）直接访问一位数组的第6个元素；2）寻访一维数组的第1、3、5个元素；3）寻访一维数组中的第4个至最后1个元素；4）寻访一维数组中大于70的元素。三、实验程序：NO.1 x=2,4x = 2 4 y=x.3+(x-0.98).2/(x+1.35).3-5*(x+1./x)y = -4.4423 42.8077NO.2 cos(pi/3)-(9-sqrt(2)(1/3)ans = -1.46493 a=3a = NO. 3 A=4A = 4 b=a2b = 9 B=b2-1B = 80 c=a+A-2*Bc =-153 C=a+2*B+cC =10NO.4A=rand(3)A = 0.5341 0.8385 0.7027 0.7271 0.5681 0.5466 0.3093 0.3704 0.4449 B=rand(3,2)B = 0.6946 0.9568 0.6213 0.5226 0.7948 0.8801 C=rand(1,5)C = 0.1730 0.9797 0.2714 0.2523 0.8757 F=A,B;CF = 0.5341 0.8385 0.7027 0.6946 0.9568 0.7271 0.5681 0.5466 0.6213 0.5226 0.3093 0.3704 0.4449 0.7948 0.8801 0.1730 0.9797 0.2714 0.2523 0.8757NO.5 A=magic(4)A = 16 2 3 13 5 11 10 8 9 7 6 12 4 14 15 1 B=rand(3,4)B = 0.7373 0.8939 0.6614 0.0648 0.1365 0.1991 0.2844 0.9883 0.0118 0.2987 0.4692 0.5828 C=A;BC = 16.0000 2.0000 3.0000 13.0000 5.0000 11.0000 10.0000 8.0000 9.0000 7.0000 6.0000 12.0000 4.0000 14.0000 15.0000 1.0000 0.7373 0.8939 0.6614 0.0648 0.1365 0.1991 0.2844 0.9883 0.0118 0.2987 0.4692 0.5828 C(:,1,2)ans =NO.6 A=-4,-2,0,2,4;-3,-1,1,3,5A = -4 -2 0 2 4 -3 -1 1 3 5 LA=abs(A)3LA = 1 0 0 0 1 0 0 0 0 1 A(LA)ans = -4 4 5NO.7 A=1,4,8,13;-3,6,-5,-9;2,-7,-12,-8A = 1 4 8 13 -3 6 -5 -9 2 -7 -12 -8 B=5,4,3,-2;6,-2,3,-8;-1,3,-9,7B = 5 4 3 -2 6 -2 3 -8 -1 3 -9 7 C=A*BC = 19 -82 30 12 27 3 -38 54 29 D=A.*BD = 5 16 24 -26 -18 -12 -15 72 -2 -21 108 -56NO.8 den=conv(1,0,2,conv(1,4,1,2)den = 1 6 10 12 16num=1,0,1,1; q,r=deconv(den,num)q = 1 6r = 0 0 9 5 10NO.9 A=10,15,8;5,16,35;16,8,26A = 10 15 8 5 16 35 16 8 26 (A=10)&(ANO.10x = 1 2 3 4 5 6 7 x(1)ans =1 x(3)ans =3 x(5)ans = 5 x(4:end)ans = 4 5 6 7 x(x70)ans = Empty matrix: 1-by-0 实验二1. 设，把区间分为125点，画出以x为横坐标，y为纵坐标的曲线。2. 设，要求在区间内画出三维曲线。3. 设，求定义域内的z值（网格取0.1见方），并绘制三维曲面。4. 设，画出的三维曲面图，并叠在上题的图中。5. 设，若，在4个子图中分别画出其曲线。6. 绘制空间曲线：7. 绘制的极坐标图。8. 在区间内，绘制曲线。三、实验程序：（1）x=linspace(0,2*pi,125);y=cos(x).*(0.5+3*sin(x)./(x.2+1);plot(x,y)（2）t=-45:0.01:45;x=t.*sin(3*t);y=t.*cos(3*t);z=t;plot3(x,y,z)（3）x=-2:0.1:2;y=x;X,Y=meshgrid(x,y);Z=(X.2)*exp(-(X.2+Y.2);surf(Z)（4）x=-2:0.1:2;y=x;X,Y=meshgrid(x,y);Z=X.2.*exp(-(X.2+Y.2);surf(Z);hold onZ2=0.05*X-0.05*Y+0.1;surf(Z2)（5）N=2;t=0:pi/20:2*pi;x=cos(t);a=0;y=sin(N*t+a);subplot(2,2,1);plot(x,y)a=pi/3;y=sin(N*t+a);subplot(2,2,2);plot(x,y)a=pi/2;y=sin(N*t+a);subplot(2,2,3);plot(x,y)a=pi;y=sin(N*t+a);subplot(2,2,4);plot(x,y)（6）t=0:0.01:2*pi;x=8*sin(t);y=4*sqrt(2).*cos(t);z=-4*sqrt(2).*cos(t);plot3(x,y,z)（7）t=0:pi/100:2*pi;y=sin(2*t).*cos(2*t);polar(t,y)（8）x=0:pi/100:2*pi;y=2*exp(-0.5*x).*sin(2*pi*x);plot(x,y)实验三MATLAB的语法规则简洁，编程效率高，作为一个完整的程序语言，MATLAB也有其各种程序流程控制、文件格式和函数调用的规则，通过函数的调用就能够组成庞大的程序。1. 一个3位整数各位数字的立方和等于该数本身则称该数为水仙花数。输出100999之间的全部水仙花数。2. 编写函数文件求半径为r的圆的面积和周长。3. 已知，当时，求y的值。4. 利用函数的递归调用，求 。5. 编制程序，计算12+n 153 370 371 407NO.2新建mj.m文件function f=mj(n)f=pi*n*n;新建zc.m文件function f=zc(n)f=2*pi*n; mj(2)ans = 12.5664 zc(2)ans = 12.5664NO.3y=0;for i=1:100 y=y+1./i.2;enddisp(y) 1.6350NO.4新建dg.m文件function f=dg(n)f=0;if n=1 f=1;elseif n=0 f=1;else f=n*dg(n-1); end dg(6)ans = 720NO.5s=0;i=1;while sn =62NO.6(1) k=0;for i=0:63 k=k+2.i;enddisp(k)1.8447e+019(2) k=0;i=1;while i t=446 714 950 1422 1634;p=7.04 4.28 3.40 2.54 2.13;P1=interp1(t,p,500)P2=interp1(t,p,1000)P3=interp1(t,p,1500)P1 = 6.4839P2 = 3.3089P3 =2.3892NO.2 x=1:0.1:1.4; y=0.25 0.2268 0.2066 0.1890 0.1736; p=polyfit(x,y,3); dx=0.001; dy=diff(p,3); df=dy./dx; xx=1.0,1.2; yy=polyval(df,xx)yy = 1.0e+003 * 4.2767 4.2767NO.3a=3,4,7,2,9,12;r=roots(a)r = -0.8612 + 1.4377i -0.8612 - 1.4377i 0.6737 + 1.0159i 0.6737 - 1.0159i -0.9583NO.4 X=3,4,7,2,9,12; poly(X)ans = 1 -37 533 -3803 14106 -25776 18144NO.5 a=2,4,9;4,2,4;9,4,18;v,d=eig(a)v = -0.8483 -0.3290 0.4149 0.4514 -0.8589 0.2419 0.2767 0.3925 0.8771d = -3.0645 0 0 0 1.7042 0 0 0 23.3603NO.6 s=inline(x-y).2.*sin(x+y).2,x,y);q=dblquad(s,pi,2*pi,0,pi)q = 54.3546NO.7 a