《MATLAB语言程序设计》练习题1(08通信)

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)。

电子信息工程系2008级《MATLAB语言》试卷A参考答案1.Fill in the blanks ( 2’ * 5 = 10’ )(1)letter / 字母(2)0(3)imag( )(4)toc(5)eye (100, 102) 或[ eye (100), zeroes (100, 2) ]2.True or False ( 2’ * 5 = 10’)(1)T(2)T(3) F(4) F(5)T3.Short Ans wers ( 5’ * 3 = 15’ )(1)a) Add semicolons to the end of the statements.在句末加分号；b) Preallocate all the arrays used in loops.对所有循环体中用到的矩阵进行预分配空间；c) V ectorization向量化编程；d) Call a function by its handle.通过函数句柄调用函数。

注：答对一条给2分，答对两条给4分，答对三条给5分。

(2)a) Clearly state the problem that you are trying to solve.清晰的描述待解决的问题；b) Define the inputs required and the outputs to be produced by theprogram.定义程序需要的输入变量和返回的输出变量；c) Design the algorithm that you intend to implement in the program.设计程序中需要用到的算法；d) Turn the algorithm into MATLAB statements.将算法转化为MATLAB语句；e) Test the resulting MATLAB program.测试生成的MATLAB程序。

注：回答时只需回答出黑体部分的关键句即可，每条1分。

matlab考试题目及答案1. 题目：编写一个MATLAB函数，实现计算并返回一个向量中所有元素的平方和。

答案：以下是一个简单的MATLAB函数，用于计算并返回一个向量中所有元素的平方和。

```matlabfunction sumOfSquares = calculateSumOfSquares(vector)sumOfSquares = sum(vector.^2);end```2. 题目：给定一个3x3的矩阵A，使用MATLAB编写代码，求出矩阵A 的转置。

答案：可以通过简单的转置操作来求得矩阵A的转置。

```matlabA = [1 2 3; 4 5 6; 7 8 9];A_transpose = A';```3. 题目：编写一个MATLAB脚本，实现对一个二维数组进行排序，并返回排序后的数组。

答案：以下是一个MATLAB脚本，用于对一个二维数组进行排序，并返回排序后的数组。

```matlabfunction sortedArray = sort2DArray(array)sortedArray = sort(array(:));end```4. 题目：给定一个向量x，使用MATLAB编写代码，计算并返回向量x的元素个数。

答案：可以通过内置函数`numel`来计算向量x的元素个数。

```matlabx = [1, 2, 3, 4, 5];numElements = numel(x);```5. 题目：编写一个MATLAB函数，实现计算并返回两个向量元素的点积。

答案：以下是一个简单的MATLAB函数，用于计算两个向量的点积。

```matlabfunction dotProduct = calculateDotProduct(vector1, vector2)dotProduct = dot(vector1, vector2);end```6. 题目：给定一个矩阵B，使用MATLAB编写代码，求出矩阵B的行列式。

