薛山matlab基础教程第三版习题解答3

MATLAB课后习题集附标准答案

第2章MATLAB概论1、与其他计算机语言相比较，MA TLAB语言突出的特点是什么？答：起点高、人机界面适合科技人员、强大而简易的作图功能、智能化程度高、功能丰富，可扩展性强.2、MA TLAB系统由那些部分组成？答：开发环境、MATLAB数学函数库、MATLAB语言、图形功能、应用程序接口3、安装MATLAB时，在选择组件窗口中哪些部分必须勾选，没有勾选的部分以后如何补安装？答：在安装MATLAB时，安装内容由选择组件窗口中各复选框是否被勾选来决定，可以根据自己的需要选择安装内容，但基本平台（即MATLAB选项）必须安装.第一次安装没有选择的内容在补安装时只需按照安装的过程进行，只是在选择组件时只勾选要补装的组件或工具箱即可.4、MATLAB操作桌面有几个窗口？如何使某个窗口脱离桌面成为独立窗口？又如何将脱离出去的窗口重新放置到桌面上？答：在MATLAB操作桌面上有五个窗口，在每个窗口的右下角有两个小按钮，一个是关闭窗口的Close 按钮，一个是可以使窗口称为独立的Undock按钮，点击Undock按钮就可以使该窗口脱离桌面称为独立窗口，在独立窗口的view菜单中选择Dock，菜单项就可以将独立的窗口重新防止的桌面上.5、如何启动M文件编辑/调试器？答：在操作桌面上选择“建立新文件”或“打开文件”操作时，M文件编辑/调试器将被启动.在命令窗口中键入edit命令时也可以启动M文件编辑/调试器.6、存储在工作空间中的数组能编辑吗？如何操作？答：存储在工作空间的数组可以通过数组编辑器进行编辑：在工作空间浏览器中双击要编辑的数组名打开数组编辑器，再选中要修改的数据单元，输入修改内容即可.7、命令历史窗口除了可以观察前面键入的命令外，还有什么用途？答：命令历史窗口除了用于查询以前键入的命令外，还可以直接执行命令历史窗口中选定的内容、将选定的内容拷贝到剪贴板中、将选定内容直接拷贝到M文件中.8、如何设置当前目录和搜索路径，在当前目录上的文件和在搜索路径上的文件有什么区别？答：当前目录可以在当前目录浏览器窗口左上方的输入栏中设置，搜索路径可以通过选择操作桌面的file菜单中的Set Path菜单项来完成.在没有特别说明的情况下，只有当前目录和搜索路径上的函数和文件能够被MATLAB运行和调用，如果在当前目录上有与搜索路径上相同文件名的文件时则优先执行当前目录上的文件，如果没有特别说明，数据文件将存储在当前目录上.9、在MA TLAB中有几种获得帮助的途径？答：（1）帮助浏览器：选择view菜单中的Help菜单项或选择Help菜单中的MATLAB Help菜单项可以打开帮助浏览器.（2）help命令：在命令窗口键入“help”命令可以列出帮助主题，键入“help 函数名”可以得到指定函数的在线帮助信息.（3）lookfor命令：在命令窗口键入“lookfor 关键词”可以搜索出一系列与给定关键词相关的命令和函数.（4）模糊查询：输入命令的前几个字母，然后按Tab键，就可以列出所有以这几个字母开始的命令和函数.注意：lookfor和模糊查询查到的不是详细信息，通常还需要在确定了具体函数名称后用help命令显示详细信息.第3章 MATLAB 数值运算3.1在MA TLAB 中如何建立矩阵⎥⎦⎤⎢⎣⎡194375，并将其赋予变量a ？答：在Command Window 窗口输入操作:>> a=[5 7 3;4 9 1]3.2有几种建立矩阵的方法？各有什么优点？答：（1）直接输入法，如a=[1 2 3 4]，优点是输入方法方便简捷；（2）通过M 文件建立矩阵，该方法适用于建立尺寸较大的矩阵，并且易于修改；（3）由函数建立，如y=sin(x)，可以由MATLAB 的内部函数建立一些特殊矩阵；（4）通过数据文件建立，该方法可以调用由其他软件产生数据.3.3在进行算术运算时，数组运算和矩阵运算各有什么要求？答：进行数组运算的两个数组必须有相同的尺寸.进行矩阵运算的两个矩阵必须满足矩阵运算规则，如矩阵a 与b 相乘（a*b ）时必须满足a 的列数等于b 的行数.渗釤呛俨匀谔鱉调硯錦。

matlab第三章课后部分答案

matlab第三章课后部分答案习题三3-2 从键盘输入一个三位整数，将它反向输出。

如输入639，输出为936程序如下：m=input('请输入一个三位整数：');m1=fix(m/100);%求m的百位整数m2=rem(fix(m/10),10);%求m的十位数字m3=rem(m,10);%求m的个位数字m=m3*100+m2*10+m1%反向输出m3-3 输入一个百分制成绩，要求输出成绩等级A,B,C,D,E。

其中90~100分为A，80~89分为B，70~79分为C，60~69分为D，60分以下为E。

要求：（1）分别用if语句和switch语句实现。

（2）输入百分制成绩后要判断该成绩的合理性，对不合理的成绩应输出出错信息。

