ex6p7.m

来自「关于数字信号处理的一些matlab例程」· M 代码 · 共 9 行

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)