%《MATLAB在電子信息課程中的應(yīng)用》第五章例5.15程序q515
% 頻率響應(yīng)：二階電路
% 電子工業(yè)出版社出版  陳懷琛 吳大正 高西全合著  2001年10月

clear, format compact
for Q=[1/3,1/2,1/sqrt(2),1,2,5]
   w=logspace(-1,1,50);			% 設(shè)定頻率數(shù)組w
   H=1./(1+j*w/Q+(j*w).^2);		% 求復(fù)頻率響應(yīng)
   figure(1)
   subplot(2,1,1),plot(w,abs(H)),hold on	% 繪制幅頻特性
   subplot(2,1,2),plot(w,angle(H)),hold on	% 繪制相頻特性
   figure(2)				% 繪制對(duì)數(shù)頻率特性
   subplot(2,1,1),semilogx(w,20*log10(abs(H))),hold on %縱坐標(biāo)為分貝
   subplot(2,1,2),semilogx(w,angle(H)),hold on	% 繪制相頻特性
end
figure(1)
subplot(2,1,1),   grid,xlabel('w'),ylabel('abs(H)')
axis([0,5,min(abs(H)),max(abs(H))])
subplot(2,1,2),   grid,xlabel('w'),ylabel('angle(H)')
axis([0,5,min(angle(H)),max(angle(H))])
figure(2)
subplot(2,1,1),   grid,xlabel('w'),ylabel('分貝')
subplot(2,1,2),   grid,xlabel('w'),ylabel('angle(H)')