%《MATLAB在電子信息課程中的應用》第五章例5.11程序q511
% 有受控源的交流電路：戴維南定理
% 電子工業出版社出版  陳懷琛 吳大正 高西全合著  2001年10月

clear, format compact
R1=4;R2=2;R3=2;			% 設定元件參數
Xl1=10;Xl2=8;XM=4;Xc=8;
Us=10;Is=10;				% 設定信號源參數
Y1=1/R1+1/(-j*Xc); Y2=1/j/(Xl1-XM); Y3=1/j/XM;
Y4=1/(R2+j*(Xl2-XM)); Y5=1/R3;
a11=Y1+Y2; a12=-Y2; a13=0;		% 設定系數矩陣A
a21=-Y2; a22=Y2+Y3+Y4; a23=-Y4;
a31=0; a32=-Y4; a33=Y4+Y5;
A=[a11,a12,a13;a21,a22,a23;a31,a32,a33];
B=[Us/R1;0;Is];					% 設定系數矩陣B
X=A\B					% X=[Ua;Ub;Uc]=A\[Us/R1;0;Is];
I1=(Us-X(1))/R1;
Pus=Us*I1',Pis=X(3)*Is'	% 求Us和Is產生的復功率