q520.m

MATLAB及在电子信息课程中的应用例题索引与程序

%《MATLAB在电子信息课程中的应用》第五章例5.20程序q520
% 桥梯形全通网络的计算
% 电子工业出版社出版  陈怀琛 吴大正 高西全合著  2001年10月

clear
for k=1:101						  % 对各个频点分别计算其频率特性
   w=0.02*(k-1)+1e-8;		  % 为避免奇异值,加一个微量
   C1=2;C3=2;C2=4/3;L2=1;L4=1;	  % 元件赋值
   Z1=1./(w*C1*j);Z3=1./(w*C3*j);	% 阻抗计算
   Z2=L2*j*w+1./(j*w*C2);
   Z4=L4*j*w;ZL=1;
   Ya=[1./Z4,-1./Z4;-1./Z4,1./Z4];	    % a网络参数Y计算
   Zb=[Z1+Z2,Z2;Z2,Z3+Z2];	Yb=inv(Zb);   	% b网络参数Y计算
   Y=Ya+Yb;	Z=inv(Y);						% a,b网络并联
   H(k)=Z(2,1)*ZL./(det(Z)+Z(1,1).*ZL);	  % 求U2/U1
   Zin(k)=(det(Z)+Z(1,1)*ZL)./(Z(2,2)+ZL);  % 求输入阻抗
end
w1=[0:100]*0.02;	  % 设定绘图横坐标数组w1,注意w只是一个标量
figure(1)
subplot(2,1,1),plot(w1,abs(H))		   % 画幅频特性
axis([0 2 0 2])      % 因abs(H)处处相同,系统无法自动取纵坐标,要人为设定
subplot(2,1,2),plot(w1,unwrap(angle(H)))
% 画相频特性,unwrap函数用以避免因取主角使小于-pi的值跳到+pi
figure(2)							% 画网络函数实频特性和虚频特性
subplot(2,1,1),plot(w1,real(H))
subplot(2,1,2),plot(w1,imag(H))
figure(3)						% 画输入阻抗的实频特性和虚频特性
subplot(2,1,1),plot(w1,real(Zin))
axis([0 2 0 2])					% 理由同前
subplot(2,1,2),plot(w1,imag(Zin),'linewidth',2)
axis([0 2 -0.1 0.2])