关于频谱混叠

fs=90;
T=1/10;
N=100*fs*T;
n=1:N;
t=n/fs;
s=cos(2*pi*10*t+pi/3)+2*cos(2*pi*50*t-pi/6);%T=1/10
F=fftshift(fft(s));
subplot(211);
plot(n*fs/N-fs/2,abs(F)/N);
xlabel('频率(Hz)');

subplot(212);
plot(n*2*pi/N-pi,abs(F)/N);
xlabel('角频率(rad)');
% 结论：
% 未发生频谱混叠时，频谱图中的角频率就是2*pi*f/fs；
% 当发生频谱混叠时，频谱图上表现的是折返后的频率，大小为：2*pi*(f/fs-1)
% 比如此例中fs=90，则折叠频率为-8*pi/9，约为2.7925
% 所以，当fs<2*fm时，折返后的频率为：fs-fm

运行结果：