%《MATLAB及在電子信息課程中的應用》第五章例5.2程序q502
% 含受控源的直流電路計算
% 電子工業出版社出版  陳懷琛 吳大正 高西全合著 2001年10月初版，2003年7月第二版
%
clear,format compact
R1=4;R2=4;R3=4;R4=2;		% 設置元件參數
is=2;k1=0.5;k2=4;
% 按A*X=B*is列寫此電路的矩陣方程，其中X=[ua;ub;i1;i2]。
a11=1/R1+1/R2;a12=-1/R2;a13=0;a14=-k1;				% 設置系數矩陣A
a21=-1/R2;a22=1/R2+1/R3+1/R4;a23=-k2/R3;a24=k1;
a31=1/R2;a32=-1/R2;a33=-1;a34=0;
a41=0;a42=1/R4;a43=0;a44=-1;
A=[a11,a12,a13,a14;a21,a22,a23,a24;a31,a32,a33,a34;a41,a42,a43,a44];
B=[1;0;0;0];										% 設置系數矩陣B
X=A\B*is;											% 解出X
i1=X(3),i2=X(4)									% 顯示要求的分量