%《MATLAB及其在理工課程中的應用指南》(十一五規劃版)第八章例8－1－5程序
% 交流電路的穩態計算
% 西安電子科技大學出版社出版  陳懷琛編著 2007年5月
%
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',3) 	% 改變相量圖線寬
set(gcf,'color','w')