hr_type3.m

西交大刘树棠写的《数字信号处理》那本书的所有源代码

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