程序如下：（1）if语句c=input('请输入成绩：');if c>=90&c<=100disp('A 成绩合理');elseif c>=80&c<=89disp('B 成绩合理');elseif c>=70&c<=79disp('C 成绩合理'); elseif c>=60&c<=69disp('D 成绩合理'); elseif c<60disp('E 成绩合理');elsedisp('成绩错误');end（2）switch语句c=input('请输入成绩：'); switch fix(c)case num2cell(90:100)disp('A 成绩合理'); case num2cell(80:89)disp('B 成绩合理'); case num2cell(70:79)disp('C 成绩合理'); case num2cell(60:69)disp('D 成绩合理'); case num2cell(0:59)disp('E 成绩合理');x=fix(rand(1,20)*89)+10;x1=fix(sum(x)/20);disp(['平均数是：',num2str(x1)])m=(rem(x,2)==0&x<x1);n=find(m);disp(['小于平均数的数是：',num2str(x(n))]); 3-6 输入20个数，求其中最大数和最小数。

Matlab基础与应用教程答案

第1章一、思考题4.( 1) B=A(2:5,1:2:5)(2)A(7)=[](3)A=A+30(4)size(A)( 5) t(find(t==0))=eps(6)t=reshape(x,3,4)(7)abs('matlab')(8)char(93)5.A=[97,67,34,10;-78,75,65,5;32,5,-23,-59;0,-12,54,7] (1)B=A(1:3,:)C=A(:,1:2)D=A(2:4,3:4)E=B*C(2)A(find(A>=50 & A<=100)) 二、实验题=-74/180*pi;y=-27/180*pi;sin(x*x+y*y)/sqrt(tan(abs(x+y)))+pi一・・,exp*a).*sin(a+=[2,4;,5];log(x+sqrt(1+x.*x))/24.A=[3,54,2;34,-45,7;87,90,15];B=[1,-2,67;2,8,74;9,3,0];(1)A*Bans =129 432 41977 -407 - 1052402 591 12489A.*Bans =3 -108 13468 -360 518783 270 0(2) A A3ans =-28917 240246 -4368137883 -259101 27669171333 252504 38673A.A3ans =27 157464 839304 -91125 343658503 729000 3375( 3) A/Bans =B/A ans =（4）[A,B]ans =3 54 2 1 -2 6734 -45 7 2 8 7487 90 15 9 3 0[&[1,3],:)砂2]ans =3 54 287 90 15600 183 -81684 282 72615 6 825=1+2i;b=3+4i; c=exp(pi*i/6) c =+ c+a*b/(a+b) ans =+ 第2章一、思考题=0; for n=0:63s=s+2A n;end disp(s)n=0:63;s=sum(2.An)二、实验题1.x=input( ' 输入一个四位整数：' );y=[fix(x/1000),mod(fix(x/100),10),mod(fix(x/10),10),mod(x,1 0)] z=mod((y+7),10)x=z(3)*1000+z(4)*100+z(1)*10+z(2)2. gh=input( ' 输入工号' ); h=input( ' 输入工时' );dj=84;if h>120gz=dj*120+*dj*(h-120);elseif h<60gz=dj*h-700;else gz=dj*h;end format bank; display([gh,gz])3. 循环结构n=input( 'input n:' );s=0;for k=1:n s=s+1/k A2;end display(sqrt(s*6)) 向量运算n=input( 'input n:' ); k=1:n;display(sqrt(sum(1./k.A2)*6))4.y=0;k=0;while y<3k=k+1; y=y+1/(2*k-1);end display([k-1,y-1/(2*k-1)])5.x0=0;x=1;k=0;a=input( 'a=' ); b=input( 'b=' );while abs(x-x0)>=1e-5 && k<500 x0=x;x=a/(b+x0);k=k+1;end display([k,x]); display([(-b+sqrt(bA2+4*a))/2,(-b-sqrt(bA2+4*a))/2]);6.y=fun(40)/(fun(30)+fun(20))(1)函数文件function f=fun(n) f=n+log(nA2+5);(2)函数文件function f=fun(n) a=1:n;f=sum(a.*(a+1));第3章一、思考题4.t=0::;y=sqrt(3)/2*exp(-4*t).*sin(4*sqrt(3)*t+pi/3);5.x=-10::10; y=linspace(-6,6,size(x,2))z=x.A3+3*x.*y.A2;plot3(x,y,z)6.x=100:100:400;y=100:100:400;z=[636,697,624,478;698,712,630,478;680,674,598,412;662,626,552,334]; [X,Y]=meshgrid(x,y);mesh(X,Y,z)二、实验题1.( 1 ) x=-10::10;plot(x,x-x.A3/3/2)( 2) plot(x,exp(-x.*x/2)/2/pi)(3)x=-8::8;plot(x,sqrt((64-x.*x)/2))( 4) t=0:pi/100:6*pi;plot(t.*sin(t),t.*cos(t))2.( 1) x1=linspace(0,1,100);y1=2*;t=linspace(0,pi,100);x=sin(3*t).*cos(t);y=sin(3*t).*sin(t);plot(x1,y1,'r-',x,y,'b:');text,1,'y='); text,,'x=sin(3t)cos(t)');text,,'y=sin(3t)sin(t)');(2)subplot(1,2,1);scatter(x1,y1,10)title('y=');subplot(1,2,2);scatter(x,y,10)3.subplot(1,2,1);x=1:1:100; y=sin(1./x);plot(x,y) subplot(1,2,2);fplot('sin(1/x)',[1,100])4.subplot(2,2,1);bar(t,y);subplot(2,2,2);stairs(t,y);subplot(2,2,3);stem(t,y) subplot(2,2,4); semilogy(t,y);3.theta=linspace(-pi,pi,100); ro=5.*cos(theta)+4; polar(theta,ro); fi=linspace(0,2*pi,100);a=1 r=a.*(1+cos(fi));polar(fi,r);4.plot3(x,y,z);t=0::1;x=t;y=t.A2;z=t.A3;plot3(x,y,z);7.x=-30::0;y=0::30; [x,y]=meshgrid(x,y);z=10.*sin(sqrt(x.A2+y.A2))./sqrt(1+x.A2+y.A2); meshc(x,y,z);8. x=linspace(-3,3,100);y=linspace(-3,3,100);[x y]=meshgrid(x,y);fxy=-5./(1+x.A2+y.A2);i=find(abs(x)<= & abs(y)<=; fxy(i)=NaN;surf(x,y,fxy)9.u=linspace(1,10,100);v=linspace(-pi,pi,100); [u,v]=meshgrid(u,v);x=3.*u.*sin(v);y=2.*u.*cos(v);z=4*u.A2;x=3*u.*sin(v);y=2*u.*cos(v);z=4*u.A2; surf(x,y,z);shading interp;light('position',[1,0,1]);10. t=0:pi/100:2*pi; y=sin(t);comet(t,y) 第4章一、思考题5.(1)A=eye(3);(2)C=100+(200-100)*rand(5,6);( 3) D=1+sqrt*randn(1,500);( 4) E=ones(size(A));( 5) A=A+30*eye(size(A));( 6) B=diag(diag(A))二、实验题1.P=pascal(5);H=hilb(5); Dp=det(P);Dh=det(H); Kp=cond(P);Kh=cond(H);P矩阵的性能更好，因为Kp较小2.A=[1,-1,2,3;0,9,3,3;7,-5,0,2;23,6,8,3] B=[3,pi/2,45;32,-76,sqrt(37);5,72,;exp(2),0,97] A1=diag(A);B1=diag(B);A2=triu(A);B2=triu(B); A3=tril(A);B3=tril(B); rA=rank(A);rB=rank(B);nA=norm(A);nb=norm(B); cA=cond(A);cB=cond(B);3.A=[31,1,0;-4,-1,0;4,-8,-2] ；[V,D]=eig(A);4.A=diag([-1,-1,-1,-1],-1)+diag([-1,-1,-1,-1],1)+diag([2,2,2,2,2])b=[1,0,0,0,0]';x1=inv(A)*b; x2=A\b; [L,U]=lu(A); x3=U\(L\b); [Q,R]=qr(a); [Q,R]=qr(A); x4=R\(Q\b) R=chol(A); x5=R\(R'\b) 5.B=sparse(A); x1=inv(B)*b; x2=B\b; [L,U]=lu(B); x3=U\(L\b); 第5章一、思考题3.A=randn(10,5); mean(A) std(A) max(max(A)) min(min(A)) sum(A,2) sum(sum(A)) sort(A,1) sort(A,2,'descend') 二、实验题1.A=rand(1,30000); mean(A) std(A) max(A) min(A) size(find(A>)/size(A)2.h=[466,715,950,1422,1635]; w=[,,,,];hh=[500,900,1500]; ww=interp1(h,w,hh,'spline')3.x=linspace(1,10,50); y=log(x);第6章一、思考题2.fx=i nlin e('1./(1+x.A 2)');[I,n]=quad(fx,-100000,100000,1e-3); [I,n]=quadl(fx,-100000,100000,1e-3); x=-100000::100000; y=1./(1+x.*x);f=polyfit(x,y,5); yy=polyval(f,x); plot(x,y,'r-',x,yy,'g.') 4. N=64; T=5;t=linspace(0,T,N); x=3*exp(-t); % dt=t(2)-t(1); f=1/dt; X=fft(x); F=X(1:N/2+1); f=f*(0:N/2)/N;plot(f,abs(F),'-*') % %% 求各采样点样本值 x% % %采样点数采样时间终点给岀N 个采样时间ti(l=1:N) 采样周期采样频率 (Hz)计算 x 的快速傅立叶变换 X% F(k)=X(k)(k=1:N/2+1)% 使频率轴 f 从零开始 % 绘制振幅 - 频率图xlabel('Frequency');ylabel('|F(k)|')5.(1)p1=[1 2 0 0 7];p2=[1 -2];p3=[1 0 0 5 1];p12=conv(p1,p2);p=p12+[zeros(1,size(p12,2)-size(p3,2)),p3]; roots(p)(2)A=[-1,4,3;2,1,5;0,5,6];Px=polyval(p,A) Pxm=polyvalm(p,A)6.(1) z=fzero('3*x-sin(x)+1',0)(2)建立函数文件 function F=myfun(X)x=X(1);y=X(2); F(1)=x*x+y*y-9;F(2)=x+y-1; 在命令窗口中输入以下命令： x=fsolve(@myfun,[3,0]',optimset('Display','of f'))trapz(x,y);3.(1)fx=inline('-2*y+2*x*x+2*x');[t,y]=ode23(fx,[0,],1)(2)fx=inline('y-exp(x)*cos(x)');[t,y]=ode23(fx,[0,3],1)二、实验题1.for x=1:3fx=[x,xA2,xA3;1,2*x,3*x;0,2,6*x]; diff(fx)end2.(1 ) x=0::1;y=x.A10+10.Ax+1./log10(x);dy=diff(y)/;(2)x=0::1;y=log(1+x);dy=diff(y,2)/;plot(x(1:99),dy)3.(1 ) fx=inline( 'x.A2.*sqrt(2*x.*x+3)');quad(fx,1,5)(2)fx=inline( 'x./sin(x).A2' ); quad(fx,pi/4,pi/3)(3)fx=inline( 'abs(cos(x+y))' ); dblquad(fx,0,pi,0,pi)(4)syms x y;fx=x*y;int(int(fx,yA2,y+2),-1,2)x的积分区间为【0, 2】时fx=inline( 'x.*y' );dblquad(fx,0,2,-1,2)4.x=::;y=[,,,,,,];trapz(x,y)5.(1)yp=i nlin e( '-+si n( 10*x))*y' );[t,y]=ode23(yp,[0,5],1);(2)令x1y, x2y ,x3 y''，则可写出原方程的状态方程形式: x1x2X2 X3cos 1x3cost 5 X3 X2 X1(t 1)2 3 sin t0 1 0 X1 0=>> x20 0 1 X2 01 / 5cos2t X3 12X3 cost3 sint (t 1)2建立函数文件fun cti on y=ztfu n( t,x)b=[0;0;cos(t)];y=[0,1,0;0,0,1;-1/(3+si n( t)),-1,5*cos(2*t)/(t+1)A2]*x+b; 解微分方程组[t,y]=ode23(@ztfu n,[0,5],[1;0;2]);6.建立函数文件fun cti on yy=ztfu n( t,y)yy=[ y( 2)*y(3);-y(1)* y(3) ;*y(1)*y(2)];解微分方程组[t,y]=ode23(@ztfu n,[0,5],[0;1;1])第7章一、思考题3.(1 )数值积分fx=i nlin e('exp(x).*(1+exp(x)).A2'); quad(fx,0,log(2)) 符号积分f=sym('exp(x)*(1+exp(x))A2');v=in t(f,0,log (2));eval(v)(2 )略二、实验题1.A=sym('[1,2,3;x,y, z; 3,2,1]')rank(A)inv(A)det(A)2.(1)y=sym('sqrt(x+sqrt(x+sqrt(x)))');y1=diff(y)y2=diff(y,'x',2)(2)syms x y;fxy=si n( x A2*y)*exp(-x A2-y);diff(diff(fxy,x),y)3.(1)syms xin t(1/(1+xA4))(2)syms x tin t((-2*x*x+1)/(2*x*x-3*x+1)A2,x,cos(t),exp(2*t))4.syms n xsymsum(1/(2* n+1)/(2*x+1)A(2* n+1), n,0,i nf)symsum(1/(2* n+1)/(2*x+1)A(2* n+1), n, 0,5)5.(1)syms xtaylor((exp(x)+exp(-x))/2,5,0)(2)syms a xtaylor(exp(-5*x)*si n(3*x+pi/3),5,a)6.(1)x=solve(sym('xA3+a*x+1=0'))(2)[x y]=solve(sym('sqrt(xA2+yA2)-100=0,3*x+5*y-8=0'))y' y17.方程转化为：’，严y1' y 1 一符号解[y1,y11]=dsolve(,Dy=y1,Dy1+y=1-t A2/pi,,,y(-2)=5,y1(-2)=5',,t,)数值解编写函数文件fun cti on yy=ztfu n( t,y)yy=[y(2);1-tA2/pi-y(1)];在命令窗口输入以下命令[t,y]=ode45(@ztfu n,[-2,7],[-5;5]);t=li nspace(-2,7,49)y2=y8.[x,y]=dsolve('Dx=3*x+4*y,Dy=-4*x+3*y','x(0)=0,y(0)=1')第9章二、实验题1.(1 )新建一个Blank GUI。

