hr_type2.m

来自「<Matlab 7辅助信号处理技术与应用》的源码」· M 代码 · 共 14 行

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

M = length(h);
L = M/2;
b = 2*h(L:-1:1);
n = [1:1:L];
n=n-0.5;
w = [0:1:500]'*2*pi/500;
Hr = cos(w*n)*b';