ex060700.m

来自「介绍了一些数字信号处理方面的matlab程序」· M 代码 · 共 10 行

% 第六章：例6.7
% 频率采样形式: 给定频率样本
%
M = 32; alpha = (M-1)/2; k = 0:M-1; wk = (2*pi/M)*k;
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)