q503.m

matlab在电子信息中的应用和实例中的源代码

%《MATLAB及在电子信息课程中的应用》第五章例5.3程序q503
% 直流电路的戴维南定理
% 电子工业出版社出版  陈怀琛 吴大正 高西全合著 2001年10月初版，2003年7月第二版
%
clear,format compact
R1=4;R2=2;R3=4;R4=8;		% 设置元件参数
is1=2;is2=0.5;
% 按A*X=B*is列写此电路的矩阵方程,其中X=[u1;u2;ua];is=[is1;is2;ia]
a11=1/R1+1/R4;a12=-1/R1;a13=-1/R4;  				% 设置系数矩阵A
a21=-1/R1;a22=1/R1+1/R2+1/R3;a23=-1/R3;
a31=-1/R4;a32=-1/R3;a33=1/R3+1/R4;
A=[a11,a12,a13;a21,a22,a23;a31,a32,a33];
B = [1,1,0;0,0,0;0,-1,1];						% 设置系数矩阵B
% 方法1：令ia=0,求uoc=u(3);再令is1=is2=0,设ia=1,求Req=ua/ia=ua.
X1=A\B*[is1;is2;0];uoc=X1(3)
X2=A\B*[0;0;1];Req=X2(3)
RL=Req;P=uoc^2*RL/(Req+RL)^2				    % 求最大负载功率
% 也可设RL为一数组,求出的负载功率也为一数组,画出曲线找极大值
RL=0:10,p=(RL*uoc./(Req+RL)).*uoc./(Req+RL),    % 设RL序列,求其功率
figure(1),plot(RL,p),grid			            % 画出功耗随RL变化的曲线如图5.3B(a)
set(gcf,'color','w')                            % 设置图形背景色为白色
% 方法2：设一个ia序列,计算一个ua序列,用线性拟合求出其等效开路电压和等效内阻
for k=1:21
   ia(k)=(k-1)*0.1;
   X=A\B*[is1;is2;ia(k)];				        % 定义X=[u1;u2;ua]
   u(k)=X(3); end
figure(2),plot(ia,u,'x'),grid	                % 线性拟合,见图5.3B(b)
set(gcf,'color','w')                            % 设置图形背景色为白色
c=polyfit(ia,u,1)	;		                    % ua=c(2)*ia+c(1),用拟合函数求c(1),c(2)
uoc=c(1),Req=c(2)