hr_type1.m

来自「the file contain many matlab signal sour」· M 代码 · 共 16 行

function [Hr,w,a,L]=hr_type1(h);
% Computes Amplitude response of Type-1 LP FIR filter
% ---------------------
% [Hr,w,a,L]=hr_type1(h)
% Hr=Amplitude Response
% w=500 frequecies between [0 pi] over which Hr is computed
% a=Type-1 LP filter coefficients
% L=Order of Hr
% h=Type-1 LP filter impulse response
% 
M=length(h);
L=(M-1)/2;
a=[h(L+1) 2*h(L:-1:1)];
n=[0:L];
w=[0:500]'*pi/500;
Hr=cos(w*n)*a';