q520.m

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

%《MATLAB及在电子信息课程中的应用》第五章例5.20程序q520
% 桥梯形全通网络的计算
% 电子工业出版社出版  陈怀琛 吴大正 高西全合著 2001年10月初版，2003年7月第二版
%
clear,close all
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))	% 画幅频特性
set(gcf,'color','w')            % 设置图形背景色为白色
axis([0 2 0 2])                 % 因abs(H)处处相同,系统无法自动取纵坐标,要人为设定
subplot(2,1,2),plot(w1,unwrap(angle(H)))
% 画相频特性,unwrap函数用以避免因取主角使小于-pi的值跳到+pi
figure(2)						% 画网络函数实频特性和虚频特性
set(gcf,'color','w')            % 设置图形背景色为白色
subplot(2,1,1),plot(w1,real(H))
subplot(2,1,2),plot(w1,imag(H))
figure(3)						% 画输入阻抗的实频特性和虚频特性
set(gcf,'color','w')            % 设置图形背景色为白色
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])