dm13504.m

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

%dm13504
%绘制滑动平均数字滤波器幅频特性曲线
M=3;                             %定义延时器数目
num=ones(1,M)/M;                 %定义系统函数分子多项式系数向量
den=[1 zeros(1,M-1)];            %定义系统函数分母多项式系数向量
[H,w]=freqz(num,den,'whole');     %计算系统频率响应
Hm=abs(H);                       %计算幅频响应
w=w/pi;                          %归一化角频率
plot(w,Hm,'b-')                   %绘制幅频特性曲线
title('离散系统幅频特性曲线')
xlabel('\theta/\pi');
ylabel('Magnitude')
hold on                          %保持原波形
M=9;
num=ones(1,M)/M;
den=[1 zeros(1,M-1)];
[H,w]=freqz(num,den,'whole');
w=w/pi;
Hm=abs(H);
plot(w,Hm,'r:')
legend('b-','M=3','r:','M=9')         %标注图例