hr_type2.m

来自「原名：数字信号处理基础及MATLAB实现,解压缩密码：lfjk2002@163.」· M 代码 · 共 9 行

function [Hr,w,b,L]=hr_type2(h);
%线性相位滤波器类型2，计算滤波器振幅响应
%Hr振幅响应；w在[0,2pi]之间计算Hr的500个频率点；
%b为2型低通滤波器系数；L为Hr的阶次；h为2型低通滤波器的脉冲响应。
M=length(h);L=M/2;
b=2*[h(L:-1:1)];            %1乘(L+1)行向量
n=[1:1:L];n=n-0.5;          %(L+1)乘1列向量
w=[0:1:500]'*2*pi/500;
Hr=cos(w*n)*b';