%《MATLAB在電子信息課程中的應(yīng)用》第五章例5.2程序q502
% 含受控源的直流電路計(jì)算
% 電子工業(yè)出版社出版  陳懷琛 吳大正 高西全合著  2001年10月

clear,format compact
%fg512
R1=4;R2=4;R3=4;R4=2;		% 設(shè)置元件參數(shù)
is=2;k1=0.5;k2=4;
% 按A*X=B*is列寫(xiě)此電路的矩陣方程，其中X=[ua;ub;i1;i2]。
a11=1/R1+1/R2;a12=-1/R2;a13=0;a14=-k1;				% 設(shè)置系數(shù)矩陣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];										% 設(shè)置系數(shù)矩陣B
X=A\B*is;											% 解出X
i1=X(3),i2=X(4)									% 顯示要求的分量