MATLAB基础教程薛山课后答案[2-12章].khda

2. (1) y e x (2) y 1
x5
பைடு நூலகம்
x 8
x sin x
>> y=@(x)exp(x)-x^5; >> x = fzero(y,8) x= 12.7132 >> fplot(y,[x-1,x+1]); >> hold on >> plot(x,y(x),'r*');
2 >> y=@(x)x*sin(x); x= 0
5 1
% find the times and places of the specified letter in the string letter = 'a'; string = 'China'; % The specified letter to be searched for % The specified tring to be searched
a 39
b 58
c 3
d
7
(1) a b (2) a c (3) a b & &b c (4) a d (5) a | b c (6) ~~ d 1 >> a=39; >> b=58;
>> c=3; >> d=7; >> a>b ans = 0 2 >> a<c ans = 0 3 >> a>b&&b>c ans = 0 4 >> a==d ans = 0 5 >> a|b>c ans = 1 6 >> ~~d ans = 1 7

matlab基础与应用部分习题答案

作业一4、写出完成下列操作的命令.(1）将矩阵A第2～5行中第1,3,5列元素赋给矩阵B。

>> M=［0：1：48]；〉〉A=reshape（M，7,7）A =0 7 14 21 28 35 421 8 15 22 29 36 432 9 16 23 30 37 443 10 17 24 31 38 454 11 18 25 32 39 465 12 19 26 33 40 476 13 20 27 34 41 48>> B=A（2:5,1：2：5）B =1 15 292 16 303 17 314 18 32（2）删除矩阵A的第七行元素。

