MATLAB习题集
MATLAB编程练习(含答案很好的)
001双峰曲线图：z=peaks(40);mesh(z);surf(z)002解方程：A=[3,4,-2;6,2,-3;45,5,4];>> B=[14;4;23];>> root=inv(A)*B003傅里叶变换load mtlb ;subplot(2,1,1);plot(mtlb);>> title('原始语音信息');>> y=fft(mtlb);>> subplot(2,1,2);>> yy=abs(y);>> plot(yy);>> title('傅里叶变换')004输入函数：a=input('How many apples\n','s')005输出函数a=[1 2 3 4 ;5 6 7 8;12 23 34 45;34 435 23 34]a =1 2 3 45 6 7 812 23 34 4534 435 23 34disp(a)a =1 2 3 45 6 7 812 23 34 4534 435 23 34b=input('how many people\n' ,'s')how many peopletwo peopleb =two people>> disp(b)two people>>006求一元二次方程的根a=1;b=2;c=3;d=sqrt(b^2-4*a*c);x1=(-b+d)/(2*a)x1 =-1.0000 + 1.4142i>> x2=(-b-d)/(2*a)x2 =-1.0000 - 1.4142i007求矩阵的相乘、转置、存盘、读入数据A=[1 3 5 ;2 4 6;-1 0 -2;-3 0 0];>> B=[-1 3;-2 2;2 1];>> C=A*BC =3 142 20-3 -53 -9>> C=C'C =3 2 -3 314 20 -5 -9>> save mydat C>> clear>> load mydat C008编写数学计算公式：A=2.1;B=-4.5;C=6;D=3.5;E=-5;K=atan(((2*pi*A)+E/(2*pi*B*C))/D) K =1.3121009A=[1 0 -1;2 4 1;-2 0 5];>> B=[0 -1 0;2 1 3;1 1 2];>> H=2*A+BH =2 -1 -26 9 5-3 1 12>> M=A^2-3*BM =3 3 -62 13 -2-15 -3 21>> Y=A*BY =-1 -2 -29 3 145 7 10>> R=B*AR =-2 -4 -1-2 4 14-1 4 10>> E=A.*BE =0 0 04 4 3-2 0 10>> W=A\BW =0.3333 -1.3333 0.66670.2500 1.0000 0.25000.3333 -0.3333 0.6667 >> P=A/BP =-2.0000 3.0000 -5.0000-5.0000 3.0000 -4.00007.0000 -9.0000 16.0000>> Z=A.\BWarning: Divide by zero.Z =0 -Inf 01.0000 0.2500 3.0000-0.5000 Inf 0.4000>> D=A./BWarning: Divide by zero.D =Inf 0 -Inf1.0000 4.0000 0.3333-2.0000 0 2.5000010a=4.96;b=8.11;>> M=exp(a+b)/log10(a+b)M =4.2507e+005011求三角形面积：a=9.6;b=13.7;c=19.4;>> s=(a+b+c)/2;>> area=sqrt(s*(s-a)*(s-b)*(s-c))area =61.1739012逻辑运算A=[-1 0 -6 8;-9 4 0 12.3;0 0 -5.1 -2;0 -23 0 -7]; >> B=A(:,1:2)B =-1 0-9 40 00 -23>> C=A(1:2,:)C =-1.0000 0 -6.0000 8.0000 -9.0000 4.0000 0 12.3000>> D=B'D =-1 -9 0 00 4 0 -23>> A*Bans =1.0000 -184.0000-27.0000 -266.90000 46.0000 207.0000 69.0000>> C<Dans =0 0 1 01 0 0 0>> C&Dans =1 0 0 00 1 0 1>> C|Dans =1 1 1 11 1 0 1>> ~C|~Dans =0 1 1 11 0 1 0013矩阵运算练习：A=[8 9 5;36 -7 11;21 -8 5]A =8 9 536 -7 1121 -8 5>> BB =-1 3 -22 0 3-3 1 9>> RT=A*BRT =-5 29 56-83 119 6-52 68 -21>> QW=A.*BQW =-8 27 -1072 0 33-63 -8 45>> ER=A^3ER =6272 3342 294415714 -856 52608142 -1906 2390 >> BF=A.^3BF =512 729 12546656 -343 13319261 -512 125 >> A/Bans =3.13414.9634 -0.4024-1.2561 12.5244 -3.2317-1.9878 6.4512 -2.0366>> EKV=B\AEKV =10.7195 -1.2683 3.52449.4756 1.5854 3.71954.8537 -1.4878 1.3171>> KDK=[A,B]KDK =8 9 5 -1 3 -236 -7 11 2 0 321 -8 5 -3 1 9 >> ERI=[A;B]ERI =8 9 536 -7 1121 -8 5-1 3 -22 0 3-3 1 9014一般函数的调用：A=[2 34 88 390 848 939];>> S=sum(A)S =2301>> min(A)ans =2>> EE=mean(A)EE =383.5000>> QQ=std(A)QQ =419.3794>> AO=sort(A)AO =2 34 88 390 848 939 >> yr=norm(A)yr =1.3273e+003>> RT=prod(A)RT =1.8583e+012>> gradient(A)ans =32.0000 43.0000 178.0000 380.0000 274.5000 91.0000 >> max(A)ans =939>> median(A)ans =239>> diff(A)ans =32 54 302 458 91>> length(A)ans =6>> sum(A)ans =2301>> cov(A)ans =1.7588e+005>>015矩阵变换：A=[34 44 23;8 34 23;34 55 2]A =34 44 238 34 2334 55 2>> tril(A)ans =34 0 08 34 034 55 2>> triu(A)ans =34 44 230 34 230 0 2>> diag(A)ans =34342norm(A)ans =94.5106>> rank(A)ans =3>> det(A)ans =-23462>> trace(A)ans =70>> null(A)ans =Empty matrix: 3-by-0>> eig(A)ans =80.158712.7671-22.9257>> poly(A)ans =1.0e+004 *0.0001 -0.0070 -0.1107 2.3462>> logm(A)Warning: Principal matrix logarithm is not defined for A with nonpositive real eigenvalues. A non-principal matrixlogarithm is returned.> In funm at 153In logm at 27ans =3.1909 + 0.1314i 1.2707 + 0.1437i 0.5011 - 0.2538i0.4648 + 0.4974i 3.3955 + 0.5438i 0.1504 - 0.9608i0.2935 - 1.2769i 0.8069 - 1.3960i 3.4768 + 2.4663i>> fumn(A)Undefined command/function 'fumn'.>> inv(A)ans =0.0510 -0.0502 -0.0098-0.0326 0.0304 0.02550.0305 0.0159 -0.0343>> cond(A)ans =8.5072>> chol(A)Error using ==> cholMatrix must be positive definite.>> lu(A)ans =34.0000 44.0000 23.00000.2353 23.6471 17.58821.0000 0.4652 -29.1816>> pinv(A)ans =0.0510 -0.0502 -0.0098-0.0326 0.0304 0.02550.0305 0.0159 -0.0343>> svd(A)ans =94.510622.345611.1095>> expm(A)ans =1.0e+034 *2.1897 4.3968 1.93821.31542.6412 1.16431.8782 3.7712 1.6625>> sqrtm(A)ans =5.2379 + 0.2003i 3.4795 + 0.2190i 1.8946 - 0.3869i0.5241 + 0.7581i 5.1429 + 0.8288i 2.0575 - 1.4644i3.0084 - 1.9461i4.7123 - 2.1276i 2.1454 + 3.7589i >>016多项式的计算：A=[34 44 23;8 34 23;34 55 2]A =34 44 238 34 2334 55 2>> P=poly(A)P =1.0e+004 *0.0001 -0.0070 -0.1107 2.3462>> PPA=poly2str(P,'X')PPA =X^3 - 70 X^2 - 1107 X + 23462017多项式的运算：p=[2 6 8 3];w=[32 56 0 2];>> m=conv(p,w)m =64 304 592 548 180 16 6 >> [q,r]=deconv(w,p)q =16r =0 -40 -128 -46>> dp=polyder(w)dp =96 112 0>> [num,den]=polyder(w,p)num =80 512 724 312 -16den =4 24 68 108 100 48 9>> b=polyfit(p,w,4)Warning: Polynomial is not unique; degree >= number of data points. > In polyfit at 74b =-0.6704 9.2037 -32.2593 0 98.1333>> r=roots(p)r =-1.2119 + 1.0652i-1.2119 - 1.0652i-0.5761018求多项式的商和余p=conv([1 0 2],conv([1 4],[1 1]))p =1 5 6 10 8>> q=[1 0 1 1]q =1 0 1 1>> [w,m]=deconv(p,q)w =1 5m =0 0 5 4 3>> cq=w;cr=m;>> disp([cr,poly2str(m,'x')])5 x^2 + 4 x + 3>> disp([cq,poly2str(w,'x')])x + 5019将分式分解a=[1 5 6];b=[1];>> [r,p,k]=residue(b,a)r =-1.00001.0000p =-3.0000-2.0000k =[]020计算多项式：a=[1 2 3;4 5 6;7 8 9];>> p=[3 0 2 3];>> q=[2 3];>> x=2;>> r=roots(p)r =0.3911 + 1.0609i0.3911 - 1.0609i-0.7822>> p1=conv(p,q)p1 =6 9 4 12 9>> p2=poly(a)p2 =1.0000 -15.0000 -18.0000 -0.0000 >> p3=polyder(p)p3 =9 0 2>> p4=polyval(p,x)p4 =31021求除式和余项：[q,r]=deconv(conv([1 0 2],[1 4]),[1 1 1])022字符串的书写格式：s='student's =student>> name='mary';>> s1=[name s]s1 =marystudent>> s3=[name blanks(3);s]s3 =marystudent>>023交换两个数：clearclca=[1 2 3 4 5];b=[6 7 8 9 10];c=a;a=b;b=c;ab24If语句n=input('enter a number,n=');if n<10nend025 if 双分支结构a=input('enter a number ,a=');b=input('enter a number ,b=');if a>bmax=a;elsemax=b;endmax026三个数按照由大到小的顺序排列：A=15;B=24;C=45;if A<BT=A;A=B;B=T;elseif A<CT=A;A=C;C=T;elseif B<CT=B;B=C;C=T;endABC027建立一个收费优惠系统：price=input('please jinput the price : price=') switch fix(price/100)case[0,1]rate =0;case[2,3,4]rate =3/100;case num2cell(5:9)rate=5/100;case num2cell(10:24)rate=8/100;case num2cell(25:49)rate=10/100;otherwiserate=14/100;endprice=price*(1-rate)028：while循环语句i=0;s=0;while i<=1212s=s+i;i=i+1;ends029，用for循环体语句：sum=0;for i=1:1.5:100;sum=sum+i;endsum030循环的嵌套s=0;for i=1:1:6;for j=1:1:8;s=s+i^j;end;end;s031continue 语句的使用：for i=100:120;if rem(i,7)~=0;continue;end;iend032x=input ('输入X的值x=')if x<1y=x^2;elseif x>1&x<2y=x^2-1;elsey=x^2-2*x+1;endy033求阶乘的累加和sum=0;temp=1;for n=1:10;temp=temp*n;sum=sum+temp;endsum034对角线元素之和sum=0;a=[1 2 3 4;5 6 7 8;9 10 11 12;13 14 15 16]; for i=1:4;sum=sum+a(i,i);endsum035用拟合点绘图A=[12 15.3 16 18 25];B=[50 80 118 125 150.8];plot(A,B)036绘制正玄曲线：x=0:0.05:4*pi;y=sin(x);plot(x,y)037绘制向量x=[1 2 3 4 5 6;7 8 9 10 11 12;13 14 15 16 17 18] plot(x)x=[0 0.2 0.5 0.7 0.6 0.7 1.2 1.5 1.6 1.9 2.3]plot(x)x=0:0.2:2*piy=sin(x)plot(x,y,'m:p')038在正弦函数上加标注：t=0:0.05:2*pi;plot(t,sin(t))set(gca,'xtick',[0 1.4 3.14 56.28])xlabel('t(deg)')ylabel('magnitude(v)')title('this is a example ()\rightarrow 2\pi')text(3.14,sin(3.14),'\leftarrow this zero for\pi')039添加线条标注x=0:0.2:12;plot(x,sin(x),'-',x,1.5*cos(x),':');legend('First','Second',1)040使用hold on 函数x=0:0.2:12;plot(x,sin(x),'-');hold onplot(x,1.5*cos(x),':');041一界面多幅图x=0:0.05:7;y1=sin(x);y2=1.5*cos(x);y3=sin(2*x);y4=5*cos(2*x);subplot(221);plot(x,y1);title('sin(x)')subplot(222);plot(x,y2);title('cos(x)')subplot(223);plot(x,y3);title('sin(2x)')subplot(224);plot(x,y4);title('cos(2x)')042染色效果图x=0:0.05:7;y1=sin(x);y2=1.5*cos(x);y3=sin(2*x);y4=5*cos(2*x);subplot(221);plot(x,y1);title('sin(x)');fill(x,y1,'r') subplot(222);plot(x,y2);title('cos(x)');fill(x,y2,'b') subplot(223);plot(x,y3);title('sin(2x)');fill(x,y3,'k') subplot(224);plot(x,y4);title('cos(2x)');fill(x,y4,'g')043特殊坐标图clcy=[0,0.55,2.5,6.1,8.5,12.1,14.6,17,20,22,22.1] subplot(221);plot(y);title('线性坐标图');subplot(222);semilogx(y);title('x轴对数坐标图');subplot(223);semilogx(y);title('y轴对数坐标图');subplot(224);loglog(y);title('双对数坐标图')t=0:0.01:2*pi;r=2*cos(2*(t-pi/8));polar(t,r)044特殊函数绘图：fplot('cos(tan(pi*x))',[-0.4,1.4])fplot('sin(exp(pi*x))',[-0.4,1.4])045饼形图与条形图：x=[8 20 36 24 12];subplot(221);pie(x,[1 0 0 0 1]);title('饼图');subplot(222);bar(x,'group');title('垂直条形图');subplot(223);bar(x,'stack');title('累加值为纵坐标的垂直条形图'); subplot(224);barh(x,'group');title('水平条形图');046梯形图与正弦函数x=0:0.1:10;y=sin(x);subplot(121);stairs(x);subplot(122);stairs(x,y);047概率图x=randn(1,1000);y=-2:0.1:2;hist(x,y)048向量图：x=[-2+3j,3+4j,1-7j];subplot(121);compass(x);rea=[-2 3 1];imag=[3 4 -7];subplot(122);feather(rea,imag);049绘制三维曲线图：z=0:pi/50:10*pi;x=sin(z);y=cos(z);plot3(x,y,z)x=-10:0.5:10;y=-8:0.5:8;[x,y]=meshgrid(x,y);z=sin(sqrt(x.^2+y.^2))./sqrt(x.^2+y.^2); subplot(221);mesh(x,y,z);title('普通一维网格曲面');subplot(222);meshc(x,y,z);title('带等高线的三维网格曲面'); subplot(223);meshz(x,y,z);title('带底座的三维网格曲面'); subplot(224);surf(x,y,z);title('充填颜色的三维网格面')050 带网格二维图x=0:pi/10:2*pi;y1=sin(x);y2=cos(x);plot(x,y1,'r+-',x,y2,'k*:')grid onxlabel('Independent Variable x') ylabel('Dependent Variable y1&y2') text(1.5,0.5,'cos(x)')051各种统计图y=[18 5 28 17;24 12 36 14;15 6 30 9]; subplot(221);bar(y)x=[4,6,8];subplot(222);bar3(x,y)subplot(223);bar(x,y,'grouped') subplot(224);bar(x,y,'stack')052曲面图x=-2:0.4:2;y=-1:0.2:1;[x,y]=meshgrid(x,y);z=sqrt(4-x.^2/9-y.^2/4); surf(x,y,z)grid on053创建符号矩阵e=[1 3 5;2 4 6;7 9 11];m=sym(e)符号表达式的计算问题因式分解：syms xf=factor(x^3-1)s=sym('sin(a+b)'); expand(s)syms x tf=x*(x*(x-8)+6)*t; collect(f)syms xf=sin(x)^2+cos(x)^2; simplify(f)syms xs=(4*x^2+8*x+3)/(2*x+1); simplify(s)通分syms x yf=x/y-y/x;[m,n]=numden(f)嵌套重写syms xf=x^4+3*x^3-7*x^2+12; horner(f)054求极限syms x a;limit(exp(-x),x,0,'left')求导数syms xdiff(x^9+x^6)diff(x^9+x^6,4)055求不定积分与定积分syms x ys=(4-3*x^2)^2;int(s)int(x/(x+y),x)int(x^2/(x+2),x,1,3) double(ans)056函数的变换：syms x ty=exp(-x^2);Ft=fourier(y,x,t)fx=ifourier(Ft,t,x)057求解方程syms a b c xs=a*x^2+b*x+c;solve(s)syms x y zs1=2*x^2+y^2-3*z-4;s2=y+z-3;s3=x-2*y-3*z;[x,y,z]=solve(s1,s2,s3)058求微分方程：y=dsolve('Dy-(t^2+y^2)/t^2/2','t')059求级数和syms x ksymsum(k)symsum(k^2-3,0,10)symsum(x^k/k,k,1,inf)060泰勒展开式syms xs=(1-x+x^2)/(1+x+x^2);taylor(s)taylor(s,9)taylor(s,x,12)taylor(s,x,12,5)061练习syms x a;s1=sin(2*x)/sin(5*x);limit(s1,x,0)s2=(1+1/x)^(2*x);limit(s2,x,inf)syms xs=x*cos(x);diff(s)diff(s,2)diff(s,12)syms xs1=x^4/(1+x^2);int(s1)s2=3*x^2-x+1int(s2,0,2)syms x y zs1=5*x+6*y+7*z-16;s2=4*x-5*y+z-7;s3=x+y+2*z-2;[x,y,z]=solve(s1,s2,s3)syms x yy=dsolve('Dy=exp(2*x-y)','x')y=dsolve('Dy=exp(2*x-y)','y(0)=0','x')n=sym('n');s=symsum(1/n^2,n,1,inf)x=sym('x');f=sqrt(1-2*x+x^3)-(1-3*x+x^2)^(1/3);taylor(f,6)062求于矩阵相关的值a=[2 2 -1 1;4 3 -1 2;8 5 -3 4;3 3 -2 2]adet=det(a)atrace=trace(a)anorm=norm(a)acond=cond(a)arank=rank(a)eiga=eig(a)063矩阵计算A=[0.1389 0.6038 0.0153 0.9318;0.2028 0.2772 0.7468 0.4660;0.1987 0.1988 0.4451 0.4186]B=var(A)C=std(A)D=range(A)E=cov(A)F=corrcoef(A)064求根及求代数式的值P=[4 -3 2 5];x=roots(P)x=[3 3.6];F=polyval(P,x)065多项式的和差积商运算：f=[1 2 -4 3 -1]g=[1 0 1]g1=[0 0 1 0 1]f+g1f-g1conv(f,g)[q,r]=deconv(f,g)polyder(f)066各种插值运算：X=0:0.1:pi/2;Y=sin(X);interp1(X,Y,pi/4)interp1(X,Y,pi/4,'nearest')interp1(X,Y,pi/4,'spline')interp1(X,Y,pi/4,'cubic')067曲线的拟合：X=0:0.1:2*pi;Y=cos(X);[p,s]=polyfit(X,Y,4)plot(X,Y,'K*',X,polyval(p,X),'r-')068求函数的最值与0点x=2:0.1:2;[x,y]=fminbnd('x.^3-2*x+1',-1,1) [x,y]=fzero('x.^3-2*x+1',1)069求多项式的表达式、值、及图像y=[1 3 5 7 19]t=poly(y)x=-4:0.5:8yx=polyval(t,x)plot(x,yx)070数据的拟合与绘图x=0:0.1:2*pi;y=sin(x);p=polyfit(x,y,5);y1=polyval(p,x)plot(x,y,'b',x,y1,'r')071求代数式的极限：syms xf=sym('log(1+2*x)/sin(3*x)');b=limit(f,x,0)072求导数与微分syms xf=sym('x/(cos(x))^2');y1=diff(f)y2=int(f,0,1)078划分网格函数[x,y]=meshgrid(-2:0.01:2,-3:0.01:5); t=x.*exp(-x.^2-y.^2);[px,py]=gradient(t,0.05,0.1);td=sqrt(px.^2+py.^2);subplot(221)imagesc(t)subplot(222)imagesc(td)colormap('gray')079求多次多项方程组的解：syms x1 x2 a ;eq1=sym('x1^2+x2=a')eq2=sym('x1-a*x2=0')[x1 x2]=solve(eq1,eq2,x1,x2)v=solve(eq1,eq2)v.x1v.x2an1=x1(1),an2=x1(2)an3=x2(1),an4=x2(2)080求解微分方程：[y]=dsolve('Dy=-y^2+6*y','y(0)=1','x')s=dsolve('Dy=-y^2+6*y','y(0)=1','x')[u]=dsolve('Du=-u^2+6*u','u(0)=1')w=dsolve('Du=-u^2+6*u','z')[u,w]=dsolve('Du=-w^2+6*w,Dw=sin(z)','u(0)=1,w(0)=0','z') v=dsolve('Du=-w^2+6*w,Dw=sin(z)','u(0)=1,w(0)=0','z')081各种显现隐含函数绘图：f=sym('x^2+1')subplot(221)ezplot(f,[-2,2])subplot(222)ezplot('y^2-x^6-1',[-2,2],[0,10])x=sym('cos(t)')y=sym('sin(t)')subplot(223)ezplot(x,y)z=sym('t^2')subplot(224)ezplot3(x,y,z,[0,8*pi])082极坐标图：r=sym('4*sin(3*x)')ezpolar(r,[0,6*pi])083多函数在一个坐标系内：x=0:0.1:8;y1=sin(x);subplot(221)plot(x,y1)subplot(222)plot(x,y1,x,y2)w=[2 3;3 1;4 6]subplot(223)plot(w)q=[4 6:3 5:1 2]subplot(224)plot(w,q)084调整刻度图像：x=0:0.1:10;y1=sin(x);y2=exp(x);y3=exp(x).*sin(x);subplot(221)plot(x,y2)subplot(222)loglog(x,y2)subplot(223)plotyy(x,y1,x,y2)085等高线等图形，三维图：t=0:pi/50:10*pi;subplot(2,3,1)plot3(t.*sin(t),t.*cos(t),t.^2) grid on[x,y]=meshgrid([-2:0.1:2])z=x.*exp(-x.^2-y.^2)subplot(2,3,2)plot3(x,y,z)box offsubplot(2,3,3)meshz(x,y,z)subplot(2,3,4)surf(x,y,z)contour(x,y,z)subplot(2,3,6)surf(x,y,z)subplot(2,3,5)contour(x,y,z)box offsubplot(2,3,6)contour3(x,y,z)axis off086统计图Y=[5 2 1;8 7 3;9 8 6;5 5 5;4 3 2]subplot(221)bar(Y)box offsubplot(222)bar3(Y)subplot(223)barh(Y)subplot(224)bar3h(Y)087面积图Y=[5 1 2;8 3 7;9 6 8;5 5 5;4 2 3];subplot(221)area(Y)grid onset(gca,'Layer','top','XTick',1:5)sales=[51.6 82.4 90.8 59.1 47.0];x=90:94;profits=[19.3 34.2 61.4 50.5 29.4];subplot(222)area(x,sales,'facecolor',[0.5 0.9 0.6], 'edgecolor','b','linewidth',2) hold onarea(x,profits,'facecolor',[0.9 0.85 0.7], 'edgecolor','y','linewidth',2) hold offset(gca,'Xtick',[90:94])set(gca,'layer','top')gtext('\leftarrow 销售量') gtext('利润')gtext('费用')xlabel('年','fontsize',14)088函数的插值：x=0:2*pi;y=sin(x);xi=0:0.1:8;yi1=interp1(x,y,xi,'linear')yi2=interp1(x,y,xi,'nearest') yi3=interp1(x,y,xi,'spline')yi4=interp1(x,y,xi,'cublic')p=polyfit(x,y,3)yy=polyval(p,xi)subplot(3,2,1)plot(x,y,'o')subplot(3,2,2)plot(x,y,'o',xi,yy)subplot(3,2,3)plot(x,y,'o',xi,yi1)subplot(3,2,4)plot(x,y,'o',xi,yi2)subplot(3,2,5)plot(x,y,'o',xi,yi3)subplot(3,2,6)plot(x,y,'o',xi,yi4)089二维插值计算：[x,y]=meshgrid(-3:0.5:3);z=peaks(x,y);[xi,yi]=meshgrid(-3:0.1:3); zi=interp2(x,y,z,xi,yi,'spline') plot3(x,y,z)hold onmesh(xi,yi,zi+15)hold offaxis tight090函数表达式;function f=exlin(x)if x<0f=-1;elseif x<1f=x;elseif x<2f=2-x;elsef=0;end091：硬循环语句:n=5;for i=1:nfor j=1:nif i==ja(i,j)=2;elsea(i,j)=0;endendendwhile 循环语句：n=1;while prod(1:n)<99^99;n=n+1endn:092 switch开关语句a=input('a=?')switch acase 1disp('It is raning') case 0disp('It do not know')case -1disp('It is not ranging')otherwisedisp('It is raning ?')end093画曲面函数：x1=linspace(-3,3,30)y1=linspace(-3,13,34)[x,y]=meshgrid(x1,y1);z=x.^4+3*x.^2-2*x+6-2*y.*x.^2+y.^2-2*y; surf(x,y,z)。
MATLAB课后习题
MATLAB课后习题第⼀部分 MATLAB 运算基础1. 先求下列表达式的值，然后显⽰MATLAB ⼯作空间的使⽤情况并保存全部变量。
(1) 0122sin 851z e=+ (2) 221ln(1)2z x x =++，其中2120.455i x +??=?- (3) 0.30.330.3sin(0.3)ln , 3.0, 2.9,,2.9,3.022a a e e az a a --+=++=--(4) 2242011122123t t z t t t t t ?≤=-≤，其中t =0:0.5:2.52. 已知：1234413134787,2033657327A B --==-求下列表达式的值：(1) A+6*B 和A-B+I （其中I 为单位矩阵） (2) A*B 和A.*B (3) A^3和A.^3 (4) A/B 及B\A (5) [A,B]和[A([1,3],:);B^2]3. 设有矩阵A 和B123453166789101769,111213141502341617181920970212223242541311A B-???(1) 求它们的乘积C 。
(2) 将矩阵C 的右下⾓3×2⼦矩阵赋给D 。
(3) 查看MATLAB ⼯作空间的使⽤情况。
4. 完成下列操作：(1) 求[100,999]之间能被21整除的数的个数。
(2) 建⽴⼀个字符串向量，删除其中的⼤写字母。
第⼆部分 MATLAB 矩阵分析与处理1. 设有分块矩阵33322322E R A O S=?，其中E 、R 、O 、S 分别为单位矩阵、随机矩阵、零矩阵和对⾓阵，试通过数值计算验证22E R RS A OS +??=?。
2. 产⽣5阶希尔伯特矩阵H 和5阶帕斯卡矩阵P ，且求其⾏列式的值Hh 和Hp 以及它们的条件数Th 和Tp ，判断哪个矩阵性能更好。
为什么？ 3. 建⽴⼀个5×5矩阵，求它的⾏列式值、迹、秩和范数。
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习题及答案2. ⽤MATLAB 语句输⼊矩阵A 和B3．假设已知矩阵A ，试给出相应的MATLAB 命令，将其全部偶数⾏提取出来，赋给B 矩阵，⽤magic(8)A =命令⽣成A 矩阵，⽤上述命令检验⼀下结果是不是正确。
4．⽤数值⽅法可以求出∑=++++++==6363622284212i i S ，试不采⽤循环的形式求出和式的数值解。
由于数值⽅法是采⽤double 形式进⾏计算的，难以保证有效位数字，所以结果不⼀定精确。
试采⽤运算的⽅法求该和式的精确值。
5．选择合适的步距绘制出下⾯的图形。
（1）)/1sin(t ，其中)1,1(-∈t ；（2）)tan(sin )sin(tan t t -，其中),(ππ-∈t6. 试绘制出⼆元函数2222)1(1)1(1),(yx yx y x f z ++++-==的三维图和三视图7. 试求出如下极限。
（1）xxxx 1)93(lim +∞→；（2）11lim0-+→→xy xy y x ；（3）22)()cos(1lim222200y x y x ey x y x +→→++-8. 已知参数⽅程-==tt t y t x sin cos cos ln ，试求出x y d d 和3/22d d π=t x y9. 假设?-=xyt t e y x f 0d ),(2，试求222222y fy x f x f y x ??+-?? 10. 试求出下⾯的极限。
（1）-++-+-+-∞→1)2(1161141121lim 2222n n ；（2）)131211(lim 2222ππππn n n n n n n ++++++++∞→ 11. 试求出以下的曲线积分。
（1）?+ls y x d )(22，l 为曲线)sin (cos t t t a x +=，)cos (sin t t t a y -=，)20(π≤≤t 。
（2）?-+++ly y y xe x e yx )dy 2(xy d )(33，其中l 为22222c y b x a =+正向上半椭圆。
matlab试题及答案
matlab试题及答案# MATLAB试题及答案一、选择题1. MATLAB的基本数据单位是：A. 矩阵B. 向量C. 标量D. 数组答案：A2. 下列哪个命令可以用来绘制函数图形？A. `plot`B. `graph`C. `draw`D. `chart`答案：A3. MATLAB中，以下哪个是正确的矩阵转置操作？A. `transpose(A)`B. `A'`C. `A^T`D. `flip(A)`答案：B二、简答题1. 简述MATLAB中矩阵的基本操作。
答案：在MATLAB中，矩阵是最基本的数据结构，可以进行加、减、乘、除等基本运算。
矩阵的创建可以使用方括号`[]`，例如`A = [1 2;3 4]`。
矩阵的转置使用单引号`'`，例如`A'`。
矩阵的求逆使用`inv`函数，例如`inv(A)`。
2. MATLAB中如何实现循环结构？答案：MATLAB中实现循环结构主要有两种方式：`for`循环和`while`循环。
`for`循环用于已知迭代次数的情况，例如：```matlabfor i = 1:5disp(i);end````while`循环用于迭代次数未知的情况，例如：```matlabi = 1;while i <= 5disp(i);i = i + 1;end```三、计算题1. 给定矩阵A和B，请计算它们的乘积C，并求C的行列式。
A = [1 2; 3 4]B = [5 6; 7 8]答案：首先计算矩阵乘积C：```matlabC = A * B;```然后计算C的行列式：```matlabdetC = det(C);```结果为：```matlabC = [19 22; 43 50]detC = -16```2. 编写一个MATLAB函数，计算并返回一个向量的范数。
答案：```matlabfunction norm_value = vector_norm(v)norm_value = norm(v);end```四、编程题1. 编写一个MATLAB脚本，实现以下功能：- 随机生成一个3x3的矩阵。
MATLAB程序设计期末复习题集
第1章 MATLAB系统环境一、选择题1．最初的MATLAB核心程序是采用（）语言编写的。
A．FORTRAN B．C C．BASIC D．PASCAL2．2016年3月发布的MATLAB版本的编号为（）。
A．MATLAB 2016Ra B．MATLAB R2016aC．MATLAB 2016Rb D．MATLAB R2016b3．下列选项中能反应MATLAB特点的是（）。
A．算法最优 B．不需要写程序C．程序执行效率高 D．编程效率高4．当在命令行窗口执行命令时，如果不想立即在命令行窗口中输出结果，可以在命令后加上（）。
A．冒号（:） B．逗号（,） C．分号（;） D．百分号（%）5．如果要重新执行以前输入的命令，可以使用（）。
A．左移光标键（←） B．右移光标键（→）C．下移光标键（↓） D．上移光标键（↑）6．MATLAB命令行窗口中提示用户输入命令的符号是（）。
A．> B．>> C．>>> D．>>>>7．plot(x,y)是一条（）命令。
A．打印 B．输出 C．绘图 D．描点8．以下两个命令行的区别是（）。
>> x=5,y=x+10>> x=5,y=x+10;A．第一个命令行同时显示x和y的值，第二个命令行只显示x的值B．第一个命令行同时显示x和y的值，第二个命令行只显示y的值C．第一个命令行只显示x的值，第二个命令行同时显示x和y的值D．第一个命令行只显示y的值，第二个命令行同时显示x和y的值9．下列命令行中，输出结果与其他3项不同的是（）。
A．>> 10+20+... B．>> ...30 10+20+30C．>> 10+20+30%5 D．>> %10+20+3010．下列选项中，不是MATLAB帮助命令的是（）。
A．lookfor B．lookfor –all C．search D．help二、填空题1．MATLAB一词来自的缩写。
精讲多练matlab习题
第一章2.设A =1.2,B =-4.6,C =8.0,D =3.5,E =-4.0计算： >> A=1.2; B=-4.6;C=8.0;D=3.5;E=-4.0; >> T=atan(2*pi*A+E/(2*pi*B*C))/D T = 0.4112>> x=pi/180*45;>> (sin(x)+sqrt(35))/72^(1/5) ans = 2.81583.设45x ︒=，计算 >> x=pi/180*45;>> (sin(x)+sqrt(35))/72^(1/5) ans = 2.81584.设5.67,7.811,a b ==计算： >> a=5.67;b=7.811; >> exp(a+b)/log10(a+b) ans =6.3351e+0055.计算2316()( 3.2)/(7.7)y x x x x=++-+在3x =时的值。
>> x=3;>> y=sqrt(x)-6*(x+1/x)+(x-3.2)^2/(x+7.7)^3 y = -18.26796.已知圆的半径为15，求其直径，周长及面积。
>> r=15;d=2*r; >> c=pi*2*r c = 94.2478 >> s=pi*r*r s = 706.85837已知某三角形的三个边的边长为8.5,14.6和18.4，求该三角形的面积。
提示：其中：()/2,s a b c =++a 、b 、c 分别为三角形三边边长。
>>a=8.5;b=14.6;c=18.4; >> s=(a+b+c)/2;>> area=sqrt(s*(s-a)*(s-b)*(s-c)) area = 60.6106第二章1.设矩阵求：>> A=[3 1 1;2 1 2;1 2 3];B=[1 1 -1;2 -1 0;1 -1 1]; （1）>> 2*A+Bans =7 3 16 1 43 3 7（2）>> 4*A*A-3*B*Bans =42 21 3840 19 4640 33 56（3）>> A*Bans =6 1 -26 -1 08 -4 2（4）>> B*Aans =4 0 04 1 02 2 2（5）>> A*B-B*A ans =2 1 -2 2 -2 0 6 -6 02.设三阶矩阵A B 、，满足16,A BA A BA -=+ 其中 求矩阵B 。
(完整版)matlab基础练习题(带答案)
Matlab 基础练习题常量、变量、表达式1、 MATLAB 中，下面哪些变量名是合法的？（ ）（A ）_num （B ）num_ （C ）num- （D ）-num 2、 在MA TLAB 中，要给出一个复数z 的模，应该使用（ ）函数。
（A ）mod(z) （B ）abs(z) （C ）double(z) （D ）angle(z) 3、 下面属于MATLAB 的预定义特殊变量的是？（ ）（A ）eps （B ）none （C ）zero （D ）exp4、 判断：在MA TLAB 的内存工作区中，存放一个英文字符 'a' 需要占用1个字节，存放一个中文字符‘啊’需要占用2个字节。
（ 错，都是2个字节 ）5、 判断：MA TLAB 中，i 和j （ 对 ）6、 判断：MA TLAB 中，pi 代表圆周率，它等于3.14。
（ 错，后面还有很多位小数 ）7、 在MA TLAB 中，若想计算的51)3.0sin(21+=πy 值，那么应该在MA TLAB 的指令窗中输入的MA TLAB 指令是__y1=2*sin(0.3*pi)/(1+sqrt(5))_。
8、 在MA TLAB 中，a = 1，b = i ，则a 占_8__个字节，b 占_16_个字节，c 占________字节。
9、 在MA TLAB 中，inf 的含义是__无穷大__，nan 的含义是__非数（结果不定）___。
数组1、 在MA TLAB 中，X 是一个一维数值数组，现在要把数组X 中的所有元素按原来次序的逆序排列输出，应该使用下面的（ ）指令。
（A ）X[end:1] （B ）X[end:-1:1] （C ）X (end:-1:1) （D ）X(end:1) 2、 在MA TLAB 中，A 是一个字二维数组，要获取A 的行数和列数，应该使用的MATLAB的命令是（ ）。
（A ）class(A) （B ）sizeof(A) （C ）size(A) （D ）isa(A)3、 在MATLAB 中，用指令x=1:9生成数组x 。
MATLAB习题及答案
一、填空题1.MATLAB于1984年由美国Mathworks公司推出，其后每年更新（两）次。
2.MATLAB是一种以（矩阵）运算为基础的交互式程序设计语言。
3.MATLAB具有卓越的数值计算能力和符号计算、文字处理、可视化建模仿真和实时控制等众多功能，其每个变量代表一个（矩阵），每个元素都看作（复数）。
4.通过命令（help）、（lookfor），可以查找所有命令或函数的使用方法。
5.执行语句a=1:2:10，得到的一维数组是（1 3 5 7 9）.6.执行语句b=linspace(1,10,10)后，一维数组b包含（10）个元素，最大值是（10）。
7.函数rem()的功能是取（余）数。
8.若p=[1 0 0;1 1 0],则p|~p=([1 1 1;1 1 1]).（注：填空时请用本题的p的方式表示结果）9.若p=[1 0 0;1 1 0],则all(p)=([1 0 0]).10.矩阵的加减运算，要求相加减的矩阵阶数相同。
若A=[1 2 3 4；2 3 1 8]，则执行语句：[n,m]=size(A)，则n=( 2 ),m=( 4 ).11.对于一维矩阵，求其长度的函数是（length() ）.12.数组和数组之间的运算，尤其是对于乘除运算和乘方运算，如果采用点方式进行计算，表明是数组的（元素）之间的运算关系。
13.求矩阵运算A*B时，要求在维度上，A的（列）数与B的（行）数相等。
二、判断题1.MATLAB只有一种数据类型，一种标准的输入输出语句，不需编译，可直接运行。
（对）2.MATLAB的特殊常量是一些预选定义好的数值变量。
（对）3.MATLAB变量名不区分大小写。
（错）4.i是特殊常量。
（对）5.NAN是非数。
（对）6.MATLAB中所有的变量都表示一个矩阵或一个向量。
（对）7.MATLAB中变量不需要先定义后使用，会自动根据实际赋值的类型对变量类型进行定义。
（对）8.clc命令可以从内存中删除一个、多个和所有变量。
matlab试题及答案
matlab试题及答案一、单项选择题（每题2分，共10分）1. MATLAB中用于创建向量的命令是：A. vectorB. arrayC. linspaceD. colon答案：D2. 在MATLAB中，以下哪个函数用于计算矩阵的行列式？A. detB. rankC. invD. eig答案：A3. MATLAB中用于进行矩阵转置的运算符是：A. 'B. .C. ^D. !答案：A4. 若A是一个3x3的矩阵，执行命令A(2,:)=0;后，矩阵A的第二行将变为：A. [0 0 0]B. [1 0 0]C. [0 1 0]D. [0 0 1]答案：A5. MATLAB中，以下哪个函数用于绘制三维曲面图？A. plotB. surfC. barD. hist答案：B二、填空题（每题3分，共15分）1. MATLAB中，使用________命令可以清除所有变量。
答案：clear2. 若要在MATLAB中创建一个从0到1的100个元素的向量，可以使用命令________。
答案：linspace(0,1,100)3. MATLAB中，使用________函数可以计算矩阵的特征值。
答案：eig4. 在MATLAB中，________函数用于计算两个矩阵的点乘。
答案：dot5. 若要在MATLAB中绘制一个圆，可以使用________函数。
答案：plot三、简答题（每题5分，共20分）1. 请解释MATLAB中矩阵索引的概念。
答案：在MATLAB中，矩阵索引指的是通过行号和列号来访问矩阵中特定元素的过程。
例如，A(2,3)表示访问矩阵A的第二行第三列的元素。
2. MATLAB中如何实现矩阵的元素乘法？答案：在MATLAB中，矩阵的元素乘法可以通过使用点乘运算符（.*）来实现。
例如，C = A .* B，其中A和B是相同大小的矩阵。
3. 请说明MATLAB中如何使用循环结构。
答案：MATLAB中可以使用for循环和while循环两种循环结构。
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习题
MATLAB习题习题一1、通常情况下，MATLAB的工作界面主要由有哪几种窗口构成？答：指令窗口、历史指令窗口、工作空间浏览器、当前目录浏览器等窗口构成。
2、在MATLAB窗口中，可以键入各种指令、函数，并显示运行结果的窗口为【 A 】A、指令窗口B、当前目录浏览器C、工作空间浏览器D、历史指令窗口3、在MATLAB窗口中，可以对内存变量直接进行操作的窗口为【 C 】A、指令窗口B、当前目录浏览器C、工作空间浏览器D、历史指令窗口4、在MATLAB窗口中，记录在MATLAB命令窗口执行过的指令操作的窗口为【 D 】A、指令窗口B、当前目录浏览器C、工作空间浏览器D、历史指令窗口5、在MATLAB窗口中，用于搜索、查看和改变MATLAB路径和文件的窗口为【 B 】A、指令窗口B、当前目录浏览器C、工作空间浏览器D、历史指令窗口6、在MATLAB语言中，如果一个命令行很长，可以在物理行之后如下续行符【 B 】A、“,”B、“…”C、“、、、”D、“，”7、在MATLAB语言中，提示其后语句为输入指令的行提示符为【 D 】A、“;”B、“：”C、“，”D、“>>”8、在MATLAB语言中，在指令窗口运行过的指令再次调出运行的按键为【 B 】A、“;”B、“↑”或“↓”C、“，”D、“>>”9、在MATLAB语言中，若希望结果不被显示，可在语句之后加上符号【 A 】A、“;”B、“：”C、“，”D、“>>”10、在MATLAB语言中，下列常用操作命令解释错误的是【BE 】A、dir为显示当前目录下的文件B、clc为清除工作区中的变量C、path为显示搜索目录D、clf为清空当前图形窗口E、clear为清除命令窗口中的内容F、path为显示搜索路径11、若只知要搜索函数的前两个字母，可用lookfor 命令搜索此关键字的帮助信息。
12、在MATLAB语言中，可采用Tab 键进行模糊查询。
Matlab 课后习题集
Matlab程序设计课程实验目录Matlab操作实验（第1章）实验指导思考练习矩阵运算实验（第2章）实验指导思考练习程序设计实验（第3章）实验指导思考练习高级绘图实验（第4章）实验指导思考练习数据分析与多项式计算实验（第5章）实验指导思考练习解方程与最优化问题求解实验（第6章）实验指导思考练习数值运算实验（第7章）实验指导思考练习符号实验（第8章）实验指导思考练习图形句柄实验（第9章）实验指导思考练习图形用户界面设计实验（第10章）实验指导思考练习Simulink动态仿真实验（第11章）实验指导思考练习外部接口技术实验（第12章）实验指导思考练习1、Matlab操作实验一、实验内容1. 第一章实验内容1：path(path,’系统盘名:\matlab文件名\要添加的文件名’)如：path(path,’c:\matlab701\my work’)实验内容-2例1-1：x=[0:0.5:360]*pi/180;plot(x,sin(x),x,cos(x));例1-2p=[3,7,9,0,-23];x=roots(p)例1-3quad(‘x.*log(1+x)’,0,1)例1-4：实验内容3Help 函数名如：help inv help plot等实验内容4（1）x=0:pi/10:2*pi;y=sin(x);二、思考练习思考题1启动：有三种常见方法，（1）、在Windows桌面，单机任务栏上的‚开始‛按钮，选择‚所有程序‛菜单项，然后选择MATLAB程序组中的MATLAB （版本号）程序选项，即可启动。
（2）、在MATLAB的安装路径中找到MATLAB系统启动程序matlab.exe，然后运行它。
（3）、利用建立快捷方式功能，将MATLAB系统启动程序以快捷方式的形式放在Windows桌面上，以后只要在桌面上双击该图标即可启动MATLAB。
关闭：常见的也有3种方法，（1）、在MATLAB主窗口file菜单中选择exit MATLAB命令。
MATLAB平时练习题(附答案)
一、填空题1、MATLAB常用操作界面包括、工作空间窗口、、、内存数组编辑器、M文件编辑/调试器、帮助导航/浏览器、图形窗口等。
2、MATLAB中Inf或inf表示、NaN或nan表示、nargout表示。
3、MATLAB中逗号主要用作；用作输入量与输入量之间的分隔符；。
4、工作空间浏览器主要用于内存变量的、和。
5、MATLAB实现将全下标转换为单下标的指令为、据单下标换算出全下标的指令为。
6、二维数组的标识有、、“逻辑1”标识。
二、简述下列命令在MATLAB中的作用1.clf %清除图对象clear %清除工作空间内的所有变量clc %清除当前屏幕上显示的所有内容，但不清除工作空间中的数据2.ceil %沿+∞方向取整factor %符号计算的因式分解3.sph2cart %球坐标变为直角坐标pow2 %2的幂4.digits(d) %设置当前的可变算术精度的位数为整数 d 位vpa(pi,100) %设置当前pi的小数位数为1005.box on %打开框状坐标轴开grid off %网格关6.logspace %对数分度向量cat %串接成高维数组7.sym2poly %符号多项式转变为双精度多项式系数向量poly2sym %双精度多项式系数转变为向量符号多项式8.plot3 %三维线图poly2str %以习惯方式显示多项式9.bar %二维直方图pie %二维饼图10.zoom on %打开图形缩放模式 edit %M文件编辑11.whos %对当前工作空间变量的信息进行列表figure %生成图形窗口12.cart2sph %直角坐标变为球坐标pol2cart % 极或柱坐标变为直角坐标13.diff %数值差分、符号微分14.ezplot3 %画三维曲线的简捷指令fix %向零取整factor %符号计算的因式分解三简答题1.叙述MATLAB的主要窗口包括哪些？其作用是什么？命令窗口（Commend Window）——供输入命令，实现计算或绘图功能。
MATLAB习题及答案
填空题1. MATLAB于1984年由美国Mathworks公司推出，其后每年更新（两次。
2. MATLAB是一种以（矩阵）运算为基础的交互式程序设计语言。
3. MATLAB具有卓越的数值计算能力和符号计算、文字处理、可视化建模仿真和实时控制等众多功能，其每个变量代表一个（矩阵），每个元素都看作（复数）。
4.通过命令（help）、（lookfor），可以查找所有命令或函数的使用方法。
5.执行语句a=1:2:10，得到的一维数组是（1 3 5 7 9）.6.执行语句b=linspace（1,10,10）后，一维数组b包含（10）个元素，最大值是10）7.函数rem（）的功能是取（余）数。
8.若p=[1 0 0;1 1 0],则p|〜p=（[1 1 1;1 1 1]）.（注：填空时请用本题的p的方式表示结果）9.若p=[1 0 0;1 1 0],则all（p）=（[1 0 0]）.10.矩阵的加减运算，要求相加减的矩阵阶数相同。
若A=[1 2 3 4；2 3 1 8]，则执行语句：[n,m]=size（A），则n=（2 ）,m=（4 ）.11.对于一维矩阵，求其长度的函数是（length（））.12.数组和数组之间的运算，尤其是对于乘除运算和乘方运算，如果采用点方式进行计算，表明是数组的（元素）之间的运算关系。
13.求矩阵运算A*B时，要求在维度上，A的（列）数与B的（行）数相等。
二、判断题1.MATLAB只有一种数据类型，一种标准的输入输出语句，不需编译，可直接运行。
（对2.MATLAB的特殊常量是一些预选定义好的数值变量。
（对3.MATLAB变量名不区分大小写。
（错4.i是特殊常量。
（对5.NAN是非数。
（对6.MATLAB中所有的变量都表示一个矩阵或一个向量。
（对7.MATLAB中变量不需要先定义后使用，会自动根据实际赋值的类型对变量类型进行定义。
（对8.clc命令可以从内存中删除一个、多个和所有变量。
matlab考试题库
matlab考试题库一、选择题1.以下哪个函数可以在MATLAB中保存变量到文件中？A. savevar()B. writevar()C. savematrix()D. writematrix()2. MATLAB的注释符号是什么？A. //B. ++C. %%D. ##3. MATLAB中，以下哪个函数可以生成一个随机数？A. random()B. rand()C. randomnum()D. randnum()4. MATLAB的关系运算符是什么？A. >>B. <=C. !!D. ^^5. 在MATLAB中，以下哪个函数可以取数组的最大值？A. maximum()B. max()C. big()D. large()二、填空题1. 在MATLAB中，使用______函数可以创建一个全零矩阵。
2. MATLAB中，________运算符可以求两个矩阵的逐元素积。
3. MATLAB中，______函数可以返回行矩阵的列数。
三、简答题1. 请简要介绍MATLAB的基本功能和应用领域。
2. 列举并说明MATLAB中常用的数据类型。
3. 在MATLAB中如何定义一个函数，并举例说明。
四、编程题1. 编写一个MATLAB程序，输入一个矩阵，然后计算该矩阵的转置并输出结果。
2. 编写一个MATLAB函数，输入一个整数n，输出一个n行n列的单位矩阵。
3. 编写一个MATLAB程序，输入一个数组，然后统计该数组中大于10的元素个数，并输出结果。
以上是MATLAB考试题库的部分题目，希望能帮助你进行练习和复习。
祝你考试顺利！。
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文件中. 謀荞抟箧飆鐸怼类蒋薔。
Matlab习题
习题 11. 执行下列指令，观察其运算结果, 理解其意义： (1) [1 2;3 4]+10-2i(2) [1 2; 3 4].*[0.1 0.2; 0.3 0.4] (3) [1 2; 3 4].\[20 10;9 2] (4) [1 2; 3 4].^2 (5) exp([1 2; 3 4]) (6)log([1 10 100]) (7)prod([1 2;3 4])(8)[a,b]=min([10 20;30 40]) (9)abs([1 2;3 4]-pi)(10) [1 2;3 4]>=[4,3;2 1](11)find([10 20;30 40]>=[40,30;20 10])(12) [a,b]=find([10 20;30 40]>=[40,30;20 10]) (提示：a 为行号，b 为列号) (13) all([1 2;3 4]>1) (14) any([1 2;3 4]>1) (15) linspace(3,4,5) (16) A=[1 2;3 4];A(:,2)2. 执行下列指令，观察其运算结果、变量类型和字节数，理解其意义： (1) clear; a=1,b=num2str(a),c=a>0, a= =b, a= =c, b= =c (2) clear; fun='abs(x)',x=-2,eval(fun),double(fun)3. 本金K 以每年n 次，每次p %的增值率(n 与p 的乘积为每年增值额的百分比)增加，当增加到rK 时所花费的时间为)01.01ln(ln p n rT +=(单位：年)用MA TLAB 表达式写出该公式并用下列数据计算：r =2, p =0.5, n =12.4．已知函数f (x )=x 4-2x 在(-2, 2)内有两个根。
取步长h =0.05, 通过计算函数值求得函数的最小值点和两个根的近似解。
Matlab练习题
Matlab练习题Matlab练习题１已知矩阵11 12 13 1421 22 23 24A= 31 32 33 3441 42 43 44(1)A(:,1) (2) A(2,:) (3) A(:,2:3)(4) A(2:3,2:3) (5) A(:,1:2:3) (6) A(2:3)(7) A(:) (8) A(:,:) (9) ones(2,2)(10) eye(2) (11) [A,[ones(2,2);eye(2)]](12) diag(A) (13) diag(A,1) (14) diag(A,-1)(15) diag(A,2)2(1)输入如下矩阵Ａ０π/3Ａ＝π/6 π/2(2) 求矩阵B1,B1中每一元素为对应矩阵Ａ中每一元素的正弦函数(3) 求矩阵B2, B2中每一元素为对应矩阵Ａ中每一元素的余弦函数(4) 求B12+B22(5) 求矩阵Ａ的特征值与特征矢量：称特征矢量为Ｍ，而特征值矩阵为Ｌ(6) 求Msin(L)M-1(7)使用funm命令求矩阵A的正弦函数（结果应该与（６）同）(8)求cosA(9) 证明sin2A+cos2A=I3 按题目要求用MATLAB命令完成下列矩阵运算(1) 使用rand命令产生5个2x2随机矩阵A,B,C,D,E(2) 求矩阵F(使用和不使用inv命令两种情况) F=A-1[B+C-1(D-1E)]4 手算和上机分别求A.*B‘ A.\B‘ A‘.\B其中A=[1;1;1] B=[2,3,4]5 已知A=[2 7 6;9 0 –10;3 0.5 6]; B=[8 0.2 0;3 2 5;4 0 7];求(1)A|B, A&B, A~B, A>B, A>=B, A<="B," a~="B</p">(2)元素值为零的元素标号(3)元素值大于6的元素标号6 某专业有三名研究生,本学期选修了四门课程,若这些研究生的姓名,学号,性别,出生年月,课程名称,考试成绩可任意假定,(1)分别用结构型变量和细胞型变量表示以上信息;(2)举例说明查阅以上任何一条信息的方法;(3)求每一个研究生的平均成绩.7 已知矩阵A=[1.2 3 5 0.9;5 1.7 5 6;3 9 0 1;1 2 3 4],求(1)A的特征多项式(2)特征多项式中未知数为20 时的值(3)特征多项式的根(4)特征多项式的导数8 已知五个数据点: (1,5.5), (2,43.1), (3,128), (4,290.7), (5,498.4)(1)用三次曲线拟合上述数据点(2)在同一图形中绘出数据点和拟合曲线(3)用适当图形表示拟合精度9 在实验中测得如下6组数据:(0,1.5), (1.5,3.4), (2.8,13), (3.8,36), (4.5,63), (4.9,78)请用三次曲线拟合以上数据并给出以下结果：(1)三次多项式的各项系数；(2)将数据点和拟合曲线以最佳效果在同一图中绘出；10 求下列函数的极限（1）lim(x2/sin2(x/3)) x 0 （2）lim((tanx-sinx)/sin3x) x 0 （3）lim(sin(a+x)-sin(a-x))/x x 0 （4）limxcos(1/x) x 0（5）lim((1+mx)n-(1+nx)m)/x2 x 0 （6）lim(1+1/n)(n+5) n∞（7）lim(1-2/x)x x ∞（8）lim(1+cosx) 3secx x л/2 11 求下列函数的积分（1）x2/sin2(x/3) （2）(tanx-sinx)/sin3x（3）(sin(a+x)-sin(a-x))/x （4）xcos(1/x)（5）((1+mx)n-(1+nx)m)/x2 (6) cos2x(7) sinaxcosbx (8)cosaxsinbx(9) arcsin(x/a) (10)1/(a+bsinx)(11) xarcsin(x/a) (12) x2arcsin(x/a)12 求下列函数的反函数（1）y=2sin3x （2）y=1+ln(x+2)（3）y=2x/(2x+1) （4）y=(ex+e(-x))/2（5）y=1/2(arccos(x/2)) （6）y=x+1/x13 求下列函数的定积分（1）(x+sinx)/(1+cosx) [0,pi/2] （2）ln(1+tanx ) [0,pi/4] （3）1/(1+cos2x ) [0,pi/2] （4）cos5xsinx [0,pi/2] （5）(3x4+3x2+1)/(x2+1) [-1,0] (6) x2+1/x4 [1,2] (7) tan2 x [0,pi/4] (8) 4cos4x [-pi/2,pi/2] (9) 1-sin3x [0,pi] (10)1/(11+5x) 3 [-2,1] (11) cosxcos2x [-pi/2 , pi/2] (12) (x 3 sin 2 x )/(x 4 + 2x 2 + 1) [-5,5] 14 求解代数方程(1) ax2+bx+c=0(2) cos(2x)+sin(x)=115 解线性方程组a1x+b1y+c1z=d1a2x+b2y+c2z=d2 其中ai, bi, ci, di为常数, x, y, z为变量a3x+b3y+d3z=d316 求方程tan(x)+sin(x)=2在区间[-2π, 2π]上的全部实数解并用图示。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
第2章 MATLAB 概论
1、 与其他计算机语言相比较，MA TLAB 语言突出的特点是什么？
2、 MA TLAB 系统由那些部分组成？
3、 安装MA TLAB 时，在选择组件窗口中哪些部分必须勾选，没有勾选的部分以后如何补安装？
4、 MA TLAB 操作桌面有几个窗口？如何使某个窗口脱离桌面成为独立窗口？又如何将脱离出去的窗口重新放置到桌面上？
5、 如何启动M 文件编辑/调试器？
6、 存储在工作空间中的数组能编辑吗？如何操作？
7、 命令历史窗口除了可以观察前面键入的命令外，还有什么用途？
8、 如何设置当前目录和搜索路径，在当前目录上的文件和在搜索路径上的文件有什么区别？
9、 在MA TLAB 中有几种获得帮助的途径？
第3章 MATLAB 数值运算
3.1 在MA TLAB 中如何建立矩阵⎥⎦
⎤
⎢
⎣⎡19
4
375，并将其赋予变量a ？
3.2 有几种建立矩阵的方法？各有什么优点？
3.3 在进行算术运算时，数组运算和矩阵运算各有什么要求？
3.4 数组运算和矩阵运算的运算符有什么区别？
3.5 计算矩阵⎥⎥⎥⎦⎤
⎢⎢⎢
⎣⎡89
7
473
535与⎥⎥⎥⎦
⎤
⎢⎢⎢⎣⎡63
8
976
242
之和。
3.6 求⎥⎦
⎤
⎢
⎣⎡+-+-+-+-++=i 44i
93i
49i
67i
23i 57i 41i 72i 53i 84x
的共轭转置。
3.7 计算⎥⎦
⎤⎢⎣⎡=57
2
396a
与⎥⎦
⎤⎢⎣⎡=86
4
142b
的数组乘积。
3.8 “左除”与“右除”有什么区别？
3.9 对于B AX =，如果⎥⎥⎥⎦
⎤
⎢⎢⎢⎣⎡=75
3
467
294A ，⎥⎥
⎥
⎦
⎤
⎢⎢⎢⎣⎡=282637B ，求解X 。
3.10 已知：⎥⎥⎥⎦
⎤
⎢⎢⎢⎣⎡=98
7
654
321
a ，分别计算a 的数组平方和矩阵平方，并观察其结果。
3.11
⎥⎦
⎤⎢
⎣⎡-=46
3
521
a ，⎥⎦
⎤⎢⎣⎡-=26
3
478b
，观察a 与b 之间的六种关系运算的结果。
3.12 []7.0802.05--=a ，在进行逻辑运算时，a 相当于什么样的逻辑量。
3.13 在sin(x )运算中，x 是角度还是弧度？
3.14 角度[]6045
30
=x ，求
x 的正弦、余弦、正切和余切。
3.15 用四舍五入的方法将数组[2.4568 6.3982 3.9375 8.5042]取整。
3.16
矩阵⎥⎥⎥⎦⎤
⎢⎢⎢
⎣⎡=72
8
365
219a ，分别对a 进行特征值分解、奇异值分解、LU 分解、QR 分解
3.17 将矩阵⎥⎦
⎤⎢⎣⎡=57
24a
、⎥⎦
⎤⎢⎣⎡=3817b
和⎥⎦
⎤⎢⎣⎡=2695c
组合成两个新矩阵：
3.18 将(x -6)(x -3)(x -8)展开为系数多项式的形式。
3.19 求解多项式x 3-7x 2+2x +40的根。
3.20 求解在x =8时多项式(x -1)(x -2) (x -3)(x -4)的值。
3.21 计算多项式乘法(x 2+2x +2)(x 2
+5x +4)。
3.22 计算多项式除法(3x 3+13x 2
+6x +8)/(x +4)。
3.23 对下式进行部分分式展开：
2
72436
45232
3
4
5
2
3
4
+++++++++x x x x x
x x
x
x
3.24 计算多项式9514124234++--x x x x 的微分和积分。
3.25 解方程组⎥⎥⎥⎦⎤
⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡661362
21143092x 。
3.26 矩阵⎥⎥⎥⎦
⎤⎢⎢⎢⎣⎡-=943
457624a ，计算a 的行列式和逆矩阵。
第4章 MATLAB 符号运算
4.1 创建符号变量有几种方法？
4.2 下面三种表示方法有什么不同的含义？ （1）f=3*x^2+5*x+2 （2）f='3*x^2+5*x+2' （3）x=sym('x') f=3*x^2+5*x+2
4.3 用符号函数法求解方程a t 2+b*t +c=0。
4.4 用符号计算验证三角等式：
sin(ϕ1)cos(ϕ2)-cos(ϕ1)sin(ϕ2) =sin(ϕ1-ϕ2)
4.5 求矩阵⎥⎦
⎤
⎢⎣⎡=22211211a a a a A 的行列式值、逆和特征根。
4.6 因式分解：6555234-++-x x x x 4.7 ⎥⎥⎦
⎤⎢⎢
⎣⎡=)sin()
log(1
2
x x e
x
x
a f ax ，用符号微分求df/dx 。
4.8 求代数方程组⎪⎩⎪⎨
⎧=+=++0
2y x c by ax 关于x,y 的解。
4.9 符号函数绘图法绘制函数x=sin(3t)cos(t)，y=sin(3t)sin(t)的图形，t 的变化范围为[0,2π]。
4.10 绘制极坐标下sin(3*t)*cos(t)的图形。
第5章MATLAB 程序设计
1 命令文件与函数文件的主要区别是什么？
2 如何定义全局变量？
3 如果x 是一个结构型数组，如何观察其中的内容？
4 if 语句有几种表现形式？
5 说明break 语句和return 语句的用法。
6 有一周期为4π的正弦波上叠加了方差为0.1的正态分布的随机噪声的信号，用循环结构编制一个三点线性滑动平均的程序。
（提示：①用0.1*randn(1,n)产生方差为0.1的正态分布的随机噪声；②三点线性滑动平均就是依次取每三个相邻数的平均值作为新的数据，如x1(2)=(x(1)+x(2)+x(3))/3，x1(3)=( x(2)+x(3)+x(4))/3……）
7 编制一个解数论问题的函数文件：取任意整数，若是偶数，则用2除，否则乘3加1，重复此过程，直到整数变为1。
8 有一组学生的考试成绩（见表），根据规定，成绩在100分时为满分，成绩在90~99之间时为优秀，成绩在80~89分之间时为良好，成绩在60~79分之间为及格，成绩在60分以下时为不及格，编制一个根据成绩划分等级的程序。
9 编写一段程序，能够把输入的摄氏温度转化成华氏温度，也能把华氏温度转换成摄氏温度。
10 keyboard 命令的作用是什么？当在函数中使用keyboard 命令时，是否可以在工作空间浏览器中观察和修改函数中的局部变量？退出keyboard 状态的命令是什么？编一个函数文件，在其中设置keyboard ，观察变量
情况。
第6章 常用图形命令
1 绘制曲线13++=x x y ，x 的取值范围为[-5,5]。
2 有一组测量数据满足-at e =y ，t 的变化范围为0~10，用不同的线型和标记点画出a=0.1、a=0.2和a=0.5三种情况下的曲线。
3 在上题结果图中添加标题-at e =y ，并用箭头线标识出各曲线a 的取值。
4 在上题结果图中添加标题-at e =y 和图例框。
5
6 x= [66 49 71 56 38]，绘制饼图，并将第五个切块分离出来。
7 2
2
y
x
xe z --=，当x 和y 的取值范围均为-2到2时，用建立子窗口的方法在同一个图形窗口中绘制出三
维线图、网线图、表面图和带渲染效果的表面图。
8 绘制peaks 函数的表面图，用colormap 函数改变预置的色图，观察色彩的分布情况。
9 用sphere 函数产生球表面坐标，绘制不通明网线图、透明网线图、表面图和带剪孔的表面图。
10 将上题中的带剪孔的球形表面图的坐标改变为正方形，以使球面看起来是圆的而不是椭圆的，然后关闭坐标轴的显示。
第7章 Simulink 基础
7.1 什么是Simulink ?
7.2 如何进行下列操作： （1）翻转模块
（2）给模型窗口加标题 （3）指定仿真时间
（4）设置示波器的显示刻度
7.3 有传递函数如下的控制系统，用Simulink 建立系统模型，并对系统的阶跃响应进行仿真。
8
41
)(2
++=
s s
s G
7.4 建立一个简单模型，用信号发生器产生一个幅度为2V 、频率为0.5Hz 的正弦波，并叠加一个0.1V 的噪声信号，将叠加后的信号显示在示波器上并传送到工作空间。
7.5 建立一个简单模型，产生一组常数（1×5），再将该常数与其5倍的结果合成一个二维数组，用数字显示器显示出来。
7.6 建立一个模拟系统，将摄氏温度转换为华氏温度（T f = 9/5T c +32）。
7.7 建立二阶系统2
2
2)(n
n n
S S
S G ωςωω++=的脉冲响应模型，设ωn =10Hz ，观察当0<ζ<1、ζ=0、ζ=1
及ζ>1时系统的响应。