dm10502.m

来自「the code of the book come form the book 」· M 代码 · 共 26 行

% dm10502
% 二阶高通滤波器的频率响应
b=[0.04 0 0];	%生成向量b
a=[0.04 0.4 2];	%生成向量a
[h,w]=freqs(b,a,100);	%求频率响应函数H(jω),设定100个频率点
h1=abs(h);							%求幅频响应
h2=angle(h);						%求相频响应 
subplot(211);
plot(w,h1);
hold on
plot([7.0711 7.0711],[0 0.707],':');  %截止频率的位置
plot([0 7.0711],[0.707 0.707],':');
axis([0 40 0 1.1]);
grid
xlabel('角频率(W)');
ylabel('幅度');
title('H(jw)的幅频特性');
subplot(212);
plot(w,h2*180/pi);
axis([0 40 0 200]);
grid
xlabel('角频率(w)');
ylabel('相位(度)');
title('H(jw)的相频特性');
%End