〉〉A（7，:）=［］A =0 7 14 21 28 35 421 8 15 22 29 36 432 9 16 23 30 37 443 10 17 24 31 38 454 11 18 25 32 39 465 12 19 26 33 40 47 (3)将矩阵A的每个元素值加30.>〉A=A+30A =30 37 44 51 58 65 7231 38 45 52 59 66 7332 39 46 53 60 67 7433 40 47 54 61 68 7534 41 48 55 62 69 7635 42 49 56 63 70 77(4求矩阵A的大小和维素。

sizeA = size（A)dA = ndims（A）sizeA =6 7dA =2（5）将向量t的0元素用机器零来代替。

〉> t=[1 2 3 4 0 5]；t =1 2 3 4 0 5〉〉find（t==0）ans =5〉> t(5)=epst =1.0000 2。

0000 3。

0000 4。

0000 0.0000 5。

0000（6）将含有12个元素的向量x转换成3＊4矩阵.〉> x=［0：11]x =0 1 2 3 4 5 6 7 8 9 10 11 〉〉y=reshape（x，3，4）y =0 3 6 91 4 7 102 5 8 11（7）求一个字符串的ASCII。

matlab程序设计第三章课后习题答案资料

m a t l a b程序设计第三章课后习题答案1. p138 第6题在同一坐标轴中绘制下列两条曲线并标注两曲线交叉点。

