ideal_hilbert_transformer_impulse_response.m

来自「basice code and example of MatLab about 」· M 代码 · 共 19 行

% dhilb.m - ideal lowpass Hilbert transformer FIR filter
%
% h = dhilb(wc, N) = row vector
%
% N = 2M+1 = filter length (odd)
% wc = cutoff frequency in [rads/sample]

function h = dhilb(wc, N)

M = (N-1)/2;

for k = -M:M,
    if k == 0,
        h(k+M+1) = 0;
    else
        h(k+M+1) = (1 - cos(wc * k)) / (pi * k);
    end
end