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中用于生成0到1之间均匀分布随机数的函数是： - A. rand- B. randi- C. randperm- D. randn答案：A2. 下列哪个命令可以用于绘制函数图像？- A. plot- B. text- C. title- D. xlabel答案：A3. MATLAB中，以下哪个不是矩阵的属性？- A. size- B. length- C. rank- D. transpose答案：D4. 以下哪个函数可以用于求解线性方程组？- A. solve- B. linsolve- C. linprog- D. fsolve答案：A5. MATLAB中，用于计算矩阵特征值的函数是：- A. eig- B. eign- C. eigen- D. eigenvalue答案：A二、填空题1. MATLAB的基本数据单位是________。
答案：矩阵2. 使用MATLAB进行数值计算时，可以利用________来存储数据。
答案：变量3. MATLAB中的向量可以是________或________。
- 答案：行向量；列向量4. 矩阵的转置可以通过________操作实现。
答案：.'5. MATLAB中，使用________函数可以计算矩阵的行列式。
答案：det三、简答题1. 简述MATLAB中使用循环结构的注意事项。
答案：在MATLAB中使用循环结构时，应注意以下几点：确保循环体内部逻辑正确，避免无限循环；使用for循环时，循环变量的初始化和步长设置要合理；使用while循环时，循环条件要明确，确保循环能够正常退出。
2. 描述MATLAB中如何实现数组的多维索引。
答案：在MATLAB中，多维数组的索引可以通过使用圆括号来实现，每个维度的索引用逗号分隔。
例如，对于一个三维数组A，可以使用A(i,j,k)来访问第i行、第j列、第k层的元素。
matlab简单编程21个题目及答案
1、设⎥⎦⎤⎢⎣⎡++=)1(sin35.0cos2xxxy，把x=0~2π间分为101点，画出以x为横坐标，y为纵坐标的曲线。
第一题的matlab源程序：①考虑cos（x）为一个整体，然后乘以中括号里面的全部x=0:2*pi/100:2*pi; %x的步长以及范围从0到2*pi y=cos(x).*(0.5+3*sin(x)./(1+x.^2)); %y的表达式plot(x,y)%画出图形图如下：②考虑对整体求解cos，先求x乘以括号中的部分x=0:2*pi/100:2*pi; %x的步长以及范围从0到2*pi y=cos(x.*(0.5+3*sin(x)./(1+x.^2))); %y的表达式plot(x,y) %画出图形图如下：2、产生8×6阶的正态分布随机数矩阵R1, 求其各列的平均值和均方差。
并求该矩阵全体数的平均值和均方差。
第二题的matlab源程序如下：R1=randn(8,6) %产生正态分布随机矩阵R1 =1.0933 -0.7697 1.5442 -0.1924 1.4193 0.21571.1093 0.3714 0.0859 0.8886 0.2916 -1.1658-0.8637 -0.2256 -1.4916 -0.7648 0.1978 -1.14800.0774 1.1174 -0.7423 -1.4023 1.5877 0.1049-1.2141 -1.0891 -1.0616 -1.4224 -0.8045 0.7223-1.1135 0.0326 2.3505 0.4882 0.6966 2.5855-0.0068 0.5525 -0.6156 -0.1774 0.8351 -0.66691.5326 1.1006 0.7481 -0.1961 -0.2437 0.1873aver=(sum(R1(1:end,1:end)))./8 %产生各行的平均值aver =0.0768 0.1363 0.1022 -0.3473 0.4975 0.1044a=std(R1(1:end,1:end)) %产生各行的均方差也就是标准差a =1.0819 0.8093 1.3456 0.8233 0.8079 1.2150aver1=(sum(R1(:)))./48 %全体数的平均值aver1 =0.0950b=std(R1(:)) %全体数的均方差即标准差b =1.01033、设x=rcost+3t,y=rsint+3,分别令r=2,3,4，画出参数t=0~10区间生成的x~y 曲线。
Matlab 考题题整理 带答案
MATLAB 考试试题(1)产生一个1x10的随机矩阵，大小位于（-5 5），并且按照从大到小的顺序排列好！（注：要程序和运行结果的截屏）答案：a=10*rand(1,10)-5;b=sort(a,'descend')1.请产生一个100*5的矩阵，矩阵的每一行都是[1 2 3 4 5]2. 已知变量：A=’ilovematlab’；B=’matlab’, 请找出：（A）B在A中的位置。
（B）把B放在A后面，形成C=‘ilovematlabmatlab’3. 请修改下面的程序，让他们没有for循环语句！A=[1 2 3; 4 5 6; 7 8 9];[r c]=size(A);for i=1:1:rfor j=1:1:cif (A(i,j)>8 | A(i,j)<2)A(i,j)=0;endendend4. 请把变量A=[1 2 3; 4 5 6; 7 8 9]写到文件里(output.xls)，写完后文件看起来是这样的1 2 3 4 5 6 7 8 95.试从Yahoo网站上获得微软公司股票的2008年9月的每日收盘价。
6.编写M文件，从Yahoo网站批量读取60000.SH至.SH在2008年9月份的每日收盘价（提示：使用字符串函数）。
7. 将金牛股份（）2005年12月14日至2006年1月10日的交易记录保存到Excel中，编写程序将数据读入MATLAB中，进一步将数据读入Access数据库文件。
8.已知资产每日回报率为0.0025，标准差为0.0208，资产现在价值为0.8亿，求5％水平下资产的10天在险价值(Var)。
9.a=[1 2 3 4 5],b=a(1)*a(5)+a(2)*a(4)+a(3)*a(3)+a(4)*a(2)+a(5)*a(1).试用MATLAB中最简单的方法计算b，注意最简单哦。
1、求下列联立方程的解3x+4y-7z-12w=45x-7y+4z+ 2w=-3x +8z- 5w=9-6x+5y-2z+10w=-8求系数矩阵的秩；求出方程组的解。
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的矩阵。
matlab20道试题及解答
试题1．“数学黑洞”：任意一个4位自然数，将组成该数的各位数字重新排列，形成一个最大数和一个最小数，之后两数相减，其差仍为一个自然数。
重复进行上述运算，最终会出现一个神秘的数，请编程输出这个神秘的数。
clear;a=input('请输入一个四位正整数：');str_a=num2str(a); %将a转化为一个字符串b_min=str2double(sort(str_a)); %形成最小数b_max=str2double(sort(str_a,'descend')); %形成最大数b=b_max-b_min; %求最大数与最小数之差while (b~=a)a=b;str_a=num2str(a); %将a转化为一个字符串b_min=str2double(sort(str_a)); %形成最小数b_max=str2double(sort(str_a,'descend')); %形成最大数b=b_max-b_min; %求最大数与最小数之差endb试题2．将数字1、2、3、4、5、6填入一个2行3列的表格中，要使得每一列右边的数字比左边的数字大，每一行下面的数字比上面的数字大。
请编写程序求出按此要求可有几种填写方法。
a(1)=1;a(6)=6;count=0; %用来计数b=perms('2345'); %产生2345的全排列[m,n]=size(b);for i=1:mtemp=b(i,:);a(2)=str2double(temp(1));a(3)=str2double(temp(2));a(4)=str2double(temp(3));a(5)=str2double(temp(4));if ((a(4)>a(2))&&(a(4)>a(3))&&(a(5)>a(3)))count=count+1;c=reshape(a,2,3); %将a向量转化为2*3矩阵输出disp(c);endenddisp(['共有',num2str(count),'种填写方法']); %输出填写方法的种数试题3．编写成绩排序程序。
matlab期末考试题目及答案
matlab期末考试题目及答案1. 题目：编写一个MATLAB函数，实现矩阵的转置操作。
答案：可以使用`transpose`函数或`.'`操作符来实现矩阵的转置。
例如，对于一个矩阵`A`，其转置可以通过`A'`或`transpose(A)`来获得。
2. 题目：使用MATLAB求解线性方程组Ax=b，其中A是一个3x3的矩阵，b是一个3x1的向量。
答案：可以使用MATLAB内置的`\`操作符来求解线性方程组。
例如，如果`A`和`b`已经定义，求解方程组的代码为`x = A\b`。
3. 题目：编写MATLAB代码，计算并绘制函数f(x) = sin(x)在区间[0, 2π]上的图像。
答案：首先定义x的范围，然后计算对应的函数值，并使用`plot`函数绘制图像。
代码示例如下：```matlabx = linspace(0, 2*pi, 100); % 定义x的范围y = sin(x); % 计算函数值plot(x, y); % 绘制图像xlabel('x'); % x轴标签ylabel('sin(x)'); % y轴标签title('Plot of sin(x)'); % 图像标题```4. 题目：使用MATLAB编写一个脚本，实现对一个给定的二维数组进行排序，并输出排序后的结果。
答案：可以使用`sort`函数对数组进行排序。
如果需要对整个数组进行排序，可以使用`sort`函数的两个输出参数来获取排序后的索引和值。
代码示例如下：```matlabA = [3, 1, 4; 1, 5, 9; 2, 6, 5]; % 给定的二维数组[sortedValues, sortedIndices] = sort(A(:)); % 对数组进行排序sortedMatrix = reshape(sortedValues, size(A)); % 将排序后的值重新构造成矩阵disp(sortedMatrix); % 显示排序后的结果```5. 题目：编写MATLAB代码，实现对一个字符串进行加密，加密规则为将每个字符的ASCII码值增加3。
Matlab编程习题
1.求方程05.01)1ln(22=---+-x x x x x 的正根。
解：syms x;2.求满足∑=+m n n 0)1ln(>100的最小m 第2题解；y=0;n=0;while y<100n=n+1;y=y+log(n);endm=n-13.已知函数f (x )=x 4-2x 在(-2,2)内有两个根。
取步长h =0.05,通过计算函数值求得函数的最小值点和两个根的近似解。
（提示：求近似根等价于求函数绝对值的最小值点）解；x=-2:0.05:2;y=x.^4-2.^x;ymin=min(y);[imin,jmin]=find(y==ymin);xmin=x(imin,jmin)ymin=y(imin,jmin)x1=-2:0.05:x(imin,jmin);y1=abs(x1.^4-2.^x1);y1min=min(y1);[i1min,j1min]=find(y1==y1min);x1min=x1(i1min,j1min)y1min=y1(i1min,j1min)x2=x(imin,jmin):0.05:2;y2=abs(x2.^4-2.^x2);y2min=min(y2);%[i2min,j2min]=find(y2==y2min);%x2min=x2(i2min,j2min)%y2min=y2(i2min,j2min)%第4题%[x,y]=solve('(x-2)^2+(y+2*x-3)^2=5,18*(x-3)^2+y^2=36'4.(椭园的交点)两个椭圆可能具有0～4个交点，求下列两个椭园的所有交点坐标(x -2)2+(y -3+2x )2=52(x -3)2+(y /3)2=4[x,y]=solve('(x-2)^2+(y+2*x-3)^2=5,18*(x-3)^2+y^2=36')%syms x y real5.考虑函数f (x,y )=y 3/9+3x 2y +9x 2+y 2+xy +9(1)作出f (x,y )在-2<x <1,-7<y <1的图，观察极值点的位置；(2)用MATLAB 函数fminsearch 求极值点和极值。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
作业一1输出x,y两个中值较大的一个值x=input(‘x’);y=input(‘y’);if x>yxelseyend2输入x，计算y的值。
计算函数的值y=x+1,x<0,y=2x-1,x≧0x=input错误！未指定书签。
(‘x);if x<0y=x+1elsey=2*x-1end3输入一学生成绩，评定其等级，方法是：90~100分为“优秀”，80~89分为“良好”，70~79分为“中等”，60~69分为“及格”，60分为“不合格”x=input(‘x’)if x>100|x<0y=’输入错误’elseif x>=90y=’优秀’elseif x>=80y=’良好’elseif x>=70y=’中等’elseif x>=60y=’及格’elsey=’不合格’emd4某超市节日期间举办购物打折的促销活动，优惠办法是：每位顾客当天一次性购物在100元以上者，按九五折优惠；在200元以上者，按九折优惠；在300元以上者，按八五折优惠；在500元以上者，按八折优惠。
x=input(‘x’);if x>=500y=x*0.8elseif x>=300y=x*0.85elseif x>=200y=x*0.9elseif x>=100y=x*0.95else y=xend 5编程计算：s=1+2+3+…+100sum=0;for i=1:100sum=sum+i;endsum引申1!-2!+3!-4!+5!- (99)sum=0;for i=1:99pdr=1;for k=1:ipdr=pdr*k;endsum=sum+pdr*(-1)^(i-1);endsum引申1*2*3*4*……*100sum=1；for i=1:100sum=sum*iendsum6计算1~100的奇数和sum=0;for i=1:2:100sum=sum+i;endsum7百元买百鸡问题。
假定小鸡每只5角，公鸡每只2元，母鸡每只3元。
编程列出所有可能的够鸡方案for x=1:33for y=1:50for z=1:200if x+y+z==100&3*x+2*y+0.5*z==100 [x,y,z]endendendend8我国有13亿人口，按人口年增长0.8%计算，多少年后我国人口超过26亿n=1;m=13;while m<26m=m*(1+0.008);n=n+1;endn9求水仙花数，水仙花数是指一个n位数（n>=3），他的每个位上的数字的n次幂之和等于它本身。
求出1000以内for i=100:999a(1)=fix(i/100);a(2)=fix((i-a(1)*100)/10);a(3)=i-a(1)*100-a(2)*10;ifa(1)^3+a(2)^3+a(3)^3==iiendend10鸡和兔子关在一个笼子里，已知共有头36个，脚100个，求笼内关了多少只兔子和多少只鸡？for i=1:50for j=1:36ifi+j==36&i*2+j*4==100[i,j]endendend作业二1.用matlab编程求算式xyz+yzz=532中x，y，z的值（其中xyz和yzz分别表示一个三位数）。
clearfor x=1:9for y=1:9for z=1:9i=100*x+10*y+z+100*y+10*z+z;if i==532[x,y,z]endendendend2.用matlab编程完成用一元人民币换成一分、两分、五分的所有兑换方案（即输出所有的组合方式）。
clearfor i=0:20for j=0:50k=100-i*5-j*2;if(k>=0)[I,j,k]endendend3.有1020个西瓜，第一天卖一半多两个，以后每天卖剩下的一半多两个，问几天后可以卖完，请用matlab编程计算。
day=0;x1=1020while(x1)x2=x1/2-2;x1=x2;day=day+1;endday4.有一堆零件（100到200之间），如果分成4个零件一组的，则多2个零件；若分成7个零件一组的，则多3个零件；若分成9个零件一组的，则多5个零件。
用matlab 编程求这堆零件总数。
for x=100:200if mod(x,4)==2if mod(x,7)==3if mod(x,9)==5xendendendend5.编写程序，求1000至9999之间的回文数．回文数是指正读与反读都一样的数，如1221.for i=1000:1000:9999for j=0:9n=i+j*100+j*10=i/1000endend6.作函数y=sinx,x∈[0,2∏]的图象，用蓝色五角星表示点。
x=linspace(0,2*pi);y=sin(x);plot(x,y,'bp');作业三1:编写M函数3222 6.30.05 3.14x x xfx x-+-=+-,计算)3()2()1(2fff+function y=fun1(x)y=(x^3-2*x^2+x-6.3)/(x^2+0.05*x-3.14) ;fun1(1)*fun1(2)+fun1(3)*fun1(3)2：编写M函数1,23,28()45,820cos sin,20x xx xf xx xx x x+<⎧⎪≤≤⎪=⎨-<≤⎪⎪+>⎩，计算(0.1)f，(1)f，(9)f，(22)f，(2)fπfunction y=fff(x)if x<2y=x+1;elseif x>=2&x<=8y=3*x;elseif x>8&x<=20y=4*x-5;elsey=cos(x)+sin(x);end3：根据pi*pi/6=1/1^2+1/2^2+1/3^2+……+1/n^2,求pi 的近似值。
当n 分别取100，1000，10000时，结果是多少？y=0;n=100;for i=1:n y=y+1/i/i;end ypi=sqrt(6*y)y=0;n=1000;for i=1:n y=y+1/i/i;end ypi=sqrt(6*y)y=0;n=10000;for i=1:n y=y+1/i/i;end ypi=sqrt(6*y)4:编程:求[100，200]之间第一个能被21整除的整数。
for n=100:200if mod(i,21)==0continue break end end5：编写函数文件求半径为r 的圆的面积和周长。
function [s,p]=fcircle(r)s=pi*r*r;p=2*pi*r;6：根据y=1+1-n 215131+⋯++,求：（1）y<3时的最大n 值。
（2）与（1）的n 值对应的y 值。
y=1;n=1;while(y<3)n=n+2;y=y+1/n;endy=y-1/n n=(n+1)/27:已知⎪⎪⎩⎪⎪⎨⎧>+-=======---3,23,12,01,1321321n n n n fff f f f f n n n n求ff 1001~中：最大值、最小值、各数之和。
f(1)=1f(i)=0f(3)=1for i=4:100f(i)=f(i-1)+2*f(i-2)+f(i-3)end min(f)max(f)sum(f)作业四1.编程计算3+33+333+3333+…..前10项之和。
m=3;n=0;s=0;for i=1:10n=n*10+m;s=s+n;end s2.编程求1-1/2+1/3-1/4+.......+1/99-1/100的值m=0;for i=1:1:100o=(1/i)*(-1)^(i+1);m=m+o;end m3.编写程序，求1-3+5-7+…-99+101的值。
s=0;j=1;for i=1:2:101s=s+(-1)^(j-1)*i;j=j+1;end s4.编写程序，输出从公元1000年至2000年所有闰年的年号。
判断公元年是否为闰年的条件是：公元年数如能被4整除，而不能被100整除，则是闰年；公元年数能被400整除也是闰年。
for year=1000:2000if(mod(year,4)==0&mod(year,100)~=0)|m od(year,400)==0yearendend5.将10个整数输入到一个数组中，然后再将其按输入顺序的逆序进行排列并输出。
for i=1:10a(i)=input('x=');endfor i=10:-1:1a(i)end6.一球从100米高度自由落下，每次落地后反跳回原高度的一半，再落下。
求它在第10次落地时，共经过多少m？第10次反弹多高？sn=100;hn=sn/2;for n=2:10sn=sn+2*hn;hn=hn/2;endsnhn7.猴子吃桃问题。
猴子第一天摘下若干个桃子，当即吃了一半，还不过瘾，又多吃了一个。
第二天早上又将剩下的桃子吃掉一半，又多吃了一个。
以后每天早上都吃了前一天剩下的一半零一个。
到第10天早上再想吃时，见只剩一个桃子了。
求第一天共摘多少桃子。
day=9;x=1;while(day>0)x=(x+1)*2;day=day-1;endx作业五上机题1.给定一个矩阵a=123456编写一个M函数，要求输入是a，输出各元素之和。
function y=f(a)[r,q]=size(a);y=0;for i=1:rfor j=1:qy=y+a(i,j);endend2.A是一个m×n的矩阵.写一M函数,算出A中有多少个零元素。
function y=f(a)[m,n]=size(a);y=0;for i=1:mfor j=1:nif a(i,j)==0y=y+1;endendend3.矩阵乘法运算要求两矩阵的维数相容，否则会出错。
先求两矩阵的乘积，若出错，则自动转去求两矩阵的点乘。
A=input(‘A’);B=input(‘B’);tryC=A*B;catchC=A.*B;endC4.一个三位整数各位数字的立方和等于该数本身则称该数为水仙花数。
输出全部水仙花数。
for i=100:999a(1)=fix(i/100);a(2)=fix((i-a(1)*100)/10);a(3)=i-a(1)*100-a(2)*10;if a(1)^3+a(2)^3+a(3)^3==iiendend5.fibonccci数组的元素满足fibonacci规则：ak+2=ak+ak+1,(k=1,2,…);且a1=a2=1。
现要求求出该数组中第一个大于10000的元素。