%《MATLAB及其在理工課程中的應用指南》第八章例8－1－2程序
% 交流電路的穩態計算
% 西安電子科技大學出版社出版  陳懷琛編著 1999年10月初版，2004年11月第二版
%
R1=2;R2=3;R3=4;XL=2; XC1=3; XC2=5; 	        % 給出原始數據
us1=8, us2=6; us3=8; us4=15;		 	    % 給出原始數據
Y1=1/R1+1/(j*XL); 			                % 用復數表示各支路導納
Y2=1/R2-1/(j*XC1);
Y3=1/R3-1/(j*XC2);
A=[Y1+Y2,-Y2;-Y2,Y2+Y3]; 		            % 按線性方程組列出ua,ub的系數矩陣
B=[us1/(j*XL)+us2/R1; us3/R3+us4/(-j*XC2)-us2/R2]; 	% 列出線性方程組右端
U=A\B;ua=U(1),ub=U(2) 				        % 求ua,ub
I1=ua*Y1,I2=(ub-ua)*Y2,I3=ub*Y3,		    % 求各支路的I 
I1R=ua/R1,I1L=ua/(j*XL),
I2R=(ub-ua)/R2,I2C=(ub-ua)/(-j*XC1),
I3R=ub/R3,I3C=ub/(-j*XC2),
hold off
H=compass([ua,ub,I1,I2,I3]);		        % 畫相量圖
set(H,'linewidth',2) 	                    % 改變相量圖線寬
set(gcf,'color','w')                        % 置圖形背景色為白色