matlab典型例题
matlab编程实例100例
1-32是：图形应用篇33-66是：界面设计篇67-84是：图形处理篇85-100是：数值分析篇实例1：三角函数曲线（1）function shili01h0=figure('toolbar','none',...'position',[198 56 350 300],...'name','实例01');h1=axes('parent',h0,...'visible','off');x=-pi:0.05:pi;y=sin(x);plot(x,y);xlabel('自变量X');ylabel('函数值Y');title('SIN( )函数曲线');grid on实例2：三角函数曲线（2）function shili02h0=figure('toolbar','none',...'position',[200 150 450 350],...'name','实例02');x=-pi:0.05:pi;y=sin(x)+cos(x);plot(x,y,'-*r','linewidth',1);grid onxlabel('自变量X');ylabel('函数值Y');title('三角函数');实例3：图形的叠加function shili03h0=figure('toolbar','none',...'position',[200 150 450 350],...'name','实例03');x=-pi:0.05:pi;y1=sin(x);y2=cos(x);plot(x,y1,...'-*r',...x,y2,...'--og');grid onxlabel('自变量X');ylabel('函数值Y');title('三角函数');实例4：双y轴图形的绘制function shili04h0=figure('toolbar','none',...'position',[200 150 450 250],...'name','实例04');x=0:900;a=1000;b=0.005;y1=2*x;y2=cos(b*x);[haxes,hline1,hline2]=plotyy(x,y1,x,y2,'semilogy','plot'); axes(haxes(1))ylabel('semilog plot');axes(haxes(2))ylabel('linear plot');实例5：单个轴窗口显示多个图形function shili05h0=figure('toolbar','none',...'position',[200 150 450 250],...'name','实例05');t=0:pi/10:2*pi;[x,y]=meshgrid(t);subplot(2,2,1)plot(sin(t),cos(t))axis equalsubplot(2,2,2)z=sin(x)-cos(y);plot(t,z)axis([0 2*pi -2 2])subplot(2,2,3)h=sin(x)+cos(y);plot(t,h)axis([0 2*pi -2 2])subplot(2,2,4)g=(sin(x).^2)-(cos(y).^2);plot(t,g)axis([0 2*pi -1 1])实例6：图形标注function shili06h0=figure('toolbar','none',...'position',[200 150 450 400],...'name','实例06');t=0:pi/10:2*pi;h=plot(t,sin(t));xlabel('t=0到2\pi','fontsize',16);ylabel('sin(t)','fontsize',16);title('\it{从0to2\pi 的正弦曲线}','fontsize',16) x=get(h,'xdata');y=get(h,'ydata');imin=find(min(y)==y);imax=find(max(y)==y);text(x(imin),y(imin),...['\leftarrow最小值=',num2str(y(imin))],...'fontsize',16)text(x(imax),y(imax),...['\leftarrow最大值=',num2str(y(imax))],...'fontsize',16)实例7：条形图形function shili07h0=figure('toolbar','none',...'position',[200 150 450 350],...'name','实例07');tiao1=[562 548 224 545 41 445 745 512];tiao2=[47 48 57 58 54 52 65 48];t=0:7;bar(t,tiao1)xlabel('X轴');ylabel('TIAO1值');h1=gca;h2=axes('position',get(h1,'position'));plot(t,tiao2,'linewidth',3)set(h2,'yaxislocation','right','color','none','xticklabel',[]) 实例8：区域图形function shili08h0=figure('toolbar','none',...'position',[200 150 450 250],...'name','实例08');x=91:95;profits1=[88 75 84 93 77];profits2=[51 64 54 56 68];profits3=[42 54 34 25 24];profits4=[26 38 18 15 4];area(x,profits1,'facecolor',[0.5 0.9 0.6],...'edgecolor','b',...'linewidth',3)hold onarea(x,profits2,'facecolor',[0.9 0.85 0.7],...'edgecolor','y',...'linewidth',3)hold onarea(x,profits3,'facecolor',[0.3 0.6 0.7],...'edgecolor','r',...'linewidth',3)hold onarea(x,profits4,'facecolor',[0.6 0.5 0.9],...'edgecolor','m',...'linewidth',3)hold offset(gca,'xtick',[91:95])set(gca,'layer','top')gtext('\leftarrow第一季度销量')gtext('\leftarrow第二季度销量')gtext('\leftarrow第三季度销量')gtext('\leftarrow第四季度销量')xlabel('年','fontsize',16);ylabel('销售量','fontsize',16);实例9：饼图的绘制function shili09h0=figure('toolbar','none',...'position',[200 150 450 250],...'name','实例09');t=[54 21 35;68 54 35;45 25 12;48 68 45;68 54 69];x=sum(t);h=pie(x);textobjs=findobj(h,'type','text');str1=get(textobjs,{'string'});val1=get(textobjs,{'extent'});oldext=cat(1,val1{:});names={'商品一：';'商品二：';'商品三：'};str2=strcat(names,str1);set(textobjs,{'string'},str2)val2=get(textobjs,{'extent'});newext=cat(1,val2{:});offset=sign(oldext(:,1)).*(newext(:,3)-oldext(:,3))/2; pos=get(textobjs,{'position'});textpos=cat(1,pos{:});textpos(:,1)=textpos(:,1)+offset;set(textobjs,{'position'},num2cell(textpos,[3,2]))实例10：阶梯图function shili10h0=figure('toolbar','none',...'position',[200 150 450 400],...'name','实例10');a=0.01;b=0.5;t=0:10;f=exp(-a*t).*sin(b*t);stairs(t,f)hold onplot(t,f,':*')hold offglabel='函数e^{-(\alpha*t)}sin\beta*t的阶梯图'; gtext(glabel,'fontsize',16)xlabel('t=0:10','fontsize',16)axis([0 10 -1.2 1.2])实例11：枝干图function shili11h0=figure('toolbar','none',...'position',[200 150 450 350],...'name','实例11');x=0:pi/20:2*pi;y1=sin(x);y2=cos(x);h1=stem(x,y1+y2);hold onh2=plot(x,y1,'^r',x,y2,'*g');hold offh3=[h1(1);h2];legend(h3,'y1+y2','y1=sin(x)','y2=cos(x)') xlabel('自变量X');ylabel('函数值Y');title('正弦函数与余弦函数的线性组合');实例12：罗盘图function shili12h0=figure('toolbar','none',...'position',[200 150 450 250],...'name','实例12');winddirection=[54 24 65 84256 12 235 62125 324 34 254];windpower=[2 5 5 36 8 12 76 14 10 8];rdirection=winddirection*pi/180;[x,y]=pol2cart(rdirection,windpower); compass(x,y);desc={'风向和风力','气象台','10月1日0：00到','10月1日12：00'};gtext(desc)实例13：轮廓图function shili13h0=figure('toolbar','none',...'position',[200 150 450 250],...'name','实例13');[th,r]=meshgrid((0:10:360)*pi/180,0:0.05:1); [x,y]=pol2cart(th,r);z=x+i*y;f=(z.^4-1).^(0.25);contour(x,y,abs(f),20)axis equalxlabel('实部','fontsize',16);ylabel('虚部','fontsize',16);h=polar([0 2*pi],[0 1]);delete(h)hold oncontour(x,y,abs(f),20)实例14：交互式图形function shili14h0=figure('toolbar','none',...'position',[200 150 450 250],...'name','实例14');axis([0 10 0 10]);hold onx=[];y=[];n=0;disp('单击鼠标左键点取需要的点'); disp('单击鼠标右键点取最后一个点'); but=1;while but==1[xi,yi,but]=ginput(1);plot(xi,yi,'bo')n=n+1;disp('单击鼠标左键点取下一个点');x(n,1)=xi;y(n,1)=yi;endt=1:n;ts=1:0.1:n;xs=spline(t,x,ts);ys=spline(t,y,ts);plot(xs,ys,'r-');hold off实例14：交互式图形function shili14h0=figure('toolbar','none',...'position',[200 150 450 250],...'name','实例14');axis([0 10 0 10]);hold onx=[];y=[];n=0;disp('单击鼠标左键点取需要的点'); disp('单击鼠标右键点取最后一个点'); but=1;while but==1[xi,yi,but]=ginput(1);plot(xi,yi,'bo')n=n+1;disp('单击鼠标左键点取下一个点');x(n,1)=xi;y(n,1)=yi;endt=1:n;ts=1:0.1:n;xs=spline(t,x,ts);ys=spline(t,y,ts);plot(xs,ys,'r-');hold off实例15：变换的傅立叶函数曲线function shili15h0=figure('toolbar','none',...'position',[200 150 450 250],...'name','实例15');axis equalm=moviein(20,gcf);set(gca,'nextplot','replacechildren') h=uicontrol('style','slider','position',...[100 10 500 20],'min',1,'max',20) for j=1:20plot(fft(eye(j+16)))set(h,'value',j)m(:,j)=getframe(gcf);endclf;axes('position',[0 0 1 1]);movie(m,30)实例16：劳伦兹非线形方程的无序活动function shili15h0=figure('toolbar','none',...'position',[200 150 450 250],...'name','实例15');axis equalm=moviein(20,gcf);set(gca,'nextplot','replacechildren') h=uicontrol('style','slider','position',...[100 10 500 20],'min',1,'max',20) for j=1:20plot(fft(eye(j+16)))set(h,'value',j)m(:,j)=getframe(gcf);endclf;axes('position',[0 0 1 1]);movie(m,30)实例17：填充图function shili17h0=figure('toolbar','none',...'position',[200 150 450 250],...'name','实例17');t=(1:2:15)*pi/8;x=sin(t);y=cos(t);fill(x,y,'r')axis square offtext(0,0,'STOP',...'color',[1 1 1],...'fontsize',50,...'horizontalalignment','center') 例18：条形图和阶梯形图function shili18h0=figure('toolbar','none',...'position',[200 150 450 250],...'name','实例18');subplot(2,2,1)x=-3:0.2:3;y=exp(-x.*x);bar(x,y)title('2-D Bar Chart')subplot(2,2,2)x=-3:0.2:3;y=exp(-x.*x);bar3(x,y,'r')title('3-D Bar Chart')subplot(2,2,3)x=-3:0.2:3;y=exp(-x.*x);stairs(x,y)title('Stair Chart')subplot(2,2,4)x=-3:0.2:3;y=exp(-x.*x);barh(x,y)title('Horizontal Bar Chart')实例19：三维曲线图function shili19h0=figure('toolbar','none',...'position',[200 150 450 400],...'name','实例19');subplot(2,1,1)x=linspace(0,2*pi);y1=sin(x);y2=cos(x);y3=sin(x)+cos(x);z1=zeros(size(x));z2=0.5*z1;z3=z1;plot3(x,y1,z1,x,y2,z2,x,y3,z3)grid onxlabel('X轴');ylabel('Y轴');zlabel('Z轴');title('Figure1:3-D Plot')subplot(2,1,2)x=linspace(0,2*pi);y1=sin(x);y2=cos(x);y3=sin(x)+cos(x);z1=zeros(size(x));z2=0.5*z1;z3=z1;plot3(x,z1,y1,x,z2,y2,x,z3,y3)grid onxlabel('X轴');ylabel('Y轴');zlabel('Z轴');title('Figure2:3-D Plot')实例20：图形的隐藏属性function shili20h0=figure('toolbar','none',...'position',[200 150 450 300],...'name','实例20');subplot(1,2,1)[x,y,z]=sphere(10);mesh(x,y,z)axis offtitle('Figure1:Opaque')hidden onsubplot(1,2,2)[x,y,z]=sphere(10);mesh(x,y,z)axis offtitle('Figure2:Transparent') hidden off实例21PEAKS函数曲线function shili21h0=figure('toolbar','none',...'position',[200 100 450 450],...'name','实例21');[x,y,z]=peaks(30);subplot(2,1,1)x=x(1,:);y=y(:,1);i=find(y>0.8&y<1.2);j=find(x>-0.6&x<0.5);z(i,j)=nan*z(i,j);surfc(x,y,z)xlabel('X轴');ylabel('Y轴');zlabel('Z轴');title('Figure1:surfc函数形成的曲面') subplot(2,1,2)x=x(1,:);y=y(:,1);i=find(y>0.8&y<1.2);j=find(x>-0.6&x<0.5);z(i,j)=nan*z(i,j);surfl(x,y,z)xlabel('X轴');ylabel('Y轴');zlabel('Z轴');title('Figure2:surfl函数形成的曲面')实例22：片状图function shili22h0=figure('toolbar','none',...'position',[200 150 550 350],...'name','实例22');subplot(1,2,1)x=rand(1,20);y=rand(1,20);z=peaks(x,y*pi);t=delaunay(x,y);trimesh(t,x,y,z)hidden offtitle('Figure1:Triangular Surface Plot'); subplot(1,2,2)x=rand(1,20);y=rand(1,20);z=peaks(x,y*pi);t=delaunay(x,y);trisurf(t,x,y,z)title('Figure1:Triangular Surface Plot'); 实例23：视角的调整function shili23h0=figure('toolbar','none',...'position',[200 150 450 350],...'name','实例23');x=-5:0.5:5;[x,y]=meshgrid(x);r=sqrt(x.^2+y.^2)+eps;z=sin(r)./r;subplot(2,2,1)surf(x,y,z)xlabel('X-axis')ylabel('Y-axis')zlabel('Z-axis')title('Figure1')view(-37.5,30)subplot(2,2,2)surf(x,y,z)xlabel('X-axis')ylabel('Y-axis')zlabel('Z-axis')title('Figure2')view(-37.5+90,30)subplot(2,2,3)surf(x,y,z)xlabel('X-axis')ylabel('Y-axis')zlabel('Z-axis')title('Figure3')view(-37.5,60)subplot(2,2,4)surf(x,y,z)xlabel('X-axis')ylabel('Y-axis')zlabel('Z-axis')title('Figure4')view(180,0)实例24：向量场的绘制function shili24h0=figure('toolbar','none',...'position',[200 150 450 350],...'name','实例24');subplot(2,2,1)z=peaks;ribbon(z)title('Figure1')subplot(2,2,2)[x,y,z]=peaks(15);[dx,dy]=gradient(z,0.5,0.5); contour(x,y,z,10)hold onquiver(x,y,dx,dy)hold offtitle('Figure2')subplot(2,2,3)[x,y,z]=peaks(15);[nx,ny,nz]=surfnorm(x,y,z);surf(x,y,z)hold onquiver3(x,y,z,nx,ny,nz)hold offtitle('Figure3')subplot(2,2,4)x=rand(3,5);y=rand(3,5);z=rand(3,5);c=rand(3,5);fill3(x,y,z,c)grid ontitle('Figure4')实例25：灯光定位function shili25h0=figure('toolbar','none',...'position',[200 150 450 250],...'name','实例25');vert=[1 1 1;1 2 1;2 2 1;2 1 1;1 1 2;12 2;2 2 2;2 1 2];fac=[1 2 3 4;2 6 7 3;4 3 7 8;15 8 4;1 2 6 5;5 6 7 8];grid offsphere(36)h=findobj('type','surface');set(h,'facelighting','phong',...'facecolor',...'interp',...'edgecolor',[0.4 0.4 0.4],...'backfacelighting',...'lit')hold onpatch('faces',fac,'vertices',vert,...'facecolor','y');light('position',[1 3 2]);light('position',[-3 -1 3]); material shinyaxis vis3d offhold off实例26：柱状图function shili26h0=figure('toolbar','none',...'position',[200 50 450 450],...'name','实例26');subplot(2,1,1)x=[5 2 18 7 39 8 65 5 54 3 2];bar(x)xlabel('X轴');ylabel('Y轴');title('第一子图');subplot(2,1,2)y=[5 2 18 7 39 8 65 5 54 3 2];barh(y)xlabel('X轴');ylabel('Y轴');title('第二子图');实例27：设置照明方式function shili27h0=figure('toolbar','none',...'position',[200 150 450 350],...'name','实例27');subplot(2,2,1)sphereshading flatcamlight leftcamlight rightlighting flatcolorbaraxis offtitle('Figure1')subplot(2,2,2)sphereshading flatcamlight leftcamlight rightlighting gouraudcolorbaraxis offtitle('Figure2')subplot(2,2,3)sphereshading interpcamlight rightcamlight leftlighting phongcolorbaraxis offtitle('Figure3')subplot(2,2,4)sphereshading flatcamlight leftcamlight rightlighting nonecolorbaraxis offtitle('Figure4')实例28：羽状图function shili28h0=figure('toolbar','none',...'position',[200 150 450 350],...'name','实例28');subplot(2,1,1)alpha=90:-10:0;r=ones(size(alpha));m=alpha*pi/180;n=r*10;[u,v]=pol2cart(m,n);feather(u,v)title('羽状图')axis([0 20 0 10])subplot(2,1,2)t=0:0.5:10;x=0.05+i;y=exp(-x*t);feather(y)title('复数矩阵的羽状图')实例29：立体透视（1）function shili29h0=figure('toolbar','none',...'position',[200 150 450 250],...'name','实例29');[x,y,z]=meshgrid(-2:0.1:2,...-2:0.1:2,...-2:0.1:2);v=x.*exp(-x.^2-y.^2-z.^2);grid onfor i=-2:0.5:2;h1=surf(linspace(-2,2,20),...linspace(-2,2,20),...zeros(20)+i);rotate(h1,[1 -1 1],30)dx=get(h1,'xdata');dy=get(h1,'ydata');dz=get(h1,'zdata');delete(h1)slice(x,y,z,v,[-2 2],2,-2)hold onslice(x,y,z,v,dx,dy,dz)hold offaxis tightview(-5,10)drawnowend实例30：立体透视（2）function shili30h0=figure('toolbar','none',...'position',[200 150 450 250],...'name','实例30');[x,y,z]=meshgrid(-2:0.1:2,...-2:0.1:2,...-2:0.1:2);v=x.*exp(-x.^2-y.^2-z.^2); [dx,dy,dz]=cylinder;slice(x,y,z,v,[-2 2],2,-2)for i=-2:0.2:2h=surface(dx+i,dy,dz);rotate(h,[1 0 0],90)xp=get(h,'xdata');yp=get(h,'ydata');zp=get(h,'zdata');delete(h)hold onhs=slice(x,y,z,v,xp,yp,zp);axis tightxlim([-3 3])view(-10,35)drawnowdelete(hs)hold offend实例31：表面图形function shili31h0=figure('toolbar','none',...'position',[200 150 550 250],...'name','实例31');subplot(1,2,1)x=rand(100,1)*16-8;y=rand(100,1)*16-8;r=sqrt(x.^2+y.^2)+eps;z=sin(r)./r;xlin=linspace(min(x),max(x),33); ylin=linspace(min(y),max(y),33); [X,Y]=meshgrid(xlin,ylin);Z=griddata(x,y,z,X,Y,'cubic'); mesh(X,Y,Z)axis tighthold onplot3(x,y,z,'.','Markersize',20)subplot(1,2,2)k=5;n=2^k-1;theta=pi*(-n:2:n)/n;phi=(pi/2)*(-n:2:n)'/n;X=cos(phi)*cos(theta);Y=cos(phi)*sin(theta);Z=sin(phi)*ones(size(theta));colormap([0 0 0;1 1 1])C=hadamard(2^k);surf(X,Y,Z,C)axis square实例32：沿曲线移动的小球h0=figure('toolbar','none',...'position',[198 56 408 468],...'name','实例32');h1=axes('parent',h0,...'position',[0.15 0.45 0.7 0.5],...'visible','on');t=0:pi/24:4*pi;y=sin(t);plot(t,y,'b')n=length(t);h=line('color',[0 0.5 0.5],...'linestyle','.',...'markersize',25,...'erasemode','xor');k1=uicontrol('parent',h0,...'style','pushbutton',...'position',[80 100 50 30],...'string','开始',...'callback',[...'i=1;',...'k=1;,',...'m=0;,',...'while 1,',...'if k==0,',...'break,',...'end,',...'if k~=0,',...'set(h,''xdata'',t(i),''ydata'',y(i)),',...'drawnow;,',...'i=i+1;,',...'if i>n,',...'m=m+1;,',...'i=1;,',...'end,',...'end,',...'end']);k2=uicontrol('parent',h0,...'style','pushbutton',...'position',[180 100 50 30],...'string','停止',...'callback',[...'k=0;,',...'set(e1,''string'',m),',...'p=get(h,''xdata'');,',...'q=get(h,''ydata'');,',...'set(e2,''string'',p);,',...'set(e3,''string'',q)']); k3=uicontrol('parent',h0,...'style','pushbutton',...'position',[280 100 50 30],...'string','关闭',...'callback','close');e1=uicontrol('parent',h0,...'style','edit',...'position',[60 30 60 20]);t1=uicontrol('parent',h0,...'style','text',...'string','循环次数',...'position',[60 50 60 20]);e2=uicontrol('parent',h0,...'style','edit',...'position',[180 30 50 20]); t2=uicontrol('parent',h0,...'style','text',...'string','终点的X坐标值',...'position',[155 50 100 20]); e3=uicontrol('parent',h0,...'style','edit',...'position',[300 30 50 20]); t3=uicontrol('parent',h0,...'style','text',...'string','终点的Y坐标值',...'position',[275 50 100 20]);实例33：曲线转换按钮h0=figure('toolbar','none',...'position',[200 150 450 250],...'name','实例33');x=0:0.5:2*pi;y=sin(x);h=plot(x,y);grid onhuidiao=[...'if i==1,',...'i=0;,',...'y=cos(x);,',...'delete(h),',...'set(hm,''string'',''正弦函数''),',...'h=plot(x,y);,',...'grid on,',...'else if i==0,',...'i=1;,',...'y=sin(x);,',...'set(hm,''string'',''余弦函数''),',...'delete(h),',...'h=plot(x,y);,',...'grid on,',...'end,',...'end'];hm=uicontrol(gcf,'style','pushbutton',...'string','余弦函数',...'callback',huidiao);i=1;set(hm,'position',[250 20 60 20]);set(gca,'position',[0.2 0.2 0.6 0.6])title('按钮的使用')hold on实例34：栅格控制按钮h0=figure('toolbar','none',...'position',[200 150 450 250],...'name','实例34');x=0:0.5:2*pi;y=sin(x);plot(x,y)huidiao1=[...'set(h_toggle2,''value'',0),',...'grid on,',...];huidiao2=[...'set(h_toggle1,''value'',0),',...'grid off,',...];h_toggle1=uicontrol(gcf,'style','togglebutton',...'string','grid on',...'value',0,...'position',[20 45 50 20],...'callback',huidiao1);h_toggle2=uicontrol(gcf,'style','togglebutton',...'string','grid off',...'value',0,...'position',[20 20 50 20],...'callback',huidiao2);set(gca,'position',[0.2 0.2 0.6 0.6])title('开关按钮的使用')实例35：编辑框的使用h0=figure('toolbar','none',...'position',[200 150 350 250],...'name','实例35');f='Please input the letter';huidiao1=[...'g=upper(f);,',...'set(h2_edit,''string'',g),',...];huidiao2=[...'g=lower(f);,',...'set(h2_edit,''string'',g),',...];h1_edit=uicontrol(gcf,'style','edit',...'position',[100 200 100 50],...'HorizontalAlignment','left',...'string','Please input the letter',...'callback','f=get(h1_edit,''string'');',...'background','w',...'max',5,...'min',1);h2_edit=uicontrol(gcf,'style','edit',...'HorizontalAlignment','left',...'position',[100 100 100 50],...'max',5,...'min',1);h1_button=uicontrol(gcf,'style','pushbutton',...'string','小写变大写',...'position',[100 45 100 20],...'callback',huidiao1);h2_button=uicontrol(gcf,'style','pushbutton',...'string','大写变小写',...'position',[100 20 100 20],...'callback',huidiao2);实例36：弹出式菜单h0=figure('toolbar','none',...'position',[200 150 450 250],...'name','实例36');x=0:0.5:2*pi;y=sin(x);h=plot(x,y);grid onhm=uicontrol(gcf,'style','popupmenu',...'string',...'sin(x)|cos(x)|sin(x)+cos(x)|exp(-sin(x))',...'position',[250 20 50 20]);set(hm,'value',1)huidiao=[...'v=get(hm,''value'');,',...'switch v,',...'case 1,',...'delete(h),',...'y=sin(x);,',...'h=plot(x,y);,',...'grid on,',...'case 2,',...'delete(h),',...'y=cos(x);,',...'h=plot(x,y);,',...'grid on,',...'case 3,',...'delete(h),',...'y=sin(x)+cos(x);,',...'h=plot(x,y);,',...'grid on,',...'case 4,',...'y=exp(-sin(x));,',...'h=plot(x,y);,',...'grid on,',...'end'];set(hm,'callback',huidiao)set(gca,'position',[0.2 0.2 0.6 0.6])title('弹出式菜单的使用')实例37：滑标的使用h0=figure('toolbar','none',...'position',[200 150 450 250],...'name','实例37');[x,y]=meshgrid(-8:0.5:8);r=sqrt(x.^2+y.^2)+eps;z=sin(r)./r;h0=mesh(x,y,z);h1=axes('position',...[0.2 0.2 0.5 0.5],...'visible','off');htext=uicontrol(gcf,...'units','points',...'position',[20 30 45 15],...'string','brightness',...'style','text');hslider=uicontrol(gcf,...'units','points',...'position',[10 10 300 15],...'min',-1,...'max',1,...'style','slider',...'callback',...'brighten(get(hslider,''value''))');实例38：多选菜单h0=figure('toolbar','none',...'position',[200 150 450 250],...'name','实例38');[x,y]=meshgrid(-8:0.5:8);r=sqrt(x.^2+y.^2)+eps;z=sin(r)./r;h0=mesh(x,y,z);hlist=uicontrol(gcf,'style','listbox',...'string','default|spring|summer|autumn|winter',...'max',5,...'min',1,...'position',[20 20 80 100],...'callback',[...'k=get(hlist,''value'');,',...'switch k,',...'case 1,',...'colormap default,',...'case 2,',...'colormap spring,',...'case 3,',...'colormap summer,',...'case 4,',...'colormap autumn,',...'case 5,',...'colormap winter,',...'end']);实例39：菜单控制的使用h0=figure('toolbar','none',...'position',[200 150 450 250],...'name','实例39');x=0:0.5:2*pi;y=cos(x);h=plot(x,y);grid onset(gcf,'toolbar','none')hm=uimenu('label','example');huidiao1=[...'set(hm_gridon,''checked'',''on''),',...'set(hm_gridoff,''checked'',''off''),',...'grid on'];huidiao2=[...'set(hm_gridoff,''checked'',''on''),',...'set(hm_gridon,''checked'',''off''),',...'grid off'];hm_gridon=uimenu(hm,'label','grid on',...'checked','on',...'callback',huidiao1);hm_gridoff=uimenu(hm,'label','grid off',...'checked','off',...'callback',huidiao2);实例40：UIMENU菜单的应用。
MATLAB习题及参考答案经典.doc
习题：1， 计算⎥⎦⎤⎢⎣⎡=572396a 与⎥⎦⎤⎢⎣⎡=864142b 的数组乘积。
2， 对于B AX =，如果⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=753467294A ，⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=282637B ，求解X 。
3， 已知：⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=987654321a ，分别计算a 的数组平方和矩阵平方，并观察其结果。
4， 角度[]604530=x ，求x 的正弦、余弦、正切和余切。
(应用sin,cos,tan.cot)5， 将矩阵⎥⎦⎤⎢⎣⎡=7524a 、⎥⎦⎤⎢⎣⎡=3817b 和⎥⎦⎤⎢⎣⎡=2695c 组合成两个新矩阵： （1）组合成一个4⨯3的矩阵，第一列为按列顺序排列的a 矩阵元素，第二列为按列顺序排列的b 矩阵元素，第三列为按列顺序排列的c 矩阵元素，即 ⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡237912685574（2）按照a 、b 、c 的列顺序组合成一个行矢量，即 []2965318772546， 将(x -6)(x -3)(x -8)展开为系数多项式的形式。
(应用poly,polyvalm)7， 求解多项式x 3-7x 2+2x +40的根。
(应用roots)8， 求解在x =8时多项式(x -1)(x -2) (x -3)(x -4)的值。
(应用poly,polyvalm)9， 计算多项式9514124234++--x x x x 的微分和积分。
(应用polyder,polyint ，poly2sym)10， 解方程组⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡66136221143092x 。
(应用x=a\b)11， 求欠定方程组⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡5865394742x 的最小范数解。
(应用pinv) 12， 矩阵⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=943457624a ，计算a 的行列式和逆矩阵。
(应用det,inv)13， y =sin(x )，x 从0到2π，∆x =0.02π，求y 的最大值、最小值、均值和标准差。
matlab程序设计例题及答案
matlab程序设计例题及答案1.编写程序：计算1/3+2/5+3/7+……+10/21法一： s=0;for i=1:10s=s+i/(2*i+1); end ss =法二：sum((1:10)./(3:2:21)) ans =2.编写程序：计算1~100中即能被3整除，又能被7整除的所有数之和。
s=0;for i=1:100if mod(i,3)==0&&mod(i,7)==0 s=s+i; end,end ss =2103.画出y=n！的图,阶乘的函数自己编写，禁用MATLAB 自带的阶乘函数。
x=1:10; for i=1:10try y(i)=y(i-1)*i; catch y(i)=1; end,end plot(x,y)106123456789104.一个数恰好等于它的因子之和，这个数就称为完数。
例如，6的因子为1，2,3，而6=1+2+3，因此6就是一个完数。
编程找出20XX以内的所有完数。
g=;for n=2:20XX s=0;for r=1:n-1if mod(n,r)==0 s=s+r; end endif s==ng=[g n]; end end gg =6 28 4965.编写一个函数，模拟numel函数的功能，函数中调用size函数。
function y=numelnumel(x) m=size(x); y=m(1)*m(2);numelnumel([1 2 3;4 5 6])ans =66. 编写一个函数，模拟length函数的功能，函数中调用size函数。
function y=lengthlength(x) m=size(x);y=max(m(1),m(2));lengthlength([1 2 3;4 5 6])ans =37.求矩阵rand的所有元素和及各行平均值，各列平均值。
s=rand(5);sum=sum(sum(s)) mean2=mean(s,2) mean1=mean(s)sum =mean2 =mean1 =8.编程判断1001,1003,1007，1009，1011为素数，若不是，输出其约数。
matlab的例题
浅议matlab 的功能和应用数值计算1、编程求满足∑=>mi i 1100002的最小m 值。
解：m=1; s=0;while s<=10000s=s+2^m; m=m+1;end m=m-1;di sp('m=');disp(m);运行结果：m=132、编写一个函数，计算下面函数的值，给出标量x 的值，调用该函数后，返回y 的值。
function [y]=myfun1(x)选择一些数据测试你编写的函数。
function y=myfun1(x) if x<=0y=sin(x); elseif x>0&x<=3 y=x; else y=6-x; end⎪⎩⎪⎨⎧>+-≤<≤=3,630,0,sin )(x x x x x x x yreturn ;3、编写一个函数，给出一个向量],,[21n x x x x=，生成如下范德蒙矩阵。
function [v]=myvander(x)⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡---112112222121111n n n n n n x x x x x x x x x例如：>>v=myvander([2 3 4 5]) 得v=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡1256427825169454321111 生成一些数据测试你写的函数。
function [v]=myvander(x) n=length(x); v(1:n)=1; for i=n+1:n*n a=i/n;r=mod(i,n); b=floor(a); if r==0 r=n; b=b-1; endv(i)=x(r)^b; endv=reshape(v,n,n)';>> myvander([1 2 3])ans =1 1 11 2 3 1 4 9>> myvander([ 1 4 6 8])ans =1 1 1 1 1 4 6 8 1 16 36 64 1 64 216 512图像的绘制1、在同一坐标系下绘制下面三个函数在t ∈[0，4π]的图象。
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是一种强大的工程计算软件，在工程领域有着广泛的应用。
对于学习Matlab程序设计来说，通过例题学习是非常有效的方法。
通过实际的例题练习，可以帮助学生更深入地理解Matlab的基本原理和应用技巧。
本文将介绍几个常见的Matlab程序设计例题，并进行详细的分析和讲解，希望能够帮助读者加深对Matlab程序设计的理解。
二、例题一：矩阵运算题目：编写一个Matlab程序，实现两个任意大小的矩阵相加的运算。
1.我们需要定义两个任意大小的矩阵A和矩阵B，可以通过rand函数生成随机矩阵，也可以手动输入矩阵的元素。
2.编写Matlab程序，使用矩阵的加法运算符“+”进行矩阵相加操作。
3.输出矩阵相加的结果，并进行验证和检查。
通过这个例题，可以加深对Matlab矩阵运算的理解，包括矩阵的定义、矩阵相加的操作和结果输出的方法。
三、例题二：函数绘图题目：编写一个Matlab程序，绘制sin函数的曲线图像。
1.我们需要定义sin函数的自变量范围和取值步长，可以使用linspace函数生成自变量的取值范围。
2.编写Matlab程序，调用sin函数计算自变量对应的函数值。
3.接下来，使用plot函数绘制sin函数的曲线图像，并设置图像的标题、坐标轴标签等其他参数。
4.输出绘制好的sin函数曲线图像，并进行观察和分析。
通过这个例题，可以了解Matlab中函数绘图的基本方法，包括自变量的定义、函数值的计算和图像的绘制等操作。
四、例题三：数值积分题目：编写一个Matlab程序，实现对指定函数在指定区间上的数值积分。
1.我们需要定义积分的目标函数和积分区间，在Matlab中可以使用function定义目标函数，也可以手动输入目标函数和积分区间。
2.编写Matlab程序，调用Matlab内置的数值积分函数（如quad、quadl等）进行数值积分的计算。
3.输出数值积分的结果，并与解析解进行比较和分析。
通过这个例题，可以了解Matlab中数值积分的计算方法，包括目标函数的定义、积分区间的设定和数值积分函数的调用。
matLAB经典例题及答案
一.对以下数据分别作二次,三次多项式拟合,并画出图形.x=1:16;y=[4,6.4,8,8.4,9.28,9.5,9.7,9.86,10,10.2,10.32,10.42,10.5, 10.55,10.58,10.6];答：程序如下（1）x=(1:16);y=erf(x);p=polyfit(x,y,2);f=polyval(p,x);plot(x,y,x,f);结果p=-0.00100.02020.9096（2）y=[4,6.4,8,8.4,9.28,9.5,9.7,9.86,10,10.2,10.32,10.42,10.5, 10.55,10.58,10.6];y=erf(x);p=polyfit(x,y,3)f=polyval(p,x);plot(x,y,x,f)结果P=0.0002-0.00710.06280.8404二.在[0,4pi]画sin(x),cos(x)(在同一个图象中);其中cos(x)图象用红色小圆圈画.并在函数图上标注“y=sin(x)”,“y=cos(x)”,x轴,y轴,标题为“正弦余弦函数图象”.答：程序如下x=[0:720]*pi/180;plot(x,sin(x),x,cos(x),'ro');x=[2.5;7];y=[0;0];s=['y=sin(x)';'y=cos(x)'];text(x,y,s);xlabel('正弦余弦函数图象'),ylabel('正弦余弦函数图象')图形如下三.选择一个单自由度线性振动系统模型，自定质量、弹簧刚度、阻尼、激振力等一组参数，分别编程（m 文件）计算自由和强迫振动时的响应，并画出振动曲线图。
（要求画出该单自由度线性振动系统模型图）其中质量为m=1000kg,弹性刚度k=48020N/m,阻尼c=1960N.s/m,激振力f(t)=0.阻尼比ζ的程序p=1960/(2*sqrt(48020*1000))求得p=0.1414而p为阻尼比ζ强迫振动时的响应程序g =tf([-101],[48020048020*1.9848020]);bode(g)图形g =tf([001],[0001]);bode(g)振动曲线图程序：函数文件function dx =rigid(t,x)dx =zeros(2,1);dx(1)=x(2);dx(2)=(-48020*x(1)-1960*x(2))/1000;命令文件options =odeset('RelTol',1e-4,'AbsTol',[1e-41e-4]);[T,X]=ode45(@rigid,[012],[11],options);plot(T,X(:,1),'-')其图形如下024681012-6-5-4-3-2-11234单自由度线性强迫振动系统模型图其中质量为m=1000kg,弹性刚度k=48020N/m,阻尼c=1960N.s/m，f(t)=cos(3*pi*t)振动曲线图程序：函数文件function dx=rigid(t,x)dx=zeros(2,1);dx(1)=x(2);dx(2)=(-48020*x(1)-1960*x(2))/1000+cos(3*pi*t);命令文件options=odeset('RelTol',1e-4,'AbsTol',[1e-41e-4]);[T,X]=ode45(@rigid,[020],[11],options);plot(T,X(:,1),'-')力等一组参数，建立Simulink仿真模型框图进行仿真分析。
matlab经典编程例题
以下各题均要求编程实现，并将程序贴在题目下方。
1．从键盘输入任意个正整数，以0结束，输出那些正整数中的素数。
clc;clear;zzs(1)=input('请输入正整数：');k=1;n=0;%素数个数while zzs(k)~=0flag=0;%是否是素数，是则为1for yz=2:sqrt(zzs(k))%因子从2至此数平方根if mod(zzs(k),yz)==0flag=1;break;%非素数跳出循环endendif flag==0&zzs(k)>1%忽略0和1的素数n=n+1;sus(n)=zzs(k);endk=k+1;zzs(k)=input('请输入正整数：');enddisp(['你共输入了' num2str(k-1) '个正整数。
它们是：'])disp(zzs(1:k-1))%不显示最后一个数0if n==0disp('这些数中没有素数！')%无素数时显示elsedisp('其中的素数是：')disp(sus)end2．若某数等于其所有因子（不含这个数本身）的和，则称其为完全数。
编程求10000以内所有的完全数。
clc;clear;wq=[];%完全数赋空数组for ii=2:10000yz=[];%ii 的因子赋空数组for jj=2:ii/2 %从2到ii/2考察是否为ii 的因子if mod(ii,jj)==0yz=[yz jj];%因子数组扩展，加上jjendendif ii==sum(yz)+1wq=[wq ii];%完全数数组扩展，加上iiendenddisp(['10000以内的完全数为：' num2str(wq)])%输出3．下列这组数据是美国1900—2000年人口的近似值（单位：百万）。
（1） 若.2c bt at y t y ++=的经验公式为与试编写程序计算出上式中的a 、b 、c;（2） 若.bt ae y t y =的经验公式为与试编写程序计算出上式中的a 、b;（3） 在一个坐标系下，画出数表中的散点图（红色五角星），c bx ax y ++=2中拟合曲线图(蓝色实心线),以及.bt ae y = (黑色点划线)。
MATLAB例题考试及答案
例1.1 分别绘制函数和的曲线。
x=-2*pi:pi/180:2*pi;plot(x,2.^(-abs(x)),':',x,sin(x));例1.2 求方程2x5-3x3 +71x2-9x+13=0的全部根。
p=[2,0,-3,71,-9,13];x=roots(p)例1.3 求解线性方程组。
a=[2,3,-1;8,2,3;45,3,9];b=[2;4;23];x=inv(a)*b例1.4 求积分quad('x.*log(1+x)',0,1)例2.2 利用M文件建立MYMAT矩阵。
(1)启动有关编辑程序或MATLAB文本编辑器(见第4章)，并输入待建矩阵：MYMAT=[101,102,103,104,105,106,107,108,109;201,202,203,204,205,206,207,208,209;301,302,303,304,305,306,307,308,309]例2.3 建立5阶方阵A，判断A的元素是否能被3整除。
A =[24,35,13,22,63;23,39,47,80,80; ...90,41,80,29,10;45,57,85,62,21;37,19,31,88,76]P=rem(A,3)==0例2.5 建立矩阵A，然后找出在[10，20]区间的元素的位置。
(1) 建立矩阵A。
A=[4,15,-45,10,6;56,0,17,-45,0](2) 找出大于4的元素的位置。
find(A>=10 & A<=20)ans =367例2.6 建立一个字符串向量，然后对该向量做如下处理：(1)取第1~5个字符组成的子字符串。
(2)将字符串倒过来重新排列。
(3)将字符串中的小写字母变成相应的大写字母，其余字符不变。
(4)统计字符串中小写字母的个数。
命令如下：ch='ABc123d4e56Fg9';subch=ch(1:5)subch =ABc12revch=ch(end:-1:1)revch =9gF65e4d321cBAk=find(ch>='a'&ch<='z');ch(k)=ch(k)-('a'-'A');char(ch)ans =ABC123D4E56FG9length(k)ans =4例3.2 建立随机矩阵：(1) 在区间[20,50]内均匀分布的5阶随机矩阵。
MATLAB数学实验100例题解
一元函数微分学实验1 一元函数的图形（基础实验）实验目的 通过图形加深对函数及其性质的认识与理解， 掌握运用函数的图形来观察和分析 函数的有关特性与变化趋势的方法，建立数形结合的思想; 掌握用Matlab 作平面曲线图性的方法与技巧。
初等函数的图形2 作出函数x y tan =和x y cot =的图形观察其周期性和变化趋势。
解：程序代码：>〉 x=linspace （0,2＊pi,600）; t=sin （x)。
/(cos （x ）+eps ）;plot(x ，t)；title （’tan （x ）')；axis (［0，2*pi ，-50,50]）; 图象:程序代码: 〉〉 x=linspace （0,2*pi,100); ct=cos （x)。
/（sin(x)+eps ）； plot(x,ct ）；title(’cot(x)')；axis (［0,2＊pi ，—50,50])； 图象：cot(x)4在区间]1,1[-画出函数xy 1sin =的图形。
解：程序代码:>> x=linspace （-1，1，10000）;y=sin(1。
/x ）; plot （x,y ）; axis （［-1，1，—2,2］） 图象：二维参数方程作图6画出参数方程⎩⎨⎧==t t t y tt t x 3cos sin )(5cos cos )(的图形：解：程序代码:>〉 t=linspace(0,2＊pi,100)； plot(cos(t ）.＊cos （5＊t ），sin(t ）。
*cos(3*t)); 图象：极坐标方程作图8 作出极坐标方程为10/t e r =的对数螺线的图形. 解：程序代码：〉〉 t=0:0.01:2*pi ； r=exp （t/10);polar(log(t+eps ）,log （r+eps)); 图象：90270分段函数作图10 作出符号函数x y sgn =的图形。
matlab十个简单案例编写
matlab十个简单案例编写1. 求解线性方程组线性方程组是数学中常见的问题之一，而MATLAB提供了用于求解线性方程组的函数。
例如，我们可以使用"linsolve"函数来求解以下线性方程组：2x + 3y = 74x - 2y = 2代码如下所示：A = [2, 3; 4, -2];B = [7; 2];X = linsolve(A, B);disp(X);解释：上述代码定义了一个2x2的矩阵A和一个2x1的矩阵B，分别表示线性方程组的系数矩阵和常数向量。
然后，使用linsolve函数求解线性方程组，结果存储在X中，并通过disp函数打印出来。
运行代码后，可以得到x=2和y=1的解。
2. 求解非线性方程除了线性方程组外，MATLAB还可以用于求解非线性方程。
例如，我们可以使用"fzero"函数求解以下非线性方程：x^2 + 2x - 3 = 0代码如下所示：fun = @(x) x^2 + 2*x - 3;x0 = 0;x = fzero(fun, x0);disp(x);解释：上述代码定义了一个匿名函数fun，表示非线性方程。
然后，使用fzero函数传入fun和初始值x0来求解非线性方程的根，并通过disp函数打印出来。
运行代码后，可以得到x=1的解。
3. 绘制函数图像MATLAB提供了强大的绘图功能，可以帮助我们可视化函数的形状和特征。
例如，我们可以使用"plot"函数绘制以下函数的图像：y = cos(x)代码如下所示：x = linspace(0, 2*pi, 100);y = cos(x);plot(x, y);解释：上述代码首先使用linspace函数生成一个从0到2π的100个等间距点的向量x，然后计算对应的cos值，并存储在向量y中。
最后，使用plot函数将x和y作为横纵坐标绘制出函数图像。
运行代码后，可以看到cos函数的周期性波动图像。
MATLAB经典例题
subplot(2,1,2);
plot(f,P(1:256));
title('频谱');
设计一个合适的FIR滤波器，100Hz-200Hz带通，要求滤波性能要好。显示其频率响应
生成一个采样频率2000，由50Hz、120Hz、200Hz正弦信号以及噪声信号组成。计算其功率谱密度并显示原始数据波形以及频谱
xlabel('频率f');ylabel('功率密度p');
title('滤波后功率谱密度');
function t3()
b=fir1(30,[0.2 0.4],'stop');
figure(1);
freqz(b,1,512);
title('频率响应');
t=0:0.001:1;
x=sin(t*2*pi*350)+sin(t*2*pi*10);
f=1000*(0:255)/512;
subplot(2,1,2);
plot(f,P(1:256));
title('频谱');
function t3()
fid=fopen('test.dat','r');
[signal,number]=fread(fid,501,'float32');
fclose(fid);
function t1()
t=[0:0.01*pi:2*pi];
y1=t.*t;y2=-1*y1;y3=y1.*sin(t);
plot(t,y1,'r',t,y2,'b',t,y3,'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编程经典例题
以下是一些经典的 MATLAB 编程例题：
1. 编写一个程序，计算一个数列的前 N 个斐波那契数（斐波那契数列是指前两个数为 1，后续每个数是前两个数之和）。
2. 编写一个程序，计算一个数的阶乘。
3. 编写一个程序，计算两个矩阵的乘积。
4. 编写一个程序，找到一个数组中的最大元素。
5. 编写一个程序，检查一个字符串是否是回文。
6. 编写一个程序，为给定的一组数据计算均值、中位数和标准差。
7. 编写一个程序，将一个字符串中的所有元音字母替换成大写字母。
8. 编写一个程序，计算一个数列的前 N 个素数（素数是指只能被 1 和自身整除的数）。
9. 编写一个程序，实现冒泡排序算法，对一个数组进行排序。
10. 编写一个程序，计算一个数的平方根。
以上是一些经典的 MATLAB 编程例题，你可以根据自己的实际需要选择其中的一个或几个进行练习和编程。
matlab典型例题
【例】水资源系统规划调度常应用系统分析方法处理，以一个水资源分配问题为例，讨论线性规划问题。
例：有甲、乙两个水库同时给A、B、C三个城市供水，甲水库的日供水量为28万m3/d，乙水库的日供水量为35万m3/d，三个城市的日需水量分别为A≥10万m3/d，B≥15万m3/d，C≥20 万m3/d。
由于水库与各城市的距离不等，输水方式不同，因此单位水费也不同。
各单位水费分别为c11＝2000元/万m3、c12＝3000元/万m3、c13＝4000元/万m3、c21＝4500元/万m3、c22＝3500元/万m3、c23＝3000元/万m3。
试作出在满足对三个城市供水的情况下，输水费用最小的方案。
设甲水库向三城市日供水量分别为x ll、x12、x13，乙水库向三城市日供水量分别为x2l、x22、x23。
建立约束条件：x11 + x21 ≥10x12 + x22 ≥15x13 + x23 ≥20x11 + x12 + x13 ≤28x21 + x22 + x23 ≤35x11，x12，x13，x21，x22，x23，≥0目标函数：fmin=c11x11+c12x12+c13x13+c21x21+c22x22+c23x23这样的问题单纯求解是非常繁琐的，而MLTLAB求解是十分简单的，只要在命令行输入：》f=[0.2 0.3 0.4 0.45 0.35 0.3]’；》A=[-1 0 0 -1 0 0；0 -1 0 0 -1 0；0 0 -1 0 0 -1；1 1 1 0 0 0；0 0 0 1 1 1]；》B=[-10 -15 -20 28 35]；》lb=zeros（6，1）；》[X，Zmin]=linprog(f，A，B，[]，[]，lb,[])最后得出x=[10 15 0 0 0 20]万m3，Zmin=12.5万元。
matlab经典例题
4 编写函数文件f17.m实现求一个向量所有元素的和（不能用sum）。
再编写一个测试脚本t17.m，随机生成一个4行5列的矩阵A，调用f17求A每列的和function sum=f17(A)a=0;b=size(A);for i=1:b(2)a=a+A(i);endsum=a;A=rand(4,5)sum=[];for j=1:5b=A(1:end,j);b=b';m=f17(b);sum=[sum,m];enddisp(sum);2 有1020个西瓜，第一天卖一半多两个，以后每天卖剩下的一半多两个，问几天后可以卖完？编写求出天数的脚本文件.day=0;x1=1020;while(x1>=0)b=round(x1/2);x2=x1-(b+2);x1=x2;day=day+1;endfprintf('%d\n',day);case {1,3,5,7,8,10,12}day_of_year = day_of_year + 31;case {4,6,9,11}day_of_year= day_of_year + 30;case 2day_of_year = day_of_year + 28 + leap_day;endendfprintf('The date %2d/%2d/%4d is day of year %d.\n', month, day, year,day_of_year);例2,商场购物，100件以下，不优惠，100~199件95折，200~399件90折，400~799件85折，800~1499件80折，1500件以上，75折。
输入所购货物的单价、件数，求实际付款数目。
编写脚本文件，用if语句switch语句实现。
%ifclc;price=input('Input price:'); num=input('Input num:'); if num>=1500d=0.75;elseif num>=800d=0.8;elseif num>=400d=0.85;elseif num>=200d=0.90;elseif num>=100d=0.95;elsed=1;endsum=num*price*d%switchclc;price=input('Input price:'); num=input('Input num:'); sel=fix(num/100);switch selcase 0d=1;case 1d=0.95;case {2,3}d=0.90;case {4,5,6,7}d=0.85;case {8,9,10,11,12,13,14}d=0.80;otherwised=0.75;endsum=num*price*d例3,用1、2、3、4共4个不同的个数字，能组成多少个互不相同且无重复数字的三位数？都是多少？clc;p=[];count=0;for i=1:4for j=1:4for k=1:4if i~=j & i~=k & j~=kn=i*100+j*10+k;p=[p,n];count=count+1;endendendenddisp('互不重复的三位数分别是：');pdisp(['这样的三位数共有',num2str(count),'个'])例4,有1020个桃子，第一天卖一半多两个，以后每天卖剩下的一半多两个，问几天后可以卖完？编写脚本文件求出天数。
MATLAB数学实验100例题解
一元函数微分学实验1 一元函数的图形（基础实验）实验目的 通过图形加深对函数及其性质的认识与理解, 掌握运用函数的图形来观察和分析 函数的有关特性与变化趋势的方法,建立数形结合的思想; 掌握用Matlab 作平面曲线图性的方法与技巧.初等函数的图形2 作出函数x y tan =和x y cot =的图形观察其周期性和变化趋势.解：程序代码：>> x=linspace(0,2*pi,600);t=sin(x)./(cos(x)+eps);plot(x,t);title('tan(x)');axis ([0,2*pi,-50,50]);图象：程序代码：>> x=linspace(0,2*pi,100);ct=cos(x)./(sin(x)+eps);plot(x,ct);title('cot(x)');axis ([0,2*pi,-50,50]);图象：4在区间]1,1[-画出函数xy 1sin =的图形. 解：程序代码：>> x=linspace(-1,1,10000);y=sin(1./x);plot(x,y);axis([-1,1,-2,2])图象：二维参数方程作图6画出参数方程⎩⎨⎧==t t t y t t t x 3cos sin )(5cos cos )(的图形： 解：程序代码：>> t=linspace(0,2*pi,100);plot(cos(t).*cos(5*t),sin(t).*cos(3*t));图象：极坐标方程作图8 作出极坐标方程为10/t e r =的对数螺线的图形.解：程序代码：>> t=0::2*pi;r=exp(t/10);polar(log(t+eps),log(r+eps));图象：分段函数作图10 作出符号函数x y sgn =的图形.解：程序代码：>> x=linspace(-100,100,10000);y=sign(x);plot(x,y);axis([-100 100 -2 2]);函数性质的研究12研究函数)3(log 3)(35x e x x f x -++=在区间]2,2[-上图形的特征.解：程序代码：>> x=linspace(-2,2,10000);y=x.^5+3*exp(x)+log(3-x)/log(3);plot(x,y);图象：实验2 极限与连续（基础实验）实验目的 通过计算与作图, 从直观上揭示极限的本质,加深对极限概念的理解. 掌握用 Matlab 画散点图, 以及计算极限的方法. 深入理解函数连续的概念,熟悉几种间断点的图形 特征,理解闭区间上连续函数的几个重要性质.作散点图14分别画出坐标为)10,,2,1(),4,(),,(3222Λ=+i i i i i i 的散点图, 并画出折线图.解：散点图程序代码：>> i=1:10;plot(i,i.^2,'.')或：>> x=1:10;y=x.^2;for i=1:10;plot(x(i),y(i),'r')hold onend折线图程序代码：>> i=1:10;plot(i,i.^2,'-x')程序代码：>> i=1:10;plot(i.^2,4*(i.^2)+i.^3,'.')>> i=1:10;plot(i.^2,4*(i.^2)+i.^3,'-x')数列极限的概念16通过动画观察当∞→n 时数列21n a n =的变化趋势.解：程序代码：>> n=1:100;an=(n.^2);n=1:100;an=1./(n.^2);n=1:100;an=1./(n.^2);for i=1:100plot(n(1:i),an(1:i)),axis([0,100,0,1])pauseend图象：函数的极限18在区间]4,4[-上作出函数xx x x x f --=339)(的图形, 并研究 )(lim x f x ∞→ 和 ).(lim 1x f x → 解：作出函数x x x x x f --=339)(在区间]4,4[-上的图形 >> x=-4::4;y=(x.^3-9*x)./(x.^3-x+eps);plot(x,y)从图上看，()f x 在x →1与x →∞时极限为0两个重要极限20计算极限 ⎪⎭⎫ ⎝⎛+→x x x x x sin 11sin lim )1(0 x x e x 2lim )2(+∞→ 30sin tan lim )3(xx x x -→ x x x 0lim )4(+→ x x x ln cot ln lim )5(0+→ x x x ln lim )6(20+→ xx x x x x sin cos sin lim )7(20-→ 125523lim )8(323+++-∞→x x x x x x x x e e x x x sin 2lim )9(0----→ x x x x cos 110sin lim )10(-→⎪⎭⎫ ⎝⎛解：（1）>> limit(x*sin(1/x)+1/x*sin(x))ans =1(2) >> limit(x^2/exp(x),inf)ans = 0(3) >> limit((tan(x)-sin(8))/x^3)ans =NaN(4) >> limit(x^x,x,0,'right')ans =1(5) >> limit(log(cot(x))/log(x),x,0,'right')ans =-1(6) >> limit(x^2*log(x),x,0,'right')ans =0(7) >> limit((sin(x)-x.*cos(x))./(x.^2.*sin(x)),x,0)ans =1/3(8) >> limit((3*x.^3-2*x.^2+5)/(5*x.^3+2*+1),x,inf)ans =3/5(9) >> limit((exp(x)-exp(-x)-2*x)./(x-sin(x)))ans =2(10) >> limit((sin(x)/x).^(1/(1-cos(x))))ans =exp(-1/3)实验3 导数（基础实验）实验目的 深入理解导数与微分的概念, 导数的几何意义. 掌握用Matlab 求导数与高 阶导数的方法. 深入理解和掌握求隐函数的导数, 以及求由参数方程定义的函数的导数的方法. 导数概念与导数的几何意义22作函数71232)(23+-+=x x x x f 的图形和在1-=x 处的切线.解：作函数71232)(23+-+=x x x x f 的图形程序代码：>> syms x;>> y=2*x^3+3*x^2-12*x+7;>> diff(y)ans =6*x^2+6*x-12>> syms x;y=2*x^3+3*x^2-12*x+7;>> f=diff(y)f =6*x^2+6*x-12>> x=-1;f1=6*x^2+6*x-12f1 =-12>> f2=2*x^3+3*x^2-12*x+7f2 =20>> x=linspace(-10,10,1000);y1=2*x.^3+3*x.^2-12*x+7;y2=-12*(x+1)+20;plot(x,y1,'r',x,y2,'g')求函数的导数与微分24求函数bx ax x f cos sin )(=的一阶导数. 并求.1⎪⎭⎫ ⎝⎛+'b a f 解：求函数bx ax x f cos sin )(=的一阶导数程序代码：>> syms a b x y;y= sin(a*x)*cos(b*x);D1=diff(y,x,1)答案：D1 =cos(a*x)*a*cos(b*x)-sin(a*x)*sin(b*x)*b 求.1⎪⎭⎫ ⎝⎛+'b a f 程序代码：>> x=1/(a+b);>> cos(a*x)*a*cos(b*x)-sin(a*x)*sin(b*x)*b答案：ans =cos(a/(a+b))*a*cos(b/(a+b))-sin(a/(a+b))*sin(b/(a+b))*b拉格朗日中值定理26对函数),2)(1()(--=x x x x f 观察罗尔定理的几何意义.(1) 画出)(x f y =与)(x f '的图形, 并求出1x 与.2x解：程序代码：>> syms x;f=x*(x-1)*(x-2);f1=diff(f)f1 =(x-1)*(x-2)+x*(x-2)+x*(x-1)>> solve(f1)ans =1+1/3*3^(1/2)1-1/3*3^(1/2)>> x=linspace(-10,10,1000);y1=x.*(x-1).*(x-2);y2 =(x-1).*(x-2)+x.*(x-2)+x.*(x-1);plot(x,y1,x,y2)(2)画出)(x f y =及其在点))(,(11x f x 与))(,(22x f x 处的切线.程序代码：>> syms x;>> f=x*(x-1)*(x-2);>> f1=diff(f)f1 =(x-1)*(x-2)+x*(x-2)+x*(x-1)>> solve(f1)ans =1+1/3*3^(1/2)1-1/3*3^(1/2)>> x=linspace(-3,3,1000);>> y1=x.*(x-1).*(x-2);>> y2 =(x-1).*(x-2)+x.*(x-2)+x.*(x-1);>> plot(x,y1,x,y2)>> hold on>> x=1+1/3*3^(1/2);>> yx1=x*(x-1)*(x-2)yx1 =>> x=1-1/3*3^(1/2);>> yx2=x*(x-1)*(x-2)yx2 =x=linspace(-3,3,1000);yx1 =*x.^0;yx2 =*x.^0;plot(x,yx1,x,yx2)28求下列函数的导数: (1) 31+=x e y ;解：程序代码：>> syms x y;y=exp((x+1)^3);D1=diff(y,1)答案：D1 =3*(x+1)^2*exp((x+1)^3) (2) )]42ln[tan(π+=x y ; 解：程序代码：>> syms x;y=log(tan(x/2+pi/4));D1=diff(y,1)答案：D1 =(1/2+1/2*tan(1/2*x+1/4*pi)^2)/tan(1/2*x+1/4*pi) (3) x x y sin ln cot 212+=; 解：程序代码：>> syms x;y=1/2*(cot(x))^2+log(sin(x));D1=diff(y,1)答案：D1 =cot(x)*(-1-cot(x)^2)+cos(x)/sin(x) (4) xy 2arctan 21=. 解：程序代码：>> y=sqrt(2)*atan(sqrt(2)/x);>> D1=diff(y,1)答案：D1 =-2/x^2/(1+2/x^2)一元函数积分学与空间图形的画法实验4 一元函数积分学(基础实验)实验目的 掌握用Matlab 计算不定积分与定积分的方法. 通过作图和观察, 深入理解 定积分的概念和思想方法. 初步了解定积分的近似计算方法. 理解变上限积分的概念. 提高应用定积分解决各种问题的能力.不定积分计算30求.)1(532⎰-dx x x解：程序代码：>> syms x y;>> y=x^2*(1-x^3)^5;>> R=int(y,x)答案：R =-1/18*x^18+1/3*x^15-5/6*x^12+10/9*x^9-5/6*x^6+1/3*x^332求.arctan 2⎰xdx x解：程序代码：>> syms x y;>> y=x^2*atan(x);>> R=int(y,x)答案：R =1/3*x^3*atan(x)-1/6*x^2+1/6*log(x^2+1)定积分计算 34 求.)(102⎰-dx x x解：程序代码：>> syms x y;>> y=x-x^2;>> R=int(y,x,0,1)答案： R =1/6变上限积分36 画出变上限函数⎰x dt t t 02sin 及其导函数的图形.解：程序代码：>> syms x y t;>> y=t*sin(t^2);>> R=int(y,x,0,x)答案：R =t*sin(t^2)*x程序代码：>> DR=diff(R,x,1)答案：DR =t*sin(t^2)实验5 空间图形的画法(基础实验)实验目的 掌握用Matlab 绘制空间曲面和曲线的方法. 熟悉常用空间曲线和空间曲面 的图形特征,通过作图和观察, 提高空间想像能力. 深入理解二次曲面方程及其图形.一般二元函数作图38作出函数2214y x z ++=的图形. 解：程序代码：>> x=linspace(-5,5,500);[x,y]=meshgrid(x);z=4./(1+x.^2+y.^2);mesh(x,y,z);xlabel('x-axis'),ylabel('y-axis'),zlabel('z-axis');title('function')40作出函数)94cos(22y x z +=的图形.解：程序代码：>> x=-10::10;[x,y]=meshgrid(x);z=cos(4*x.^2+9*y.^2);mesh(x,y,z);xlabel('x-axis'),ylabel('y-axis'),zlabel('z-axis');title('function')讨论：坐标轴选取范围不同时，图形差异很大，对本题尤为明显，如右图为坐标轴[-1，1]二次曲面42作出单叶双曲面1941222=-+z y x 的图形.(曲面的参数方程为 ,tan 3,cos sec 2,sin sec u z v u y v u x === (.20,2/2/πππ≤≤<<-v u ))解：程序代码：>> v=0:pi/100:2*pi;>> u=-pi/2:pi/100:pi/2;>> [U,V]=meshgrid(u,v);>> x=sec(U).*sin(V);>> y=2*sec(U).*cos(V);>> z=3*tan(U);>> surf(x,y,z)44 可以证明: 函数xy z =的图形是双曲抛物面. 在区域22,22≤≤-≤≤-y x 上作出它的图形.解：程序代码：>> x=-2::2;[x,y]=meshgrid(x);>> z=x.*y;>> mesh(x,y,z);46 画出参数曲面]2,001.0[],4,0[)5/2/ln(tan cos sin sin sin cos ∈∈⎪⎩⎪⎨⎧++===v u u v v z vu y v u x π 的图形.解：程序代码： >> v=::2;>> u=0:pi/100:4*pi;>> [U,V]=meshgrid(u,v);>> x=cos(U).*sin(V);>> y=sin(U).*sin(V);>> z=cos(V)+log(tan(V/2)+U/5);>> mesh(x,y,z);空间曲线48 作出空间曲线)60(2,sin ,cos π≤≤===t t z t t y t t x 的图形.解：程序代码：>> syms t;ezplot3(t*cos(t),t*sin(t),2*t,[0,6*pi])50绘制参数曲线 ⎪⎪⎩⎪⎪⎨⎧=+==t z t y t x arctan 211cos 2的图形. 解：程序代码：>> t=-2*pi:pi/100:2*pi;x=cos(t).*cos(t);y=1./(1+2*t);z=atan(t);plot3(x,y,z);grid;xlabel('x'),ylabel('y'),zlabel('z')多元函数微积分实验6 多元函数微分学（基础实验）实验目的 掌握利用Matlab 计算多元函数偏导数和全微分的方法, 掌握计算二元 函数极值和条件极值的方法. 理解和掌握曲面的切平面的作法. 通过作图和观察, 理解二元 函数的性质、方向导数、梯度和等高线的概念.求多元函数的偏导数与全微分52设),(cos )sin(2xy xy z +=求.,,,222yx z x z y z x z ∂∂∂∂∂∂∂∂∂ 解：程序代码：>> syms x y;S=sin(x*y)+(cos(x*y))^2;D1=diff(S,'x',1);D2=diff(S,'y',1);D3=diff(S,'x',2);D4=diff(S,'y',2);D1,D2,D3,D4答案： D1 = cos(x*y)*y-2*cos(x*y)*sin(x*y)*yD2 = cos(x*y)*x-2*cos(x*y)*sin(x*y)*xD3 =-sin(x*y)*y^2+2*sin(x*y)^2*y^2-2*cos(x*y)^2*y^2D4 = -sin(x*y)*x^2+2*sin(x*y)^2*x^2-2*cos(x*y)^2*x^2实验7 多元函数积分学（基础实验）实验目的掌握用Matlab 计算二重积分与三重积分的方法; 深入理解曲线积分、曲面积分的 概念和计算方法. 提高应用重积分和曲线、曲面积分解决各种问题的能力.计算重积分54计算,2dxdy xyD ⎰⎰ 其中D 为由,,2y x y x ==+ 2=y 所围成的有界区域.解：程序代码：>> syms x y;int(int(x*y^2,x,2-y,sqrt(y)),y,1,2)答案：ans =193/120重积分的应用56求旋转抛物面224y x z --=在Oxy 平面上部的面积.S解：程序代码：>> int(2*pi*r,r,0,2)答案： ans =4*pi无穷级数与微分方程实验8 无穷级数（基础实验）实验目的观察无穷级数部分和的变化趋势,进一步理解级数的审敛法以及幂级数部分和对函数的 逼近. 掌握用Matlab 求无穷级数的和, 求幂级数的收敛域, 展开函数为幂级数以及展 开周期函数为傅里叶级数的方法.数项级数58(1) 观察级数∑∞=121n n 的部分和序列的变化趋势. 解：程序代码：for i=1:100 s=0;for n=1:i s=s+1/n^2;endplot(i,s,'.');hold on;end(2) 观察级数∑∞=11n n 的部分和序列的变化趋势. >> for i=1:100 s=0;for n=1:i s=s+1/n;endplot(i,s,'.'); hold on;end60 求∑∞=++123841n n n 的值. 解：程序代码：>> syms n;score=symsum(1/(4*n^2+8*n+3),1,inf)答案： score =1/6函数的幂级数展开62求x arctan 的5阶泰勒展开式.>> syms x;>> T5=taylor(atan(x),6)答案：T5 =x-1/3*x^3+1/5*x^5实验9 微分方程（基础实验）实验目的 理解常微分方程解的概念以及积分曲线和方向场的概念,掌握利用Matlab 求微分方程及方程组解的常用命令和方法.求解微分方程64求微分方程 22x xe xy y -=+'的通解.解：程序代码：>> y=dsolve('Dy+2*x*y=x*exp(-x^2)','x')答案：y =(1/2*x^2+C1)*exp(-x^2)66求微分方程x e y y y x 2cos 52=+'-''的通解.解：程序代码：>> y=dsolve('D2y-2*Dy+5*y=exp(x)*cos(2*x)','x')答案： y =exp(x)*sin(2*x)*C2+exp(x)*cos(2*x)*C1+1/4*exp(x)*sin(2*x)*x68求微分方程组⎪⎪⎩⎪⎪⎨⎧=--=++02y x dtdy e y x dt dx t 在初始条件0,100====t t y x 下的特解. 解：程序代码：>> [x,y]=dsolve('Dx+x+2*y-exp(t)','Dy-x-y','x(0)=1','y(0)=0','t')答案： x = cos(t)y = 1/2*sin(t)-1/2*cos(t)+1/2*exp(t)70求解微分方程,)1(122/5+=+-x x y dx dy 并作出积分曲线. 解：程序代码：>> syms x yy=dsolve('Dy-2*y/(x+1)-(x+1)^(5/2)','x')答案：y =(2/3*(x+1)^(3/2)+C1)*(x+1)^2做积分曲线由>> syms x yx=linspace(-5,5,100);C=input('请输入C 的值:');y=(2/3*(x+1).^(3/2)+C).*(x+1).^2;plot(x,y)例如对应有： 请输入C 的值:2 请输入C 的值:20矩阵运算与方程组求解实验10 行列式与矩阵实验目的掌握矩阵的输入方法. 掌握利用Matlab 对矩阵进行转置、加、减、数乘、相乘、乘方等运算, 并能求矩阵的逆矩阵和计算方阵的行列式.矩阵A 的转置函数Transpose[A]72 求矩阵⎪⎪⎪⎪⎪⎭⎫ ⎝⎛411365243271的转置. 解：程序代码： >> A=[1,7,2;3,4,2;5,6,3;1,1,4];>> Sove=A'答案：Sove =1 3 5 17 4 6 12 23 4矩阵线性运算73设,291724,624543⎪⎪⎭⎫ ⎝⎛=⎪⎪⎭⎫ ⎝⎛=B A 求.24,A B B A -+ 解：程序代码：>> A=[3,4,5;4,2,6];B=[4,2,7;1,9,2];S1=A+BS2=4*B-2*A答案：S1 =7 6 125 11 8S2 =10 0 18-4 32 -474设,148530291724,36242543⎪⎪⎪⎪⎪⎭⎫ ⎝⎛=⎪⎪⎭⎫ ⎝⎛=mb ma 求矩阵ma 与mb 的乘积. 解：程序代码：>> ma=[3,4,5,2;4,2,6,3];>> mb=[4,2,7;1,9,2;0,3,5;8,4,1];>> Sove=ma*mb答案：Sove =32 65 5642 56 65矩阵的乘法运算75设,101,530291724⎪⎪⎪⎭⎫ ⎝⎛=⎪⎪⎪⎭⎫ ⎝⎛=B A 求AB 与,A B T 并求.3A解：程序代码：>> A=[4 2 7;1 9 2;0 3 5];B=[1;0;1];>> AB=A*BAB =1135>> BTA=B'*ABTA =4 5 12>> A3=A^3A3 =119 660 555141 932 44454 477 260求方阵的逆76 设,5123641033252312⎪⎪⎪⎪⎪⎭⎫ ⎝⎛=A 求.1-A 解：程序代码：>> A=[2,1,3,2;5,2,3,3;0,1,4,6;3,2,1,5];Y=inv(A)答案：Y =77 设,221331317230,5121435133124403⎪⎪⎪⎪⎪⎭⎫ ⎝⎛=⎪⎪⎪⎪⎪⎭⎫ ⎝⎛=B A 求.1B A - 解：程序代码：>> A=[3 0 4 4 ;2 1 3 3 ;1 5 3 4;1 2 1 5];B=[0 3 2 ;7 1 3;1 3 3 ;1 2 2];Solve=A'*B答案：Solve =16 16 1714 20 2225 26 2830 37 3978 解方程组⎪⎩⎪⎨⎧-=-+=+-=++.2442,63,723z y x z y x z y x解：程序代码：>> A=[3 2 1;1 -1 3;2 4 -4];b=[7 6 -2];>> A\b'答案：ans =求方阵的行列式79 求行列式 .3351110243152113------=D 解：程序代码：>> A=[3,1,-1,2;-5,1,3,-4;2,0,1,-1;1,-5,3,-3];D=det(A)答案：D =4080求.11111111111122222222d d d d c c c c b b b b a a a a D ++++= 解：程序代码：>> syms a b c d;D=[a^2+1/a^2 a 1/a 1;b^2+1/b^2 b 1/b 1;c^2+1/c^2 c 1/c 1;d^2+1/d^2 d 1/d 1]; det(D)答案：ans =-(-c*d^2*b^3+c^2*d*b^3-c^3*d^2*a+c^3*d*a^2*b^4+c*d^2*a^3-c^3*d^2*a*b^4-c^2*d*a^3-c*d^2*b^3*a^4+c^2*d*b^3*a^4+c^3*d^2*b*a^4-c^3*d*b^2*a^4-c^2*d^3*b*a^4+c*d^3*b^2*a ^4+c*d^2*a^3*b^4-c^2*d*a^3*b^4+c^3*d^2*b-c^3*d*b^2-c^2*d^3*b+c*d^3*b^2+c^3*d*a^2+c^2*d^3*a-c*d^3*a^2-b*d^2*a^3+b^2*d*a^3+b^3*d^2*a-b^3*d*a^2-b^2*d^3*a+b*d^3*a^2+b*c^2*a^3-b^2*c*a^3-b^3*c^2*a+b^3*c*a^2+b^2*c^3*a-b*c^3*a^2+c^2*d^3*a*b^4-c*d^3*a^2*b^4-b*d^2*a^3*c^4+b^2*d*a^3*c^4+b^3*d^2*a*c^4-b^3*d*a^2*c^4-b^2*d^3*a*c^4+b*d^3*a^2*c^4+b*c^2*a^3*d^4-b^2*c*a^3*d^4-b^3*c^2*a*d^4+b^3*c*a^2*d^4+b^2*c^3*a*d^4-b*c^3*a^2*d^4)/a^2/c^2/d^2/b^281 计算范德蒙行列式.1111145444342413534333231252423222154321x x x x x x x x x x x x x x x x x x x x 解：程序代码：>> syms x1 x2 x3 x4 x5;>> A=[1,1,1,1,1;x1,x2,x3,x4,x5;x1^2,x2^2,x3^2,x4^2,x5^2;x1^3,x2^3,x3^3,x4^3,x5^3;x1^4,x2^4,x3^4,x4^4,x5^4];>> DC=det(A);>> DS=simple(DC)答案：DS =(-x5+x4)*(x3-x5)*(x3-x4)*(-x5+x2)*(x2-x4)*(x2-x3)*(-x5+x1)*(x1-x4)*(x1-x3)*(x1-x2)82 设矩阵 ,60975738723965110249746273⎪⎪⎪⎪⎪⎪⎭⎫ ⎝⎛----=A 求.),(|,|3A A tr A 解：程序代码：>> A=[3,7,2,6,-4;7,9,4,2,0;11,5,-6,9,3;2,7,-8,3,7;5,7,9,0,-6];>> D=det(A),T=trace(A),A3=A^3答案：D =11592T =3A3=726 2062 944 294 -3581848 3150 26 1516 2281713 2218 31 1006 4041743 984 -451 1222 384801 2666 477 745 -125向量的内积83 求向量}3,2,1{=u 与}0,1,1{-=v 的内积.解：程序代码：>> u=[1 2 3];v=[1 -1 0];solve=dot(u,v)答案：solve =-184设,001001⎪⎪⎪⎭⎫⎝⎛=λλλA 求.10A 一般地?=k A (k 是正整数).解：程序代码：>> syms r;>> A=[r,1,0;0,r,1;0,0,r];>> A^10答案：ans =[ r^10, 10*r^9, 45*r^8][ 0, r^10, 10*r^9][ 0, 0, r^10]85.求⎪⎪⎪⎪⎪⎪⎭⎫ ⎝⎛+++++a a a aa 1111111111111111111111111的逆. 解：程序代码：>> syms aA=[1+a,1,1,1,1;1,1+a,1,1,1;1,1,1+a,1,1;1,1,1,1+a,1;1,1,1,1,1+a];solve=inv(A)答案：solve =[ 1/a*(a+4)/(a+5), -1/a/(a+5), -1/a/(a+5), -1/a/(a+5), -1/a/(a+5)][ -1/a/(a+5), 1/a*(a+4)/(a+5), -1/a/(a+5), -1/a/(a+5), -1/a/(a+5)][ -1/a/(a+5), -1/a/(a+5), 1/a*(a+4)/(a+5), -1/a/(a+5), -1/a/(a+5)][ -1/a/(a+5), -1/a/(a+5), -1/a/(a+5), 1/a*(a+4)/(a+5), -1/a/(a+5)][ -1/a/(a+5), -1/a/(a+5), -1/a/(a+5), -1/a/(a+5), 1/a*(a+4)/(a+5)]实验11 矩阵的秩与向量组的极大无关组实验目的 学习利用Matlab 求矩阵的秩,作矩阵的初等行变换; 求向量组的秩与极大无关组.求矩阵的秩86 设,815073*********⎪⎪⎪⎭⎫ ⎝⎛-------=M 求矩阵M 的秩.解：程序代码：>> M=[3,2,-1,-3,-2;2,-1,3,1,-3;7,0,5,-1,-8];R=rank(M)答案：R=2向量组的秩87求向量组)0,3,0,2(),2,5,4,0(),1,1,2,1(231=--=-=ααα的秩.解：程序代码：>> A=[1,2,-1,1;0,-4,5,-2;2,0,3,0];R=rank(A)答案：R =288向量组)7,5,1,3(),5,4,3,1(),1,1,1,1(),3,2,1,1(4321==-==αααα是否线性相关? 解：由>> A=[1 1 2 3;1 -1 1 1;1 3 4 5;3 1 5 7];rank(A)ans = 3即rank(A)=3 小于阶数489向量组)3,1,1(),2,1,3(),7,2,2(321=-==ααα是否线性相关?解：由>> A3=[2,2,7;3,-1,2;1,1,3];R=rank(A3)得 R = 3即rank(A3)=3 等于阶数3故向量组线性无关。
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. 请编写一个程序，求出一个函数的积分。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
例2.1>> muw0=1.785;>> a=0.03368;>> b=0.000221;>> t=0:20:80;>> muw=muw0./(1+a*t+b*t.^2)例2.2 数值数组和字符串的转换>> a=[1:5];>> b=num2str(a);>> a*2ans =2 4 6 8 10>> b*2ans =98 64 64 100 64 64 102 64 64 104 64 64 106例2.9比较左除和右除求解恰定方程>> rand('seed',12);>> a=rand(100)+1.e8;>> x=ones(100,1);>> b=a*x;>> cond(a)ans =5.0482e+011>> tic;x1=b'/a;t1=toct1 =0.4711>> er1=norm(x-x1')er1 =139.8326>> re1=norm(a*x1'-b)/norm(b)re1 =4.3095e-009>> tic;x1=a\b;t1=toct1 =0.0231>> tic;x1=a\b;t1=toct1 =0.0011>> er2=norm(x-x1)er2 =1.5893e-004>> re1=norm(a*x1-b)/norm(b)re1 =4.5257e-016例2.14：计算矩阵的指数>> b=magic(3);>> expm(b)ans =1.0e+006 *1.0898 1.0896 1.08971.0896 1.0897 1.08971.0896 1.0897 1.0897 例2.18：特征值条件数>> a=[-149 -50 -154;537 180 546; -27 -9 -25]a =-149 -50 -154537 180 546-27 -9 -25>> [V,D,s]=condeig(a)V =0.3162 -0.4041 -0.1391-0.9487 0.9091 0.9740-0.0000 0.1010 -0.1789D =1.0000 0 00 2.0000 00 0 3.0000例2.41 5阶多项式在【0,2pi】最小二乘拟合>> x=0:pi/20:pi/2;>> y=sin(x);>> a=polyfit(x,y,5);>> x1=0:pi/30:pi*2;>> y1=sin(x1);>> y2=a(1)*x1.^5+a(2)*x1.^4+a(3)*x1.^3+a(4)*x1.^2+a(5)*x1+a(6); >> plot(x1,y1,'b-',x1,y2,'r*')>> legend('原曲线','拟合曲线')>> axis([0,7,-1.2,4])例3.7 gradient绘制矢量图>> x=0:pi/20:pi/2;>> y=sin(x);>> a=polyfit(x,y,5);>> x1=0:pi/30:pi*2;>> y1=sin(x1);>> y2=a(1)*x1.^5+a(2)*x1.^4+a(3)*x1.^3+a(4)*x1.^2+a(5)*x1+a(6); >> plot(x1,y1,'b-',x1,y2,'r*')>> legend('原曲线','拟合曲线')>> axis([0,7,-1.2,4])>>>> [x,y]=meshgrid(-2:.2:2,-2:.2:2); >> z=x.*exp(-x.^2-y.^2);>> [px,py]=gradient(z,.2,.2);>> contour(z),>> hold on>> quiver(px,py)>> hold off例基本绘图命令rand(100,1);plot(y)例4.1 绘制如图>> x=1:0.1*pi:2*pi;>> y=sin(x);>> z=cos(x);>> plot(x,y,'--k',x,z,'-.rd')例4.5 绘制如图>> x=1:10;>> y=rand(10,1);>> bar(x,y);>> x=0:0.1*pi:2*pi;>> y=x.*sin(x);>> feather(x,y)例 4.6 绘制如图>> lim=[0,2*pi,-1,1];>> fplot('[sin(x),cos(x)]',lim)例4.7绘图如下>> x=[2,4,6,8];>> pie(x,{'math','english','chinese','music'}) 例4.9 绘图如下三维螺旋线>> x=0:pi/50:10*pi;>> y=sin(x);>> x=0:pi/50:10*pi;>> y=sin(x);>> z=cos(x);>> plot3(x,y,z);例4.10 绘图如下。
矩阵三维图>> [x,y]=meshgrid(-2:0.1:2,-2:0.1:2); >> z=x.*exp(-x.^2-y.^2);>> plot3(x,y,z)例4.13绘图如下>> [X,Y]=meshgrid([-4:0.5:4]);>> Z=sqrt(X.^2+Y.^2);>> meshc(Z)例4.19 绘制柱面图>> x=0:pi/20:pi*3;>> r=5+cos(x);>> [a,b,c]=cylinder(r,30);>> mesh(a,b,c)例4.20 地球表面气温分布示意图>> [a,b,c]=sphere(40);>> t=abs(c);>> surf(a,b,c,t);>> axis('equal')>> axis('square')>> colormap('hot')例4.24坐标标注函数应用示意图>> x=1:0.1*pi:2*pi;>> y=sin(x);>> plot(x,y)>> xlabel('x(0-2\pi)','fontweight','bold');>> ylabel('y=sin(x)','fontweight','bold');>> title('正弦函数','fontsize',12,'fontweight','bold','fontname','隶书') >>例4.30 同一张图绘制几个三角函数>> x=0:0.1*pi:2*pi;>> y=sin(x);>> z=cos(x);>> plot(x,y,'-*')>> hold on>> plot(x,z,'-o')>> plot(x,y+z,'-h')>> legend('sin(x)','cos(x)','sin(x)+cos(x)',0)>> hold off例4.31 4个子图中绘制不同的三角函数图>> x=0:0.1*pi:2*pi;>> subplot(2,2,1);>> plot(x,sin(x),'-*');>> title('sin(x)');>> subplot(2,2,2);>> plot(x,cos(x),'-o');>> title('cos(x)');>> subplot(2,2,3);>> plot(x,sin(x).*cos(x),'-x');>> title('sin(x)*cos(x)');>> subplot(2,2,4);>> plot(x,sin(x)+cos(x),'-h');>> title('sin(x)+cos(x)');例7.3正弦曲线插值示例>> x=0:0.1:10;>> y=sin(x);>> xi=0:.25:10;>> yi=interp1(x,y,xi);>> plot(x,y,'o',xi,yi)例7.7 x 0.5 1.0 2.0 2.5 3.0y 1.75 2.45 3.81 4.80 8.00 8.60 y=span{1,x,x^2}，最小二乘法拟合>> x=[0.5 1.0 1.5 2.0 2.5 3.0];>> y=[1.75 2.45 3.81 4.80 8.00 8.60];>> a=polyfit(x,y,2)a =0.4900 1.2501 0.8560>> x1=[0.5:0.05:3.0];>> y1=a(3)+a(2)*x1+a(1)*x1.^2;>> plot(x,y,'*')>> hold on>> plot(x1,y1,'-r')例7.8最小二乘法求y=a+b*x^2的经验公式Xi 19 25 31 38 44Yi 19.0 32.3 49.0 73.3 98.8>> x=[19 25 31 38 44];>> y=[19.0 32.3 49.0 73.3 98.8];>> x1=x.^2x1 =361 625 961 1444 1936 >> x1=[ones(5,1),x1']x1 =1 3611 6251 9611 14441 1936>> ab=x1\y'ab =0.59370.0506>> x0=[19:0.2:44];>> y0=ab(1)+ab(2)*x0.^2;>>>> clf>> plot(x,y,'o')>> hold on>> plot(x0,y0,'-r')例7.10求积分function y=fun(t)y=exp(-0.5*t).*sin(t+pi/6);>> d=pi/1000;>> t=0:d:3*pi;>>>> nt=length(t);>> y=fun(t);>> sc=cumsum(y)*d;>> scf=sc(nt)scf =0.9016>> z=trapz(y)*dz =0.9008例7.12用Newton-cotes公式求积分Fun.mfunction f=fun(x)f=exp(-x/2);quad8('fun',1,3,1e-10)例微分函数>> x=sym('x');>> diff(sin(x^2))ans =2*x*cos(x^2)例题7-44 273-274页fun.mfunction f=fun(x,y)f=-2*y+2*x.^2+2*x;>> [x,y]=ode23('fun',[0,0.5],1);>> x'ans =Columns 1 through 70 0.0400 0.0900 0.1400 0.1900 0.2400 0.2900 Columns 8 through 120.3400 0.3900 0.4400 0.4900 0.5000>> y'ans =Columns 1 through 71.0000 0.9247 0.8434 0.7754 0.7199 0.6764 0.6440 Columns 8 through 120.6222 0.6105 0.6084 0.6154 0.6179例题7-45tic;p1=flops;[x,y]=ode23('fun',[0,0.5],1);p2=flops;t=toc;p=p2-p1;>> bj例题7-46function f=f(x,y)f=[-2 1;988 -999]*y+[2*sin(x);999*(cos(x)-sin(x))];>> ode23('f',[0,10],[2,3]);>> a=[-2 1;998 -999]; %求方程的刚性比>> b1=max(abs(real(eig(a))));>> b2=min(abs(real(eig(a))));>> s=b1/b2s =1000例7-17/18 246页>> a=[0.4096, 0.1234, 0.3678, 0.2943;0.2246, 0.3872, 0.4015, 0.1129;0.3645, 0.1920, 0.3781, 0.0643;0.1784, 0.4002, 0.2786, 0.3927]; >> aa =0.4096 0.1234 0.3678 0.29430.2246 0.3872 0.4015 0.11290.3645 0.1920 0.3781 0.06430.1784 0.4002 0.2786 0.3927>> b=[0.4043 0.1550 0.4240 -0.2557]';>> x=a\bx =-0.1819-1.66302.2172-0.4467265页，例7-39 （非线性方程组的符号解法）g.mfunction y=g(x)y(1)=0.7*sin(x(1))+0.2*cos(x(2));y(2)=0.7*cos(x(1))-0.2*sin(x(2));>> x0=[0.5 0.5];>> fsolve('g',x0)No solution found.fsolve stopped because the problem appears regular as measured by the gradient, but the vector of function values is not near zero as measured by thedefault value of the function tolerance.<stopping criteria details>ans =-0.0493 1.5215307页，例9-21>> x=[0.236 0.238 0.248 0.245 0.243;0.257 0.253 0.255 0.254 0.261;0.258 0.264 0.259 0.267 0.262];>> anova1(x')ans =1.3431e-005308页，例9-22>> a=[58.2000 56.2000 65.3000;52.6000 41.2000 60.8000;49.1000 54.1000 51.6000;42.8000 50.5000 48.4000;60.1000 70.9000 39.2000;58.3000 73.2000 40.7000;75.8000 58.2000 48.7000;71.5000 51.0000 41.4000];>> anova2(a,2)ans =0.0035 0.0260 0.0001例9.23 （309页）。