samp4_1.m

来自「MATLAB7.x数字信号处理 光盘内容」· M 代码 · 共 19 行

%Samp4_1
clf;
b=[1 2 1]; a=[1 0.5 0.25];   %传递函数分子和分母多项式系数
m=0:length(b)-1;     n=0:length(a)-1;  %序列序号
K=512;   %采用512个点绘制幅频特性
k=1:K;w=pi*k/K;     %转换为归一化频率，Nyquist频率对应于pi。
num=b*exp(-j*m'*w);    %分子多项式值序列
den=a*exp(-j*n'*w);   %分母多项式值序列
H=num./den;        %计算复数频率响应
magH=20*log10(abs(H));   %转换为分贝
phaH=angle(H);   %提取复数频率响应的相位
subplot(2,1,1),plot(w/pi,magH);   %绘制幅频响应
xlabel('w/pi');ylabel('|H(jw)|/dB');title('幅频图')
grid on;
subplot(2,1,2)
plot(w/pi,phaH*180/pi);  %绘制相频响应
xlabel('w/pi');ylabel('相位角/^o');
title('相频图');grid on;