matlab动画演示17个实例

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 1

t=4;

dt=0.001;

w=0;

dw=0.001;

w=0;%设置球弹起的初始位置

%设置球弹起的高度

while t<=4.12

t=dt+t;

if w<=1

w=dw+w;

else

w=-1;

end

y=(-w*w)+1;

set(head,'xdata',t,'ydata',y);%设置球的运动

%set(body,'xdata',t,'ydata',y);%描绘轨迹线

drawnow;

end %结束程序

w=0;%设置球弹起的初始位置

%设置球弹起的高度

while t>=2.11

t=t-dt;

if w<=1

w=dw+w;

else

w=-1;

end

y=(-w*w)+1;

set(head,'xdata',t,'ydata',y);%设置球的运动%set(body,'xdata',t,'ydata',y);%描绘轨迹线drawnow;

end %结束程序

w=0;%设置球弹起的初始位置

%设置球弹起的高度

while t>=1.11

t=t-dt;

if w<=1

w=dw+w;

else

w=-1;

end

y=(-w*w)/4+1;

set(head,'xdata',t,'ydata',y);%设置球的运动%set(body,'xdata',t,'ydata',y);%描绘轨迹线drawnow;

end %结束程序

w=-0.71;%设置球弹起的初始位置

%设置球弹起的高度

while t>=-0.62

t=t-dt;

if w<=1

w=dw+w;

else

w=-1;

end

y=(-w*w)/2+1;

set(head,'xdata',t,'ydata',y);%设置球的运动%set(body,'xdata',t,'ydata',y);%描绘轨迹线drawnow;

end %结束程序

w=-0.71;%设置球弹起的初始位置

%设置球弹起的高度

while t>=-2.31

t=t-dt;

if w<=1

w=dw+w;

else

w=-1;

end

y=(-w*w)/2+0.75;

set(head,'xdata',t,'ydata',y);%设置球的运动

%set(body,'xdata',t,'ydata',y);%描绘轨迹线

drawnow;

end %结束程序

w=-0.71;%设置球弹起的初始位置

%设置球弹起的高度

while t>=-4

t=t-dt;

if w<=1

w=dw+w;

else

w=-1;

end

y=(-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)

end

end

end

2.嫦娥奔月演示程序

figure('name','嫦娥一号与月亮、地球关系');%设置标题名字s1=[0:.01:2*pi];

hold on;axis equal;%建立坐标系

axis off % 除掉Axes

r1=10;%月亮到地球的平均距离

r2=3;%嫦娥一号到月亮的平均距离

w1=1;%设置月亮公转角速度

w2=12%设置嫦娥一号绕月亮公转角速度

t=0;%初始时刻为０

pausetime=.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','.','markersize',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 1

set(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); %暂停一会

drawnow

end

3.曲柄连杆机构充气球动画演示程序

hf=figure('name','打气筒吹气球');

set(hf,'color','g');

axis([-10,10,-4,4]);

hold on

axis 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.7

c=0.7

a1=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','-','linewidth',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','-','linewidth',5);%设置打气筒

a14=line([-2.7;-2.7],[0.1;-0.1],'color','b','linestyle','-','linewidth',9);%设置气筒嘴

a16=line([-3.2;-3.2],[0.1;-0.1],'color','r','linestyle','-','linewidth',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.6 0.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=.0001

s=0;

ds=1;

t=0;

dt=0.015*pi;

while t<=15.68

t=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>200

s=0;

d=0;

end

set(gcf,'doublebuffer','on');%消除震动

plot(d,0,'color','r','marker','.','markersize',s);%画气球

set(gcf,'doublebuffer','on') %消除抖动

pause(pausetime); %暂停一会

drawnow;

end

end

x=[-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);

end

pausetime2=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);

end

pause(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<4

a=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);

end

a=0;

while a<.4

a=a+da;

set(b5,'xdata',[5.5;5.5],'ydata',[6.1-s-a,5.7-s-a]);

pause(pausetime4);

end

s=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);

end

5.电路演示程序

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,'erasemode','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','-','linewidth',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,'erasemode','xor');%绘制箭头横线

line(6.7,9.4,'color','b','linestyle','>','erasemode','xor','markersize',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,'erasemode','xor');

