%《MATLAB及在電子信息課程中的應(yīng)用》第五章例5.3程序q503
% 直流電路的戴維南定理
% 電子工業(yè)出版社出版  陳懷琛 吳大正 高西全合著 2001年10月初版，2003年7月第二版
%
clear,format compact
R1=4;R2=2;R3=4;R4=8;		% 設(shè)置元件參數(shù)
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;  				% 設(shè)置系數(shù)矩陣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];						% 設(shè)置系數(shù)矩陣B
% 方法1：令ia=0,求uoc=u(3);再令is1=is2=0,設(shè)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				    % 求最大負(fù)載功率
% 也可設(shè)RL為一數(shù)組,求出的負(fù)載功率也為一數(shù)組,畫出曲線找極大值
RL=0:10,p=(RL*uoc./(Req+RL)).*uoc./(Req+RL),    % 設(shè)RL序列,求其功率
figure(1),plot(RL,p),grid			            % 畫出功耗隨RL變化的曲線如圖5.3B(a)
set(gcf,'color','w')                            % 設(shè)置圖形背景色為白色
% 方法2：設(shè)一個(gè)ia序列,計(jì)算一個(gè)ua序列,用線性擬合求出其等效開路電壓和等效內(nèi)阻
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	                % 線性擬合,見(jiàn)圖5.3B(b)
set(gcf,'color','w')                            % 設(shè)置圖形背景色為白色
c=polyfit(ia,u,1)	;		                    % ua=c(2)*ia+c(1),用擬合函數(shù)求c(1),c(2)
uoc=c(1),Req=c(2)