sfreq1.m

来自「里面囊括了基于matlab滤波器设计的各种.m文件」· M 代码 · 共 17 行

%f = sfreq1(h,omega)
%Computation of the frequency response of filter h
%at the point omega by using the simple method
%but the order of h must be even
%See also SFREQ

function f=sfreq1(h,omega)
N = size(h,2);
for I = 0:(N-1)/2
 c(I+1) = cos(I*omega);
 if I==0
  s(I+1) = h((N+1)/2);
 else
  s(I+1) = 2*h((N+1)/2-I); 
 end
end
f = abs(s*c');