q515.m

MATLAB及在电子信息课程中的应用,电子工业出版社出版 陈怀琛 吴大正 高西全合著

%《MATLAB及在电子信息课程中的应用》第五章例5.15程序q515
% 频率响应：二阶电路
% 电子工业出版社出版  陈怀琛 吴大正 高西全合著 2001年10月初版，2003年7月第二版
%
clear, format compact
for Q=[1/3,1/2,1/sqrt(2),1,2,5]
   w=logspace(-1,1,50);			            % 设定频率数组w
   H=1./(1+j*w/Q+(j*w).^2);		            % 求复频率响应
   figure(1)
   subplot(2,1,1),plot(w,abs(H)),hold on	% 绘制幅频特性
   subplot(2,1,2),plot(w,angle(H)),hold on	% 绘制相频特性
   figure(2),subplot(2,1,1),				% 绘制对数频率特性
   semilogx(w,20*log10(abs(H))),hold on     % 纵坐标为分贝
   subplot(2,1,2),
   semilogx(w,angle(H)),hold on	            % 绘制相频特性
end
figure(1)
subplot(2,1,1),   
set(gcf,'color','w')                        % 设置图形背景色为白色
grid on,xlabel('w'),ylabel('abs(H)')
axis([0,5,min(abs(H)),max(abs(H))])
subplot(2,1,2),   
grid on,xlabel('w'),ylabel('angle(H)')
axis([0,5,min(angle(H)),max(angle(H))])
figure(2)
subplot(2,1,1),   
set(gcf,'color','w')                        % 设置图形背景色为白色
grid on,xlabel('w'),ylabel('分贝')
subplot(2,1,2),   
grid on,xlabel('w'),ylabel('angle(H)')