squirrelmovie1.m

来自「关于电机的.m程序」· M 代码 · 共 44 行

%    squirrel-cage machine picture
clc,clear,close all
figure('units','normalized','numbertitle','off','name','AC MACHINE','menubar','none','units','normalized',...
    'position',[0.1 .1 .8 .8]);
%
set(gcf,'color','w')
t=0;e=0;
dir=[0 0 1];

warning off
avisqui=avifile('squimovie5.avi');

for p=1:39
[X,Y,Z]=cylinder(1.44,160);h=surf(X,Y,Z,'facecolor',[.7 .8 .8],'edgecolor',[.7 .8 .8 ]);hold on    
for k=1:24;
 or0=[cos(t+e) sin(t+e) 0];  
 [X,Y,Z]=cylinder(.04,50);h0=mesh(X,Y,Z,'facecolor','b','edgecolor','b');%hold on
 rotate(h0,[0 90],90,or0);
 t=t+pi/12;
end 

q=0:pi/20:2*pi;
h1=plot(sqrt(2)*cos(q) ,sqrt(2)*sin(q),'b','linewidth',8);
rotate(h1,dir,90,[0 0 -1])
h2=plot(sqrt(2)*cos(q) ,sqrt(2)*sin(q),'b','linewidth',8);
rotate(h2,dir,90,[0 0 0])
h3=plot(2*cos(q) ,2*sin(q),'color',[.8 .8 .8],'linewidth',44);
%h3=fill(2*cos(q) ,2*sin(q),[.8 .8 .8]);
rotate(h3,dir,90,[0 0 -1])
%fill(.9*cos(q),0.9*sin(q),[.6 .6 .6],'LineWidth',1.5)
[X,Y,Z]=cylinder(1.7,160);h=surf(X,Y,Z,'facecolor',[.8 .8 .8],'edgecolor',[.9 .9 .9 ]);
    [X,Y,Z]=cylinder(2.3,160);h=surf(X,Y,Z,'facecolor',[.8 .8 .8],'edgecolor',[.9 .9 .9]);

set(gca,'Xlim',[-2.5 2.5],'Ylim',[-2.5 2.5],'visible','off')
axis equal;

 hold off
      %          
      f=getframe(gca);
      avisqui=addframe(avisqui,f);
 e=e+pi/(12*40);
end
  avisqui=close(avisqui);   
rotate3d