q501.m

MATLAB及在电子信息课程中的应用,电子工业出版社出版 陈怀琛 吴大正 高西全合著

%《MATLAB及在电子信息课程中的应用》第五章例5.1程序q501
% 直流电路的稳态计算
% 电子工业出版社出版  陈怀琛 吴大正 高西全合著 2001年10月初版，2003年7月第二版
%
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=input('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=input('给定u42= ');
k1=i3/us;k2=u4/us;k3=u7/us;			% 由问题（1）得出待求量与us的比例系数
us2=u42/k2,i32=k1/k2*u42,u72=k3/k2*u42	% 按比例方法求出所需变量