MATLAB实验题答案

1、-0.33330.83331求以下变量的值，并在 MATLAB验证。( 1 ) a = 1 : 2 : 5a =1 3 5 ( 2 ) b = a' , a' , a' ；a b =1 1 1333555135( 3 ) c = a + b ( 2 , : )c =4 6 8>> result8=a.cAttempt to reference field of non-structure array.>> result9=a.b result9 =2.00002.00000.25000.6000>> result10=a.A2resu

2、lt10 =1 4 9162536>> result11=2.Aaresult11 =2、下列运算是否合法，为什么？如合法，结果是多少？>> result2=a*bError using *Inner matrix dimensions must agree.>> result3=a+bresult3 =3625>> result4=b*d result4 =811312222404913>> result5=b;c'*dresult5 =312222404913-5-87>> result6=a.*bresult

3、6 =2 4 816 32 643、用MATLAB求解下面的的方程组。7212x1491532x27(1)22115x3113213x40>> A=7 2 1 -2;9 15 3 -2;-2 -2 11 5;1 3 2 13>> B=4 7 -1 0>> B=B'>> x=inv(A)*Bxyz1x2yzw8(2)2xy3w 33x3y5z 6w 5>> A1=1 1 1 0;1 2 1 -1;2 -1 0 -3;3 3 5 -6>> B2=1;8;3;5>> x2=inv(A1)*B22 8-34 1

4、530>> result7=a./bresult7 =0.50000.5000-3.00004.00001.66671.2000721291532A221154、已知13213(1) 求矩阵A 的秩(rank)求矩阵 A的行列式(determinant)分为A, 8089分为B, 7079分为C, 6069求矩阵A的逆(inverse)分为D, 60分以下为E。对超出百分制范围(4) 求矩阵 A 的特征值及特征向量 (eigenvalue 的成绩，给出错误提示信息。and eigenvector)if 结构程序：>> A3=7 2 1 -2;9 15 3 -2;-2 -

5、2 11 5;1 3 2 13x=input('please enter score=');>> r=rank(A3)>> b=inv(A3)>> a=det(A3)>> V,D=eig(A3)if x>=90&&x<=100disp('A')elseif x<90&&x>=80 disp('B')elseif x<80&&x>=702n2 102910n 10,求 y= ?(运行 format long g 查看

6、y 的值)m1=0;for m=-10:10m仁 m1+2m;endm1m1 =2047.99902343756、求分段函数的值。2 xx6x0y2 x5x 6 0x52 xx1x5、 disp('C') elseif x<70&&x>=609102d2isp('D')elseif x<60&&x>=0命令后，disp('E')elsedisp('error')endswitch 结构程序： x=input('please enter score=');swi

7、tch fix(x/10) case10,9 if x>100 disp('error') else disp('A') endcase8disp('B')case7disp('D') case0,1,2,3,4,5 disp('E') otherwisedisp('C')用if语句实现，算出下列表中x对应的y值。case6 x=input('enter x='); if x<0y=xA2+x-6;elseif x>=0&&x<5y=xA2-5

