hr_type3.m

来自「DSP2 using matlab v4」· M 代码 · 共 33 行

function [Hr,w,c,L] = Hr_Type3(h);

% Computes Amplitude response Hr(w) of a Type-3 LP FIR filter

% -----------------------------------------------------------

% [Hr,w,c,L] = Hr_Type3(h)

% Hr = Amplitude Response

%  w = frequencies between [0 pi] over which Hr is computed

%  c = Type-3 LP filter coefficients

%  L = Order of Hr

%  h = Type-3 LP impulse response

%

 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';