>> t=0:0.01:pi;>> x1=t;>> y1=2*x1-0.5;>> x2=sin(3*t).*cos(t);>> y2=sin(3*t).*sin(t);>> plot(x1,y1,'r-',x2,y2,'g-')>> axis([-1,2,-1.5,1])>> hold on>> s=solve('y=2*x-0.5','x=sin(3*t)*cos(t)','y=sin(3*t)*sin(t)');>> plot(double(s.x),double(s.y),'*');截图：p366 第4题绘制极坐标曲线，并分析对曲线形状的影响。

function [ output_args ] = Untitled2( input_args )%UNTITLED2 Summary of this function goes here% Detailed explanation goes heretheta=0:0.01:2*pi;a=input('请输入a的值：');b=input('请输入b的值：');n=input('请输入n的值：');rho=a*sin(b+n*theta);polar(theta,rho,'k');end下面以a=1，b=1，n=1的极坐标图形为基础来分析a、b、n的影响。

对a的值进行改变：对比发现a只影响半径值的整倍变化对b的值进行改变：对比发现b的值使这个圆转换了一定的角度对n的值进行改变：对比发现当n>=2时有如下规律1、当n为整数时，图形变为2n个花瓣状的图形2、当n为奇数时，图形变为n个花瓣状的图形分别让n为2、3、4、5同时改变b和n的值，对比发现图形变为2n个花瓣状，并且还旋转了一定角度1 绘制sin(x)曲线，通过句柄修改曲线的颜色，要求：a）生成曲线的同时，获取句柄，再修改颜色；>> x=0:pi/50:2*pi;>> y=sin(x);>> h=plot(x,y)h =174.0101>> set(h,'color','r');截图如下：b）先生成曲线，然后再获取句柄，再完成颜色修改。

matlab基础与应用教程课后答案

matlab基础与应用教程课后答案matlab基础与应用教程课后答案【篇一：matlab教程基本应用练习题及解答】txt>要求：将每题的答案（命令行和运行结果、图片或m文件的文件名及具体内容）直接拷贝插入到各题的下方：（1）若为命令行，要求将提示符“”一起拷入，并在右侧用“％”注明命令行的每条命令的作用；（2）若为多个运行结果，拷入后要求解释每个结果具体对应题目中的哪个要求，也在右侧用“％”注明；（3）如果为m 文件，除了将文件名和此文件的具体内容全部拷入外，再将所有原始m文件和本练习题电子版放在一个以“专业班级＋本人姓名＋学号”命名的文件夹内一起上传。

ans =1 0 00 1 00 0 1ans =0 0 00 0 0ans =1 1 11 1 11 1 11 1 120+40*rand(1,10) %生成10个在区间[20,60]上均匀分布的随机数。

ans =44.617351.677556.872549.528327.050636.228257.418856.6 76236 .410855.7460二．1）计算向量（2，4，6，8）的最大值，最小值，平均值，中值，排序，总和值； max([2,4,6,8])%最大值ans =8min([2,4,6,8])% 最小值ans =2mean([2,4,6,8])% 平均值ans =5median([2,4,6,8])% 中值ans =5sort([2,4,6,8])% 排序ans =2 4 6 8sum([2,4,6,8])% 总和值ans =202）在行向量（2，4，6，8）和（1，3，5，7）之间实施加减乘除及幂运算； a=[2 4 6 8];b=[1 3 5 7];a+b%加ans =3 7 11 15a-b%减ans =a.*b%乘ans =2 12 30 56a./b%除ans =2.0000 1.3333 1.2000 1.1429a.^b%幂ans =2 64 7776 20971523）在向量（1，3，5）与标量2之间实施加减乘除及幂运算。