8、*x+6;elsey=xA2-x-1;endy7、分别用 if 和 switch 语句实现，将百分制成绩转换为成绩等级A、B、C、D E。其中90100disp('error')end8、思考题 设计程序，完成成两位数的加、减、乘、除四则运算， 即产生两个两位随机整数， 再输入一 个运算符号， 做相应的运算， 并显示相应的 结果。x=input('请输入运算符') a=n um2str(floor(ra nd(1)*90+10); a值分别是多少？a=input('请输入数值') n=0;b=nu m2str(floor(ra nd(1)*90+

9、10); bif x='+'y=a+b;elseif x='-'y=a-b;elseif x='*'y=a*b;elseif x='/'y=a/b;elsefor m=1:100n=n+1/(m*m); endn=6* n;x=sqrt (n);x11、编程求满足值。m kk12k10000的最小mdisp('error') end9、启动 MATLAB，点击 File|New|M-File动 MATLAB的程序编辑及调(Editor/Debugger )，编辑以下程序，启试器点击sum=0;m=2;a=1;whi

10、le afor i=1:mFile|Save保存程序，注意文件名最好用英 文字符。点击Debug|Run运行程序，在命令 窗口查看运行结果，程序如有错误则改正。 注：数论中一个有趣的题目：任意一个正整数，若为偶数，则用2除之，若为奇数，则与 相乘再加上1。重复此过程，最终得到的结 果为1。n=input('请输入 n 值:');a=n;while n>1if rem(n, 2)=0n=n/2;elsesum=sum+2Ai;if sum>10000a=0;endendm=m+1;end3nd12、已思考题知 y的函数y(t)t2 /2!t3/3!t20 / 20!n

11、=3*n+1;enda=a, n;x2根据当n分别取100、1 丄2232.2.4.6.8.0求下面表格中与对应的y值tt=input('请输入 t 值：) sum=0;for i=1:20b=factorial(i); sum=sum+tAi/b;endsum=sum+1;1000、10000 时，求 x 的，得到的x为200个(0, 1)15、求13、编写一个函数，计算下面函数的值， 量X的值，调用该函数后，返回fun ctio n y=myfu n1(x)给出标y的值。sin x,X0y(x)x,0x 3x 6,X3选择一些数据测试你编写的函数。fun ctio ny=myfu

12、n1(x)if x<=0y=si n( x);elseif x>0& x<=3y=x;elseif x>3y=_x+6;end14、编写一个函数求向量x中元素的平均值、最大值、最小值、均方根值。fun ctio n m_x,max_x,min_x,rms_x=myfu n2(x)方均根值(Root Mean Square) 的计算公式为：用下面数据测试你写的函数：(1) x=sin(0:0.01:6*pi)(2) x=rand(1,200)之间均匀分布的随机数。fun ctio n m_x,max_x,min_x,rms_x=myfu n2(x) % 平均值su

13、m_x=sum(x); %向量元素求和m,n=size(x);%最好用 n=length(x);m_x=sum_x/ n;%求最大值采用逐个比较方式if x(1)>x(2)max_x=x(1);elsemax_x=x(2);endfor k=3: nif max_x<x(k) max_x=x(k);elsemax_x=max_x; %可省略endend%求最小值if x(1)<x(2)min_x=x(1);elsemin_x=x(2);endfor k=3:nif min_x>x(k) min_x=x(k);elsemin_x=min_x;%可省略endend%求均方根

14、值sum_x2=0;for k=1: nsum_x2=sum_x2+x(k).A2;rms_x=sqrt(sum_x2/n);endm_x;max_x;min_x;rms_x;%按照函数值行参顺序输出结果编写一个函数，给出一个向量 XxX，生成如下范德蒙矩阵。function v=myva nder(x)111Xix2xn222X1X2Xnn 1n1n1X1X2Xn例如：>>v=myva nder(2 3 4 5)得v=1111234549162582764125function v=myva nder(x) v1=va nder(x);%生成范德蒙矩阵v2=v1：v=flipud

15、(v2);%实现矩阵上下翻转编写程序，用如下迭代公式求；a，a的值分别18、编写程序，选择合适的步距，绘制下面函数16、思考题17、在同一坐标系下绘制下面三个函数在t 0 ,4的图象。y 1 ty 2、ty 34 e 0.1t sin(t)>> t=O:O.1:4*pi;>> y1=t;>> y2=tA(0.5);>> y2=sqrt(t);>> y3=4*pi*eps(-0.1*t).*si n(t);>> plot(t,y1,t,y2,t,y3)为：3，17113。迭代的终止条件为在区间-6，6中的图象。Xn 1Xn1

16、0，迭代初值x°1.0 ,sin x,迭代次数不超过100次。分别对迭代结果和y(x)x,准确值进行比较，并统计迭代次数。X 00x3Xn 1Xna22xn注：x2a2x2 x2 ax 6,function y=fun 1(x)2=; axn ax2or x2xxn1 7 2xnif x0<=0fun ctio nx, n=sqrt_a(a)x=1.0;for k=1:100m=x;x=x/2+a/(2*x);if abs(x-m)<=10A(-5) break endendx;n=k; s=(x-sqrt(a);if s<=10A(-5) disp('正确

17、'); elsedisp('错误');endy=y,s in (x0);elseif x0>0&& x0<=3y=y,x0；elseif x0>3 y=y,6-x0;endend>> x=-6:0.1:6;>> y=fun 1(x);>> plot(x,y)19、用compass函数画下面相量图ua = 1 ; ub = cos(-2*pi/3)+sin(-2*pi/3)*i;uc=cos(2*pii/3)+si n(2*pi/3)*i;compass(ua,ub,uc,ua-ub,ub-uc,uc-

18、ua)>> ua=1;>> ob=cos(-2*pi/3)>> ob=cos(-2*pi/3)+sin(-2*pi/ 3)*i;>> ub=cos(-2*pi/3)+sin(-2*pi/ 3)*i;>> uc=cos(2*pi/3)+si n( 2*pi/ 3)*i;>> compass(ua,ub,uc,ua-ub,ub-uc,uc-ua)20、三维空间曲线绘制z=0:0.1:4*pi;x=cos(z);y=s in(z);plot3(x,y,z)>> z=0:0.1:4*pi;>> x=cos(

19、z);>> y=si n( z);>> plot3(x,y,z)21、用mesh或surf函数，绘制下面方程所表示的三维空间曲面，x和y的取值范围设为 -3，3。>> x=-3:0.01:3;>> y=-3:0.01:3;>> x,y=meshgrid(x,y);>> z=-x.*x/10+y.*y/10;>> mesh(x,y,z)>> z=x.*x/10+y.*y/10;>> mesh(x,y,z)22、思考题用对分法求解方程2e %sinx 在01内的解，并验证，在程序中统计出对分

20、次数。 提示：先将原f (x) 2e对分法的基本思想是：f(x1)*f(x2)<0方程转化x sin x0的形式。解。取该区间的中点f(x1)和 f(x2)若 f(x1)*f(xm)<0x1,xm一个一元方程f(x)=0 ，若 ，则在x1,x2区间内有实数xm=(x1+x2)/2，判定 者中哪一个与f(xm)异i=0;while y=0i=i+1;x3=(x0+x1)/2;y=2*eps(-x3)-si n(x3);y0=2*eps(-x0)-si n(x0);y仁 2*eps(-x1)-si n(x1);if y*y0<0 x0=x3;elseif y*y1<0x仁

21、x3;elseif y=0x3=x3;endif i=1000000y=0;endEnd23、 求 代 数 方 程3x54x4 7x32x29x 12 0 的 5个根，并将其用星号(*)标记在复平面图上。(用roots 和plot函数)>> a=3,4,7,2,9,12;>> y=roots(a)> plot(y,'*')求代数方程X510的5个根，并将其用星号(* )标记在复平面图上。(用roots和plot 函数) 程序：>> a仁1,0,0,0,0,-1;，则解存在的区间缩小为 否则解存在的区间缩小为号，>> y1=r

22、oots(a1)>> plot(y1,'*')xm,x2。重复这样的步骤，直到区间的长 24、求下面函数在0.5,4区间内的过零点。(用度小于一个可以接受的小数(比如1e-10), fzero 函)则认为中点即是原方程的解。>> x0=0;32f (x) x 2x sin(x) 5x cos(x)x1=1; y=1;>> x=0.5:0.01:4;>> y2=x.A3-2*x.A2.*si n(x)+5*x.*cos(x)+1./x;>> plot(x,y2)函数：fun ctio n y2=fz(x)y2=x.A3-2*x.A2.*si n(x)+5*x.*cos(x)+1./x;分析：根据图像可知零点在 x=1和x=