%《MATLAB及在電子信息課程中的應(yīng)用》第五章例5.11程序q511
% 有受控源的交流電路：戴維南定理
% 電子工業(yè)出版社出版  陳懷琛 吳大正 高西全合著 2001年10月初版，2003年7月第二版
%
clear, format compact
R1=4;R2=2;R3=2;			% 設(shè)定元件參數(shù)
Xl1=10;Xl2=8;XM=4;Xc=8;
Us=10;Is=10;				% 設(shè)定信號(hào)源參數(shù)
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;		% 設(shè)定系數(shù)矩陣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];					% 設(shè)定系數(shù)矩陣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產(chǎn)生的復(fù)功率