📄 im3phase.m
字号:
% MOVIE OF INDUCTION MACHINE CONFIGURATION
clc,clear,close all
disp(date)
info= ...
[' WOUND ROTOR INDUCTION MACHINE '
' The 3-phase stator is represented by 3 concentrated coils displaced 120o from each other. '
' The wound rotor is also represented by 3 concentrated coils and it is rotating at the speed '
' n=0.8 pu. With balanced 3-phase voltage excitation of the stator, a rotating magnetic field '
' at 1.00 pu speed is established (represented by the yellow space vector) which overtakes the '
' rotor at slip speed (1-n), thereby inducing in the rotor coils slip-frequency voltages and '
' corresponding currents. ' ];
disp(info)
warning off
aviwave=avifile('im3phase.avi','compression','cinepak');
phis=pi/2;phim=0;g=pi/3;alpha=2*pi/3;rs=.75;rr=.45;n=0.8;b=1.25;
%
figure('units','normalized','Position',[.36 .09 .63 .72],'Name','INDUCTION MOTOR',...
'NumberTitle','off');
set(gcf,'color','w')
k=0;
for an=0:6:1800
k=k+1;
theta=an*pi/180;
fa1=rr*cos(phim-b);fb1=rr*sin(phim-b);
Fa1=rs*cos(phis);Fb1=rs*sin(phis);
fa2=rr*cos(phim+g-b);fb2=rr*sin(phim+g-b);
Fa2=rs*cos(phis+g);Fb2=rs*sin(phis+g);
fa3=rr*cos(phim+2*g-b);fb3=rr*sin(phim+2*g-b);
Fa3=rs*cos(phis+2*g);Fb3=rs*sin(phis+2*g);
fa4=rr*cos(phim+3*g-b);fb4=rr*sin(phim+3*g-b);
Fa4=rs*cos(phis+3*g);Fb4=rs*sin(phis+3*g);
fa5=rr*cos(phim+4*g-b);fb5=rr*sin(phim+4*g-b);
Fa5=rs*cos(phis+4*g);Fb5=rs*sin(phis+4*g);
fa6=rr*cos(phim+5*g-b);fb6=rr*sin(phim+5*g-b);
Fa6=rs*cos(phis+5*g);Fb6=rs*sin(phis+5*g);
fAo=cos(theta);
fBo=cos(theta-alpha);
fCo=cos(theta+alpha);
tA=abs(fAo);
tB=abs(fBo);
tC=abs(fCo);
fao=cos((1-n)*theta+pi/2-1.4);
fbo=cos((1-n)*theta-alpha+pi/2-1.4);
fco=cos((1-n)*theta+alpha+pi/2-1.4);
ta=abs(fao);
tb=abs(fbo);
tc=abs(fco);
q=0:pi/20:2*pi;
fill(.9*cos(q),0.9*sin(q),[.6 .6 .6],'LineWidth',1.5),hold on
fill(.62*cos(q),0.62*sin(q),'w','LineWidth',1.5)
fill(.58*cos(q),0.58*sin(q),[.6 .6 .6],'LineWidth',1.5)
plot(Fa1,Fb1,'ro','MarkerSize',18,'MarkerEdgeColor','k','MarkerFacecolor','w')
plot(Fa1,Fb1,'ro','MarkerSize',tA*16,'MarkerFacecolor','r')
if fAo>0
plot(Fa1,Fb1,'k.','MarkerSize',10,'LineWidth',4)
else
plot(Fa1,Fb1,'kx','MarkerSize',10,'LineWidth',2)
end
plot(Fa2,Fb2,'bo','MarkerSize',18,'MarkerEdgeColor','k','MarkerFacecolor','w')
plot(Fa2,Fb2,'bo','MarkerSize',tC*16,'MarkerFacecolor',[.5 1 .8])
if fCo>0
plot(Fa2,Fb2,'kx','MarkerSize',10,'LineWidth',2)
else
plot(Fa2,Fb2,'k.','MarkerSize',10,'LineWidth',4)
end
plot(Fa3,Fb3,'go','MarkerSize',18,'MarkerEdgeColor','k','MarkerFacecolor','w')
plot(Fa3,Fb3,'go','MarkerSize',tB*16,'MarkerFacecolor','g')
if fBo>0
plot(Fa3,Fb3,'k.','MarkerSize',10,'LineWidth',4)
else
plot(Fa3,Fb3,'kx','MarkerSize',10,'LineWidth',2)
end
plot(Fa4,Fb4,'ro','MarkerSize',18,'MarkerEdgeColor','k','MarkerFacecolor','w')
plot(Fa4,Fb4,'ro','MarkerSize',tA*16,'MarkerFacecolor','r')
if fAo<0
plot(Fa4,Fb4,'k.','MarkerSize',10,'LineWidth',4)
else
plot(Fa4,Fb4,'kx','MarkerSize',10,'LineWidth',2)
end
plot(Fa5,Fb5,'bo','MarkerSize',18,'MarkerEdgeColor','k','MarkerFacecolor','w')
plot(Fa5,Fb5,'bo','MarkerSize',tC*16,'MarkerFacecolor',[.5 1 .8])
if fCo<0
plot(Fa5,Fb5,'kx','MarkerSize',10,'LineWidth',2)
else
plot(Fa5,Fb5,'k.','MarkerSize',10,'LineWidth',4)
end
plot(Fa6,Fb6,'go','MarkerSize',18,'MarkerEdgeColor','k','MarkerFacecolor','w')
plot(Fa6,Fb6,'go','MarkerSize',tB*16,'MarkerFacecolor','g')
if fBo<0
plot(Fa6,Fb6,'k.','MarkerSize',10,'LineWidth',4)
else
plot(Fa6,Fb6,'kx','MarkerSize',10,'LineWidth',2)
end
plot(fa1,fb1,'ro','MarkerSize',18,'MarkerEdgeColor','k','MarkerFacecolor','w')
plot(fa1,fb1,'ro','MarkerSize',ta*16,'MarkerFacecolor','r')
if fao>0
plot(fa1,fb1,'k.','MarkerSize',10,'LineWidth',4)
else
plot(fa1,fb1,'kx','MarkerSize',10,'LineWidth',2)
end
plot(fa2,fb2,'bo','MarkerSize',18,'MarkerEdgeColor','k','MarkerFacecolor','w')
plot(fa2,fb2,'bo','MarkerSize',tc*16,'MarkerFacecolor',[.5 1 .8])
if fco>0
plot(fa2,fb2,'k.','MarkerSize',10,'LineWidth',4)
else
plot(fa2,fb2,'kx','MarkerSize',10,'LineWidth',2)
end
plot(fa3,fb3,'go','MarkerSize',18,'MarkerEdgeColor','k','MarkerFacecolor','w')
plot(fa3,fb3,'go','MarkerSize',tb*16,'MarkerFacecolor','g')
if fbo>0
plot(fa3,fb3,'k.','MarkerSize',10,'LineWidth',4)
else
plot(fa3,fb3,'kx','MarkerSize',10,'LineWidth',2)
end
plot(fa4,fb4,'ro','MarkerSize',18,'MarkerEdgeColor','k','MarkerFacecolor','w')
plot(fa4,fb4,'ro','MarkerSize',ta*16,'MarkerFacecolor','r')
if fao<0
plot(fa4,fb4,'k.','MarkerSize',10,'LineWidth',4)
else
plot(fa4,fb4,'kx','MarkerSize',10,'LineWidth',2)
end
plot(fa5,fb5,'bo','MarkerSize',18,'MarkerEdgeColor','k','MarkerFacecolor','w')
plot(fa5,fb5,'bo','MarkerSize',tc*16,'MarkerFacecolor',[.5 1 .8])
if fco<0
plot(fa5,fb5,'k.','MarkerSize',10,'LineWidth',4)
else
plot(fa5,fb5,'kx','MarkerSize',10,'LineWidth',2)
end
plot(fa6,fb6,'go','MarkerSize',18,'MarkerEdgeColor','k','MarkerFacecolor','w')
plot(fa6,fb6,'go','MarkerSize',tb*16,'MarkerFacecolor','g')
if fbo<0
plot(fa6,fb6,'k.','MarkerSize',10,'LineWidth',4)
else
plot(fa6,fb6,'kx','MarkerSize',10,'LineWidth',2)
end
g3=plot(0,0,'ok',[0 .6 .55 .6 .55],[0 0 .06 0 -.06],'-y','LineWidth',8);
rotate(g3,[0 0 1],an-70,[0 0 0])
text(-1.06,1.13,'WOUND-ROTOR INDUCTION MOTOR','FontSize',16,'FontWeight','bold','color','blue')
set(gca,'Xlim',[-1.2 1.2],'Ylim',[-1.2 1.2],'visible','off')
axis square
%
hold off
%
f=getframe(gca);
aviwave=addframe(aviwave,f);
phim=n*theta;
M(k)=getframe;
end
aviwave=close(aviwave);
uicontrol('Units','normalized', ...
'Style','PushButton',...
'BackgroundColor',[.2 .3 .8], ...
'Callback','movie(M)', ...
'FontSize',10, ...
'FontWeight','bold', ...
'ForegroundColor',[1 1 1], ...
'Position',[0.31 0.15 0.09 .06], ...
'String','Repeat', ...
'TooltipString','Repeat');
uicontrol('Units','normalized', ...
'Style','PushButton',...
'BackgroundColor',[.2 .3 .8], ...
'Callback','close,clear,close all', ...
'FontSize',10, ...
'FontWeight','bold', ...
'ForegroundColor',[1 1 1], ...
'Position',[0.67 0.15 0.09 .06], ...
'String','Quit', ...
'TooltipString','Quit');
%%%%%%%
⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -