%《MATLAB及其在理工課程中的應用指南》第八章例8－3－2程序
% 三相感應電機旋轉磁場的形成
% 西安電子科技大學出版社出版  陳懷琛編著 1999年10月初版，2004年11月第二版
%
%感應電動機定子三相合成磁場矢量圖
clear, clf, whitebg('w')
I=10; freq = 50; w = 2*pi*freq;               % 50 Hz角速度 (rad/s)
t = 0:1/5000:1.2/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相電流
% 建立三個分量磁場的表達式：
kmag = 1/I;                                         % 選適當的繞組常數，把最大磁場歸一化為1 
Baa = kmag * Ia * (cos(0) + j*sin(0)); 	            % 括號中的項表示A磁場空間方向為0度
Bbb = kmag * Ib * (cos(2*pi/3) + j*sin(2*pi/3)); 	% B磁場空間方向為120度
Bcc = kmag * Ic * (cos(-2*pi/3) + j*sin(-2*pi/3)); 	% C磁場空間方向為-120度
Bnet = Baa + Bbb + Bcc; 			                % 計算合成磁場 
% 分別畫出合成磁場Bnet和三相磁場Baa,Bbb,Bcc的矢量幅值和方向，
% Bnet為紅色，Baa為黑色，Bbb為藍色，Bcc為磁灰色。
for ii = 1:length(t)
    plot(Bnet,'k');	                                % 畫出合成磁場向量端點的軌跡作參考，它是一個園
    set(gcf,'color','w')                            % 置圖形背景色為白色
    hold on;
    % 畫出四個磁場相量，前三個方向固定，大小隨時間變化，其合成磁場則成為旋轉向量
    plot([0 real(Baa(ii))],[0 imag(Baa(ii))],'k','LineWidth',2);	
    plot([0 real(Bbb(ii))],[0 imag(Bbb(ii))],'b','LineWidth',2);
    plot([0 real(Bcc(ii))],[0 imag(Bcc(ii))],'m','LineWidth',2);
    plot([0 real(Bnet(ii))],[0 imag(Bnet(ii))],'r','LineWidth',3);
    axis square;   axis([-2 2 -2 2]);drawnow;       % 立即繪圖
    hold off;
end