安徽农业大学2008―2009学年第二学期《MATLAB 程序设计》试卷（A 卷）考试形式: 闭卷笔试，2小时 适用专业： 07级通信工程一、填空题(本题满分30分，每空3分)1.设有程序A=[2,0,0,9;0,5,3,1;0, 3 ,1, 2];B=[1,2,-1;1,3,6];D=A([1 3] ,1:3); E=size(D)+ length(B(:,3));F= D.*B 将下列命令的运行结果填在横线上D= ;E = ; F= .2. A=[1,2,3,4;2,3,4,5;3,4,5,6]; A1=sum(A>2,2);A(:,4)=[],A2=A.^2A1= ; A2= ;3.P=[1,2,5,4;5,4,6,1];max(P,[ ],2)= ;mean(P,2)= . 院： 专业班级： 姓名： 学号：装 订 线二、(本题满分12分)试编写计算程序画出曲线与曲面图形.22,22,(2).;),12sin((1).222.0≤≤-≤≤-++=≤≤+=-yxyxyxzxxey xπ三、(本题满分8分)V.与特征向量U 的特征值A 求(4).;A )(解矩阵.(3));A B)((A*R R 秩(2).|;B))(A B 2(A |D (1).,480200192,32112032111E B B X A A B A T T T T -=--=-+=⎪⎪⎪⎭⎫⎝⎛=⎪⎪⎪⎭⎫ ⎝⎛--=--方程计算矩阵的计算行列式设四、(本题满分16分)运用符号运算功能写出下列各题Matlab 程序解求方程1sin 3.12=-x xMdx dx dM x xx e e M x x ⎰⎥⎦⎤⎢⎣⎡=--π01.0,,3cos 2sin sin .2并计算生成符号矩阵3. 解求微分方程t e x x x -=-'-''654. ).()(,)(2s F Laplace t f e t t f t 变换函数的求设-=五、（本题满分12分）根据药学的蛋白质结合理论，蛋白的某种单分子层具有吸附功能。

matlab语言练习题MATLAB是一种面向科学计算和工程应用的高级编程语言和环境，广泛应用于各个领域。

通过进行一些MATLAB语言的练习题，可以帮助我们更好地理解和熟练应用这门语言。

本文将介绍一些常见的MATLAB语言练习题，并逐步给出解答过程。

一、编写函数1. 编写一个函数，求解给定数列的和。

数列以向量的形式输入，函数应返回数列的和。

解答：```matlabfunction sum_res = sum_vector(vector)sum_res = sum(vector);end```2. 编写一个函数，判断给定的数是否是素数。

素数是指只能被1和自身整除的正整数。

解答：```matlabfunction is_prime = check_prime(num)is_prime = true;if num < 2is_prime = false;elsefor i = 2:sqrt(num)if rem(num, i) == 0is_prime = false;break;endendendend```二、图形绘制1. 绘制正弦曲线。

使用`sin`函数和绘图命令`plot`绘制一个在0到2π范围内的正弦曲线。

解答：```matlabx = linspace(0, 2*pi, 100); % 生成0到2π的100个等间距点y = sin(x);plot(x, y);xlabel('x');ylabel('sin(x)');title('Sine Curve');```2. 绘制柱状图。

使用`bar`函数绘制一个柱状图，展示不同城市的降雨量。

解答：```matlabcities = {'City A', 'City B', 'City C', 'City D'};rainfall = [80, 60, 120, 90];bar(rainfall);xticks(1:4);xticklabels(cities);ylabel('Rainfall (mm)');title('Rainfall in Different Cities');```三、矩阵运算1. 矩阵相乘。

课程号：1003025《Matlab语言基础》期末考试试卷（A卷）考试形式：闭卷考试考试时间：120分钟班号学号姓名得分一、（5分）简答：分别说明Matlab命令clc、clear、close的作用及调用方式。

二、（5分）简答：Matlab与一般程序设计语言相比数据结构有何特殊性，为什么说Matlab有很强的扩充性？三、（5分）简答：简述Matlab的矩阵运算、矩阵的数组运算和符号运算的本质及不同之处。

四、（10分）简答：输入一个3行2列的数值矩阵、一个2行3列的字符串矩阵，最后生成一个2行1列的元胞矩阵，要求元胞矩阵的两个元素分别是刚输入的数值矩阵和字符串矩阵；比较这三类矩阵输入方式的异同。

五、（15分）编程：公历闰年判定遵循的规律为：四年一闰，百年不闰，四百年再闰；即能被4整除而不能被100整除或能被400整除的年份为闰年。

分别使用for循环和while循环编写一段Matlab程序，判断并输出公元2000年至3000年间（包含2000年和3000年）闰年的个数，并屏幕输出闰年年份。

求余函数的Matlab帮助：modModulus after divisionSyntaxM = mod(X,Y)Examplesmod(13,5)ans =3六、（20分）编程：单窗口多曲线分幅绘图。

要求如下：①（2分）新建一个绘图窗口，分两行一列进行绘图；②（8分）第一幅绘制sin曲线，要求x轴取对数，使用红色实线绘制，并对图形使用所有加注（xlabel、ylabel、title、text、legend、grid、gtext）；③（10分）第二幅绘制cos曲线，要求使用直角坐标系，并使用句柄绘图功能改变窗口名称（name）；窗口左下角坐标为[0.1 0.1]，右上角坐标为[0.7 0.8]（units，position）；把y轴设置为一个[R G B]的自定义色（ycolor）、把x轴翻转（xdir）、y轴方向绘制网格线（ygrid）；要求把点型改为圆形（linestyle）、把点符号大小改为12（markersize）、把点符号边缘色改为红色（markeredgecolor）、把点符号填充颜色改为灰色（markerfacecolor）、把线对象的y值改为原先的1/2（ydata）。

《MATLAB 语言与应用技术》2008年试题说明： 1，从3题中任选2题；2，交电子文档（包括程序、示例及作业）及打印稿各一份； 3，完成时间：2009年7月1日前。

1， 如图1所示，有三个小车（重物）用弹簧连接，开始处于平衡状态，当分别施加外力P 1, P 2, P 3时，小车状态将改变，三个小车从其初始平衡位置（x i =0，i 值取1，2，3）处开始振动，其位移可通过求系统的势能（PE ）的极小值求得：P X KX X PE ''5.0-=其中，⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡++---++---++=654545532343321k k k k k k k k k k k k k k k K []'=321P P P P[]'=321x x x X输入数据为：k 1=4500N/m k 4=2250N/m P 1=1100Nk 2=1650N/m k 5=550N/m P 2=1800Nk 3=1100N/m k 6=9300N/m P 2=3300N试求出三个小车的平衡位置（提示：用fminunc 函数）。

图12，图2是一个电动自行车的方框图。

其中，倾斜仪是用来检测车体偏离垂直方向角度y (t )的。

倾斜仪的输出与期望垂直角度r (t )相比较，将误差送回控制器并产生操纵信号。

系统的任何扰动d (t )都作为车子的输入。

方框图中的传递函数为：29()9G s s =+，222()2F s s s ωζωω=++，()s z C s k s p +=+ 其中，C (s )为超前控制器。

①倾斜仪的调整时间为0.2s ，带宽125Hz 。

传感器的参数250ωπ=，20/ζω=。

当单位阶跃响应的超调小于20％，调节时间小于0.4s 时，求k 、z 、p 的值。

②假设有另一类倾斜仪可选，它基于钟摆原理，其谐振频率为7.4Hz( 1.8ωπ=),阻尼系数0.4ζ=。

试分析该倾斜仪是否可用。

第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))二、实验题1.x=-74/180*pi;y=-27/180*pi;sin(x*x+y*y)/sqrt(tan(abs(x+y)))+pi2.a=-3.0:0.1:3.0;exp(-0.3*a).*sin(a+0.3)3.x=[2,4;-0.45,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^3ans =-28917 240246 -4368137883 -259101 27669171333 252504 38673A.^3ans =27 157464 839304 -91125 343658503 729000 3375（3）A/Bans =-5.8845 5.3549 -0.20286.3554 -5.6596 4.3293-6.6325 6.2078 9.0241B/Aans =30.2855 19.2643 -8.561433.4394 21.1547 -9.3974-0.7443 -0.3938 0.2830（4）[A,B]ans =3 54 2 1 -2 6734 -45 7 2 8 7487 90 15 9 3 0[A([1,3],:);B^2]ans =3 54 287 90 15600 183 -81684 282 72615 6 8255.a=1+2i;b=3+4i;c=exp(pi*i/6)c =0.8660 + 0.5000ic+a*b/(a+b)ans =1.6353 + 1.8462i第2章一、思考题3.s=0;for n=0:63s=s+2^n;enddisp(s)n=0:63;s=sum(2.^n)二、实验题1.x=input('输入一个四位整数：');y=[fix(x/1000),mod(fix(x/100),10),mod(fix(x/10),10),mod(x,10)] 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+1.15*dj*(h-120);elseif h<60gz=dj*h-700;elsegz=dj*h;endformat bank;display([gh,gz])3.循环结构n=input('input n:');s=0;for k=1:ns=s+1/k^2;enddisplay(sqrt(s*6))向量运算n=input('input n:');k=1:n;display(sqrt(sum(1./k.^2)*6))4.y=0;k=0;while y<3k=k+1;y=y+1/(2*k-1);enddisplay([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<500x0=x;x=a/(b+x0);k=k+1;enddisplay([k,x]);display([(-b+sqrt(b^2+4*a))/2,(-b-sqrt(b^2+4*a))/2]);6.y=fun(40)/(fun(30)+fun(20))（1）函数文件fun.mfunction f=fun(n)f=n+log(n^2+5);（2）函数文件fun.mfunction f=fun(n)a=1:n;f=sum(a.*(a+1));第3章一、思考题4.t=0:0.01:1.5;y=sqrt(3)/2*exp(-4*t).*sin(4*sqrt(3)*t+pi/3);5.x=-10:0.01:10;y=linspace(-6,6,size(x,2))z=x.^3+3*x.*y.^2;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:0.01:10;plot(x,x-x.^3/3/2)（2）plot(x,exp(-x.*x/2)/2/pi)（3）x=-8:0.01: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*x1-0.5;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(0.4,1,'y=2x-0.5');text(-0.4,0.5,'x=sin(3t)cos(t)');text(-0.4,0.3,'y=sin(3t)sin(t)');（2）subplot(1,2,1);scatter(x1,y1,10)title('y=2x-0.5');subplot(1,2,2);scatter(x,y,10)3.subplot(1,21);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);5.theta=linspace(-pi,pi,100);ro=5.*cos(theta)+4;polar(theta,ro);fi=linspace(0,2*pi,100);a=1r=a.*(1+cos(fi));polar(fi,r);6.t=0:pi/20:2*pi;x=exp(-t./20).*cos(t);y=exp(-t./20).*sin(t);z=t; plot3(x,y,z);t=0:0.01:1;x=t;y=t.^2;z=t.^3;plot3(x,y,z);7.x=-30:0.1:0;y=0:0.1:30;[x,y]=meshgrid(x,y);z=10.*sin(sqrt(x.^2+y.^2))./sqrt(1+x.^2+y.^2); meshc(x,y,z);8.x=linspace(-3,3,100);y=linspace(-3,3,100);[x y]=meshgrid(x,y);fxy=-5./(1+x.^2+y.^2);i=find(abs(x)<=0.8 & abs(y)<=0.5);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.^2;x=3*u.*sin(v);y=2*u.*cos(v);z=4*u.^2;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(0.2)*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,4.5e-4;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>0.5))/size(A)2.h=[466,715,950,1422,1635];w=[7.04,4.28,3.40,2.52,2.13];hh=[500,900,1500];ww=interp1(h,w,hh,'spline')3.x=linspace(1,10,50);y=log(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); % 给出N个采样时间ti(I=1:N) x=3*exp(-t); % 求各采样点样本值xdt=t(2)-t(1); % 采样周期f=1/dt; % 采样频率(Hz)X=fft(x); % 计算x的快速傅立叶变换XF=X(1:N/2+1); % F(k)=X(k)(k=1:N/2+1) f=f*(0:N/2)/N; % 使频率轴f从零开始plot(f,abs(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）建立函数文件myfun.mfunction 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','off')) 第6章一、思考题2.fx=inline('1./(1+x.^2)');[I,n]=quad(fx,-100000,100000,1e-3);[I,n]=quadl(fx,-100000,100000,1e-3);x=-100000:0.01:100000;y=1./(1+x.*x);trapz(x,y);3.(1)fx=inline('-2*y+2*x*x+2*x');[t,y]=ode23(fx,[0,0.5],1)(2)fx=inline('y-exp(x)*cos(x)');[t,y]=ode23(fx,[0,3],1)二、实验题1.for x=1:3fx=[x,x^2,x^3;1,2*x,3*x;0,2,6*x];diff(fx)end2.（1）x=0:0.01:1;y=x.^10+10.^x+1./log10(x);dy=diff(y)/0.01;（2）x=0:0.01:1;y=log(1+x);dy=diff(y,2)/0.01;plot(x(1:99),dy)3.（1）fx=inline('x.^2.*sqrt(2*x.*x+3)');quad(fx,1,5)（2）fx=inline('x./sin(x).^2');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,y^2,y+2),-1,2)x 的积分区间为【0，2】时fx=inline('x.*y');dblquad(fx,0,2,-1,2)4.x=0.3:0.2:1.5;y=[0.3895,0.6598,0.9147,1.1611,1.3971,1.6212,1.8325]; trapz(x,y)5.（1）yp=inline('-(1.2+sin(10*x))*y');[t,y]=ode23(yp,[0,5],1);（2）令''3,,21y x y x y x ='==，则可写出原方程的状态方程形式：123233221sin 31)1(2cos 5cos x tx x t t t x x x x x +--++='='='⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡+⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡+-+-=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡'''t x x x t t t x x x cos 00)1(2cos 51sin 311000103212321 建立函数文件ztfun.mfunction y=ztfun(t,x)b=[0;0;cos(t)];y=[0,1,0;0,0,1;-1/(3+sin(t)),-1,5*cos(2*t)/(t+1)^2]*x+b; 解微分方程组[t,y]=ode23(@ztfun,[0,5],[1;0;2]);6.建立函数文件ztfun.mfunction yy=ztfun(t,y)yy=[y(2)*y(3);-y(1)*y(3);-0.51*y(1)*y(2)];解微分方程组[t,y]=ode23(@ztfun,[0,5],[0;1;1])第7章一、思考题3.（1）数值积分fx=inline('exp(x).*(1+exp(x)).^2');quad(fx,0,log(2))符号积分f=sym('exp(x)*(1+exp(x))^2');v=int(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=sin(x^2*y)*exp(-x^2-y);diff(diff(fxy,x),y)3.(1) syms xint(1/(1+x^4))(2) syms x tint((-2*x*x+1)/(2*x*x-3*x+1)^2,x,cos(t),exp(2*t)) 4.syms n xsymsum(1/(2*n+1)/(2*x+1)^(2*n+1),n,0,inf)symsum(1/(2*n+1)/(2*x+1)^(2*n+1),n,0,5)5.(1) syms xtaylor((exp(x)+exp(-x))/2,5,0)(2) syms a xtaylor(exp(-5*x)*sin(3*x+pi/3),5,a)6.(1)x=solve(sym('x^3+a*x+1=0'))(2) [x y]=solve(sym('sqrt(x^2+y^2)-100=0,3*x+5*y-8=0'))7. 方程转化为： π21'11't y y y y -=+=⎩⎨⎧符号解[y1,y11]=dsolve('Dy=y1,Dy1+y=1-t^2/pi','y(-2)=5,y1(-2)=5','t') 数值解编写函数文件ztfun.mfunction yy=ztfun(t,y)yy=[y(2);1-t^2/pi-y(1)];在命令窗口输入以下命令[t,y]=ode45(@ztfun,[-2,7],[-5;5]);t=linspace(-2,7,49)y2=y8.[x,y]=dsolve('Dx=3*x+4*y,Dy=-4*x+3*y','x(0)=0,y(0)=1')。

matlab考试题及答案百度网盘MATLAB考试题及答案百度网盘1. MATLAB的基本数据类型有哪些？- 答案：MATLAB的基本数据类型包括数值型（如整数、浮点数）、字符型、逻辑型、结构体、单元数组和对象。

2. 如何在MATLAB中创建一个3x3的单位矩阵？- 答案：可以使用`eye(3)`命令来创建一个3x3的单位矩阵。

3. MATLAB中如何实现矩阵的转置？- 答案：可以使用`.'`操作符或者`transpose()`函数来实现矩阵的转置。

4. 在MATLAB中，如何计算两个矩阵的点积？- 答案：可以使用`dot()`函数或者`*`操作符（当两个矩阵都是向量时）来计算两个矩阵的点积。

5. MATLAB中如何绘制函数f(x)=x^2的图像？- 答案：可以使用`plot()`函数结合`fplot()`函数或者直接使用`plot()`函数配合`linspace()`函数来绘制函数f(x)=x^2的图像。

6. MATLAB中如何求解线性方程组Ax=b？- 答案：可以使用`\`操作符或者`linsolve()`函数来求解线性方程组Ax=b。

7. MATLAB中如何实现循环结构？- 答案：可以使用`for`循环或者`while`循环来实现循环结构。

8. MATLAB中如何定义一个函数？- 答案：可以使用`function`关键字来定义一个函数，例如：`function y = myFunction(x) y = x^2; end`。

9. MATLAB中如何使用条件语句？- 答案：可以使用`if`、`elseif`和`else`关键字来实现条件语句。

10. MATLAB中如何实现文件的读写操作？- 答案：可以使用`fopen()`、`fread()`、`fwrite()`和`fclose()`等函数来实现文件的读写操作。

以上为MATLAB考试题及答案的示例，具体内容可根据实际考试要求进行调整。

MATLAB作业11 、在你的机器上安装MATLAB 语言环境，并键入demo 命令，由给出的菜单系统和对话框原型演示程序，领略MATLAB 语言在求解数学问题方面的能力与方法。

>> demo>>左击MATLAB前‘+’，Programming前‘+’，左击Manipulating Multidimensional ArraysMATLAＢ具备高性能的数值计算和可视化的科学工程计算功能，支持解释性语言输入，编程实现简单，具有丰富的数学函数功能支持。

MATLAB允许与Ｃ、Fortran语言接口，其部件Simulink可以采用图形输入的方式来搭构所研究的系统。

包含丰富的工具包，在系统仿真、数字信号处理、图形图像分析、数理统计、通信及自动控制领域得到广泛的应用。

总的来说，该软件有三大特点。

一是功能强大。

具有数值计算和符号计算、计算结果和编程可视化、数学和文字统一处理、离线和在线计算等功能；二是界面友善、语言自然。

MATLAB以复数处理作为计算单元，指令表达与标准教科书的数学表达式相近；三是开放性强。

该公司本身就推出了30多个应用工具箱，而世界上超过200家公司开发出与MATLAB兼容的第三方产品，这些产品向用户提供更多的工具箱、模块集、与其他商业产品的接口等。

2、启动MATLAB 环境，并给出语句tic, A=rand(500); B=inv(A); norm(A*B-eye(500)),toc，试运行该语句，观察得出的结果，并利用help 命令对你不熟悉的语句进行帮助信息查询，逐条给出上述程序段与结果的解释。

>> tic, A=rand(500); B=inv(A); norm(A*B-eye(500)),tocans =1.2333e-012Elapsed time is 0. seconds.tic, A=rand(500);B=inv(A);norm(A*B-eye(500)),toc3、试用符号元素工具箱支持的方式表达多项式5432()34236f x x x x x x =+++++，并令11s x s -=+，将f(x) 替换成s 的函数。

MatLab考试题题库（必做题）（带答案）一，1.请登陆美国 MathWorks 公司的网站 ()，查看看现在大概有多少本 MATLAB-based books （以 MATLAB 为基本软件，来说明各个专业领域的教科书或工具书）。

哪一个领域的 MATLAB-based books 最多？中文书共有几本？答：1612本，数学方面的最多，中文书共有37本。

2.请在 MATLAB中直接输入下列常数，看它们的值是多少：a.ib.jc.epsd.infe.nanf.pig.realmaxh.realmin依次解为：ans = 0 + 1.0000i ans = 0 + 1.0000i ans =2.2204e-016 ans =Inf ans = NaN ans =3.1416 ans =1.7977e+308 ans =2.2251e-3083.试写一函数 regPolygon(n)，其功能为画出一个圆心在 (0, 0)、半径为 1 的圆，并在圆内画出一个内接正 n 边形，其中一顶点位于 (0, 1)。

例如 regPolygon(8) 可以画出如下之正八边型：解：新建regPolygon.m文件如下：function y=regPolyfon(n)n=8;%要画的n边形R=1; %圆的半径t=0:0.01:2*pi;x=R*cos(t);y=R*sin(t);m=linspace(pi/2,5/2*pi,n+1);xz=R*cos(m);yz=R*sin(m);hold onplot(x,y,xz,yz);axis 'equal';4.一条参数式的曲线可由下列方程式表示：x = sin(t), y = 1 - cos(t) + t/10当 t 由 0 变化到 4*pi 时，请写一个 MATLAB 的脚本 plotParam.m，画出此曲线在 XY 平面的轨迹。

解：新建plotParam.m ：t = linspace(0, 4*pi);x = sin(t);y = 1-cos(t)+t/10;plot(x, y, '-o'); -1-0.8-0.6-0.4-0.200.20.40.60.8100.511.522.535. 当一个小圆轮沿着一条曲线行进时，轮缘任一点的轨迹就会产生变化丰富的摆线。

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 已知A=[2 5 6；1 7 8]；B=［1 4 3;4 7 0］；写出下列各指令运行的结果。

A ＆B 的运行结果为［1 1 1;1 1 0] 。

A== B的运行结果为[0 0 0 ;0 1 0] 。

2 产生3阶单位阵的命令为eye（3) ;产生4阶全1方阵的命令为ones(4）。

3 A=rand（5，8）;b=size（A）;c=length(A);则b和c的值分别为[5,8]和8.4 求矩阵的特征值的指令为eig .5符号; 可以使命令行不显示运算结果，% 用来表示该行为注释行。

6在线形图型格式的设置中，字符r表示什么颜色（红色）7fix函数四舍五入对数组［3。

17，5.34，8.90，2.52］取整，结果为［3 5 8 2］。

8 P，Q分别是个多项式的系数矢量，求P对应的多项式的根,使用的命令是roots(P）9使2个plot的图形在同一个坐标显示，使用hold on 命令进行图形保持，使用grid on 命令为图形添加网格。

10 数据拟合采用___polyfit___函数;求不定积分采用___ int ___函数。

11 A=rand(2,5）；b=size（A)；c=length(A）；则b和c的值分别为_2 5和5__ 12已知A是矩阵,求A的对角矩阵函数是___diag(A）____,求A的下三角矩阵函数是__tril（A)____。

13条形图由函数__bar___来实现。

14已知s=‘显示“hey”’，则s的元素个数是9二解答题1、A=［71，3，—8; 2,-9，8; 0，4，5］B=A’问B=（）答案：B=［71，2，0；3，—9,4 ; —8，8,5］2、A=[1,2；3,4］；B[1，0；0,1];请问A*B=？,A.＊B=?答案：A＊B=［1,2;3，4],A.*B=[1,0;0,4]3、A=［57，19，38；-2，31，8；0，84，5］B=rot90（A）B=( )答案：B=［38,8，5；19，31，84；57,-2,0］4、已知:A=［2,2,—1,1；4,3,—1,2; 8，5，-3，4; 3，3，-2,2］R=rank(A）R=（）A=［2，2,3; 4，5，—6；7，8，9］trace（A)=( ）答案：4; 165、A=［0，0,0; 0,1，0; 0,0，1］pinv(A）=（）答案：[0，0，0; 0，1，0；0,0，1］三下列命令执行后，L1，L2，L3，L4的值分别是多少？A=1：9; B=10—A；L1=A～=B;(3分）L2=A>=5; (3分）L3=A>=2＆A〈=6；（4分）L4=find(A〉4&A<=7）; (5分）答案：L1 =1 1 1 1 0 1 1 1 1 L2 =0 0 0 0 1 1 1 1 1 L3 =0 1 1 1 1 1 0 0 0 L4 =5 6 7四写出执行以下代码后B，C，D的值。

MATLAB练习题和答案MATLAB练习题和答案控制系统仿真实验Matlab部分实验结果⽬录实验⼀MATLAB基本操作 (1)实验⼆Matlab编程 (8)实验三Matlab底层图形控制 (10)实验四控制系统古典分析 (23)实验五控制系统现代分析 (28)实验六PID控制器的设计 (35)实验七系统状态空间设计 (40)实验九直流双闭环调速系统仿真 (44)实验⼀ MATLAB 基本操作1 ⽤MATLAB 可以识别的格式输⼊下⾯两个矩阵12332357135732391894A ??=??144367823355422675342189543ii B i +??+?=+?再求出它们的乘积矩阵C ，并将C 矩阵的右下⾓2×3⼦矩阵赋给D 矩阵。

赋值完成后，调⽤相应的命令查看MATLAB ⼯作空间的占⽤情况。

A=[1,2,3,3;2,3,5,7;1,3,5,7;3,2,3,9;1,8,9,4];B=[1+4i,4,3,6,7,8;2,3,3,5,5,4+2i;2,6+7i,5,3,4,2;1,8,9,5,4,3]; C=A*B;D=C(4:5,4:6); whosName Size Bytes Class AttributesA 5x4 160 doubleB 4x6 384 doublecomplexC 5x6 480 double complexD 2x3 96 double complex2 选择合适的步距绘制出下⾯的图形(,)t∈-sin(/)t，其中111t=[-1:0.1:1];y=sin(1./t);plot(t,y)3 对下⾯给出的各个矩阵求取矩阵的⾏列式、秩、特征多项式、范数、特征根、特征向量和逆矩阵。

75350 083341009103150037193......A =??-????，5765710876810957910B ??=12345678910111213141516C=，33245518118575131D --??-?=----A=[7.5,3.5,0,0;8,33,4.1,0;0,9,103,-1.5;0,0,3.7,19.3];B=[5,7,6,5;7,10,8,7;6,8,10,9;5,7,9,10];C=[1:4;5:8;9:12;13:1rtf6];D=[3,-3,-2,4;5,-5,1,8;11,8,5,-7;5,-1,-3,-1];det(A);det(B);det(C);det(D); rank(A); rank(B); rank(C); rank(D); a=poly(A); b=poly(B); c=poly(C); d=poly(D);norm(A);norm(B);norm(C);norm(D);[v,d]=eig(A,'nobalance');[v,d]=eig(B,'nobalance');[v,d]=eig(C,'nobalance');[v,d]=eig(D,'nobalance');m=inv(A);n=inv(B);p=inv(C);q=inv(D);4 求解下⾯的线性代数⽅程，并验证得出的解真正满⾜原⽅程。

MATLAB语言2008年笔试考试答案及评分标准一、填空题(每空1分，共20分)1、help who whos clear close2、A*B=[19,22;43,50] A.*B=[5,12;21,32]3、size(a)=[3 3] sum(a)=[15 27 39]4、plot(x,y) title(‘正弦波’) xlabel(‘时间’) ylabel(‘幅度’)5、while end break input6、0 1 1 signal二、选择题(每题4分，共20分)C AD D B三、（共10分)ch = menu('select a num','1','2','3','4'); %3分switch ch % switch 语句格式写对3分case 1y=10-20*rand(1,20); %1分case 2y=zeros(4,4); %1分case 3y=5*ones(3,4); %1分case 4y=randn(3,3); %1分end四、A=[2 -1 2 6;0 4 3 0;0 6 7 3;3 1 0 5]; %3分b=[8 7 4 2]'; %3分x=A\b; 或x=inv(A)*b; %4分另一种答案:A=[2 0 0 3;-1 4 6 1;2 3 7 0;6 0 3 5]; %3分b=[8 7 4 2]; %3分x=b/A;或x=b*inv(A); %4分五、函数文件：function y=myfunval(x) %2分if x>=-3 & x<-1y=(-x^2-4*x-3)/2; %1分elseif x>=-1 & x<1y=-x^2+1; %1分elsey=(-x^2+4*x-3)/2 %1分end脚本文件：x=[-3:0.01:3]; %1分y=zeros(1,length(x)); %1分for i=1:length(x) %2分y(i)=myfunval(x(i));endplot(x,y); %1分六、%构建结构阵列,5分,也可用赋值语句逐个赋值student=struct('name',{'liming' 'wangwei' 'liujun' 'zhangxiao'},...'eng',{85 83 80 90},'phy',{75 84 69 75},'math',{78 92 77 81});%加权求平均，并排序，7分w=[3 4 5];w=w/sum(w);student_num=length(student);syn_score=zeros(1,student_num);for i=1:student_numsyn_score(i)=student(i).eng*w(1)+student(i).phy*w(2)...+student(i).math*w(3);end[syn_score,index]=sort(syn_score); %从小到大排序%显示结果，3分，用sprintf按格式输出或用字符串连接strcat等均可disp(['排名',' ','姓名',' ','综合测评成绩']);for i=1:student_numdisp([num2str(i),' ',student(index(student_num+1-i)).name,…' ',num2str(syn_score(student_num+1-i))]);end七、f1=2; f3=8; f2=5; fs=20; n=0:255; h=[2.5 7 9 7 2.5]; %3分s1=10*sin(2*pi*f1*n/fs); %1分s2=5*sin(2*pi*f2*n/fs); %1分s3=2*sin(2*pi*f3*n/fs); %1分s=s1+s2+s3; %1分output=conv(s,h); %2分s_f=abs(fft(s)); %1分output_f=abs(fft(output)); %1分figure;subplot(2,2,1); plot(s); %1分subplot(2,2,3); plot(s_f); %1分subplot(2,2,2); plot(output); %1分subplot(2,2,4); plot(output_f); %1分。

Matlab考试题库及答案(教师出卷参考专用)一、选择题1.以下哪个函数用于在Matlab中创建一个图形窗口？A. figureB. plotC. graphD. window答案：A2.在Matlab中，以下哪个选项可以用来定义一个矩阵？A. A = [1 2 3; 4 5 6]B. A = (1, 2, 3, 4, 5, 6)C. A = {1, 2, 3, 4, 5, 6}D. A = 1 2 3; 4 5 6答案：A3.以下哪个函数用于求解线性方程组Ax=b？A. solveB. linsolveC. solve(A, b)D. linsolve(A, b)答案：D4.在Matlab中，如何计算矩阵A和矩阵B的乘积？A. A BB. A \ BC. A . BD. A .\ B答案：A5.以下哪个函数用于在Matlab中绘制三维散点图？A. scatterB. scatter3C. plot3D. bar3答案：B二、填空题1.在Matlab中，要创建一个名为"myfig"的图形窗口，可以使用______函数。

答案：figure('Name', 'myfig')2.在Matlab中，要计算矩阵A的行列式，可以使用______函数。

答案：det(A)3.在Matlab中，若要计算变量x的平方，可以使用______运算符。

答案：.^24.在Matlab中，若要计算矩阵A的逆矩阵，可以使用______函数。

答案：inv(A)5.在Matlab中，要绘制一个正弦波形，可以使用______函数。

答案：plot(sin(x))三、判断题1.在Matlab中，矩阵的索引从1开始计数。

（）答案：正确2.在Matlab中，可以使用逻辑运算符"&&"和"||"。

（）答案：错误3.在Matlab中，矩阵乘法满足交换律。

西安电子科技大学考试时间 90 分钟MATLAB语言试题1.考试形式：闭卷；2.本试卷共七大题，满分100分。

班级学号姓名任课教师一、填空题(每空1分，共20分)1.MATLAB的通用命令中若对函数功能和用法不熟悉，可用命令，列出内存中的变量目录可用或命令，清除工作空间的变量可用命令，关闭图形窗口可用命令。

2.设A=[1,2;3,4];B=[5,6;7,8]，则A*B= ，A.*B= 。

3.A=[1,2,3;10 20 30;4 5 6]，则size(A)= ，sum(A)= 。

4.x=0:pi/20:pi;y=sin(x);以x为横坐标以y为纵坐标的曲线图绘制命令为：，给该图形加上“正弦波”的标题的命令为：，给该图形的横坐标标注为“时间”，纵坐标标注为“幅度”的命令分别为：和。

5.MATLAB中用于控制不定次重复循环语句为，若在循环执行过程中需要终止该循环时采用的语句为，当用户需要从键盘输入参数时采用的语句为。

6.设：str1＝’matlab’, str2＝’matche’, str3＝’Matlab’, str4＝’signalprocessing using MATLAB’;则strcmp(str1,str2)= ，strncmp(str1,str2,3)= ，strcmpi(str1,str3)= ，strtok(str4) = 。

二、选择题(每题4分，共20分)1.在Command Window窗口中分别输入下列命令，对应输出结果正确的是( )A. x=[2:-2] x=[-2 1 0 1 2]B. x=zeros(1,2);x>0 ans=1C. y=diag(eye(3),2).’y=[0 0]D. 5-10*rand(1,2) ans=[-5.0501 1.2311]2.下列关于脚本文件和函数文件的描述中不正确的是( )A. 函数文件可以在命令窗口直接运行;B. 去掉函数文件第一行的定义行可转变成脚本文件;C. 脚本文件可以调用函数文件;D. 函数文件中的第一行必须以function开始;第2页共4页3. 在MATLAB的Editor中输入下列语句：x=[1 2 3 4];y=polyval(x,1);运行后，在Command Window 中对应y的输出结果为( )A.5B.8C.24D.104. 在MATLAB中定义：A=randn(5,4,3,2),则下列关于A的操作中正确的是：（）A y=eig(A) B. y=reshape(A,[4 3 6 7])C. y＝permute(A,[2 3 4 5])D. y=sin(A)5.设str1=’hello’,str2=’help!’,str3=’xidian university’,有：k1=sum(str1==str2), k2=findstr(str3,’an’)，则k1,k2的结果分别为：（）A 3, 4 B.3, 5 C.4,4 D. 4,5三、(10分) 利用menu函数输入选择参数ch，输入菜单的名称为“select an num”，根据ch的取值分别给变量y赋不同的值。

Matlab考试题库+答案(教师出卷参考专用)[1 2；3 4][1 0；0 4]21. A=[1 2 3;6 5 4]; reshape(A,3,2)=______________.[1 5；6 3；2 4]22. A=[1,2;3,1]; B=[1,0;0,1]; A~=B=___________.[0 1；1 0]23. A=[1,2,3;1,4,9;1,7,27]; triu(A,-1)=_ ____.提取矩阵上三角矩阵的函数[1 2 3；1 4 9；0 7 27] 24. A=[4,15,-45,10,6;56,0,17,-45,0]; find(A>=10 & A<=20)=___ __________.3 6 725.已知A=[0 9 6;1 3 0];B=[1 4 3;1 5 0];写出下列各指令运行的结果。

A &B 的运行结果_______；A ./B 的运行结果_______。

ans=[0 1 1;1 1 0]ans= [0 2.25 2;1 0.6 NaN]26. 产生4阶全0方阵的命令为_______；产生3阶全1方阵的命令为。

zeros(4)ones(3)27. A=rand(2,5);b=size(A);c=length(A)；则b和c的值分别为_______ 。

2 5和 528. 用if判断语句判断80≥ x >60,在MATLAB中if语言后的判断应写为_______。

if x<=80 & x>6029.为了使两个pl ot的图形在同一个坐标显示，可以使用_______命令进行图形保持；可以使用_______命令为图形添加网格。

hold ongrid on30. MATLAB的工作空间中有三个变量v1, v2, v3，写出把它们保存到文件my_data.mat中的指令_______；save my_date31. MATLAB提供了两种运行方式，是按_______和把_______文件copy到matlab里，回车。

《MATLAB语言程序设计》练习题1（2010.5）一、填空题1．设有程序A=[1,2,3,4;6,-5,6,3];B=[2,4,5,7];A1=A(:,[1 3 4]);A2=A(1,[2,3]).将下列命令的运行结果填在横线上size(A)= ;A(1,: ).*B= ;length(B)= ;A1= ; A2= ;A+3= ;B.^2=2. A=[1,2,3,4;6,-5,6,3];B=[-2,4,-6,7,9];A1=(A<3);A2=all(A(:,1)<5);A3=any(A(2,:)<10);A4=find(abs(B)<3| abs(B)>7); A5=find(abs(B)<7& abs(B)>3)A1= ; A2= ; A3= ; A4= ;A5= .;).4(;2(3).);BA R(A R 2).A (|;)(|).1(,214154101,321214153111V U A E B ABAX AB D B A ：、TT与特征向量的特征值求解矩阵方程的秩设题试编写程序计算下列各二+=+==⎪⎪⎪⎭⎫⎝⎛--=⎪⎪⎪⎭⎫⎝⎛-=---⎪⎪⎪⎩⎪⎪⎪⎨⎧=+-=++-=++=-+=++⎪⎪⎩⎪⎪⎨⎧=-++-=-++=-++=+++⎩⎨⎧=-++=+++3512332232.3123322352:.2221432:.1321321321321321432143214321432143214321x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x 、超定方程组恰定方程组欠定方程组程组的程序写出求解下列各线性方三.)(,cos 3)sin()sin(cos 2sin )(.3;)(,,1sin )ln(.2;02.1:,432122223变换的计算并计算生成符号矩阵的根求方程下列各题运用符号运算功能计算四Laplace t A wt ewt t t tett t te tt A dx x A dxAdx x e s x x A xx、utttx⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎦⎤⎢⎢⎣⎡+==+---⎰的解函数求解微分方程用1)0(')0(,sin '2".4===++x x t x x x dsolve（1） 若.2c bt aty t y ++=的经验公式为与试编写程序计算出上式中的a 、b 、c;（2） 若.btaey t y =的经验公式为与试编写程序计算出上式中的a 、b;（3） 画出数表中的散点图（红色点）, c bx axy ++=2中拟合曲线图(蓝色实心线),以及.btaey = (绿色点划线_.).（4） 图形标注要求:无网格线,横标注“时间t ”，纵标注“人口数(百万)”，图形标题“美国1900—2000年的人口数据”。