4-9.m

来自「是《MATLAB数字信号处理与应用》一书的源代码,该书由李正周编著 清华大学出版」· M 代码 · 共 16 行

%例程4－9   Ⅲ型滤波器的幅频函数

function [Hr,w,c,L]=hr_type3(h);         
%计算所设计的3型滤波器的振幅响应
%Hr=振幅响应
%b=3型滤波器的系数
%L=Hr的阶次
%h=3型滤波器的单位冲击响应

M=length(h);
L=(M-1)/2;
c= [2*h(L+1:-1:1)];
n=[0:1:L];
w=[0:1:500]'*2*pi/500;
Hr=sin(w*n)*c';