分别计算16点序列x(n)=cos 5π/16 n,0≤n≤15的16点和32点fft,绘出幅度谱图形，并绘出该序列的DTFT幅频。matlab程序

以下是MATLAB程序实现：

% 16点FFT
x = cos(5*pi/16*(0:15));
X = fft(x,16);
subplot(2,2,1);
stem(abs(X));
title('16-point FFT');

% 32点FFT
X2 = fft(x,32);
subplot(2,2,2);
stem(abs(X2));
title('32-point FFT');

% DTFT幅频
w = linspace(0,2*pi,1024);
H = zeros(size(w));
for i = 1:length(w)
    H(i) = sum(x.*exp(-1j*w(i).*(0:length(x)-1)));
end
subplot(2,2,3:4);
plot(w,abs(H));
title('DTFT magnitude');

该程序将生成一个包含四个子图的图形。第一个子图显示16点FFT的幅度谱，第二个子图显示32点FFT的幅度谱，第三个和第四个子图显示x(n)的DTFT幅度谱。