4-22.m

来自「是《MATLAB数字信号处理与应用》一书的源代码,该书由李正周编著 清华大学出版」· M 代码 · 共 14 行

%例程4－22   Ⅰ、Ⅱ型理想带阻滤波器的单位冲激响应 计算

function hd=ideal_bs(Wcl,Wch,N)
%compute the ideal bandstop fiter unit pulse respondence hd(n)
%wcl: low cutoff frequency
%wch: high cutoff frequency
%N: window length
%hd: unit pulse respondence

alpha= (N-1)/2;
n=0:1:N-1;
m=n-alpha+eps;
hd=[sin(Wcl*m)+sin(pi*m)-sin(Wch*m)]./(pi*m);