Matlab作业3(数值分析)答案

Matlab作业3（数值分析）机电工程学院（院、系）专业班组学号姓名实验日期教师评定1.计算多项式乘法(x2+2x+2)(x2+5x+4)。

答：2. （1）将(x-6)(x-3)(x-8)展开为系数多项式的形式。

（2）求解在x=8时多项式(x-1)(x-2) (x-3)(x-4)的值。

答：（1）（2）3. y=sin(x)，x从0到2π，∆x=0.02π，求y的最大值、最小值、均值和标准差。

4.设ｘ＝[0.00.30.8 1.1 1.6 2.3]'，y=[0.500.82 1.14 1.25 1.351.40]'，试求二次多项式拟合系数，并据此计算x1=[0.9 1.2]时对应的y1。

解：x=[0.0 0.3 0.8 1.1 1.62.3]'; %输入变量数据xy=[0.50 0.82 1.14 1.25 1.35 1.40]'; %输入变量数据yp=polyfit(x,y,2) %对x,y用二次多项式拟合,得到系数px1=[0.9 1.2]; %输入点x1y1=polyval(p,x1) %估计x1处对应的y1p =-0.2387 0.9191 0.5318y1 =a) 1.29095.实验数据处理：已知某压力传感器的测试数据如下表p为压力值，u为电压值，试用多项式dcpbpappu+++=23)(来拟合其特性函数，求出a,b,c,d，并把拟合曲线和各个测试数据点画在同一幅图上。

解：>> p=[0.0,1.1,2.1,2.8,4.2,5.0,6.1,6.9,8.1,9.0,9.9];u=[10,11,13,14,17,18,22,24,29,34,39];x=polyfit(p,u,3) %得多项式系数t=linspace(0,10,100);y=polyval(x,t); %求多项式得值plot(p,u,'*',t,y,'r') %画拟和曲线x =0.0195 -0.0412 1.4469 9.8267。

matlab基础作业参考答案.doc

Matlab基础作业参考答案一、选择或填空，每个4分。

1、标点符号一_可以使命令行不显示运算结果，―用来表示该行为注释行。

2、在MATLAB语言中变量的命名应遵循的规则有：(1) 变量名区分大小写(2) 变量名不超过31位，31位后被忽略________ (3) 变量名以字母开头，由下划线、字母、数字组成，不能用标点3、兀为°〜4兀，步长为°」兀的向量，使用命令__________ 创建。

4、在命令行后加_________ 表示该命令执行但不显示执行结果_1323 --575、输入矩阵_569,使用全下标方式用A(2,2)取出元素“ 一5”，使用单下标方式用A(5) 取出元素“一5”。

6、M atlab中常利用load函数读入外部数据，其调用方法为__load+文件名［参数］________7、M脚本文件和M函数文件的主要区别是M函数文件有函数声明行和____ O8、在循环结构中跳出循环，但继续下次循环的命令为 c °(A) return; (B) break; (C) continue; (D) keyboad二、按要求在Matlab中作图，并将图形粘贴到试卷中。

1、(8分)在同一个坐标系屮画出yl=cosx和y2=0.5sin2x的图像，要求(1) x的范围为0.5Pi到4pi (2分)(2) yl用红色虚线，y2用黑色波折线(3分)(3) 在图中指明曲线对应的函数，并加上网格(3分)程序：clearclcx=0 ・5*pi:pi/100:4*pi;yl=cos (x);y2=0.5*sin(2*x);plot(x z yl, 'r-. ',x f y2f *k-- *)gridgtext(1cosx1)gtext(10.5sin2x 1)2] ) ,title(1 sin(x) !); 2]),title('sin (2x) *);2])z title(1 sin (3x) 1);-2 2])八itle(「sin (5x) 1);2、（10分）在一个画面中分别画出yl 二sinx 、y2=sin2x> y3=sin3x> y4=sin5x 的图像要求（1）每个函数图像占用不同的坐标系（4分）（2）控制x 范围为［到5, y 的范围为-2到2,且x 轴和y 轴的单位长度相同（3分）（3）给每个图像对应的函数作为标题加上（3分）程序：clear clc x=0:0.1:15; yl = sin(x); y2 = sin(2*x); y3 = sin(3*x); y4 = sin(5*x);subplot(2,2,1),plot(x,yl), axis([1 subplot(2,2,2)9 plot(x z y2)9 axis([1subplot(2,2r 3)9 plot(x z y3) , axis([1 subplot (2f 2r 4)r plot(x,y4)faxis([1 5Illa5Vorb osnrb orbsOHB •nrb ocn1,^x 啟K0 sinx'xs^s專syms x-imit(sin(x)、xxo)总遅一2,^f(xll5x>3+10x>2,33x+2134[n^®^“syms x7r5*x>3+10*x>2,33*x+21diff(f)fis汩“ 15*x>2 + 20*x — 33程序：syms x a bf=a*x A2-b*x+23solve(f,x)结果：(b + (b A2 - 92*a)A(l/2))/(2*a)(b・(bH・ 92*a)A(l/2))/(2*a)四、编程题，每个15分。

MATLAB习题及参考答案

2 9 0 1310，解方程组3 4 11 x 6。

(应用x=a\b)2 2 6 611,求欠定方程组294 73 54x68的最小范数解。

（应用pinv）5习题:6 9 3 2 4 1与b 的数组乘积。

2 7 5 4 6 81,计算a2, 对于AX B，如果A 3726，求解X。

283,已知：a 1 2 34 5 6，分别计算a的数组平方和矩阵平方，并观察其结果。

7 8 94，角度x 30 45 60，求x的正弦、余弦、正切和余切。

（应用sin,cos,4 2 7 15 95，将矩阵a 、b 和c 组合成两个新矩阵：5 7 8 36 2（1）组合成一个4 3的矩阵，第一列为按列顺序排列的a矩阵元素，第二列为按列顺序排列的元素，第三列为按列顺序排列的c矩阵元素，即b矩阵4 7 55 8 62 1 97 3 2（2）按照a、b、c的列顺序组合成一个行矢量，即4527781356926，将（x-6）（x-3）（x-8）展开为系数多项式的形式。

（应用poly,polyvalm）7，求解多项式X3-7X2+2X+40的根。

（应用roots）8，求解在x=8 时多项式(x-1)( x-2) ( x-3)( x-4)的值。

(应用poly,polyvalm)9, 计算多项式4x412x314x25x 9的微分和积分。

(应用polyder,polyint ，poly2sym)2 2计算表达式z 10 x 3 y 5 e x y 的梯度并绘图。

（应用meshgrid, gradient, con tour, holdon, quiver)15，用符号函数法求解方程a t 2+b*t +c=0。

（应用solve ）16，用符号计算验证三角等式：（应用syms,simple ）用 syms,ezplot)用 plot,title,text,legend)24, x= [66 49 71 56 38] ，绘制饼图，并将第五个切块分离出来。

