关于FFT频率确定的问题

twel2e

求fft高手为我解答！
我现在遇到的问题是：
& Z( g6 v3 r; t' R在matlab中我的输入信号为：
$ {9 q) `# K  o. N9 o) y, f! s% GS=Adc+A1*cos(2*pi*50*t+pi*P1/180)+A2*cos(2*pi*75*t+pi*P2/180);
取256点，通过FFT后，结果如图一所示：通常看图形的左半面，可以容易的看出50赫兹和75赫兹的信号
# k* r% {) S  v2 m7 o1 x而我把输入信号改为：
( T, s& Y5 a- Q% T$ Q8 YS=Adc+A1*cos(2*pi*50*t+pi*P1/180)+A2*cos(2*pi*150*t+pi*P2/180);
. R5 n4 U! p" k# M此时看图二，还是看左半面，我看到的是50赫兹和130赫兹的信号。。。真正的信号跑到右半面去了。。。。
如果我们的信号是未知的，那到底如何通过频谱图分析信号呢？
  ?/ e! ?4 A7 ?5 t' L5 e# Z7 }我就没有多少财富。。。。就50啦。。。希望高手看下。。。
我matlab的程序为：
close all; %先关闭所有图片
, w. e) }4 l0 g/ VAdc=2; %直流分量幅度
A1=3; %频率F1信号的幅度
A2=1.5; %频率F2信号的幅度
9 \* U5 [. k  T5 {% L; R. MF1=50; %信号1频率(Hz)
# \8 U/ z+ {& _3 i# L% d& `F2=75; %信号2频率(Hz)
2 v" b5 P3 c+ L" q5 t* SFs=256; %采样频率(Hz)
P1=-30; %信号1相位(度)
P2=90; %信号相位(度)
N=256; %采样点数
t=[0:1/Fs:N/Fs]; %采样时刻
' M. I9 y' C( z%信号
) P' B3 a  U* j, Y8 ~  @# g5 V8 AS=Adc+A1*cos(2*pi*F1*t+pi*P1/180)+A2*cos(2*pi*F2*t+pi*P2/180);
%显示原始信号
$ F4 P! b% T5 }+ Y  eplot(S);
title('原始信号');
7 W' q7 o4 D9 c& y4 M9 P( [! Lfigure;
Y = fft(S,N); %做FFT变换
Ayy = (abs(Y)); %取模
1 X: r4 B( A, }; a( R* A5 mplot(Ayy(1:N)); %显示原始的FFT模值结果
6 ^4 p; n9 z6 r! Ktitle('图1');
figure;
Ayy=Ayy/(N/2); %换算成实际的幅度
: g; h+ K5 Z  I& ?9 r! ]Ayy(1)=Ayy(1)/2;
. F- X* p& `: p: ?. e( ~F=([1:N]-1)*Fs/N; %换算成实际的频率值
plot(F(1:N/2),Ayy(1:N/2)); %显示换算后的FFT模值结果
title('幅度-频率曲线图');
5 |* v% H5 _7 e! P; h0 Qfigure;
Pyy=[1:N/2];
for i=1:N/2
Pyy(i)=phase(Y(i)); %计算相位
7 N6 _* n) e, \- m. }Pyy(i)=Pyy(i)*180/pi; %换算为角度
, _$ O$ |. e4 \5 h, b+ n  |0 i+ eend;
plot(F(1:N/2),Pyy(1:N/2)); %显示相位图
! X1 E. @) P1 M& J' btitle('相位-频率曲线图');

land

再看看别人是怎么说的

ESCAPE

采样率不对，会导致混叠

Memory00

256Hz的采样频率只能分析小于128Hz的信号，所以在采样以前应该设计一个低通滤波器，避免大于128Hz的信号混叠到带内。

芦根苏木

首先ADC采样要满足Nyquist采样定义，最简单就是FS>=2F0。你的FS是256Hz，而第二个仿真中出现了一个150Hz的信号分量，也就是不能满足采样定理。在仿真中的表现就是150Hz的信号混叠到采样带宽内了。
9 B7 e* s/ I+ R! K9 x