有趣的MATLAB动画演示程序汇总

MATLAB 动画演示程序汇总1.弹性蹦球演示程序figure(1);%定义函数axis([-5.1,5,-0.05,1.05]);%绘制二维图形hold on;%保持当前图形及轴系所有的特性axis('off');%覆盖坐标刻度,并填充背景%通过填充绘出台阶及两边的挡板fill([4.12,4.22,4.22,4.12],[-0.05,-0.05,1.05,1.05],'y');fill([-5,-3.2,-3.2,-5],[-0.05,-0.05,0,0],'g');fill([-3.2,-2.8,-2.8,-3.2],[-0.05,-0.05,0.2,0.2],'g');fill([-3.2,-1.4,-1.4,-3.2],[0.2,0.2,0.25,0.25],'g');fill([-1.4,-1,-1,-1.4],[0.2,0.2,0.45,0.45],'g');fill([-1.4,0.4,0.4,-1.4],[0.45,0.45,0.5,0.5],'g');fill([0.4,0.8,0.8,0.4],[0.45,0.45,0.7,0.7],'g');fill([0.4,2.0,2.0,0.4],[0.7,0.7,0.75,0.75],'g');fill([2.0,2.3,2.3,2.0],[-0.05,-0.05,0.75,0.75],'g');fill([2.3,4.12,4.12,2.3],[-0.05,-0.05,0,0],'g');%x2=line([-5,5],[0.25,0.25],'color','g','linestyle','-','markersize',50)%设置台阶边框线,颜色,擦试方式%line([-5,5],[0.5,0.5],'color','b','linestyle','-', 'markersize',50)%设置球与地面接触面的颜色,擦试方式%line([-5,5],[0.75,0.75],'color','b','linestyle','-','markersize',50)%设置球与地面接触面的颜色,擦试方式head=line(-5,1,'color','r','linestyle','.','erasemode','xor','markersize',60);%设置小球颜色,大小,线条和擦试方式%body=line(-5,1,'color','b','linestyle','-','erasemode','none'); %描绘轨迹线%设置初始条件while 1t=4;dt=0.001;w=0;dw=0.001;w=0;%设置球弹起的初始位置%设置球弹起的高度while t<=4.12t=dt+t;if w<=1w=dw+w;elsew=-1;endy=(-w*w)+1;set(head,'xdata',t,'ydata',y);%设置球的运动%set(body,'xdata',t,'ydata',y);%描绘轨迹线drawnow;end %结束程序w=0;%设置球弹起的初始位置%设置球弹起的高度while t>=2.11t=t-dt;if w<=1w=dw+w;elsew=-1;endy=(-w*w)+1;set(head,'xdata',t,'ydata',y);%设置球的运动%set(body,'xdata',t,'ydata',y);%描绘轨迹线drawnow;end %结束程序w=0;%设置球弹起的初始位置%设置球弹起的高度while t>=1.11t=t-dt;if w<=1w=dw+w;elsew=-1;endy=(-w*w)/4+1;set(head,'xdata',t,'ydata',y);%设置球的运动%set(body,'xdata',t,'ydata',y);%描绘轨迹线drawnow;end %结束程序w=-0.71;%设置球弹起的初始位置%设置球弹起的高度while t>=-0.62t=t-dt;if w<=1w=dw+w;elsew=-1;endy=(-w*w)/2+1;set(head,'xdata',t,'ydata',y);%设置球的运动%set(body,'xdata',t,'ydata',y);%描绘轨迹线drawnow;end %结束程序w=-0.71;%设置球弹起的初始位置%设置球弹起的高度while t>=-2.31t=t-dt;if w<=1w=dw+w;elsew=-1;endy=(-w*w)/2+0.75;set(head,'xdata',t,'ydata',y);%设置球的运动%set(body,'xdata',t,'ydata',y);%描绘轨迹线drawnow;end %结束程序w=-0.71;%设置球弹起的初始位置%设置球弹起的高度while t>=-4t=t-dt;if w<=1w=dw+w;elsew=-1;endy=(-w*w)/2+0.5;set(head,'xdata',t,'ydata',y);%设置球的运动%set(body,'xdata',t,'ydata',y);%描绘轨迹线drawnow;end %结束程序for i=0:0.01:200%设置延时y=i+2;if(y>199)endendend2.嫦娥奔月演示程序figure('name','嫦娥一号与月亮、地球关系');%设置标题名字s1=[0:.01:2*pi];hold on;axis equal;%建立坐标系axis off % 除掉Axesr1=10;%月亮到地球的平均距离r2=3;%嫦娥一号到月亮的平均距离w1=1;%设置月亮公转角速度w2=12%设置嫦娥一号绕月亮公转角速度t=0;%初始时刻为0pausetime=.002;%设置暂停时间sita1=0;sita2=0;%设置开始它们都在水平线上set(gcf,'doublebuffer','on') %消除抖动plot(-20,18,'color','r','marker','.','markersize',40);text(-17,18,'地球');%对地球进行标识p1=plot(-20,16,'color','b','marker','.','markersize',20);text(-17,16,'月亮');%对月亮进行标识p1=plot(-20,14,'color','w','marker','.','markersize',13);text(-17,14,'嫦娥一号');%对嫦娥一号进行标识plot(0,0,'color','r','marker','.','markersize',60);%画地球plot(r1*cos(s1),r1*sin(s1));%画月亮公转轨道set(gca,'xlim',[-20 20],'ylim',[-20 20]);p1=plot(r1*cos(sita1),r1*sin(sita1),'color','b','marker','.','markers ize',30);%画月亮初始位置l1=plot(r1*cos(sita1)+r2*cos(s1),r1*sin(sita1)+r2*sin(s1));%画嫦娥一号绕月亮公转轨道p2x=r1*cos(sita1)+r2*cos(sita2);p2y=r1*sin(sita1)+r2*sin(sita2);p2=plot(p2x,p2y,'w','marker','.','markersize',20);%画嫦娥一号的初始位置orbit=line('xdata',p2x,'ydata',p2y,'color','r');%画嫦娥一号的运动轨迹while 1set(p1,'xdata',r1*cos(sita1),'ydata',r1*sin(sita1));%设置月亮的运动过程set(l1,'xdata',r1*cos(sita1)+r2*cos(s1),'ydata',r1*sin(sita1)+r2*sin(s1));%设置嫦娥一号绕月亮的公转轨道的运动过程ptempx=r1*cos(sita1)+r2*cos(sita2);ptempy=r1*sin(sita1)+r2*sin(sita2); set(p2,'xdata',ptempx,'ydata',ptempy);%设置嫦娥一号的运动过程p2x=[p2x ptempx];p2y=[p2y ptempy];set(orbit,'xdata',p2x,'ydata',p2y);%设置嫦娥一号运动轨迹的显示过程sita1=sita1+w1*pausetime;%月亮相对地球转过的角度sita2=sita2+w2*pausetime;%嫦娥一号相对月亮转过的角度pause(pausetime); %暂停一会drawnowend3.曲柄连杆机构充气球动画演示程序hf=figure('name','打气筒吹气球');set(hf,'color','g');axis([-10,10,-4,4]);hold onaxis off; %除掉坐标xa0=-2.5;%活塞左顶点坐标xa1=-1.8;%活塞右顶点坐标xb0=-2;%连杆左顶点坐标xb1=5;%连杆右顶点坐标x3=5.6;%转轮坐标y3=0;%转轮坐标x4=xb1;%设置连杆头的初始位置横坐标y4=0;%设置连杆头的初始位置纵坐标x5=xa1;y5=0;x6=x3;%设置连轴初始横坐标y6=0;%设置连轴初始纵坐标a=0.7;b=0.7c=0.7a1=line([xa0;xa1],[0;0],'color','m','linestyle','-','linewidth',20); %设置活塞a8=line([-2.7;2.3],[0.3;0.3],'color','b','linestyle','-','linewidth', 5);%设置打气筒a9=line([-2.6;2.3],[-0.3;-0.3],'color','b','linestyle','-','linewidth ',5);%设置打气筒a10=line([-2.6;-2.6],[-0.1;-0.37],'color','b','linestyle','-','linewi dth',5);%设置打气筒a11=line([2.0;2.0],[0.1;0.37],'color','b','linewidth',5);%设置打气筒a12=line([2.0;2.0],[-0.1;-0.37],'color','b','linewidth',5);%设置打气筒a13=line([-2.6;-2.6],[0.1;0.37],'color','b','linestyle','-','linewidt h',5);%设置打气筒a14=line([-2.7;-2.7],[0.1;-0.1],'color','b','linestyle','-','linewidt h',9);%设置气筒嘴a16=line([-3.2;-3.2],[0.1;-0.1],'color','r','linestyle','-','linewidt h',25);%设置气筒嘴a2=line([xb0;xb1],[0;0],'color','m','linewidth',5);%设置连杆a5=line(x5,y5,'color','black','linestyle','.','markersize',25);%设置连杆活塞连接头a4=line(x4,y4,'color','black','linestyle','.','markersize',25);%设置连杆连接头a6=line([xb1;x3],[0;0],'color','b','linestyle','-','linewidth',7);%设置连杆连接轴a7=line(x3,0,'color','m','linestyle','.','markersize',50);%设置运动中心a3=line(x3,y3,'color' ,[0.5 0.60.3],'linestyle','.','markersize',85);%设置手轮len1=6.8;%连杆长len2=0.7;%活塞长r=1.3;%运动半径dd=0.01;d=-4;plot(d,0,'color','r','marker','.','markersize',10);pausetime=.0001s=0;ds=1;t=0;dt=0.015*pi;while t<=15.68t=t+dt;drawnow;lena1=sqrt((len1)^2-(r*sin(2*t))^2);%连杆在运动过程中横轴上的有效长度rr1=r*cos(t);%半径在运动过程中横轴上的有效长度xaa1=x3-sqrt(len1^2-(sin(2*t)*r)^2)-(r*cos(2*t));%活塞在运动过程中的右顶点坐标位置xaa0=xaa1-2;%%活塞在运动过程中的左顶点坐标位置x55=x3-cos(2*t)*r;%连杆在运动过程中横坐标位置y55=y3-sin(2*t)*r*0.32;%连杆在运动过程中纵坐标位置set(a4,'xdata',x55,'ydata',y55);%设置连杆顶点运动set(a1,'xdata',[xaa1-0.2;xaa1],'ydata',[0;0]);%设置活塞运动set(a2,'xdata',[xaa1;x55],'ydata',[0;y55]);set(a5,'xdata',xaa1);%设置活塞与连杆连接头的运动set(a6,'xdata',[x55;x3],'ydata',[y55;0]);if (sign(y55-y3)>0)s=s+ds;d=d-dd;if s>200s=0;d=0;endset(gcf,'doublebuffer','on');%消除震动plot(d,0,'color','r','marker','.','markersize',s);%画气球set(gcf,'doublebuffer','on') %消除抖动pause(pausetime); %暂停一会drawnow;endendx=[-8.0 -7.0 -6.0 -5.5 -4.5 -3.5 -3.7 -3 -4 -5.5 -6.5 -7.5 -9.1 -7.3 ]; y=[-1.8 -1.2 -2.8 -1.5 -2.8 -0.5 -0.7 0.7 0.2 1.5 0.8 1.2 1.9 0.5]; fill(x,y,'r');text(-6,0,'pa','fontsize',22);4.打桩机动画演示程序clear; %清除工作区clc; %清除命令区figure('name','武工院打桩机); %设置标题axis ([0 ,10,0,10]); %建立坐标系hold on;axis off; %除掉坐标text(3,9.8,'武工院打桩机','fontsize',20,'color','r');%画打桩机支架c1=line([1;5.5],[ 0.15;0.15],'color','k','linewidth',8);c2=line([1.7;5.2],[ 1.5;0.15],'color','k','linewidth',4);c3=line([1.8;5.2],[ 0.1;8],'color','k','linewidth',4);c4=line([1.7;5.2],[ 1.5;8],'color','k','linewidth',2);c5=line([1.2;1.8],[ 0.5;0.5],'color','k','linewidth',13);fill([4.9,5.3,5.3,4.9],[8.0,8.0,0.3,0.3],[1,0.1,0.5]);fill([4.4,4.6,5.6,5.8,5.6,4.6],[8.4,8.3,8.3,8.4,8.0,8.0],[1,0.1,0.5]) %画打桩机运动部分和水泥桩b1=line([8;8],[1;6],'color','b','linewidth',6);b2=line([5.5,5.5],[7.9,8],'color','k','linewidth',1);b3=line([5.3,5.55],[7.9,7.9],'color','k','linewidth',3);b4=line([5.5,5.5],[7.9,6.3],'color','k','linewidth',3);b5=line([5.5,5.5],[7,6.6],'color','k','linewidth',10);b6=line([5.5,5.5],[6.6,6.3],'color','k','linewidth',12);pausetime=1.6; %设置暂停时间pause(pausetime);%吊装水泥桩s=0;ds=0.01;pausetime1=.002;while s<2.5 %水泥桩向左移动s=s+ds;set(b1,'xdata',[8-s;8-s],'ydata',[1;6]);pause(pausetime1);endpausetime2=1;pause(pausetime2);s=0;ds=0.01;while s<1 %水泥桩向下移动s=s+ds;set(b1,'xdata',[5.5;5.5],'ydata',[1-s;6-s]);pause(pausetime1);endpause(pausetime2);s=0;ds=0.01;while s<1.3 %打桩机运动部分下移与水泥桩顶端相接触s=s+ds;set(b2,'xdata',[5.5;5.5],'ydata',[7.9-s,8]);set(b3,'xdata',[5.3;5.55],'ydata',[7.9-s,7.9-s]); set(b4,'xdata',[5.5;5.5],'ydata',[7.9-s,6.3-s]);set(b5,'xdata',[5.5;5.5],'ydata',[7-s,6.6-s]);set(b6,'xdata',[5.5;5.5],'ydata',[6.6-s,6.3-s]); pause(pausetime1);end%打桩机开始打桩pausetime3=.4;pause(pausetime2);s=0;ds=0.06; %设定打桩进度while s<4a=0;da=0.01;pausetime4=.0002;while a<.4 %两个子循环设置打桩机铁锺上下反复运动a=a+da;set(b5,'xdata',[5.5;5.5],'ydata',[5.7-s+a,5.3-s+a]); pause(pausetime4);enda=0;while a<.4a=a+da;set(b5,'xdata',[5.5;5.5],'ydata',[6.1-s-a,5.7-s-a]); pause(pausetime4);ends=s+ds;%打桩机运动部分与水泥桩同时向下移动set(b1,'xdata',[5.5;5.5],'ydata',[0;5-s]);set(b2,'xdata',[5.5;5.5],'ydata',[6.6-s,8]);set(b3,'xdata',[5.3;5.55],'ydata',[6.6-s,6.6-s]);set(b4,'xdata',[5.5;5.5],'ydata',[6.6-s,5-s]);set(b5,'xdata',[5.5;5.5],'ydata',[5.7-s,5.3-s]);set(b6,'xdata',[5.5;5.5],'ydata',[5.3-s,5-s]);pause(pausetime3);end5.电路演示程序figure('name','基本电路的模拟');axis([-3,12,0,10]);%建立坐标系hold on %保持当前图形的所有特性axis('off'); %关闭所有轴标注和控制%下面是画电池的过程fill([-1.5,-1.5,1.5,1.5],[1,5,5,1],[0.5,1,1]);%确定坐标轴范围并填充fill([-0.5,-0.5,0.5,0.5],[5,5.5,5.5,5],[0,0,0]); %确定坐标轴范围并填充text(-0.5,1.5,'负极');%在坐标上标注说明文字text(-0.5,3,'电池'); %在坐标上标注说明文字text(-0.5,4.5,'正极'); %在坐标上标注说明文字%下面是画导电线路的过程plot([0;0],[5.5;6.7],'color','r','linestyle','-','linewidth',4);%绘制二维图形线竖实心红色plot([0;4],[6.7;6.7],'color','r','linestyle','-','linewidth',4); %绘制二维图形线实心红色为导线a=line([4;5],[6.7;7.7],'color','b','linestyle','-','linewidth',4,'era semode','xor');%画开关蓝色plot([5.2;9.2],[6.7;6.7],'color','r','linestyle','-','linewidth',4);%绘制图导线为红色plot([9.2;9.2],[6.7;3.7],'color','r','linestyle','-','linewidth',4);% 绘制图导线竖线为红线plot([9.2;9.7],[3.7;3.7],'color','r','linestyle','-','linewidth',4); % 绘制图导线横线为红色plot([0;0],[1;0],'color','r','linestyle','-','linewidth',4);%如上画红色竖线plot([0;10],[0;0],'color','r','linestyle','-','linewidth',4);%如上画横线plot([10;10],[0;3],'color','r','linestyle','-','linewidth',4);%画竖线%下面是画灯泡的过程fill([9.8,10.2,9.7,10.3],[3,3,3.3,3.3],[0 0 0]);%确定填充范围plot([9.7,9.7],[3.3,4.3],'color','b','linestyle','-','linewidth',0.5) ;%绘制灯泡外形线为蓝色plot([10.3,10.3],[3.3,4.45],'color','b','linestyle','-','linewidth',0 .5); %绘制灯泡外形线为蓝色%以下为绘圆x=9.7:pi/50:10.3;%绘圆plot(x,4.3+0.1*sin(40*pi*(x-9.7)),'color','b','linestyle','-','linewi dth',0.5); %绘圆t=0:pi/60:2*pi; %绘圆plot(10+0.7*cos(t),4.3+0.6*sin(t),'color','b'); %绘圆%下面是箭头及注释的显示text(4.5,10,'电流运动方向'); %在坐标上标注说明文字line([4.5;6.6],[9.4;9.4],'color','r','linestyle','-','linewidth',4,'e rasemode','xor');%绘制箭头横线line(6.7,9.4,'color','b','linestyle','>','erasemode','xor','markersiz e',10);% %绘制箭头三角形pause(1);%下面是开关闭合的过程t=0;y=7.7;while y>6.7 %电路总循环控制开关动作条件x=4+sqrt(2)*cos(pi/4*(1-t));y=6.7+sqrt(2)*sin(pi/4*(1-t));set(a,'xdata',[4;x],'ydata',[6.7;y]);drawnow;t=t+0.1;end%下面是开关闭合后模拟大致电流流向的过程pause(1);light=line(10,4.3,'color','y','marker','.','markersize',40,'erasemode ','xor');%画灯丝发出的光:黄色%画电流的各部分h=line([1;1],[5.2;5.6],'color','r','linestyle','-','linewidth',4,'era semode','xor');g=line(1,5.7,'color','b','linestyle','^','erasemode','xor','markersiz e',10);%给循环初值t=0;m2=5.7;n=5.7;while n<6.3;%确定电流竖向循环范围m=1;n=0.05*t+5.7;set(h,'xdata',[m;m],'ydata',[n-0.5;n-0.1]);set(g,'xdata',m,'ydata',n);t=t+0.01;drawnow;endt=0;while t<2;%在转角处的停顿时间m=1.2-0.2*cos((pi/4)*t);n=6.3+0.2*sin((pi/4)*t);set(h,'xdata',[m-0.5;m-0.1],'ydata',[n;n]);set(g,'xdata',m,'ydata',n);t=t+0.05;drawnow;endt=0;while t<0.5 %在转角后的停顿时间t=t+0.5;g=line(1.2,6.5,'color','b','linestyle','^','markersize',10,'erasemode ','xor');%绘制第二个箭头g=line(1.2,6.5,'color','b','linestyle','>','markersize',10,'erasemode ','xor'); %绘制第二个箭头set(g,'xdata',1.2,'ydata',6.5);drawnow;endpause(0.5);t=0;while m<8 % 确定第二个箭头的循环范围m=1.1+0.05*t;n=6.5;set(g,'xdata',m+0.1,'ydata',6.5);set(h,'xdata',[m-0.4;m],'ydata',[6.5;6.5]);t=t+0.05;drawnow;endt=0;while t<2 %%在转角后的停顿时间m=8.1+0.2*cos(pi/2-pi/4*t);n=6.3+0.2*sin(pi/2-pi/4*t);set(g,'xdata',m,'ydata',n);set(h,'xdata',[m;m],'ydata',[n+0.1;n+0.5]);t=t+0.05;drawnow;endt=0;while t<0.5 %在转角后的停顿时间t=t+0.5;%绘制第三个箭头g=line(8.3,6.3,'color','b','linestyle','>','markersize',10,'erasemode ','xor');g=line(8.3,6.3,'color','b','linestyle','v','markersize',10,'erasemode ','xor');set(g,'xdata',8.3,'ydata',6.3);drawnow;endpause(0.5);t=0;while n>1 %确定箭头的运动范围m=8.3;n=6.3-0.05*t;set(g,'xdata',m,'ydata',n);set(h,'xdata',[m;m],'ydata',[n+0.1;n+0.5]);t=t+0.04;drawnow;endt=0;while t<2%箭头的起始时间m=8.1+0.2*cos(pi/4*t);n=1-0.2*sin(pi/4*t);set(g,'xdata',m,'ydata',n);set(h,'xdata',[m+0.1;m+0.5],'ydata',[n;n]);t=t+0.05;drawnow;endt=0;while t<0.5t=t+0.5;%绘制第四个箭头g=line(8.1,0.8,'color','b','linestyle','v','markersize',10,'erasemode ','xor');g=line(8.1,0.8,'color','b','linestyle','<','markersize',10,'erasemode ','xor');set(g,'xdata',8.1,'ydata',0.8);drawnow;endpause(0.5);t=0;while m>1.2 %箭头的运动范围m=8.1-0.05*t;n=0.8;set(g,'xdata',m,'ydata',n);set(h,'xdata',[m+0.1;m+0.5],'ydata',[n;n]);t=t+0.04;drawnow;endt=0;while t<2 %停顿时间m=1.2-0.2*sin(pi/4*t);n=1+0.2*cos(pi/4*t);set(g,'xdata',m,'ydata',n);set(h,'xdata',[m;m+0.5],'ydata',[n-0.1;n-0.5]);t=t+0.05;drawnow;endt=0;while t<0.5 %画第五个箭头t=t+0.5;g=line(1,1,'color','b','linestyle','<','markersize',10,'erasemode','x or');g=line(1,1,'color','b','linestyle','^','markersize',10,'erasemode','x or');set(g,'xdata',1,'ydata',1);drawnow;endt=0;while n<6.3 %循环范围m=1;n=1+0.05*t;set(g,'xdata',m,'ydata',n);set(h,'xdata',[m;m],'ydata',[n-0.5;n-0.1]);t=t+0.04;drawnow;end%下面是开关断开后的情况t=0;y=6.7;while y<7.7 %开关的断开x=4+sqrt(2)*cos(pi/4*t);y=6.7+sqrt(2)*sin(pi/4*t);set(a,'xdata',[4;x],'ydata',[6.7;y]);drawnow;t=t+0.1;endpause(0.5);%开关延时作用nolight=line(10,4.3,'color','y','marker','.','markersize',40,'erasemode','xor');end6.电梯动画演示程序figure('name','自控电梯');axis([-2.0,15.0,-2.0,15.0])hold onfill([-2,15,15,-2],[-2,-2,15,15],[0.5,0.2,0.3]);x1=[8 8 10 10];y1=[6 0 0 6];x2=[10 10 12 12];text(-1,8,'上升','fontsize',10,'color','c');text(4,8,'下降','fontsize',10,'color','c');text(6,12.5,'控制电机','fontsize',10,'color','c');text(12.5,3,'电梯','fontsize',10,'color','c');text(5,4.5,'关闭','fontsize',10,'color','c');text(5,1.5,'打开','fontsize',10,'color','c');l1=line([2;5.5],[11;11],'color','c','linestyle','-','linewidth',2); l2=line([2;2],[9;11],'color','c','linestyle','-','linewidth',2);l3=line([1;1],[7;9],'color','c','linestyle','-','linewidth',2);l4=line([3;3],[7;9],'color','c','linestyle','-','linewidth',2);l5=line([1;3],[9;9],'color','c','linestyle','-','linewidth',2);%绘制导线k1=line([2;1],[6;7],'color','r','linestyle','-','linewidth',2);%单刀双掷开关k2=line([4;5],[3;4],'color','b','linestyle','-','linewidth',2);%单刀双掷开关g1=line([7.7;8],[3;3],'color','b','linestyle','-','linewidth',2);g2=line([7.7;7.7],[3;4],'color','b','linestyle','-','linewidth',2); g3=line([7.4;7.7],[2;4],'color','b','linestyle','-','linewidth',2); g4=line([7.4;7.4],[2;4],'color','b','linestyle','-','linewidth',2); g5=line([7.1;7.4],[2;4],'color','b','linestyle','-','linewidth',2); g6=line([7.1;7.1],[2;4],'color','b','linestyle','-','linewidth',2); g7=line([6.8;7.1],[2;4],'color','b','linestyle','-','linewidth',2); g8=line([6.8;6.8],[2;4],'color','b','linestyle','-','linewidth',2); g9=line([6.5;6.8],[2;4],'color','b','linestyle','-','linewidth',2); g10=line([6.5;6.5],[2;3],'color','b','linestyle','-','linewidth',2); g11=line([6.0;6.5],[3;3],'color','b','linestyle','-','linewidth',2); g12=line([6;6],[2;4],'color','c','linestyle','-','linewidth',2);g13=line([5;6],[2;2],'color','c','linestyle','-','linewidth',2);g14=line([5;6],[4;4],'color','c','linestyle','-','linewidth',2);%绘制电梯门伸缩控制开关door1=patch(x1,y1,[0 1 1]);door2=patch(x2,y1,[0 1 1]);%画电梯的两面门t=0:pi/100:2*pi;fill(6+0.5*sin(t),11+cos(t),[0.7,0.85,0.9]);%电机左端fill(8.5+0.5*sin(t),11+cos(t),[0.7,0.85,0.9]);%电机右端e0=line([9;10],[11;11],'color','r','linewidth',2);%e1=line([10;10],[6;11],'color','b','linewidth',2);%连接电机中轴和电梯的线%画电机的表面(用八根不同颜色的线代替,每根之间相差pi/4)%为简便起见,初始条件下可将八根线分成两组放在电机的顶端和底端s1=line([6;8.5],[12;12],'color','c','linestyle','-','linewidth',2); s2=line([6;8.5],[10;10],'color','m','linestyle','-','linewidth',2); s3=line([6;8.5],[12;12],'color','b','linestyle','-','linewidth',2); s4=line([6;8.5],[10;10],'color','w','linestyle','-','linewidth',2); s5=line([6;8.5],[12;12],'color','k','linestyle','-','linewidth',2); s6=line([6;8.5],[10;10],'color','g','linestyle','-','linewidth',2); s7=line([6;8.5],[12;12],'color','r','linestyle','-','linewidth',2); s8=line([6;8.5],[10;10],'color','b','linestyle','-','linewidth',2); a=0; %设定电机运转的初始角度da=0.05;%设定电机正转的条件s=0; %设定门运动的初始条件ds=0.05;%设定门运动的周期while s<5 %条件表达式 (当0<s<5时,电机正转,门上升)a=a+da;xa1=6+abs(0.5*sin(a)); %当线运动到电机背面时会覆盖电机左端,用abs可解决这一问题增强逼真感)xa2=8.5+0.5*sin(a);ya1=11+cos(a);ya2=11+cos(a);%设定s1的两端点坐标(s1是对应0的线)xb1=6+0.5*abs(sin(a+pi));xb2=8.5+0.5*sin(a+pi);yb1=11+cos(a+pi);yb2=11+cos(a+pi); %设定s2的两端点坐标(s2是对应pi的线)xc1=6+abs(0.5*sin(a+pi/2));xc2=8.5+0.5*sin(a+pi/2);yc1=11+cos(a+pi/2);yc2=11+cos(a+pi/2);%设定s3的两端点坐标(s3是对应pi/2的线)xd1=6+0.5*abs(sin(a-pi/2));xd2=8.5+0.5*sin(a-pi/2);yd1=11+cos(a-pi/2);yd2=11+cos(a-pi/2);%设定s4的两端点坐标(s4是对应-pi/2的线)xe1=6+abs(0.5*sin(a+pi/4));xe2=8.5+0.5*sin(a+pi/4);ye1=11+cos(a+pi/4);ye2=11+cos(a+pi/4);%设定s5的两端点坐标(s5是对应pi/4的线)xf1=6+0.5*abs(sin(a+pi*3/4));xf2=8.5+0.5*sin(a+pi*3/4);yf1=11+cos(a+pi*3/4);yf2=11+cos(a+pi*3/4);%设定s6的两端点坐标(s6是对应pi*3/4的线)xg1=6+abs(0.5*sin(a-pi*3/4));xg2=8.5+0.5*sin(a-3*pi/4);yg1=11+cos(a-3*pi/4);yg2=11+cos(a-3*pi/4);%设定s7的两端点坐标(s7是对应-3*pi/4的线)xh1=6+0.5*abs(sin(a-pi/4));xh2=8.5+0.5*sin(a-pi/4);yh1=11+cos(a-pi/4);yh2=11+cos(a-pi/4); %设定s8的两端点坐标(s8是对应-pi/4的线)set(s1,'xdata',[xa1;xa2],'ydata',[ya1;ya2]);set(s2,'xdata',[xb1;xb2],'ydata',[yb1;yb2]);set(s3,'xdata',[xc1;xc2],'ydata',[yc1;yc2]);set(s4,'xdata',[xd1;xd2],'ydata',[yd1;yd2]);set(s5,'xdata',[xe1;xe2],'ydata',[ye1;ye2]);set(s6,'xdata',[xf1;xf2],'ydata',[yf1;yf2]);set(s7,'xdata',[xg1;xg2],'ydata',[yg1;yg2]);set(s8,'xdata',[xh1;xh2],'ydata',[yh1;yh2]); %绘制电机表面各线条的运动s=s+ds;set(door1,'xdata',x1,'ydata',[6+0.5*s 0+0.5*s 0+0.5*s 6+0.5*s]);set(door2,'xdata',x2,'ydata',[6+0.5*s 0+0.5*s 0+0.5*s 6+0.5*s]); %绘制门的向上运动set(e1,'xdata',[10;10],'ydata',[6+0.5*s;11]); %绘制门顶的绳索的向上运动set(gcf,'doublebuffer','on');%消除振动drawnow;endb=0;%设定电机反转的条件db=0.05;while s<10 %条件表达式 (当5<s<10时,电机反转,门下降)b=b-db;xa1=6+abs(0.5*sin(a+b));xa2=8.5+0.5*sin(a+b);ya1=11+cos(a+b);ya2=11+cos(a+b);%设定s1的两端点坐标(s1是对应0的线)xb1=6+0.5*abs(sin(a+pi+b));xb2=8.5+0.5*sin(a+pi+b);yb1=11+cos(a+pi+b);yb2=11+cos(a+pi+b);%设定s2的两端点坐标(s2是对应pi的线)xc1=6+abs(0.5*sin(a+pi/2+b));xc2=8.5+0.5*sin(a+pi/2+b);yc1=11+cos(a+pi/2+b);yc2=11+cos(a+pi/2+b);%设定s3的两端点坐标(s3是对应pi/2的线)xd1=6+2*abs(sin(a-pi/2+b));xd2=8.5+0.5*sin(a-pi/2+b);yd1=11+cos(a-pi/2+b);yd2=11+cos(a-pi/2+b);%设定s4的两端点坐标(s4是对应-pi/2的线)xe1=6+abs(0.5*sin(a+pi/4+b));xe2=8.5+0.5*sin(a+pi/4+b);ye1=11+cos(a+pi/4+b);ye2=11+cos(a+pi/4+b);%设定s5的两端点坐标(s5是对应pi/4的线)xf1=6+0.5*abs(sin(a+pi*3/4+b));xf2=8.5+0.5*sin(a+pi*3/4+b);yf1=11+cos(a+pi*3/4+b);yf2=11+cos(a+pi*3/4+b);%设定s6的两端点坐标(s6是对应pi*3/4的线)xg1=6+abs(0.5*sin(a-pi*3/4+b));xg2=8.5+0.5*sin(a-3*pi/4+b);yg1=11+cos(a-3*pi/4+b);yg2=11+cos(a-3*pi/4+b);%设定s7的两端点坐标(s7是对应-3*pi/4的线)xh1=6+0.5*abs(sin(a-pi/4+b));xh2=8.5+0.5*sin(a-pi/4+b);yh1=11+cos(a-pi/4+b);yh2=11+cos(a-pi/4+b);%设定s8的两端点坐标(s8是对应-pi/4的线)%绘制电机表面各线条的运动set(s1,'xdata',[xa1;xa2],'ydata',[ya1;ya2]);set(s2,'xdata',[xb1;xb2],'ydata',[yb1;yb2]);set(s3,'xdata',[xc1;xc2],'ydata',[yc1;yc2]);set(s4,'xdata',[xd1;xd2],'ydata',[yd1;yd2]);set(s5,'xdata',[xe1;xe2],'ydata',[ye1;ye2]);set(s6,'xdata',[xf1;xf2],'ydata',[yf1;yf2]);set(s7,'xdata',[xg1;xg2],'ydata',[yg1;yg2]);set(s8,'xdata',[xh1;xh2],'ydata',[yh1;yh2]);s=s+ds;set(k1,'xdata',[2;3],'ydata',[6;7]); %绘制闸刀的换向运动set(door1,'xdata',x1,'ydata',[11-0.5*s 5-0.5*s 5-0.5*s 11-0.5*s]); set(door2,'xdata',x2,'ydata',[11-0.5*s 5-0.5*s 5-0.5*s 11-0.5*s]);%绘制门的向下运动set(e1,'xdata',[10;10],'ydata',[11-0.5*s;11]); %绘制门顶绳索的向下运动set(gcf,'doublebuffer','on');%消除振动drawnow;endfor i=1:400a1=10-0.005*i;a2=10+0.005*i;x1=[8 8 a1 a1];x2=[a2 a2 12 12];set(door1,'xdata',x1);set(door2,'xdata',x2);set(k2,'xdata',[4;5],'ydata',[3;2]); %绘制闸刀的换向运动set(gcf,'doublebuffer','on');%消除振动drawnow;end %绘制电梯门的打开运动for i=1:400a3=8+0.005*i;a4=12-0.005*i;x1=[8 8 a3 a3];x2=[a4 a4 12 12];set(door1,'xdata',x1);set(door2,'xdata',x2);set(k2,'xdata',[4;5],'ydata',[3;4]); %绘制闸刀的换向运动set(gcf,'doublebuffer','on');%消除振动drawnow;end %绘制电梯门的关闭运动7.防汛检测系统动画演示程序for j=0:10axis([-1 1 -1 1]);%设置x,y的坐标范围axis('off');%覆盖坐标刻度x1=[0 0 0.8 0.8];y1=[-0.6 -0.8 -0.8 -0.6];%对水槽中的水进行初设置line([0;0],[0.2;-0.8],'color','k','linewidth',3);%水槽左壁的颜色和宽度line([0;0.8],[-0.8;-0.8],'color','k','linewidth',3);%水槽底部的颜色和宽度line([0.8;0.8],[-0.7;-0.8],'color','k','linewidth',3);%水槽右边出水口的下面的颜色和宽度line([0.8;0.8],[0.2;-0.6],'color','k','linewidth',3);%水槽右边出水口的上面的颜色和宽度line([0.8;0.85],[-0.7;-0.7],'color','k','linewidth',3);%出水口的下壁的颜色和宽度line([0.8;0.85],[-0.6;-0.6],'color','k','linewidth',3);%出水口的上壁的颜色和宽度line(-0.35,0,'Color','r','linestyle','.', 'markersize',20);%给水线处。

合集下载

matlab动画-流星雨

matlab动画-流星雨

clc;%清屏ysy=figure('name','一起去看流星雨');%标题axis([0 80 0 80]);hold on;axis off;set(ysy,'color','k');%设置背景为黑色%%%%%以下是画月亮t6=-pi/3:0.01:0;x6=60+10*cos(t6);y6=60+10*sin(t6);plot(x6,y6,'color',[1 1 0]);%画第一条弧线t7=-pi*2/3:0.01:pi/3;x7=67.5+5*cos(t7);y7=60-sqrt(3)*5/2+5*sin(t7);plot(x7,y7,'color',[1 1 0]);%画第二条弧线fill([x6,x7],[y6,y7],'y');%将月亮填充为黄色%%%%%以下为画五角星line([30,33,34,35,38,36,38,34,30,32,30],[60,60,62.5,60,60,59,56.5,58,56.5,59,60],...'color','b','MarkerEdgeColor','b','linewidth',2);%画五角星fill([30,33,34,35,38,36,38,34,30,32,30],...[60,60,62.5,60,60,59,56.5,58,56.5,59,60],'b');%填充五角星%%%%%以下是画人%%%%%%%男t8=0:0.01:2*pi;x8=40+2*cos(t8);y8=24+4*sin(t8);plot(x8,y8,'b','linewidth',2);%%%%%画头line([40,40.8],[26,26],'color','y');%左眼line([41.4,41.8],[26,26],'color','y');%右眼line(41.3,24.1,'color','y','linestyle','.','erasemode','xor','markersize',5);%鼻line([41,41.8],[22,22],'color','y');%嘴line([39.8,39.8,38.3,38.3,41.7,41.7,40.8,40.8,39.8],[20,19,18.7,0,0.5,18.7,19,20,20],...'color','b','linewidth',2);%躯干line([41.7,41.7],[13.4,16.2],'color','k','linewidth',2)%清除被手臂挡住的线line([40.7,40.7],[0.35,8.7],'color','b','linewidth',2);%用一条线将矩形分开,就是腿line([40,41,44,45,41,40],[18.2,16,18.8,17.8,13.1,15.3],'color','b','linewidth',2);%%胳膊%%%%%%%%女t9=0:0.01:2*pi;x9=50+2*cos(t8);y9=22+4*sin(t8);plot(x9,y9,'color','m','linewidth',2);%%%%%画头line([50,49.2],[24,24],'color','y');%左眼line([48.2,48.6],[24,24],'color','y');%右眼line(48.7,22.1,'color','y','linestyle','.','erasemode','xor','markersize',5);%鼻line([48.2,49.2],[20,20],'color','y');%嘴line([49.8,49.8,48.3,48.3,46.3,48.3,48.3,51.7,51.7,53.7,51.7,51.7,50.8,50.8],...[18,17,16.7,8,4,4,0.5,0,4,4,8,16.7,17,18],'color','m','linewidth',2);%躯干line([48.3,48.3],[14.1,16.1],'color','k','linewidth',2)%清除挡住的线line([49.3,49.3],[0.35,4],'color','m','linewidth',2);%腿line([50,47.8,45.8,45,47.8,50],[16.4,16,18.4,17.8,14,14.4],'color','m','linewidth',2);%%胳膊b1=text(20,70,'哇塞!流星雨来了!','fontsize',16,'color','r');%输入字pause(1);%停顿一秒delete(b1);%删除b1b2=text(22,68,'听说...','fontsize',16,'color','r');pause(1);delete(b2);b3=text(24,66,'在流星雨下许的愿,一定会实现!','fontsize',16,'color','r');pause(2);delete(b3);b4=text(54,28,'愿我们相伴一生!','fontsize',16,'color','r');b5=text(10,30,'愿执子之手,','fontsize',16,'color','r');b6=text(12,22,'与子偕老!','fontsize',16,'color','r');a1=line(5,70,'color','w','linestyle','.','erasemode','xor','markersize',10);%流星1a2=line(5,65,'color','w','linestyle','.','erasemode','xor','markersize',10);%流星2a3=line(5,60,'color','w','linestyle','.','erasemode','xor','markersize',10);%流星3a4=line(5,55,'color','w','linestyle','.','erasemode','xor','markersize',10);%流星4a5=line(5,50,'color','w','linestyle','.','erasemode','xor','markersize',10);%流星5x1=5;%流星1的初始横坐标y1=70;;%流星1的初始纵坐标x2=5;y2=65;x3=5;y3=60;x4=5;y4=55;x5=5;y5=50;dt1=0.4375;%递增量dt2=0.1;%递增量t1=x1;t2=x2;t3=x3;t4=x4;t5=x5;r1=1;%画心形中的小圆半径r2=4;%画心形中的大圆半径c1=39;%小圆1的初始横坐标d1=33;%小圆1的初始纵坐标c2=41;%小圆2的初始横坐标d2=33;%小圆2的初始纵坐标c3=38;%大圆1的初始横坐标d3=33;%大圆1的初始纵坐标c4=42;%大圆2的初始横坐标d4=33;%大圆2的初始纵坐标f1=c1+10;%第二个心形中小圆1的横坐标f2=c2+10;%第二个心形中小圆2的横坐标f3=c3+10;%第二个心形中大圆1的横坐标f4=c4+10;%第二个心形中大圆2的横坐标%第二个心形中的纵坐标跟第一个共用t=0.09;%递增量t10=0:0.01:pi;%设置弧度t11=-pi*11/30:0.01:0;t12=-pi:0.01:-pi*19/30;%流星雨开始及心形移动while t5<80 %循环开始t1=t1+dt1;t2=t2+dt1;t3=t3+dt1;t4=t4+dt1;t5=t5+dt1;x1=t1;x2=t2;x3=t3;x4=t4;x5=t5;y1=y1+dt2;y2=y2+dt2;y3=y3+dt2;y4=y4+dt2;y5=y5+dt2;set(a1,'xdata',[x1,x1],'ydata',[y1,y1]);%流星1的移动set(a2,'xdata',[x2,x2],'ydata',[y2,y2]);%流星2的移动set(a3,'xdata',[x3,x3],'ydata',[y3,y3]);%流星3的移动set(a4,'xdata',[x4,x4],'ydata',[y4,y4]);%流星4的移动set(a5,'xdata',[x5,x5],'ydata',[y5,y5]);%流星5的移动c1=c1+t/3.2;c2=c2+t/3.2;c3=c3+t/3.2;c4=c4+t/3.2;f1=f1-t/3.2;f2=f2-t/3.2;f3=f3-t/3.2;f4=f4-t/3.2;d1=d1+t;d2=d2+t;d3=d3+t;d4=d4+t;c5=c1+r1*cos(t10);d5=d1+r1*sin(t10);f5=f1+r1*cos(t10);c6=c2+r1*cos(t10);d6=d2+r1*sin(t10);f6=f2+r1*cos(t10);c7=c3+r2*cos(t11);d7=d3+r2*sin(t11);f7=f3+r2*cos(t11);c8=c4+r2*cos(t12);d8=d4+r2*sin(t12);f8=f4+r2*cos(t12);plot(c5,d5,'r');%小圆弧1plot(c6,d6,'r');%小圆弧2plot(c7,d7,'r');%大圆弧1plot(c8,d8,'r');%大圆弧2plot(f5,d5,'r');%第二个心形小圆弧1plot(f6,d6,'r');%小圆弧2plot(f7,d7,'r');%大圆弧1plot(f8,d8,'r');%大圆弧2pause(0.01);%停顿0.01秒drawnow;fill([c1-1.1,c2+1.1,c2+1.1,c1-1.1],[d1-4,d1-4,d2+1,d2+1],'k');%除去显示出来的心形1,以便移动fill([f1-1.1,f2+1.1,f2+1.1,f1-1.1],[d1-4,d1-4,d2+1,d2+1],'k');%除去显示出来的心形2,以便移动end;%循环结束plot(c5,d5,'r');plot(c6,d6,'r');plot(c7,d7,'r');plot(c8,d8,'r');%显示最后的心形。

matlab动态绘图Animation

matlab动态绘图Animation

matlab动态绘图Animation 主要是查看matlab帮助的Animation动画页⾯1、函数播放录制的电影桢将轴或图形捕获为影⽚帧返回与电影帧相关联的图像数据将图像转换为电影帧创建动画线⼆维彗星图三维彗星图更新数据并处理回调指定数据源后刷新图中的数据线动画2、官⽅例⼦(1)画出sin(x)图像,并且追踪⼀个点%%%创建sinxx = linspace(0,10,1000);y = sin(x);plot(x,y)hold onp = plot(x(1),y(1),'o','MarkerFaceColor','red'); %标记第⼀个点,设置为红⾊hold offaxis manual %将轴限制模式设置为⼿动,以避免在整个动画循环中重新计算限制。

%%%使⽤drawnow动态画图for k = 2:length(x)p.XData = x(k);p.YData = y(k);drawnowend(2)画出sin(x)图像,追踪⼀个点,显⽰坐标变换%%%画出sin(x)图像x = linspace(-6,6,1000);y = sin(x);plot(x,y)axis manual%%%设置图⽚格式ax = gca;h = hgtransform('Parent',ax);hold onplot(x(1),y(1),'o','Parent',h);hold offt = text(x(1),y(1),num2str(y(1)),'Parent',h,...'VerticalAlignment','top','FontSize',14);%%%使⽤drawnow动态画图for k = 2:length(x)m = makehgtform('translate',x(k)-x(1),y(k)-y(1),0);h.Matrix = m;t.String = num2str(y(k));drawnowend(3)这个例⼦展⽰了如何通过更新三⾓形的数据属性来让⼀个三⾓形在圆的内部循环。

matlab综合题:采用动画的形式书写本人的名字,具体字体自定

matlab综合题:采用动画的形式书写本人的名字,具体字体自定

matlab综合题:采用动画的形式书写本人的名字,具体字体自定为了在MATLAB中创建一个动画来书写您的名字,我们可以使用plot函数来绘制每个字符,并使用pause函数来逐帧显示动画。

以下是一个简单的示例,演示如何创建这样的动画。

假设您的名字是"张三",我们将使用宋体字体。

matlab% 初始化参数fontsize = 20; % 字体大小x = 0:0.1:10; % x轴范围y = sin(x); % 示例数据% 获取您的名字name = input('请输入您的名字: ','s');% 创建动画帧for i = 1:length(name)% 绘制当前字符plot(x, y, 'fontsize', fontsize);text(5, 0, name(i), 'fontsize', fontsize);% 暂停并显示当前帧pause(0.1); % 控制动画速度clear x y name; % 清空变量以准备下一个字符end在上述代码中,我们首先定义了一些参数,如字体大小和x轴范围。

然后,我们使用input函数从用户那里获取名字。

接下来,我们使用一个for循环来逐个绘制名字中的每个字符。

对于每个字符,我们使用plot函数绘制背景,并使用text函数在适当的位置绘制字符。

然后,我们使用pause函数暂停程序,以便用户可以看到当前帧。

最后,我们使用clear函数清空变量,以便绘制下一个字符。

注意:这只是一个基本的示例,您可能需要根据自己的需求进行修改和优化。

matlab动态爱心代码

matlab动态爱心代码

matlab动态爱心代码Matlab动态爱心代码介绍Matlab是一种广泛使用的数学软件,它具有强大的计算能力和丰富的绘图功能。

本文将介绍如何使用Matlab编写一个动态爱心代码。

步骤1. 创建一个新的Matlab脚本文件。

2. 定义变量t,从0到2*pi生成100个点。

```t = linspace(0, 2*pi, 100);```3. 计算x和y坐标。

```x = 16*sin(t).^3;y = 13*cos(t)-5*cos(2*t)-2*cos(3*t)-cos(4*t); ```4. 绘制爱心图形。

```plot(x, y, 'r', 'LineWidth', 2);axis equal off;hold on;fill(-x, y, 'r');```5. 添加动画效果。

```for i = 1:50fill(-x, y, 'w');x = x + randn(size(x))/10;y = y + randn(size(y))/10;plot(x, y, 'r', 'LineWidth', 2);axis equal off;hold on;fill(-x, y, 'r');pause(0.1);end```6. 运行脚本文件,观察动态爱心效果。

完整代码```matlabt = linspace(0, 2*pi, 100);x = 16*sin(t).^3;y = 13*cos(t)-5*cos(2*t)-2*cos(3*t)-cos(4*t);plot(x, y, 'r', 'LineWidth', 2);axis equal off;hold on;fill(-x, y, 'r');for i = 1:50fill(-x, y, 'w');x = x + randn(size(x))/10;y = y + randn(size(y))/10;plot(x, y, 'r', 'LineWidth', 2);axis equal off;hold on;fill(-x, y, 'r');pause(0.1);end```总结本文介绍了如何使用Matlab编写一个动态爱心代码。

Matlab-制作动画

Matlab-制作动画
11
保存动画的示例1
• 三维动画
[x,y,z]=peaks(30); %peaks产生一个凹凸有致的曲面 surf(x,y,z) %用X和Y定义x-y坐标网格,Z定义网格上每一点的高度,来生成三维曲面 axis([-3 3 -3 3 -10 10]) m=moviein(15); for i=1:15 view(-37.5+24*(i-1),30) m(:,i)=getframe; end movie(m)
surf(sin(2*pi*j/20)*Z,Z) F(j) = getframe; end movie(F,20)
15
其它动画示例程序
clear;clc; [x,y]=meshgrid(-8:.1:8); for j=1:20 f=@(x,y)(sin(sqrt((11-j)*(x.^2+y.^2)))./sqrt((11-j)*(x.^2+y.^2)+eps)); z=f(x,y); surf(x,y,z);shading interp; M(j) = getframe; if j==1 [I,map]=rgb2ind(M(j).cdata,256); imwrite(I,map,'out.gif','DelayTime',.1) else imwrite(rgb2ind(M(j).cdata,map),map,'out.gif','WriteMode','append','DelayTime',.1) end end movie2avi(M,'out.avi')
• 一般格式:
for i=1:n figure(i) F(i)=getframe; end; movie2avi(F,'...filename.avi')

matlab画立体蝴蝶的代码

matlab画立体蝴蝶的代码

matlab画立体蝴蝶的代码Matlab是一款强大的工具,可以帮助我们完成许多数据分析、数学建模和图像处理等工作。

在Matlab中,我们可以用简单的代码绘制出许多精美的图像,例如立体蝴蝶。

下面就来介绍一下如何使用Matlab绘制立体蝴蝶的代码。

1.准备工作我们需要在Matlab中打开一个新的文件,然后输入下面这句代码创建一个3D画布:figure('units','normalized','outerposition',[0 0 1 1])这句代码将创建一个全屏的3D画布,供我们后续用来绘制立体蝴蝶。

2.绘制蝴蝶的翅膀接下来,我们需要定义蝴蝶的两个翅膀,具体代码如下:t = 0:pi/10:2*pi;x = sin(t).*(exp(cos(t))-2*cos(4*t)-sin(t/12).^5);y = cos(t).*(exp(cos(t))-2*cos(4*t)-sin(t/12).^5);z = sin(t/2).^2.*cos(5*t+pi/2);plot3(x,y,z,'r','LineWidth',2)这段代码的意思是先设置一个参数t,然后用该参数计算出x、y、z坐标分别对应蝴蝶翅膀上的点的位置。

最后使用plot3函数将这些点连接起来,形成一个立体的蝴蝶翅膀。

这时我们已经完成了蝴蝶的一只翅膀的绘制。

接下来我们需要将另一只翅膀按照一定的规律旋转、平移后绘制出来。

代码如下:hold onfor i=1:3xx = x*cos(pi*i/3) + y*sin(pi*i/3) - 5*(cos(pi*i/3)+1)/2; yy = -x*sin(pi*i/3) + y*cos(pi*i/3);zz = z+2*i;plot3(xx,yy,zz,'r','LineWidth',2)end这段代码首先使用for循环来遍历翅膀的三个位置。

Matlab动画演示教学

y s ( A: =i x n) x a: a(; m x m xx ) y i  ̄n ) mn = (; y a i Ox Ixy n1 ; odOl x ̄ HB mi D hl l po x , ; l (y t , ' h l e'oo'10O, re'.' res e,0"r ̄mo e,o3 =i ( lr[ lmak r , kri ' 。 ae d 'xr; n e , " ,"ma  ̄ 4 e ' n lnt x; -e g ) h(
维普资讯
22 1

建 电

20 0 8年第 6期
Ma a 画演 示教 学 t b动 l

【 摘
萸 .张敬华 z
(1 . 福建农林 大学 计算机 与信 息学院 福建 福州 3 0 0 . 50 2 2 福建农林大学 计算机与信息 学院 福 建 福州 30 0 2、 50 2 要】 :介 绍了 Ma a l fb的功能与特点 , 阐述 了Mal t b在 高等数 学教 学中的作 用. 以高等数 学教 学中遇到的 问题为 a 并
例 , 细讨 论 了 MA L B应 用 于 高 等数 学 的教 学模 式 。 详 TA
【 关键词】 :Ma a ; t b 高等数 学教 学; l 可视化
0 Ma a 、 n b的动 画功 能与 特 点 ; i m了 1 MAIAB是 Ma1 rs 司 开发 的新 一 代 的科 学 与 工 程 计 例 1 l , 1 t l k公 Wo 算 软件 。 已经 成 为全 球 应 用 最 广 泛 最 流行 的软 件 之 一 . a a 自 m tb l 在 课 本 中 . 于 关 产 生之 日起 就具 有 方便 的数 据 可 视 化 功能 . 以将 向量 和 距 阵 用 l —lX i Sl a r i

MATLAB动画演示效果

MATLAB动画演示效果一%曲柄滑块机构hf=figure('name','曲柄滑块机构');set(hf,'color','g');hold onaxis([-6,6,-4,4]);grid onaxis('off');xa0=-5;%活塞左顶点坐标xa1=-2.5;%活塞右顶点坐标xb0=-2.5;%连杆左顶点坐标xb1=2.2;%连杆右顶点坐标x3=3.5;%转轮坐标y3=0;%转轮坐标x4=xb1;%设置连杆头的初始位置横坐标y4=0;%设置连杆头的初始位置纵坐标x5=xa1;y5=0;x6=x3;%设置连轴初始横坐标y6=0;%设置连轴初始纵坐标a=0.7;b=0.7c=0.7a1=line([xa0;xa1],[0;0],'color','b','linestyle','-','linewidth',40); %设置活塞a3=line(x3,y3,'color',[0.5 0.6 0.3],'linestyle','.','markersize',300);%设置转轮a2=line([xb0;xb1],[0;0],'color','black','linewidth',10);%设置连杆a5=line(x5,y5,'color','black','linestyle','.','markersize',40);%设置连杆活塞连接头a4=line(x4,y4,'color','black','linestyle','.','markersize',50);%设置连杆连接头a6=line([xb1;x3],[0;0],'color','black','linestyle','-','linewidth',10);a7=line(x3,0,'color','black','linestyle','.','markersize',50);%设置运动中心a8=line([-5.1;-0.2],[0.7;0.7],'color','y','linestyle','-','linewidth',5);%设置汽缸壁a9=line([-5.1;-0.2],[-0.72;-0.72],'color','y','linestyle','-','linewidth',5);%设置汽缸壁a10=line([-5.1;-5.1],[-0.8;0.75],'color','y','linestyle','-','linewidth',5);%设置汽缸壁a11=fill([-5,-5,-5,-5],[0.61,0.61,-0.61,-0.61],[a,b,c]);%设置汽缸气体len1=4.8;%连杆长len2=2.5;%活塞长r=1.3;%运动半径dt=0.015*pi;t=0;while 1t=t+dt;if t>2*pit=0;endlena1=sqrt((len1)^2-(r*sin(t))^2);%连杆在运动过程中横轴上的有效长度rr1=r*cos(t);%半径在运动过程中横轴上的有效长度xaa1=x3-sqrt(len1^2-(sin(t)*r)^2)-(r*cos(t));%活塞在运动过程中的右顶点坐标位置xaa0=xaa1-2.5;%%活塞在运动过程中的左顶点坐标位置x55=x3-cos(t)*r;%连杆在运动过程中横坐标位置y55=y3-sin(t)*r;%连杆在运动过程中纵坐标位置set(a4,'xdata',x55,'ydata',y55);%设置连杆顶点运动set(a1,'xdata',[xaa1-2.5;xaa1],'ydata',[0;0]);%设置活塞运动set(a2,'xdata',[xaa1;x55],'ydata',[0;y55]);set(a5,'xdata',xaa1);%设置活塞与连杆连接头的运动set(a6,'xdata',[x55;x3],'ydata',[y55;0]);set(a11,'xdata',[-5,xaa0,xaa0,-5]);%设置气体的填充set(gcf,'doublebuffer','on');%消除震动drawnow;end文献出处:/blog/static/13485835420091124584320/MATLAB动画效果演示二MATLAB 2009-12-24 17:15:48 阅读135 评论0 字号:大中小订阅%理想弹簧阵子简谐运动%Clearrectangle('position',[12,8.5,2,0.3],'FaceColor',[0.5,0.3,0.4]);axis([0,15,-1,10]);%画顶板hold onplot([13,13],[7,8.5],'r','linewidth',2);%画直线y=2:.2:7;M=length(y);x=12+mod(1:M,2)*2;x(1)=13;x(end-3:end)=13;D=plot(x,y);%弹簧C=0:.1:2*pi;r=0.35;t1=r*sin(C);F1=fill(13+r*cos(C),2+t1,'r');% 球set(gca,'ytick',[0:2:9]);set(gca,'yticklabels',num2str([-1:3]'));plot([0,15],[3.3,3.3],'black');H1=plot([0,13],[3.3,3.3],'y');% 句柄[黄线]Q=plot(0,3.8,'color','r');% 运动曲线;td=[];yd=[];T=0;text(2,9,'理想中的弹簧振子简谐振动','fontsize',16);set(gcf,'doublebuffer','on');while T<12;pause(0.2);Dy=(3/2-1/2*sin(pi*T))*1/2;Y=-(y-2)*Dy+7;Yf=Y(end)+t1;td=[td,T];yd=[yd,Y(end)];set(D,'ydata',Y);set(F1,'ydata',Yf,'facecolor',rand(1,3));set(H1,'xdata',[T,13],'ydata',[Y(end),Y(end)]);set(Q,'xdata',td,'ydata',yd) ;T=T+0.1;end文献出处:/blog/static/1348583542009112451548632/%台球完全碰撞模拟示例axis([-1.6,12.6,-1.6,10.7])%确定坐标轴参数范围hold on %保持当前图形及轴系的所有特性fill([-2,13,13,-2],[-2,-2,11,11],[0,1,0]);%填充底座背景fill([-1,12,12,-1],[-1,-1,10,10],[0,0.5,0]);%填充底座背景ball1=line(0,5,'color','r','marker','.','erasemode','xor','markersize',60);%设置小球颜色,大小,线条的擦拭方式ball2=line(8,9,'color','g','marker','.','erasemode','xor','markersize',60);%设置小球颜色,大小,线条的擦拭方式ball3=line(-1,-1,'color','g','marker','.','erasemode','xor','markersize',80);%设置左下角圆的颜色,大小,线条的擦拭方式ball4=line(12,-1,'color','g','marker','.','erasemode','xor','markersize',80);%设置右下角圆的颜色,大小,线条的擦拭方式ball3=line(-1,10,'color','g','marker','.','erasemode','xor','markersize',80);%设置左上角圆的颜色,大小,线条的擦拭方式ball4=line(12,10,'color','g','marker','.','erasemode','xor','markersize',80);%设置右上角圆的颜色,大小,线条的擦拭方式title('完全非弹碰在模拟台球比赛的应用', 'color','r','fontsize',15);%图形标题pause(1)%设定暂停时间的长度t=0;dt=0.005;%设制初始数值while t<7.2%设定横轴范围t=t+dt;%设制横轴计算公式y=1/2*t+5;%设制纵轴计算公式set(ball1,'xdata',t,'ydata',y)%设制球的运动drawnow;%刷新屏幕end %结束while t<8.8%设定横轴范围t=t+dt;%设制横轴计算公式y=1/2*t+5;%设制纵轴计算公式set(ball2,'xdata',t,'ydata',y)%设制球的运动drawnow;%刷新屏幕end %结束while t<11.5%设定横轴范围t=t+dt;%设制横轴计算公式y=-1/2*t+14.3;%设制纵轴计算公式set(ball2,'xdata',t,'ydata',y)%设制球的运动drawnow;%刷新屏幕end%结束while t>-0.5%设制横轴范围t=t-dt;%设制横轴计算公式y=1/2*t+2.90;%设制纵轴计算公式set(ball2,'xdata',t,'ydata',y)%设制球的运动drawnow;%刷新屏幕end%结束while t<6%设制横轴范围t=t+dt;%设制横轴计算公式y=-1/2*t+2.40;%设制纵轴计算公式set(ball2,'xdata',t,'ydata',y)%设制球的运动drawnow;%刷新屏幕end%结束while t<11.5%设制横轴范围t=t+dt;%设制横轴计算公式y=1/2*t-3.0;%设制纵轴计算公式set(ball2,'xdata',t,'ydata',y)%设制球的运动drawnow;%刷新屏幕end%结束while t>-2%设制横轴范围t=t-dt;%设制横轴计算公式y=-t*7.65/12.9+9.57;%设制纵轴计算公式set(ball2,'xdata',t,'ydata',y)%设制球的运动drawnow;%刷新屏幕end%结束text(2,5,'好!进球了啊!!!恭喜!','fontsize',16,'color','r'); %显示字幕的颜色和大小文献出处:/blog/static/134858354200911245185379/MATLAB动画效果演示四MATLAB 2009-12-24 17:22:18 阅读50 评论2 字号:大中小订阅%能量守恒验证示例fill([6,7,7,6],[5,5,0,0],[0,0.5,0]);%右边竖条的填充hold on; %保持当前图形及轴系的所有特性fill([2,6,6,2],[3,3,0,0],[0,0.5,0]);%左边竖条的填充hold on;% 保持当前图形及轴系的所有特性t1=0:pi/60:pi;plot(4-2*sin(t1-pi/2),5-2*cos(t1-pi/2));%绘制中间的凹弧图形grid;%添加网格线axis([0,9,0,9]);%定义坐标轴的比例%axis('off');%关闭所有轴标注,标记,背景fill([1,2,2,1],[5,5,0,0],[0,0.5,0]);%中间长方形的填充hold on;% 保持当前图形及轴系的所有特性title('31608118');%定义图题x0=6;y0=5;head1=line(x0,y0,'color','r','linestyle','.','erasemode','xor','markersize',30);head2=line(x0,y0,'color','r','linestyle','.','erasemode','xor','markersize',50); %设置小球颜色,大小,线条的擦拭方式t=0;%设置小球的初始值dt=0.001;%设置运动周期t1=0;%设置大球的初始值dt1=0.001;while 1%条件表达式t=t+dt;x1=9-1*t;y1=5;x3=6;y3=5;if t>0x2=6;y2=5;%设置小球的运动轨迹endif t>2.8t=t+dt;a=sin(t-3);x1=6.1;y1=5.1;x3=4-2*sin(1.5*a);y3=5-2*cos(1.5*a);%设置大球的运动轨迹endset(head1,'xdata',x1,'ydata',y1);%设置球的运动set(head2,'xdata',x3,'ydata',y3);drawnow;end文献出处:/blog/static/1348583542009112452218157/。

MATLAB经典教程(全)PPT课件


THANKS FOR WATCHING
感谢您的观看
信号时域分析和频域分析
时域分析
研究信号随时间变化的规律,包括波形、幅度、频率、相位等。
频域分析
将信号转换为频域表示,研究信号的频谱结构和频率特性,包括幅 度谱、相位谱、功率谱等。
时域与频域关系
时域和频域是信号分析的两个方面,它们之间存在对应关系,可以 通过傅里叶变换相互转换。
数字信号处理基础
数字信号表示
MATLAB工作环境与界面
MATLAB工作环境
包括命令窗口、工作空间、命令历史窗口、当 前文件夹窗口等。
界面介绍
详细讲解MATLAB界面的各个组成部分,如菜 单栏、工具栏、编辑器窗口等。
基本操作
介绍如何在MATLAB环境中创建、保存、运行脚本和函数,以及如何进行基本 的文件操作。
基本数据类型与运算
数据统计描述性分析
描述性统计量
介绍均值、中位数、众数、方差、标准差等常见 描述性统计量的计算方法和意义。
数据分布形态
通过直方图、箱线图等图形展示数据的分布形态 ,帮助用户了解数据的整体特征。
数据间关系
探讨协方差、相关系数等统计量在揭示数据间关 系方面的应用。
数据可视化方法
二维图形绘制
详细讲解MATLAB中二维图形的绘制方法,包括线图、散点图、 柱状图等。
特征值与特征向量
特征值与特征向量的定义
设A为n阶方阵,若存在数λ和n维非零向量x,使得Ax=λx ,则称λ为A的特征值,x为A的对应于特征值λ的特征向量 。
特征值与特征向量的性质
包括特征值的和等于方阵对角线元素之和、特征值的积等 于方阵的行列式等性质。
MATLAB求解
使用MATLAB内置函数`eig`求解方阵的特征值和特征向量 。

matlab的有趣程序

曲线之美——塞尔曲线计算机科学与技术一班周星在我们的电脑中有这样一个屏幕保护程序叫做贝塞尔曲线,不知道大家留意过没有?因为这个学期我既学了c语言,又学了matlab,所以对编程很感兴趣。

我认为自己编出屏幕保护程序来既是一种学习知识的方式,也是一种开拓智力、放松心情的休闲方式。

上网查了贝塞尔曲线的数学表示式之后,我就开始用自己可怜的一点matlab知识解决它(其实本学期matlab学的挺多的,只是本人掌握的少,辜负了老师的期望)。

网上的给出的贝塞尔曲线的定义如下:在图形图像编程时,我们常常需要根据一系列已知点坐标来确定一条光滑曲线。

其中有些曲线需要严格地通过所有的已知点,而有些曲线却不一定需要。

在后者中,比较有代表性的一类曲线是贝塞尔曲线(Bézier Splines)。

在历史上,研究贝塞尔曲线的人最初是按照已知曲线参数方程来确定四个点的思路设计出这种矢量曲线绘制法。

涕淌为了向大家介绍贝塞尔曲线的公式,也故意把问题的已知和所求颠倒了一下位置:如果已知一条曲线的参数方程,系数都已知,并且两个方程里都含有一个参数t,它的值介于0、1之间,表现形式如下所示:x(t) = ax * t ^ 3 + bx * t ^ 2 + cx * t + x0y(t) = ay * t ^ 3 + by * t ^ 2 + cy * t + y0由于这条曲线的起点(x0,y0)是已知的,我们可以用以下的公式来求得剩余三个点的坐标:x1 = x0 + cx / 3x2 = x1 + ( cx + bx ) / 3x3 = x0 + cx + bx + axy1 = y0 + cy / 3y2 = y1 + ( cy + by ) / 3y3 = y0 + cy + by + ay你细细观察一下就知道了,无论方程的已知和所求是什么,总是有六个未知数,并且我们总能找到六个等式(记住(x0,y0)总是已知的),也就是说,上面的方法是完全可逆的,因此我们可以根据四个已知点坐标来反求曲线参数公式的系数。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
相关文档
最新文档