%《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])