matlab课后习题答案

第2章 MATLAB 矩阵运算基础2.1 在MA TLAB 中如何建立矩阵⎥⎦⎤⎢⎣⎡194375，并将其赋予变量a ？ >> a=[5 7 3;4 9 1]2.5 计算矩阵⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡897473535与⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡638976242之和。

>> a=[5 3 5;3 7 4;7 9 8];>> b=[2 4 2;6 7 9;8 3 6]; >> a+b ans =7 7 7 9 14 13 15 12 142.6 求⎥⎦⎤⎢⎣⎡+-+-+-+-++=i 44i 93i 49i 67i 23i 57i 41i 72i 53i 84x 的共轭转置。

>> x=[4+8i 3+5i 2-7i 1+4i 7-5i;3+2i 7-6i 9+4i 3-9i 4+4i];>> x’ ans =4.0000 - 8.0000i 3.0000 - 2.0000i 3.0000 -5.0000i 7.0000 +6.0000i 2.0000 +7.0000i 9.0000 - 4.0000i 1.0000 - 4.0000i 3.0000 + 9.0000i 7.0000 + 5.0000i 4.0000 - 4.0000i2.7 计算⎥⎦⎤⎢⎣⎡=572396a 与⎥⎦⎤⎢⎣⎡=864142b 的数组乘积。

>> a=[6 9 3;2 7 5];>> b=[2 4 1;4 6 8]; >> a.*b ans =12 36 3 8 42 402.9 对于B AX =，如果⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=753467294A ，⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=282637B ，求解X 。

>> A=[4 9 2;7 6 4;3 5 7];>> B=[37 26 28]’;-0.5118 4.0427 1.33182.10 已知：⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=987654321a ，分别计算a 的数组平方和矩阵平方，并观察其结果。

MATLAB 实用教程课后习题答案

第二章1.计算复数3+4i与5-6i的乘积。

a=3+4ib=5-6ic=a*b2.构建结构体Students，属性包含Name、age和Email，数据包括｛’Zhang’，18，[‘*************’，’*************’]｝、｛’Wang’，21，[]｝和｛’Li’，[],[]｝，构建后读取所有Name属性值，并且修改’Zhang’的Age属性值为19。

Students(1).Age=18Students(1).Email='*************','*************'Students(2).Name='Wang'Students(2).Age=21Students(2).Email=[]Students(3).Name='Li'Students(3).Age=[]Students(3).Email=[]Student(1).Age(1)=19Student.Age3.用满矩阵和稀疏矩阵存储方式分别构造下属矩阵：A=[0 1 0 0 0;1 0 0 0 0;0 0 0 0 0;0 0 0 1 0]A=[0 1 0 0 0;1 0 0 0 0;0 0 0 0 0;0 0 0 1 0]S=sparse(A)S=sparse([2,1,4],[1,2,4],[1,1,1],4,5)4.采用向量构造符得到向量[1,5,9....,41].A=1:4:415.按水平和竖直方向分别合并下述两个矩阵：A=[1 0 0;1 1 0;0 0 1],B=[2 3 4;5 6 7;8 9 10]A=[1 0 0;1 1 0;0 0 1]B=[2 3 4;5 6 7;8 9 10]C=[A B]D=[A;B]6.分别删除第五题两个结果的第2行。

