%《MATLAB及其在理工課程中的應(yīng)用指南》第八章例8－3－2動(dòng)畫改進(jìn)程序ex832a
% 三相感應(yīng)電機(jī)旋轉(zhuǎn)磁場(chǎng)的形成（坐標(biāo)系不更新并畫出磁極）
% 西安電子科技大學(xué)出版社出版  陳懷琛編著 1999年10月初版，2004年11月第二版
%
%感應(yīng)電動(dòng)機(jī)定子三相合成磁場(chǎng)矢量圖
clear, clf, whitebg('w')
I=10; freq = 50; w = 2*pi*freq;               % 50 Hz角速度 (rad/s)
t = 0:1/5000:12.3/50;
Ia=I*sin(w*t); Ib=I*sin(w*t-2*pi/3); Ic=I*sin(w*t+2*pi/3); % A,B,C相電流
% 建立三個(gè)分量磁場(chǎng)的表達(dá)式：
kmag = 1/I;                % 選適當(dāng)?shù)睦@組常數(shù)，把最大磁場(chǎng)歸一化為1 
Baa = kmag * Ia * (cos(0) + j*sin(0)); 	% 括號(hào)中的項(xiàng)表示A磁場(chǎng)空間方向?yàn)?度
Bbb = kmag * Ib * (cos(2*pi/3) + j*sin(2*pi/3)); 	% B磁場(chǎng)空間方向?yàn)?20度
Bcc = kmag * Ic * (cos(-2*pi/3) + j*sin(-2*pi/3)); 	% C磁場(chǎng)空間方向?yàn)?120度
Bnet = Baa + Bbb + Bcc; 			% 計(jì)算合成磁場(chǎng) 
% 分別畫出合成磁場(chǎng)Bnet和三相磁場(chǎng)Baa,Bbb,Bcc的矢量幅值和方向，
% Bnet為紅色，Baa為黑色，Bbb為藍(lán)色，Bcc為磁灰色。
plot(Bnet,'k');	 % 畫出合成磁場(chǎng)向量端點(diǎn)的軌跡作參考，它是一個(gè)園
set(gcf,'color','w')                        % 置圖形背景色為白色
axis square;   axis([-2 2 -2 2]);hold on;
% 分別畫出三相磁極對(duì)
theta=linspace(-pi/8,pi/8,15);				% 給出單個(gè)磁極的張角數(shù)組
th1=theta;th2=theta+pi/3;th3=theta+2*pi/3;	% 給出六個(gè)磁極的幅角數(shù)組
th4=theta+pi;th5=theta+4*pi/3;th6=theta+5*pi/3;
r=1.1*abs(Bnet(1));						% 給出磁極的半徑，使之略大于Bnet的長(zhǎng)度
% 依次用三種顏色的寬線條畫出三對(duì)磁極
plot(r*cos(th1),r*sin(th1),'k',r*cos(th4),r*sin(th4),'k','LineWidth',6);
plot(r*cos(th2),r*sin(th2),'m',r*cos(th5),r*sin(th5),'m','LineWidth',6);
plot(r*cos(th3),r*sin(th3),'b',r*cos(th6),r*sin(th6),'b','LineWidth',6);
% 畫出四個(gè)磁場(chǎng)相量，前三個(gè)方向固定，大小隨時(shí)間變化，其合成磁場(chǎng)則成為旋轉(zhuǎn)向量
% 用改變圖形對(duì)象(直線)參數(shù)的方法做動(dòng)畫,坐標(biāo)不必更新,速度快,動(dòng)畫連續(xù)
% 先設(shè)置幾根直線線的初始位置,設(shè)置其圖柄名,圖柄參數(shù)'Erasemode'設(shè)為'xor'表示
% 當(dāng)該對(duì)象被賦予新值后,原來(lái)的對(duì)象圖形就不再保持
l1=plot([0 real(Baa(1))],[0 imag(Baa(1))]);	
set(l1,'erasemode','xor','LineWidth',2,'color','k')
l2=plot([0 real(Bbb(1))],[0 imag(Bbb(1))]);
set(l2,'erasemode','xor','LineWidth',2,'color','b')
l3=plot([0 real(Bcc(1))],[0 imag(Bcc(1))]);
set(l3,'erasemode','xor','LineWidth',2,'color','m')
l4=plot([0 real(Bnet(1))],[0 imag(Bnet(1))]);
set(l4,'erasemode','xor','LineWidth',2,'color','r')
% 不斷給對(duì)象的數(shù)據(jù)參數(shù)xdata,ydata賦予新值
for ii = 1:length(t)
   set(l1,'xdata', [0 real(Baa(ii))],'ydata',[0 imag(Baa(ii))])
   set(l2,'xdata', [0 real(Bbb(ii))],'ydata',[0 imag(Bbb(ii))])
   set(l3,'xdata', [0 real(Bcc(ii))],'ydata',[0 imag(Bcc(ii))])
   set(l4,'xdata', [0 real(Bnet(ii))],'ydata',[0 imag(Bnet(ii))])
   fft(randn(2^8));                             % 減慢速度
   drawnow
end
hold off;