exa3_4.m

来自「matlab应用实例。如：FIR、IIR数字滤波器设计等。」· M 代码 · 共 14 行

%  exa3-4_dir2fs.m, for example 3-4
% to test dir2fs.m and to obtain the step response.

clear;
M=32;
alpha=(M-1)/2;
magHk=[1,1,1,0.5,zeros(1,25),0.5,1,1];
k1=0:15;
k2=16:M-1;
angHk=[-alpha*(2*pi)/M*k1,alpha*(2*pi)/M*(M-k2)];
H=magHk.*exp(j*angHk);
h=real(ifft(H,M));
[C,B,A]=dir2fs(h)