A=[1 0 0;1 1 0;0 0 1]B=[2 3 4;5 6 7;8 9 10]C=[A B]D=[A;B]C(2,:)=[]D(2,:)=[]7.分别将第5题两个结果的第2行最后3列的数值改为[11 12 13]。

matlab课后习题答案刘

matlab课后习题答案刘MATLAB课后习题答案刘在本学期的MATLAB课程中，我们学习了MATLAB的基础知识和一些高级功能。

为了帮助同学们更好地复习和巩固所学知识，以下是针对课后习题的一些参考答案。

1. 向量和矩阵的基本操作- 创建一个3x3的单位矩阵。

```matlabA = eye(3);```- 计算矩阵A和B的乘积，其中B是A的转置。

```matlabB = A';C = A * B;```2. 条件语句和循环- 编写一个循环，打印出1到10的所有奇数。

```matlabfor i = 1:2:10disp(i);end```- 编写一个条件语句，检查一个数是否为素数。

```matlabfunction isPrime(n)if n <= 1isPrime = false;elsefor i = 2:sqrt(n)if mod(n, i) == 0isPrime = false;return;endendisPrime = true;endend```3. 函数和脚本- 编写一个函数，计算给定半径的圆的面积。

```matlabfunction area = circleArea(radius)area = pi * radius^2;end```- 编写一个脚本来调用上述函数，并打印出半径为5的圆的面积。

```matlabradius = 5;disp(['Area of circle with radius ', num2str(radius), ' is ', num2str(circleArea(radius))]);```4. 数据结构和文件I/O- 读取一个文本文件中的数据到一个cell数组。

```matlabfilename = 'data.txt';fileID = fopen(filename, 'r');data = textscan(fileID, '%s');fclose(fileID);```- 将一个矩阵写入到一个CSV文件中。

MATLAB习题答案(清华大学)

高等应用数学问题MATLAB 求解习题参考解答（薛定宇著）目录第1 章计算机数学语言概述2第2 章MATLAB 语言程序设计基础5第3 章微积分问题的计算机求解17第4 章线性代数问题的计算机求解29第5 章积分变换与复变函数问题的计算机求解43第6 章代数方程与最优化问题的计算机求解53第7 章微分方程问题的计算机求解71第8 章数据插值、函数逼近问题的计算机求解93第9 章概率论与数理统计问题的计算机求解114第10 章数学问题的非传统解法127第A章自由数学语言Scilab 简介136第1 章计算机数学语言概述1 在你的机器上安装MATLAB 语言环境，并键入demo 命令，由给出的菜单系统和对话框原型演示程序，领略MATLAB 语言在求解数学问题方面的能力与方法。

【求解】在MATLAB 提示符>> 下键入demo 命令，则将打开如图1-1 所示的窗口，窗口左侧的列表框可以选择各种不同组合的演示内容。

图1-1 MATLAB 演示程序界面例如，用户选择MATLAB ! Graphics ! Volume Vlsulization 演示，则将得出如图1-2 所示的演示说明，单击其中的Run this demo 栏目，则将得出如图1-3 所示的演示界面。

用户可以在该界面下按按钮，逐步演示相关内容，而实现这样演示的语句将在该程序界面的下部窗口中给出。

2 作者用MATLAB 语言编写了给出例子的源程序，读者可以自己用type 语句阅读一下源程序，对照数学问题初步理解语句的含义，编写的源程序说明由下表列出。

第1 章计算机数学语言概述3图1-2 MATLAB 演示程序界面举例序号文件名程序说明例1.1 c1ex1.m 利用MATLAB 的符号运算工具箱求解微分问题例1.2 c1ex2.m 分别利用MATLAB 的符号运算工具箱和数值运算功能求解多项式方程，其中用数值方法得出的结果有误差例1.3 c1ex3.m 分别利用MATLAB 的符号运算工具箱和数值运算功能计算Hilbert 矩阵的行列式，其中用数值方法得出的结果有很大误差例1.4 c1ex4.m 令x1 = y; x2 = y_，则可以将原来的二阶微分方程转换成一阶微分方程组，然后就可以求解微分方程的数值解了，原方程是非线性微分方程，故不存在解析解。

matlab练习题和答案

matlab 练习题和答案控制系统仿真切验Matlab部分实验结果目录实验一MATLAB基本操作..................................................................... ....................... 1实验二Matlab编程..................................................................... ............................... 5实验三Matlab基层图形控制..................................................................... ..................... 6实验四控制系统古典分析 ...................................................................... .......................12实验五控制系统现代分析..................................................................... . (15)实验六PID控制器的设计..................................................................... ......................19实验七系统状态空间设计 ...................................................................... .......................23实验九直流双闭环调速系统仿真..................................................................... . (25)实验一MATLAB基本操作1用 MATLAB能够识其他格式输入下边两个矩阵1233，， 1443678， i ，， ,,2357,,,,2335542，i,,,, A,1357B,,,2675342，i,,3239,,,,189543，， ,,1894 ，，再求出它们的乘积矩阵C，并将 C矩阵的右下角 2×3子矩阵赋给 D矩阵。