你好 和我做的以前做的实验一样 代码如下:
n=0:31;Fs=64;T=1/Fs;
f1=10;xn1=cos(2*pi*f1*n*T); % 由f=10Hz的余弦序列Fs=64Hz采样得序列xn1
f2=11;xn2=cos(2*pi*f2*n*T); % 由f=11Hz的余弦序列Fs=64Hz采样得序列xn2
yn1=fft(xn1,32); % 序列的DFT计算
yn2=fft(xn2,32);
Fs=32;T=1/fs;
xn3=cos(2*pi*f2*n*T); % 由f=11Hz的余弦序列Fs=32Hz采样得序列xn2
yn3=fft(xn3,32)
subplot(311);stem(n,abs(yn1),'b','.'); % 绘图部分
xlabel('');ylabel('|X(k)|');
title('f=10Hz,Fs=64Hz,N=32时余弦信号频谱');
subplot(312);stem(n,abs(yn2),'r','*');
xlabel('k');ylabel('|X(k)|');
title('f=11Hz,Fs=64Hz,N=32时余弦信号频谱');
subplot(313);stem(abs(yn3),'g','*');
xlabel('k');ylabel('|X(k)|');
title('f=11Hz,Fs=32Hz,N=32时余弦信号频谱');