《DSPusingMATLAB》示例Example6.7

代码&＃xff1a;

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

 运行结果&＃xff1a;