%《MATLAB在電子信息課程中的應用》第五章例5.1程序q501
% 直流電路的穩態計算
% 電子工業出版社出版  陳懷琛 吳大正 高西全合著  2001年10月

clear,format compact
R1=2;R2=4;R3=12;R4=4;R5=12;R6=4;R7=2;		%為給定元件賦值
% 解問題（1）
display('解問題（1）')
a11=R1+R2+R3;a12=-R3;a13=0;	%將系數矩陣各元素賦值
a21=-R3;a22=R3+R4+R5; a23=-R5; 
a31=0;a32=-R5;a33=R5+R6+R7;
b1=1;b2=0;b3=0;
us=sym('us')                  % 輸入解(1)的已知條件
A=[a11,a12,a13;a21,a22,a23;a31,a32,a33] 				% 列出系數矩陣A
B=[b1;0;0]; I=A\B*us;						% I=[ia;ib;ic]
ia=I(1);ib=I(2);ic=I(3);
i3=ia-ib,u4=R4*ib,u7=R7*ic			% 解出所需變量
% 利用電路的線性性質及問題（1）的解
display('解問題（2）')
u42=sym('u42');
k1=i3/us;k2=u4/us;k3=u7/us;			% 由問題（1）得出待求量與us的比例系數
us2=u42/k2,i32=k1/k2*u42,u72=k3/k2*u42	% 按比例方法求出所需變量