matlab简单编程21个题目及答案
MATLAB40道题与答案修正版
操作方法：首先点击File —New —M-File ，输入相应代码，然后点击运行按钮即可。
P1．已知两序列[]0.8{[][5]}nx n u n u n =--，[]{1,1,1,1,1}h n =计算两序列的卷积并绘制其波形。
%function a1 for i=1:5x(i)=0.8^(i-1); endh=[1 1 1 1 1]; y=conv(x,h) m=0:8;stem(m,y,'filled')P2。
已知复指数序列(1.52)[] 1.2j nx n e +=，绘制20点该序列的实部和虚部。
n=0:19;x=1.2*exp((1.5+j*2)*n); subplot(211); stem(n,real(x)); ylabel('real part'); xlabel('n'); legend('x[n]'); subplot(212); stem(n,imag(x)); ylabel('imag part'); xlabel('n'); legend('x[n]');P3．编写长度为5的中值滤波器程序。
原始未受干扰的序列为：s[n]=3[n(0.5) ]n ，加性噪声信号d[n]为随机序列，幅度0.4，分别绘制长度为40的受干扰序列，以及中值滤波器的输出。
% Program 2_5% Illustration of Median Filtering N = 5; R = 40;a = 0.8*rand(1,R)-0.4;b = round(a); % Generate impulse noise m = 0:R-1;s =3*m.*(0.5.^m); % Generate signalx = s + b; % Impulse noise corrupted signal y = medfilt1(x,N); % Median filtering subplot(2,1,1) stem(m,x);xlabel('n');ylabel('Amplitude');title('Impulse Noise Corrupted Signal'); subplot(2,1,2) stem(m,y);xlabel('n');ylabel('Amplitude'); title('Output of Median Filter'); P4. 已知序列x1[n]={2.2,3,1.5,4.2,1.8}, x2[n]= {0.8,1,1.6,0.8}，x[n]=x1[n]∗x2[n] （卷积），分别绘制序列x1[n] ,x2[n]和x [n]的波形。
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 实现方法（程序段）和运行结果。
1．求出满足∑=>mi i 1100000的最小m 值。
clc clear all sum=0; i=0;while sum<100000 i=i+1; sum=sum+i; end sum i Result ： sum =100128 i = 4472. 分别用循环结构和向量化方式求解级数求和问题i i i 21110012+∑=。
循环结构 clc clear allfor i=1:100 j=1/i^2+1/2^i; s=s+j; end 向量化： clc clear all a=zeros(1,100); for i=1:100a(i)=1/i^2+1/2^i; end s=sum(a) 结果s =2.63503. 试用函数plot 和函数fplot 绘制)tan(sin )sin(tan x x y -=在[]ππ,-∈x 区间内的曲线。
Plot ： clc clear all pi=3.14; for i=1:2*pi*100 u(i)=-pi+0.01*(i-1);y(i)=sin(tan(u(i)))-tan(sin(u(i))); end plot(u,y); Fplot ： clc clear all pi=3.14;y='sin(tan(x))-tan(sin(x))';-4-3-2-101234-3-2-11234. 若函数)cos(5)sin(),(2222y x e y y x x y x f y x ++++=+，试求其偏导数),(y x f y '，并绘制它们的曲线。
clear all clc syms x y pi=3.14;f=x^2*sin(x+y^2)+y^2*exp((x+y)+5*cos(x^2+y)); df=diff(f,y) 结果 df =2*x^2*cos(x+y^2)*y+2*y*exp(x+y+5*cos(x^2+y))+y^2*(1-5*sin(x^2+y))*exp(x+y+5*cos(x^2+y))x=0:pi/50:pi;y=x;f=x.^2.*sin(x+y .^2)+y .^2.*exp((x+y)+5.*cos(x.^2+y));df=2*x.^2.*cos(x+y .^2).*y+2*y .*exp(x+y+5*cos(x.^2+y))+y .^2.*(1-5.*sin(x.^2+y)).*exp(x+subplot(1,2,1) plot3(f,x,y) grid on title('f') subplot(1,2,2) plot3(df,x,y) grid on title('df')-5510x 10502401234-112x 106024012345.假设某概率密度函数有以下分段函数表示：⎪⎩⎪⎨⎧-≤+≤+<->+=+-------1,5457.011,7575.01,5457.0),(215.175.575.0216215.175.575.02112122212212122x x e x x e x x e x x p x x x x x x x x ，试以三维曲面的形式来表示这一函数，并在同一窗口显示该三维曲面的三视图。
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简单编程21个题目及答案复习课程
1、 设 ⎥⎦⎤⎢⎣⎡++=)1(sin 35.0cos 2x x x y ，把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、一个三位整数各位数字的立方和等于该数本身，则称该数为水仙花数。
试编写程序找到并输出全部水仙花数。
2、若一个数等于它的各个真因子之和，则称该数为完数，如6=1+2+3，所以6 是完数。
求[1 500]之间的全部完数。
3、写一个函数pifun 计算下列级数：F(n)=4*(1-1/3+1/5-1/7+...)并使用tic 和toc 指令来测量pifun(1000000)的计算时间。
4、输入一个百分制成绩，要求输出成绩等级为A~E，即90~100 为A，80~89 为B，70~79 为C，60~69 为D，60 以下为E。
5、输入x,y 的值，并将它们的值互换后输出。
6、某商场对顾客所购买的商品实行打折销售，标准如下(商品价格用price 来表示)：price<200 没有折扣；200≤price<500 3%折扣；500≤price<1000 5%折扣；1000≤price<2500 8%折扣；2500≤price<5000 10%折扣；5000≤price 14%折扣输入所售商品的价格，求其实际销售价格。
7、由指令A=rand(3,5)生成二维数组A，试求该数组中所有大于0.5 的元素的位置。
8、写一程序求1！+2！+ (10)9、编写一个函数，找出矩阵A 中最大元素和最小元素所在的行列号。
10、从键盘输入的10 个整数中，找出第一个能被7 整除的数。
若找到，打印此数后退出循环；若未找到，打印“not exist”。
function shuixian(n) %n为用户输入的一个整数n=3;W=zeros(1,n);ss=0;for m=10^(n-1):10^n-1for i=1:nW(i)=fix(rem(m,10^(n+1-i))/10^(n-i));endW=W.^n;for j=1:nss=ss+W(j);endif m==ssdisp(m);endss=0;W=zeros(1,n);end2.for m=1:500sum=0;for n=1:m/2if rem(m,n)==0sum=sum+n;endendif sum==mmendend3.n=input('请输入n值：');tic;f=0;for i=1:nf=f+4*((-1)^(i+1))/(2*i-1) ;endftoc4.1m=input('请输入一个百分制成绩：');m=ceil(m/10);switch mcase {10,9}disp('A');case 8disp('B');case 7disp('C');case 6disp('D');case {5,4,3,2,1,0}disp('E');otherwisedisp('输入有误！')end4.2a=input('请输入一个百分制成绩：');if a>=90 & a<=100disp('A');else if a>=80 & a<90disp('B'); else if a>=70 & a<80 disp('C'); else if a>=60 & a<70 disp('D'); else if a>=0 & a<60 disp('E'); elsedisp('输入数据有误，请核对......');endendendendend5．x = input('x=');y = input('y=');disp ([x,y]);t=x;x=y;y=t;disp ([x,y]);x=input('price:')if x<200f=x;else if x<500f=x*0.97;else if x<1000f=x*0.95;else if x<2500f=x*0.92;else if x<5000 f=x*0.90;elsef=x*0.86;endendendendendf7.A=rand(3,5)for i=1:3for j=1:5if A(i,j)>0.5Ai=i;Aj=j;Aij=[i,j]endendendsum=0;for m=1:10i=1;for n=1:mi=i*n;endsum=sum+i;endsum10.a=zeros(1,10);for i=1:10a(i)=input('please enter number'); endfor i=1:10if (rem(a(i),7)==0)a(i)breakendendif (i==0&&rem(a(i),7)~=0)disp('not exist')end。
matlab习题及答案
matlab习题及答案《Matlab习题及答案：提升编程技能，解决实际问题》Matlab是一种强大的数值计算和数据分析工具，它被广泛应用于工程、科学和其他领域的计算和模拟。
为了帮助大家提升编程技能，解决实际问题，我们为大家准备了一些Matlab习题及答案，希望能够帮助大家更好地掌握Matlab的使用。
1. 习题一：编写一个Matlab程序，计算斐波那契数列的前20个数字，并将结果打印出来。
答案：下面是一个简单的Matlab程序，用于计算斐波那契数列的前20个数字。
```matlaba = 0;b = 1;fib = zeros(1, 20);fib(1) = a;fib(2) = b;for i = 3:20fib(i) = fib(i-1) + fib(i-2);enddisp(fib);```2. 习题二：编写一个Matlab程序，求解一个二次方程ax^2 + bx + c = 0的根。
答案：下面是一个简单的Matlab程序，用于求解二次方程的根。
```matlaba = 1;b = -3;c = 2;delta = b^2 - 4*a*c;if delta > 0x1 = (-b + sqrt(delta))/(2*a);x2 = (-b - sqrt(delta))/(2*a);disp(['The roots are ', num2str(x1), ' and ', num2str(x2)]);elseif delta == 0x = -b/(2*a);disp(['The root is ', num2str(x)]);elsedisp('The equation has no real roots');end```通过以上两个习题及答案的示例，我们可以看到Matlab的强大功能和灵活性。
通过练习这些习题，我们可以更好地掌握Matlab的基本语法和常用函数，从而在实际问题中更快更准确地解决数值计算和数据分析的挑战。
matlab习题与答案
matlab习题与答案
MATLAB习题与答案
MATLAB是一种强大的数学软件，广泛应用于工程、科学和金融等领域。
通过MATLAB，用户可以进行数据分析、图像处理、模拟建模等多种操作。
为了帮
助大家更好地掌握MATLAB的应用，我们为大家准备了一些习题与答案，希望
能够帮助大家更好地理解和掌握MATLAB的使用。
习题一：编写一个MATLAB程序，实现对给定矩阵的转置操作。
解答：可以使用MATLAB中的transpose函数来实现矩阵的转置操作。
例如，
对于一个3x3的矩阵A，可以使用以下代码实现转置操作：
A = [1 2 3; 4 5 6; 7 8 9];
B = transpose(A);
习题二：编写一个MATLAB程序，实现对给定矩阵的逆矩阵计算。
解答：可以使用MATLAB中的inv函数来实现对矩阵的逆矩阵计算。
例如，对
于一个3x3的矩阵A，可以使用以下代码实现逆矩阵计算：
A = [1 2 3; 4 5 6; 7 8 9];
B = inv(A);
习题三：编写一个MATLAB程序，实现对给定矩阵的特征值和特征向量计算。
解答：可以使用MATLAB中的eig函数来实现对矩阵的特征值和特征向量计算。
例如，对于一个3x3的矩阵A，可以使用以下代码实现特征值和特征向量计算：A = [1 2 3; 4 5 6; 7 8 9];
[V, D] = eig(A);
通过以上习题与答案的学习，相信大家对MATLAB的应用已经有了更深入的了
解。
希望大家能够多加练习，不断提升自己的MATLAB技能，为今后的工作和学习打下坚实的基础。
matlab练习题及答案详解(附程序)
Matlab练习题及答案详解一、基础题1、下表为1)2) 把任意两点距离的表中，自己到自己的距离，转变成无穷大clcclearclose alldata=xlsread('ti1.xls');%从excel里导入数据x=data(1,:);%将data里第一行的所有列赋值给xy=data(2,:);%将data里第二行的所有列赋值给xn=length(x);%表x的长度dis=zeros(n,n);%对距离表进行初始化for i=1:nfor j=1:ndis(i,j)=sqrt((x(i)-x(j))^2+(y(i)-y(j))^2);%平面坐标距离公式endenddis(dis==0)=inf;%将距离表中为0的数值，变成无穷大redis=dis;2、按要求步骤，操作下表要求：1) 对每列进行[0,1]区间化'(min())/(max min())=--；x x x x2) 把区间化后的表的右边3列，连接在左边3列下边，构成16行3列的表；3) 把16行3列进行转置，得到3行16列的表。
然后对列进行操作，把3行16列的表变成一列，后转置成一行，得到最终的结果。
clcclearclose allx=load('ti2.txt');%从txt里导入数据big=max(x);small=min(x);%算出每列的最大最小值[m,n]=size(x);%表x的行和列的大小newx1=zeros(m,n);%对newx1进行初始化，大小和初始x一样大for i=1:nnewx1(:,i)=(x(:,i)-small(i))./(big(i)-small(i));endnewx2=zeros(2*m,n/2);%对最新的表进行初始化大小为，行为newx1的2倍，列为一半newx2(1:m,1:n/2)=newx1(1:m,1:n/2);%把左边3列不操作的值，直接赋值给newx2上边newx2(m+1:end,1:n/2)=newx1(1:m,n/2+1:end);%把右边的3列，赋值给newx2下边[p,q]=size(newx2);%表newx2的行和列的大小 newx3=zeros(p*q,1);%对newx3进行初始化 for j=1:qnewx3((j-1)*p+1:j*p,1)=newx2(:,q);%对newx3进行操作 end% newx3=reshape(newx2,[],1); renewx=newx3';%最终的结果3、数值运算1)用matlab 求下列式子极限(a) lim x → (b) 0x →2)用matlab 求下列函数的导数(a)22sin (1)y x =+ (b) arccos()(,)x f x y ye =3)用matlab 求下列式子的定积分(a) 4⎰(b) 312-⎰⎰clc clear %===求极限 %====3.1.a=== syms xy=(x^(1/2)-8)/(x^(1/3)-4);re1a=limit(y,x,64)%====3.1.b===syms xy=((x^2-2*x+4)^(1/2)-2)/x;re1b=limit(y,x,0)%===求导数%====3.2.a===syms xf=(sin(x^2+1))^2;re2a=diff(f,x)%对x求导%====3.2.b===syms x yf=y*exp(acos(x));re2bx=diff(f,x)%对x求偏导re2by=diff(f,y)%对y求偏导%===求定积分===%====3.3.a===syms xre3a=int((x+2)/(2*x+1)^(1/2),x,0,4)%一重积分%====3.3.b===syms x yre3b=int(int(y/(x^2-1),x,-2,-sqrt(2)),y,1,3)%二重积分二、提高题1、画一个正五边形，边长自定。
matlab试题及答案开卷
matlab试题及答案开卷1. MATLAB基础操作在MATLAB中，如何创建一个名为“myMatrix”的3x3矩阵，其元素从1到9按顺序排列？答案：在MATLAB中，可以使用以下命令创建名为“myMatrix”的3x3矩阵：```matlabmyMatrix = [1 2 3; 4 5 6; 7 8 9];```或者使用内置函数`magic`：```matlabmyMatrix = magic(3);```2. 矩阵运算给定两个矩阵A和B，其中A = [2 3; 4 5]，B = [6 7; 8 9]，计算A+B的结果。
答案：在MATLAB中，可以使用加号`+`来计算两个矩阵的和：```matlabA = [2 3; 4 5];B = [6 7; 8 9];C = A + B;```计算结果C为：```matlabC = [8 10; 12 14];```3. 向量操作创建一个从0到1的等差数列，步长为0.1。
答案：在MATLAB中，可以使用`linspace`函数或`:`操作符来创建等差数列：```matlabv = 0:0.1:1;```或者使用`linspace`函数：```matlabv = linspace(0, 1, 11);```两种方法都将得到一个包含11个元素的向量，从0开始，到1结束，步长为0.1。
4. 条件语句编写一个MATLAB脚本，判断一个给定的数n是否为素数。
答案：可以通过以下MATLAB脚本来判断一个数n是否为素数：```matlabfunction isPrime = isPrimeNumber(n)if n <= 1isPrime = false;elsefor i = 2:sqrt(n)if mod(n, i) == 0isPrime = false;return;endendisPrime = true;endend```调用此函数时，传入一个数值n，函数将返回一个布尔值，表示n 是否为素数。
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习题及答案 (入门)
习题1， 计算⎥⎦⎤⎢⎣⎡=572396a 与⎥⎦⎤⎢⎣⎡=864142b 的数组乘积。
a.*bans =12 36 3 8 42 402， 对于B AX =，如果⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=753467294A ，⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=282637B ，求解X 。
inv(a)*bans =-0.5118 4.0427 1.33183， 已知：⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=987654321a ，分别计算a 的数组平方和矩阵平方，并观察其结果。
a.*aans =1 4 9 16 25 36 49 64 81>> a^2ans =30 36 42 66 81 96 102 126 1504， 角度[]604530=x ，求x 的正弦、余弦、正切和余切。
(应用sin,cos,tan.cot) >> sin(x./pi)ans =-0.1242 0.9826 0.2465>> cos(x./pi)ans =-0.9923 -0.1857 0.9692>> tan(x./pi)ans =0.1252 -5.2915 0.2543>> cot(x./pi)ans =7.9894 -0.1890 3.93215， 将矩阵⎥⎦⎤⎢⎣⎡=7524a 、⎥⎦⎤⎢⎣⎡=3817b 和⎥⎦⎤⎢⎣⎡=2695c 组合成两个新矩阵： （1）组合成一个4⨯3的矩阵，第一列为按列顺序排列的a 矩阵元素，第二列为按列顺序排列的b 矩阵元素，第三列为按列顺序排列的c 矩阵元素，即 ⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡237912685574 [reshape(a,4,1),reshape(b,4,1),reshape(c,4,1)]ans =4 75 5 86 2 1 97 3 2（2）按照a 、b 、c 的列顺序组合成一个行矢量，即[]272548716953[reshape(a,1,4),reshape(b,1,4),reshape(c,1,4)]ans =4 5 2 7 7 8 1 3 5 6 9 26，将(x-6)(x-3)(x-8)展开为系数多项式的形式。
Matlab考试题库及答案(教师出卷参考专用)
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中，矩阵乘法满足交换律。
MATLAB考试试题及答案
MATLAB考试试题及答案一、选择题（每题5分，共25分）1. 在MATLAB中，下列哪个命令用于创建一个行向量？A. v = [1; 2; 3]B. v = [1 2 3]C. v = [1, 2, 3]D. v = (1, 2, 3)答案：B2. 在MATLAB中，下列哪个命令用于计算矩阵A的行列式？A. det(A)B. det(A')C. det(inv(A))D. det(A^2)答案：A3. 在MATLAB中，下列哪个命令用于计算矩阵A的逆？A. inv(A)B. A^(-1)C. pinv(A)D. A\B答案：A4. 在MATLAB中，下列哪个命令用于求解线性方程组Ax= b？A. A\bB. A/BC. B/AD. A^-1b答案：A5. 在MATLAB中，下列哪个命令用于绘制二维图形？A. plot(x, y)B. scatter(x, y)C. bar(x, y)D. pie(x, y)答案：A二、填空题（每题5分，共25分）6. 在MATLAB中，可以使用______命令创建一个等差数列。
答案：linspace7. 在MATLAB中，可以使用______命令创建一个等比数列。
答案：logspace8. 在MATLAB中，可以使用______命令计算矩阵A的特征值。
答案：eig(A)9. 在MATLAB中，可以使用______命令计算矩阵A的特征向量。
答案：eigenvector(A)10. 在MATLAB中，可以使用______命令计算矩阵A的奇异值。
答案：svd(A)三、解答题（每题25分，共75分）11. 编写MATLAB程序，求解以下线性方程组：2x + 3y - z = 1x - y + 2z = 03x + 2y - 4z = -3答案：```A = [2 3 -1; 1 -1 2; 3 2 -4];b = [1; 0; -3];x = A\b;disp('解为：');disp(x);```12. 编写MATLAB程序，绘制以下函数的图形：y = sin(x) + cos(x)，x ∈ [0, 2π]答案：```x = linspace(0, 2pi, 100);y = sin(x) + cos(x);plot(x, y);title('y = sin(x) + cos(x)');xlabel('x');ylabel('y');grid on;```13. 编写MATLAB程序，计算以下矩阵的特征值和特征向量：A = [1 2 3; 4 5 6; 7 8 9]答案：```A = [1 2 3; 4 5 6; 7 8 9];[V, D] = eig(A);disp('特征值：');disp(diag(D));disp('特征向量：');disp(V);```14. 编写MATLAB程序，使用牛顿迭代法求解方程f(x) = x^3 - 4x + 2 = 0在x = 1附近的根。
matlab经典编程例题30道
MATLAB是一款功能强大的数学软件，其编程功能也受到越来越多人的关注。
下面介绍的是30个matlab经典编程例题，可以帮助大家熟悉matlab的编程语法，提高matlab 编程技能。
1. 请编写一个程序，计算出两个数的和。
2. 请编写一个程序，计算出两个数的最大值和最小值。
3. 请编写一个程序，计算出一组数据的平均值和标准差。
4. 请编写一个程序，将一个矩阵转置。
5. 请编写一个程序，求出两个矩阵的乘积。
6. 请编写一个程序，求出一个矩阵的逆矩阵。
7. 请编写一个程序，求出一个矩阵的行列式。
8. 请编写一个程序，计算出一元二次方程的解。
9. 请编写一个程序，计算出两个数组的相似度。
10. 请编写一个程序，计算出一个矩阵的特征值和特征向量。
11. 请编写一个程序，求出两个矩阵的秩。
12. 请编写一个程序，求出一个矩阵的特定元素。
13. 请编写一个程序，求出一组数据的最高值和最低值。
14. 请编写一个程序，求出两个数组的交集。
15. 请编写一个程序，求出一个矩阵的行和列之和。
16. 请编写一个程序，使用循环语句计算出100以内所有奇数的和。
17. 请编写一个程序，使用循环语句计算出1到1000以内的和。
18. 请编写一个程序，使用递归函数计算出斐波那契数列的第n项。
19. 请编写一个程序，求出一个多项式的导数。
20. 请编写一个程序，求出一个函数的极值点。
21. 请编写一个程序，求出一个数组的非零元素个数。
22. 请编写一个程序，计算出函数的不定积分。
23. 请编写一个程序，计算出函数的定积分。
24. 请编写一个程序，求出一个矩阵的秩。
25. 请编写一个程序，求出函数的极限值。
26. 请编写一个程序，求出一个矩阵的特征值分解。
27. 请编写一个程序，求出一个矩阵的LU分解。
28. 请编写一个程序，求出一个矩阵的QR分解。
29. 请编写一个程序，求出三次多项式的根。
30. 请编写一个程序，求出一个函数的积分。
（完整版）matlab经典习题及解答
（完整版）matlab经典习题及解答第1章 MATLAB 概论1.1 与其他计算机语⾔相⽐较，MATLAB 语⾔突出的特点是什么？MATLAB 具有功能强⼤、使⽤⽅便、输⼊简捷、库函数丰富、开放性强等特点。
1.2 MATLAB 系统由那些部分组成？MATLAB 系统主要由开发环境、MATLAB 数学函数库、MATLAB 语⾔、图形功能和应⽤程序接⼝五个部分组成。
1.4 MATLAB 操作桌⾯有⼏个窗⼝？如何使某个窗⼝脱离桌⾯成为独⽴窗⼝？⼜如何将脱离出去的窗⼝重新放置到桌⾯上？在MATLAB 操作桌⾯上有五个窗⼝，在每个窗⼝的右上⾓有两个⼩按钮，⼀个是关闭窗⼝的Close 按钮，⼀个是可以使窗⼝成为独⽴窗⼝的Undock 按钮，点击Undock 按钮就可以使该窗⼝脱离桌⾯成为独⽴窗⼝，在独⽴窗⼝的view 菜单中选择Dock ……菜单项就可以将独⽴的窗⼝重新防⽌的桌⾯上。
1.5 如何启动M ⽂件编辑/调试器？在操作桌⾯上选择“建⽴新⽂件”或“打开⽂件”操作时，M ⽂件编辑/调试器将被启动。
在命令窗⼝中键⼊edit 命令时也可以启动M ⽂件编辑/调试器。
1.6 存储在⼯作空间中的数组能编辑吗？如何操作？存储在⼯作空间的数组可以通过数组编辑器进⾏编辑：在⼯作空间浏览器中双击要编辑的数组名打开数组编辑器，再选中要修改的数据单元，输⼊修改内容即可。
1.7 命令历史窗⼝除了可以观察前⾯键⼊的命令外，还有什么⽤途？命令历史窗⼝除了⽤于查询以前键⼊的命令外，还可以直接执⾏命令历史窗⼝中选定的内容、将选定的内容拷贝到剪贴板中、将选定内容直接拷贝到M ⽂件中。
1.8 如何设置当前⽬录和搜索路径，在当前⽬录上的⽂件和在搜索路径上的⽂件有什么区别？当前⽬录可以在当前⽬录浏览器窗⼝左上⽅的输⼊栏中设置，搜索路径可以通过选择操作桌⾯的file 菜单中的Set Path 菜单项来完成。
在没有特别说明的情况下，只有当前⽬录和搜索路径上的函数和⽂件能够被MATLAB 运⾏和调⽤，如果在当前⽬录上有与搜索路径上相同⽂件名的⽂件时则优先执⾏当前⽬录上的⽂件，如果没有特别说明，数据⽂件将存储在当前⽬录上。
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中，如何创建一个名为"myMatrix"的3x3矩阵，其元素为1到9的连续整数？答案：可以使用以下命令创建矩阵：```matlabmyMatrix = [1 2 3; 4 5 6; 7 8 9];```或者使用`magic`函数：```matlabmyMatrix = magic(3);```2. 在MATLAB中，如何计算一个向量的元素平方？答案：可以使用`.^`运算符来计算向量元素的平方：```matlabv = [1 2 3 4];v_squared = v.^2;```3. MATLAB中，如何将一个向量逆序？答案：可以使用`flip`函数或者`end`索引来逆序一个向量：```matlabv = [1 2 3 4];v_reversed = flip(v);```或者```matlabv_reversed = v(end:-1:1);```4. 在MATLAB中，如何计算矩阵的行列式？答案：可以使用`det`函数来计算矩阵的行列式：```matlabA = [1 2; 3 4];detA = det(A);```5. MATLAB中，如何使用循环结构来计算1到100的自然数之和？答案：可以使用`for`循环来计算总和：```matlabsum = 0;for i = 1:100sum = sum + i;end```或者使用`sum`函数直接计算：```matlabsum = sum(1:100);```6. 在MATLAB中，如何创建一个复数？答案：可以通过指定实部和虚部来创建一个复数：```matlabcomplexNumber = 3 + 4i;```7. MATLAB中，如何绘制一个正弦波？答案：可以使用`plot`函数和`sin`函数来绘制正弦波：```matlabx = linspace(0, 2*pi, 100);y = sin(x);plot(x, y);```8. 如何在MATLAB中计算两个矩阵的点积？答案：可以使用`dot`函数或`*`运算符来计算两个矩阵的点积：```matlabA = [1 2; 3 4];B = [5 6; 7 8];dotProduct = dot(A, B);```或者```matlabdotProduct = A * B;```9. MATLAB中，如何使用条件语句来检查一个数是否为偶数？答案：可以使用`if`语句和模运算符来检查一个数是否为偶数：```matlabnum = 10;if mod(num, 2) == 0disp('The number is even.');elsedisp('The number is odd.');end```10. 在MATLAB中，如何将一个字符串转换为小写？答案：可以使用`lower`函数来将字符串转换为小写：```matlabstr = 'Hello World';lowerStr = lower(str);```结束语：以上是MATLAB基础考试题及答案，希望这些题目能帮助你更好地理解和掌握MATLAB的基础知识。
MATLAB程序设计例题及答案
1.编写程序：计算1/3+2/5+3/7+……+10/21法一：s=0;for i=1:10s=s+i/(2*i+1);endss=4.4096法二：sum((1:10)./(3:2:21))ans=4.40962.编写程序：计算1~100中即能被3整除，又能被7整除的所有数之和。
s=0;for i=1:100if mod(i,3)==0&&mod(i,7)==0s=s+i;end,endss=2103.画出y=n！的图（1<=n<=10）,阶乘的函数自己编写，禁用MATLAB自带的阶乘函数。
x=1:10;for i=1:10try y(i)=y(i-1)*i;catch y(i)=1;end,endplot(x,y)1234567891000.511.522.533.5464.一个数恰好等于它的因子之和，这个数就称为完数。
例如，6的因子为1，2,3，而6=1+2+3，因此6就是一个完数。
编程找出2000以内的所有完数。
g=[];for n=2:2000s=0;for r=1:n-1if mod(n,r)==0s=s+r;end end if s==n g=[g n];end end g g =6284965.编写一个函数，模拟numel函数的功能，函数中调用size函数。
function y=numelnumel(x)m=size(x);y=m(1)*m(2);numelnumel([123;456])ans=66.编写一个函数，模拟length函数的功能，函数中调用size函数。
function y=lengthlength(x)m=size(x);y=max(m(1),m(2));lengthlength([123;456])ans=37.求矩阵rand（5）的所有元素和及各行平均值，各列平均值。
s=rand(5);sum=sum(sum(s))mean2=mean(s,2)mean1=mean(s)sum=13.8469mean2=0.55830.55360.65540.49310.5090mean1=0.63310.50060.64870.71240.27458.编程判断1001,1003,1007，1009，1011是否为素数，若不是，输出其约数。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
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源程序如下：t=0:0.1:10; %t的取值以0.1的步长从0到10r=2; %r=2x=r*cos(t)+3*t; %x的表达式y=r*sin(t)+3; %y的表达式plot(x,y,'r-') %用红色线绘制曲线hold on %保持窗口figure窗口不关闭t=0:0.1:10; %t的取值以0.1的步长从0到10r=3; %r=3x=r*cos(t)+3*t; %x的表达式y=r*sin(t)+3; %y的表达式plot(x,y,'b-.') %用蓝色点划线绘制曲线hold on %保持窗口figure窗口不关闭t=0:0.1:10; %t的取值以0.1的步长从0到10r=4; %r=4x=r*cos(t)+3*t; %x的表达式y=r*sin(t)+3; %y的表达式plot(x,y,'m-') %用红紫色线绘制曲线4、设f(x)=x5- 4x4 +3x2- 2x+ 6(1) 在x=[-2,8]之间取100个点，画出曲线，看它有几个过零点。
（提示：用polyval 函数）(2) 用roots函数求此多项式的根。
(1)matalab的源程序如下：（利用的是在同一图中展现两条线，通过放大看交点）x=-2:0.1:8; %x的取值范围y=0; %y=0plot(x,y,'m') %用红紫色表示绘制出的y=0的图像hold on %保持figure窗口不关闭x=-2:0.1:8; %x的取值范围p=[1,-4,0,3,-2,6]; %p表示函数表达式的系数矩阵y=polyval(p,x),plot(x,y) %y=polyval表示求多项式x处的y值，绘制图像将图形进行放大（如下），发现在实数域上有3个零点。
（2）matlab源程序如下：roots(p) %求函数的根ans =3.7999-1.26071.34790.0564 + 0.9623i0.0564 - 0.9623i可见实数解有三个，复数解有两个。
5、建立一个字符串向量，然后对该向量做如下处理：(1) 取第1～5个字符组成的子字符串。
(2) 将字符串倒过来重新排列。
(3) 统计字符串中小写字母的个数。
整个题的源程序如下所示：ch='abcj245DFGK'; %建立一个字符串subch=ch(1:5) %取1~5个字符subch =abcj2revch=ch(end:-1:1) %倒序排列revch =KGFD542jcbaa=find(ch>='a'&ch<='z'); %找到小写字母的位置length(a) %输出小写字母的个数ans =46、设x=sint, y=sin(nt+a)，(1)若a=1，令n =1,2,3,4,在四个子图中分别画出其曲线。
(2)若n=2,取a=0，π/3，π/2，及π，在四个子图中分别画出其曲线。
（1）matlab源程序如下;t=-pi:0.1:pi; %t的取值范围a=1;n=1; %a=1，n=1x=sin(t); %x的表达式y=sin(n*t+a); %y的表达式subplot(221),plot(x,y) %在子图第一张，画出图像hold on %保持figure不关闭t=-pi:0.1:pi; %t的取值范围a=1;n=2; %a=1，n=2x=sin(t); %x的表达式y=sin(n*t+a); %y的表达式subplot(222),plot(x,y) %在子图的第二张画出图像hold on %保持figure不关闭t=-pi:0.1:pi; %t的取值范围a=1;n=3; %a=1，n=3x=sin(t); %x的表达式y=sin(n*t+a); %y的表达式subplot(223),plot(x,y) %在子图的第三张，画出图像hold on %保持figure不关闭t=-pi:0.1:pi; %t的取值范围a=1;n=4; %a=1，n=4x=sin(t); %x的表达式subplot(224),plot(x,y) %在子图的第四张上，绘制图像图形如下所示：（2）matlab源程序如下所示：t=-pi:0.1:pi; %t的取值a=0;n=2; %a与n的取值x=sin(t); %x的表达式y=sin(n*t+a); %y的表达式subplot(221),plot(x,y) %在子图的第一张绘制图像hold on %保持figure窗口不关闭t=-pi:0.1:pi; %t的取值a=pi./3;n=2; %a与n的取值x=sin(t); %x的表达式y=sin(n*t+a); %y的表达式subplot(222),plot(x,y) %在子图的第二张绘制图像hold on %保持figure窗口不关闭t=-pi:0.1:pi; %t的取值a=pi./2;n=2; %a与n的取值x=sin(t); %x的表达式y=sin(n*t+a); %y的表达式subplot(223),plot(x,y) %在子图的第三张绘制图像hold on %保持figure窗口不关闭t=-pi:0.1:pi; %t的取值范围a=pi;n=2; %a与n的取值y=sin(n*t+a); %y 的表达式subplot(224),plot(x,y) %在四个字图的第四张画出图像图形如下所示：7、绘制)(222y x e x z +-=在定义域x=[-2,2],y=[-2,2]内的曲面。
Matlab 源程序如下：[x,y]=meshgrid([-2:0.1:2]); %此处描述x ，y 的范围z=x.^2.*exp(-(x.^2+y.^2)); %z 的表达式surf(x,y,z) %绘制曲面图形如下所示：经过旋转之后的图形：8、⎪⎩⎪⎨⎧≤≤-+-<≤-+--<≤----=312/)34(111132/)34(222x x x x x x x x y ，画出在[]3,3-上的曲线所建立的M 文件如下所示：function y=A(x); %文件名为Aif x>=-3&x<-1 %在-3《x<-1的条件下y=(-x.^2-4*x-3)./2; %y 所对应上述条件下的表达式 elseif x>=-1&x<1 %-1《x<1y=-x.^2+1; %y 所对应的上述条件下的表达式 else x>=1&x<=3 %1《x 《3y=(-x.^2+4*x-3)./2; %y 所对应的上述条件下的表达式 end %结束在command windows 中输入：fplot('A',[-3,3]) %绘制在-3到3上的曲线所得图形如下所示：9、（1）设,求syms x %定义变量xy=x.*exp(2.*x); %y 的表达式s=diff(y,x,20) %对y 求x 的20阶导数10485760*exp(2*x) + 1048576*x*exp(2*x)（2）设,求syms x %定义变量xy=sin(x).^4+cos(x).^4; %y 的表达式s=diff(y,x,10) %对y 求x 的10阶导数s =- 262144*cos(x)^4 + 1572864*cos(x)^2*sin(x)^2 - 262144*sin(x)^4（3）设,求syms x y %定义自变量x ，yz=y.^4.*(sin(x).^4)+x.^4.*(cos(y).^4); %z 的表达式s=diff(z,x) %对z 求x 的偏导数s =4*x^3*cos(y)^4 + 4*y^4*cos(x)*sin(x)^3（4）⎰++dx x x x ))(1(122 syms x %定义自变量xi=int(1/((x.^2+1).*(x.^2+x))); %不定积分表达式pretty(i) %求解不定积分/ i \ log(x + 1) / i \log(x) + log(x + i) | - - - 1/4 | - ---------- + log(x - i) | - - 1/4 |\ 4 / 2 \ 4 /（5）dx x x ⎰-632)9(syms x %定义自变量 xi=int(sqrt((9-x.^2).^3)./(x.^6)); %不定积分表达式pretty(i) %求解不定积分/ 9 \ 2 3 1/2| -- - 1 | (- (x - 9) )| 2 |\ x /- ---------------------------345 x（6）⎰+213x x dx syms x %定义自变量xi=int(1./(x+x.^3),2,1) %定积分表达式i = %定积分的解log(10^(1/2)/4)（7）求级数的和: ∑∞=+122n n n syms n %定义自变量nf=(n+2)./(2.^n); %级数的表达式symsum(f,n,n-1,inf) %对级数求和ans =2*2^(1 - n)*(n + 2)（8）利用函数int 计算二重不定积分syms x y %定义自变量x ，yi=int(int((x+y).*exp(-x.*y))); %二重不定积分表达式pretty(i) %求解二重不定积分2y + x y + 2------------3y exp(x y)10、编制m文件，输入n（正整数），显示所有小于n的质数。