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

clear, format compact
Z1=-j*250;Z2=250;ki=0.5;Is=2;	% 設定元件參數
a11=1/Z1+1/Z2;a12=-1/Z2;a13=0;	% 設定系數矩陣A
a21=-1/Z2;a22=1/Z2;a23=-ki;
a31=1/Z1;a32=0;a33=-1;
A=[a11,a12,a13;a21,a22,a23;a31,a32,a33];
B=[1,0;0,1;0,0];					% 設定系數矩陣B
X0=A\B*[Is;0];		% X=[Ua;Ub;I1]=A\B*[Is;Ib];
Uoc=X0(2),				% Uoc等于Ib=0，Is=2∠0°時的Ub
abs(Uoc),angle(Uoc)*180/pi
X1=A\B*[0;1]; Zeq=X1(2),	% Zeq等于Is=0,Ib=1時的Ub
% 最大負載功率發生在ZLpmax=Zeq'時
ZLpmax=Zeq',PLmax=(abs(Uoc))^2/4/real(Zeq)