g=line(1,5.7,'color','b','linestyle','^','erasemode','xor','markersize',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;

end

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;

end

t=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;

end

pause(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;

end

t=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;

end

t=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;

pause(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;

end

t=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;

end

t=0;

while t<0.5

t=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;

end

pause(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;

end

t=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;

end

t=0;

while t<0.5 %画第五个箭头

t=t+0.5;

g=line(1,1,'color','b','linestyle','<','markersize',10,'erasemode','xor');

g=line(1,1,'color','b','linestyle','^','markersize',10,'erasemode','xor');

set(g,'xdata',1,'ydata',1);

drawnow;

end

t=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;

end

pause(0.5);%开关延时作用

nolight=line(10,4.3,'color','y','marker','.','markersize',40,'erasemode','xor'); end

6.电梯动画演示程序

figure('name','自控电梯');

axis([-2.0,15.0,-2.0,15.0])

hold on

fill([-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

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;

end

b=0;%设定电机反转的条件

db=0.05;

while s<10 %条件表达式 (当5

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;

end

for i=1:400

a1=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:400

a3=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:10

axis([-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);%给水线处小圆的颜色和尺寸

line(-0.35,-0.6,'Color','r','linestyle','.', 'markersize',20);%警戒线出小圆的颜色和尺寸

line([-0.45;-0.35],[0;0],'color','k','linewidth',2);%给水线处线条的颜色和宽度

line([-0.45;-0.35],[-0.6;-0.6],'color','k','linewidth',2);%警戒线处线条的颜色和宽度

line([-0.5;-0.5],[0.2,-1],'color','b','linewidth',15);%标杆的颜色和宽度

text(-0.8,0,'给水线');%文字标注“给水线”

text(-0.8,-0.6,'警戒线');%文字标注“警戒线”

text(-0.4,0.6,'防汛检测系统');%文字标注“防汛检测系统”

text(0.6,-0.9,'与江河连接');%文字标注“与江河连接”

water=patch(x1,y1,[0 1 1]);%设置水的颜色及运动路径

ball1=line(0.4,-0.6,'EraseMode','xor','Color','b','linestyle','.', 'markersize',100);%设置水槽中小球的颜色、大小和擦除方式

ball2=line(-0.3,-0,'EraseMode','xor','Color','r','linestyle','.', 'markersize',50);%设置标杆处小球的颜色、大小和擦除方式

色、大小和擦除方式

for i=1:120

a=-0.6+0.005*i;%设置系统的运动规律

y1=[a -0.8 -0.8 a];%设置水的上升运动过程

yy1=a;%设置水槽中小球的上升运动过程

yy2=-a-0.6%设置标杆处小球的上升运动过程

set(water,'ydata',y1);%设置水的上升运动

set(ball1,'ydata',yy1);%设置水槽中小球的上升运动

set(ball2,'ydata',yy2);%设置标杆处小球的上升运动

set(gan,'ydata',[yy2 yy1]);%设置两球之间的杆的运动

drawnow;

end%水的上升过程

for i=1:120

a=-0.005*i;%设置系统运动规律

y1=[a -0.8 -0.8 a];%设置水的下降运动过程

yy1=a;%设置水槽中小球的下降运动过程

yy2=-a-0.6%设置标杆处小球的下降运动过程

set(water,'ydata',y1);%设置水的下降运动

set(ball1,'ydata',yy1);%设置水槽中小球下降的运动

set(ball2,'ydata',yy2);%设置标杆处小球的下降运动

set(gan,'ydata',[yy2 yy1]);%设置两球之间的杆的下降运动

drawnow;

end%水的下降过程

water=patch(x1,y1,[0 1 1]);%设置水的颜色及运动路径

ball1=line(0.4,-0.6,'EraseMode','xor','Color','b','linestyle','.', 'markersize',100);%设置水槽中小球的颜色、大小和擦除方式

ball2=line(-0.3,-0,'EraseMode','xor','Color','r','linestyle','.', 'markersize',50);%设置标杆处小球的颜色、大小和擦除方式

一个简单的Matlab_GUI编程实例

Matlab GUI编程教程（适用于初学者） 1．首先我们新建一个GUI文件：如下图所示； 选择Blank GUI(Default) 2．进入GUI开发环境以后添加两个编辑文本框，6个静态文本框，和一个按钮，布置如下

图所示； 布置好各控件以后，我们就可以来为这些控件编写程序来实现两数相加的功能了。3．我们先为数据1文本框添加代码； 点击上图所示红色方框，选择edit1_Callback，光标便立刻移到下面这段代码的位置。 1. 2. 3.function edit1_Callback(hObject, eventdata, handles) 4.% hObject handle to edit1 (see GCBO) 5.% eventdata reserved - to be defined in a future version of MATLAB

6.% handles structure with handles and user data (see GUIDATA) 7.% Hints: get(hObject,'String') returns contents of edit1 as text 8.% str2double(get(hObject,'String')) returns contents of edit1 as a double 复制代码 然后在上面这段代码的下面插入如下代码： 1. 2.%以字符串的形式来存储数据文本框1的内容. 如果字符串不是数字，则现实空白内容input = str2num(get(hObject,'String')); %检查输入是否为空. 如果为空,则默认显示为0if (isempty(input)) set(hObject,'String','0')endguidata(hObject, handles); 复制代码 这段代码使得输入被严格限制，我们不能试图输入一个非数字。 4．为edit2_Callback添加同样一段代码 5 现在我们为计算按钮添加代码来实现把数据1和数据2相加的目的。 用3中同样的方法在m文件中找到pushbutton1_Callback代码段 如下； 1.function pushbutton1_Callback(hObject, eventdata, handles) 2.% hObject handle to pushbutton1 (see GCBO) 3.% eventdata reserved - to be defined in a future version of MATLAB 4.% handles structure with handles and user data (see GUIDATA) 复制代码

实例matlab-非线性规划-作业

实例matlab-非线性规划-作业

现代设计方法-工程优化理论、方法与设计 姓名 学号 班级 研 问题 ： 某厂向用户提供发动机，合同规定，第一、二、三季度末分别交货40台、60台、80台。每季度的生产费用为 （元），其中x 是该季生产的台数。若交货后有剩余，可用于下季度交货，但需支付存储费，每台每季度c 元。已知工厂每季度最大生产能力为100台，第一季度开始时无存货，设a=50、b=0.2、c=4，问工厂应如何安排生产计划，才能既满足合同又使总费用最低。讨论a 、b 、c 变化对计划的影响，并作出合理的解释。 问题的分析和假设： 问题分析：本题是一个有约束条件的二次规划问题。决策变量是工厂每季度生产的台数，目标函数是总费用（包括生产费用和存储费）。约束条件是生产合同，生产能力的限制。在这些条件下需要如何安排生产计划，才能既满足合同又使总费用最低。 问题假设： 1、工厂最大生产能力不会发生变化； 2、合同不会发生变更； 3、第一季度开始时工厂无存货； 4、生产总量达到180台时，不在进行生产； 5、工厂生产处的发动机质量有保证，不考虑退货等因素； 6、不考虑产品运输费用是否有厂家承担等和生产无关的因素。 符号规定： x1——第一季度生产的台数； x2——第二季度生产的台数； 180-x1-x2——第三季度生产的台数； y1——第一季度总费用； y2——第二季度总费用； y3——第三季度总费用； y ——总费用（包括生产费用和存储费）。 ()2bx ax x f +=

建模： 1、第一、二、三季度末分别交货40台、60台、80台； 2、每季度的生产费用为 （元）； 3、每季度生产数量满足40 ≤x1≤100，0≤x2≤100，100≤x1+x2 ≤180； 4、要求总费用最低，这是一个目标规划模型。 目标函数: y1 2111x b x a Z ?+?= y2()4012222-?+?+?=x c x b x a Z y3()()()10018018021221213 -+?+--?+--?=x x c x x b x x a Z y x x x x x x Z Z Z Z 68644.04.04.0149201 212221321--+++=++= 40≤x1≤100 0≤x2≤100 100≤x1+x2≤180 ()2 bx ax x f +=

MATLAB简单程序大全

MATLAB简单程序大全 求特征值特征向量 A=[2 3 4;1 5 9;8 5 2] det(A) A' rank(A) inv(A) rref(A) eig(A)%求特征值和特征向量 卫星运行问题 h=200,H=51000,R=6378; a=(h+H+2*R)/2; c=(H-h)/2; b=(a^2-c^2)^(1/2); e=c/a; f=sqrt(1-exp(2).*cos(t)^2); l=int(f,t,0,pi/2) L=4*a.*l 动态玫瑰线 n=3;N=10000; theta=2*pi*(0:N)/N; r=cos(n*theta); x=r.*cos(theta); y=r.*sin(theta); comet(x,y) 二重积分 syms x y f=x^2*sin(y); int(int(f,x,0,1),y,0,pi) ezmesh(f,[0,1,0,pi]) 函数画图 syms x;f=exp(-0.2*x)*sin(0.5*x); ezplot(f,[0,8*pi])

玫瑰线 theta=0:0.01:2*pi; r=cos(3*theta); polar(theta,r,'r') 求x^2+y^2=1和x^2+z^2=1所围成的体积 syms x y z R r=1; Z=sqrt(1-x^2); y0=Z; V=8*int(int(Z,y,0,y0),x,0,1) 求导数及图像 f='1/(5+4*cos(x))'; subplot(1,2,1);ezplot(f) f1=diff(f) subplot(1,2,2);ezplot(f1) 绕x轴旋转 t=(0:20)*pi/10; r=exp(-.2*t).*sin(.5*t); theta=t; x=t'*ones(size(t)); y=r'*cos(theta); z=r'*sin(theta); mesh(x,y,z) colormap([0 0 0]) 某年是否闰年 year=input('input year:='); n1=year/4; n2=year/100; n3=year/400; if n1==fix(n1)&n2~=fix(n2) disp('是闰年') elseif n1==fix(n1)&n3==fix(n3) disp('是闰年') else

用matlab制作简单仿真动画

用matlab制作简单仿真动画，并生成.avi格式的电影文件 MATALB知识点2008-05-01 13:17:54 阅读152 评论0 字号：大中小订阅 第一种形式：利用for循环，在一定时间内控制图形窗口图像的显示，产生一段动态的演示过程： 如下： set(gcf,'color','green'); grid on; set(gca,'zlim',[-10,10]'); set(gca,'xlim',[-10,10]); set(gca,'ylim',[-10,10]); for i=1:10 set(gca,'view',[-i*3.75,3*i] ); pause(0.2); end 第二种形式：利用moviein和movie函数，现将生成的动画存入一个由movien 函数定义的数组中，每一帧为数组的一个元素，最后用movie重复演示，movie后面的数字代表演示次数。 x=[-30:0.2:30]; y=[-30:0.2:30]; [x,y]=meshgrid(x,y); n=5; M = moviein(n); for i=1:n z=sin(sqrt(2*(x).^2+2*(y).^2)-2*pi*i/10); zz=plot3(x,y,z,'parent',gca); mesh(x,y,z); grid on; colormap([0,0.9,0.5]); light('position',[1,1,2],'style','local','color','white'); material([0.5,0.4,0.3,10,0.3]); set(gca,'zlim',[-10,10]'); M(i)=getframe(gca); end movie(M,20) 创建电影剪辑文件，并存储起来，如下： aviobj=avifile('文件名.avi','fps',3);%定义一个avi文件， %AVIOBJ = AVIFILE(FILENAME,'PropertyName',VALUE,'PropertyName',VALUE,...) %各属性详细说明见matlab帮助 for i=1:n %在当前窗体上生成一帧图像

BP神经网络matlab实例(简单而经典)

p=p1';t=t1'; [pn,minp,maxp,tn,mint,maxt]=premnmx(p,t); %原始数据归一化 net=newff(minmax(pn),[5,1],{'tansig','purelin'},'traingdx'); %设置网络,建立相应的BP网络 net.trainParam.show=2000; % 训练网络 net.trainParam.lr=0.01; net.trainParam.epochs=100000; net.trainParam.goal=1e-5; [net,tr]=train(net ,pn,tn); %调用TRAINGDM算法训练BP网络 pnew=pnew1'; pnewn=tramnmx(pnew,minp,maxp); anewn=sim(net,pnewn); %对BP网络进行仿真 anew=postmnmx(anewn,mint,maxt); %还原数据 y=anew'; 1、BP网络构建 （1）生成BP网络 = net newff PR S S SNl TF TF TFNl BTF BLF PF (,[1 2...],{ 1 2...},,,) R?维矩阵。 PR：由R维的输入样本最小最大值构成的2 S S SNl：各层的神经元个数。 [ 1 2...] { 1 2...} TF TF TFNl：各层的神经元传递函数。 BTF：训练用函数的名称。 （2）网络训练 [,,,,,] (,,,,,,) = net tr Y E Pf Af train net P T Pi Ai VV TV （3）网络仿真 = [,,,,] (,,,,) Y Pf Af E perf sim net P Pi Ai T {'tansig','purelin'},'trainrp'

有趣的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 1 t=4; dt=0.001; w=0; dw=0.001; w=0;%设置球弹起的初始位置 %设置球弹起的高度 while t<=4.12 t=dt+t; if w<=1 w=dw+w; else w=-1;

MATLAB程序设计报告--基于MATLAB动画播放及音乐播放

《MATLAB程序设计》课程设计报告 设计题目：基于MATLAB的动画演示 及背景音乐插入 专业：2011级通信工程 姓名（学号）：储兆雄1162310213 邓少林1162310214 徐凯越1162310223 指导教师：倪建军（博士/副教授） 时间：2013年12月20日

目录 1、设计目的 2、总体设计 3、具体设计（功能实现） 4、结果分析 5、改进方向 6、心得体会 文献 附录

1、设计目的 学会运用matlab工具箱实现matlab GUI设计，处理动画运行，以及添加背景音乐，并实现其动态操作，如继续、暂停等功能。 2、总体设计

主要包括：动画模块，音乐模块，动画显示模块 3、具体设计（功能实现） 1) 动画模块 （1）打开动画文件：从文件打开对话框选择动画程序，实现动画播放的可选择性 程序实现代码如下： function btnvopen_Callback(hObject, eventdata, handles) % hObject handle to btnvopen (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) global vfname %动画文件名 [vfname vpname vindex]=uigetfile('*.m','choose moive file'); len=length(vfname); if vindex set(handles.txtvname,'string',vfname(1:len-2)) end 打开对话框效果如下：

MATLAB动画演示效果

MATLAB动画演示效果一 %曲柄滑块机构 hf=figure('name','曲柄滑块机构'); set(hf,'color','g'); hold on axis([-6,6,-4,4]); grid on axis('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.7 c=0.7 a1=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 1 t=t+dt; if t>2*pi t=0; end lena1=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');%消除震动

Matlab简单实例学习

Matlab 程序代码 绘制y = 10e-1.5t sin( 7.75t ) 的函数图象 7.75 fv clear; t=0:0.02:10; f1=10/sqrt(7.75).*exp(-1.5*t); f2=sin(sqrt(7.75).*t); y=f1.*f2; plot(t,y,'-k',t,y,'ok'); xlabel('t');ylabel('y(t) ');title('函数图像') axis([-2 10 -0.5 2]) 拉氏变换 clear; clc; syms s t fs1 fs2 fs3 ft1 ft2 ft3; L=1,C=0.1,R=[1.5 3 5]; h1=1/(L*C*s^2+R(1)*C*s+1); h2=1/(L*C*s^2+R(2)*C*s+1); h3=1/(L*C*s^2+R(3)*C*s+1); fs1=h1*(1/s); fs2=h2*(1/s); fs3=h3*(1/s);

ft1=ilaplace(fs1,s,t); ft2=ilaplace(fs2,s,t); ft3=ilaplace(fs3,s,t); ezplot(t,ft1); hold on; ezplot(t,ft2); hold on; ezplot(t,ft3); 信号编码 对[1 1 0 1 1 1 0 1 0 0 1]进行编码。 clear; clc; c=[1 1 0 1 1 1 0 1 0 0 1] for i=1:length(c) if i==1 d1(i)=0;d2(i)=0; elseif i==2 d1(i)=c(i-1);d2(i)=c(i-1); elseif i==3 d1(i)=mod(c(i-1)+c(i-2),2); d2(i)=c(i-1); else d1(i)=mod(c(i-1)+c(i-2),2); d2(i)=mod(c(i-1)+c(i-3),2); end

Matlab动画及其在物理中的应用

实验报告实验课程名称计算物理 实验项目名称Matlab动画及其在物理中的应用 年级 09级 专业物理学 学生姓名 学号 理学院 实验时间：2012 年 4 月 4 日

学生实验室守则 一、按教学安排准时到实验室上实验课，不得迟到、早退和旷课。 二、进入实验室必须遵守实验室的各项规章制度，保持室内安静、整洁，不准在室内打闹、喧哗、吸烟、吃食物、随地吐痰、乱扔杂物，不准做与实验内容无关的事，非实验用品一律不准带进实验室。 三、实验前必须做好预习（或按要求写好预习报告），未做预习者不准参加实验。 四、实验必须服从教师的安排和指导，认真按规程操作，未经教师允许不得擅自动用仪器设备，特别是与本实验无关的仪器设备和设施，如擅自动用或违反操作规程造成损坏，应按规定赔偿，严重者给予纪律处分。 五、实验中要节约水、电、气及其它消耗材料。 六、细心观察、如实记录实验现象和结果，不得抄袭或随意更改原始记录和数据，不得擅离操作岗位和干扰他人实验。 七、使用易燃、易爆、腐蚀性、有毒有害物品或接触带电设备进行实验，应特别注意规范操作，注意防护；若发生意外，要保持冷静，并及时向指导教师和管理人员报告，不得自行处理。仪器设备发生故障和损坏，应立即停止实验，并主动向指导教师报告，不得自行拆卸查看和拼装。 八、实验完毕，应清理好实验仪器设备并放回原位，清扫好实验现场，经指导教师检查认可并将实验记录交指导教师检查签字后方可离去。 九、无故不参加实验者，应写出检查，提出申请并缴纳相应的实验费及材料消耗费，经批准后，方可补做。 十、自选实验，应事先预约，拟订出实验方案，经实验室主任同意后，在指导教师或实验技术人员的指导下进行。 十一、实验室内一切物品未经允许严禁带出室外，确需带出，必须经过批准并办理手续。

matlab函数计算的一些简单例子1

MATLAB作业一1、试求出如下极限。 （1） 23 25 (2)(3) lim (5) x x x x x x x ++ + →∞ ++ + ，（2） 23 3 1 2 lim () x y x y xy x y →- → + + ，（3） 22 22 22 1cos() lim ()x y x y x y x y e+ → → -+ + 解：（1）syms x； f=((x+2)^(x+2))*((x+3)^(x+3))/((x+5)^(2*x+5)) limit（f，x，inf） =exp（-5） （2）syms x y； f=(x^2*y+x*y^3)/(x+y)^3； limit(limit(f,x,-1),y,2) =-6； （3）syms x y； f=(1-cos(x^2+y^2))/(x^2+y^2)*exp(x^2+y^2)； limit(limit(f,x,0),y,0) =0 2、试求出下面函数的导数。 （1 ）() y x=, (2)22 atan ln() y x y x =+ 解; （1）syms x; f=sqrt(x*sin(x)*sqrt(1-exp(x))); g= diff(f,x); g== (sin(x)*(1 - exp(x))^(1/2) + x*cos(x)*(1 - exp(x))^(1/2) - (x*exp(x)*sin(x))/(2*(1 - exp(x))^(1/2)))/(2*(x*sin(x)*(1 - exp(x))^(1/2))^(1/2)) pretty(g)= (2)syms x y; f=atan(y/x)-log(x^2+y^2) pretty(-simple(diff(f,x)/diff(f,y)))= 2 x + y =------- x - 2 y (3) 假设1 cos u- =，试验证 22 u u x y y x ?? = ???? 。 解：syms x y; u=1/cos(sqrt(x/y)); diff(diff(u,x),y)-diff(diff(u,y),x)=0; 所以： 22 u u x y y x ?? = ????

matlab仿真实例

matlab 仿真实例 实验五MATLAB 及仿真实验一、控制系统的时域分析 (一)稳定性 1、系统传递函数为G(s)，试判断其稳定性。 程序： >> nu m=[3,2,5,4,6]; >> den=[1,3,4,2,7,2]; >> sys=tf( nu m,de n); >> figure(1); >> pzmap(sys); >> title(' 零极点图') 由图可知：在S 右半平面有极点，因此可知系统是不稳定的。 2、用MATLA 求 出 G(s)=(s A 2+2*s+2)/(s A 4+7*s A 3+5*s+2) 的极点。 程序及结果： >> sys=tf([1,2,2],[1,7,3,5,2]); >> p=pole(sys) 矿'. 赳 _ ■ —

-6.6553 0.0327 + 0.8555i 0.0327 - 0.8555i -0.4100 (二)阶跃响应 1、二阶系统G(s)=10/s A2+2*s+10 1)键入程序，观察并记录单位阶跃响应曲线： 程序： >> sys=tf(10,[1,2,10]); >> step(sys); >> title('G(s)=10/sA2+2*s+10 单位阶跃响应曲线') 2）计算系统闭环跟、阻尼比、无阻尼振荡频率，并记录程序及结果： >> sys=tf(10,[1,2,10]); >> p=pole(sys)

p = -1.0000 + 3.0000i -1.0000 - 3.0000i >> [wn,z]=damp(sys) wn = 3.1623 3.1623 z = 0.3162 0.3162 3)记录实际测取的峰值大小，峰值时间和过渡过程时间，并填表实际值理论值峰值Cmax 1.35s 峰值时间tp 1.05s 过渡时间+5% 3.54s ts +2% 3.18s 程序: >> sys=tf(10,[1,2,10]); >> step(sys); >> title('G(s)=10/sA2+2*s+10 单位阶跃响应曲线')

matlab函数计算的一些简单例子2

MATLAB 作业二 1、请将下面给出的矩阵A 和B 输入到MATLAB 环境中，并将它们转换成符号矩阵。若某一 矩阵为数值矩阵，另以矩阵为符号矩阵，两矩阵相乘是符号矩阵还是数值矩阵。 57651653 5501232310014325462564206441211346,3 9636623 51521210760077410120172440773 473 78 124867217110 7 681 5A B ???? ?????????? ????? ?==?????? ????? ?---????????--??? ? 解：A 转换为符号矩阵；a=sym（A） a=[5,7,6,5,1,6,5] [2,3,1,0,0,1,4][6,4,2,0,6,4,4][3,9,6,3,6,6,2][10,7,6,0,0,7,7][7,2,4,4,0,7,0][4,8,6,7,2,1,7]B 转换为符号矩阵；b=sym（B）b = [3,5,5,0,1,2,3][3,2,5,4,6,2,5][1,2,1,1,3,4,6][3,5,1,5,2,1,2][4,1,0,1,2,0,1][-3,-4,-7,3,7,8,12][1,-10,7,-6,8,1,5] 若某一矩阵为数值矩阵，另以矩阵为符号矩阵，两矩阵相乘是符号矩阵例；a*B= [48,3,64,48,159,106,194][17,-26,47,-8,62,26,59][48,-8,52,12,108,64,124][59,22,41,69,151,101,184][43,-22,91,13,175,121,220][22,39,4,53,88,94,147][75,11,115,36,151,70,151] 2、利用MATLAB 语言提供的现成函数对习题1中给出的两个矩阵进行分析，判定它们是否 为奇异矩阵，得出矩阵的秩、行列式、迹和逆矩阵，检验得出的逆矩阵是否正确。 解：由于a=det（A）=3.7396e+04；故A 是非奇异矩阵。B=det（B）=0，故B 是奇异矩阵； 由于a=rank(A)=7，故A 的秩为7；由于b=rank(B)=5，故B 的秩为5；由于a=trace(A)=27，b=trace(B)=26，故A，B 的迹为27,26；由a=inv(A)得A 的逆矩阵如下；

Matlab简单实例学习

Matlab 程序代码 绘制 1.5 10sin(7.75)7.75 t y e t -= 的函数图象。 fv clear; t=0:0.02:10; f1=10/sqrt(7.75).*exp(-1.5*t); f2=sin(sqrt(7.75).*t); y=f1.*f2; plot(t,y,'-k',t,y,'ok'); xlabel('t');ylabel('y(t) ');title('函数图像') axis([-2 10 -0.5 2]) 拉氏变换 clear; clc; syms s t fs1 fs2 fs3 ft1 ft2 ft3; L=1,C=0.1,R=[1.5 3 5]; h1=1/(L*C*s^2+R(1)*C*s+1); h2=1/(L*C*s^2+R(2)*C*s+1);

h3=1/(L*C*s^2+R(3)*C*s+1); fs1=h1*(1/s); fs2=h2*(1/s); fs3=h3*(1/s); ft1=ilaplace(fs1,s,t); ft2=ilaplace(fs2,s,t); ft3=ilaplace(fs3,s,t); ezplot(t,ft1); hold on; ezplot(t,ft2); hold on; ezplot(t,ft3); 信号编码 对[1 1 0 1 1 1 0 1 0 0 1]进行编码。 clear; clc; c=[1 1 0 1 1 1 0 1 0 0 1] for i=1:length(c) if i==1 d1(i)=0;d2(i)=0; elseif i==2 d1(i)=c(i-1);d2(i)=c(i-1); elseif i==3

Matlab动画程序 弹性蹦球演示过程

Matlab动画程序弹性蹦球演示过程 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)%设置球与地面接触面的颜

Matlab作业

《Matlab/Simulink电力系统建模与仿真》 上机实验报告 班级：15电气工程及其自动化二班 学号：154139240096 姓名：汤嘉旺 实验一：Powergui在简单电力系统潮流计算中的应用实例 1.1实验内容与要求 完成2机5节点电力系统的潮流计算，以2机5节点电力系统为模型进行Matlab/Simulink电力系统建模与仿真。并完成电力系统元件的模型选择、模型参数的计算及设置、计算结果及比较。 2机5节点电力系统图 1.2 Simulink建模原理图和主要模块参数设置 1、发电机模型 在该系统中的两台发电机均选用p.u.标准同步电机模块

“Synchronous Machine pu Standard”,该模块使用标幺值参数，以转子dq 轴建立的坐标系为参数，定子绕组为星形连接。 2、变压器模块 系统中的两台变压器均选用三相双绕组变压器模块 “Three-phase Transformer(Two Windings)”,采用Y-Y连接方式。 3、线路模块 系统中带有地导纳的线路选用三相“II”形等值模块 “Three Phase PI Section Line”,没有对地导纳的线路选用三相串联RLC支路模块“Three Phase Series RLC Branch”。 4、负荷模块 在SimPowerSystems库中，利用R、L、C的串联或并联组合，提供了两个静态三相负荷模块。这两种模块是用恒阻抗支路模拟负荷，在仿真时，在给定的频率下负荷阻抗为常数。 5、母线模块 选择带有测量元件的母线模型，及三相电压电流测量元 件“Three-Phase V-I Measurement”来模拟系统中的母线，同时方便测量流过线路的潮流，在线路元件的两端也设置了该元件。 2机5节点电力系统潮流计算仿真模型图

Matlab简单实例学习

Matlab 程序代码 绘 制 1.5sin(7.75)7.75t y e t -=的函数图象。 fv clear; t=0:0.02:10; f1=10/sqrt(7.75).*exp(-1.5*t); f2=sin(sqrt(7.75).*t); y=f1.*f2; plot(t,y,'-k',t,y,'ok'); xlabel('t');ylabel('y(t) ');title('函数图像') axis([-2 10 -0.5 2]) 拉氏变换 clear; clc; syms s t fs1 fs2 fs3 ft1 ft2 ft3; L=1,C=0.1,R=[1.5 3 5]; h1=1/(L*C*s^2+R(1)*C*s+1); h2=1/(L*C*s^2+R(2)*C*s+1);

h3=1/(L*C*s^2+R(3)*C*s+1); fs1=h1*(1/s); fs2=h2*(1/s); fs3=h3*(1/s); ft1=ilaplace(fs1,s,t); ft2=ilaplace(fs2,s,t); ft3=ilaplace(fs3,s,t); ezplot(t,ft1); hold on; ezplot(t,ft2); hold on; ezplot(t,ft3); 信号编码 对[1 1 0 1 1 1 0 1 0 0 1]进行编码。clear; clc; c=[1 1 0 1 1 1 0 1 0 0 1] for i=1:length(c) if i==1 d1(i)=0;d2(i)=0; elseif i==2 d1(i)=c(i-1);d2(i)=c(i-1); elseif i==3

matlab解题实例

Matlab 综合作业：运用matlab 建模解题 导弹追踪问题 1. 位于坐标原点的甲舰向位于x 轴上点A (1, 0)处的乙舰发射导弹， 导弹头始终对准乙舰．如果乙舰以最大的速度v 0(常数)沿平行于y 轴的直线行驶，导弹的速度是5v 0，求导弹运行的曲线方程．乙舰行驶多远时，导弹将它击中？ 解：假设t 时刻导弹的位置为P (x (t ), y (t ))，乙舰位于),1(0t v Q 由于导弹头 始终对准乙舰，故此时直线PQ 就是导弹的轨迹曲线弧OP 在点P 处的切线， 即有 x y t v y --=1'0 即 y y x t v +-=')1(0 (1) 又根据题意，弧OP 的长度为AQ 的5倍， 即 0d 5x x v t =? (2) 由(1),(2)消去t, 整理得模型: (3) '151 ")1(2y y x +=- 值条件为: 0)0(=y 0)0('=y 令y 1=y , y 2=y 1`，将方程（3）化为一阶微分方程组． 2 151 '')1(y y x +=- ? ?????-+==)1/(151''21221x y y y y 1.建立M 文件a1.m function dy=a1(x,y) dy=zeros(2,1); dy(1)=y(2); dy(2)=1/5*sqrt(1+y(1)^2)/(1-x); 2. 取x0=0，xf=0．9999，建立主程序ff6．m 如下： x0=0,xf=0.9999 [x,y]=ode15s('a1',[x0 xf],[0 0]);

plot(x,y(:,1),'b.') hold on y=0:0.01:2; plot(1,y,'b*') 运行得图： 结论: 导弹大致在（1，0．2）处击中乙舰.

【谷速软件】matlab源码-电梯动画演示程序

figure('name','自控电梯'); axis([-2.0,15.0,-2.0,15.0]) hold on fill([-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','-','linewi dth',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);

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','-','linewidt h',2); g2=line([7.7;7.7],[3;4],'color','b','linestyle','-','linewi dth',2); g3=line([7.4;7.7],[2;4],'color','b','linestyle','-','linewi dth',2); g4=line([7.4;7.4],[2;4],'color','b','linestyle','-','linewi dth',2); g5=line([7.1;7.4],[2;4],'color','b','linestyle','-','linewi dth',2); g6=line([7.1;7.1],[2;4],'color','b','linestyle','-','linewi dth',2); g7=line([6.8;7.1],[2;4],'color','b','linestyle','-','linewi dth',2); g8=line([6.8;6.8],[2;4],'color','b','linestyle','-','linewi dth',2);