%《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('解問題（3）')
syms ia ib ic us
s1='(R1+R2+R3)*ia-R3*ib=us'	%將系數矩陣各元素賦值
s2='-R3*ia+(R3+R4+R5)*ib-R5*ic=0' 
s3='-R5*ia+(R5+R6+R7)*ib=0'
[ia,ib,ic]=solve(s1,s2,s3)						% 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	% 按比例方法求出所需變量