关于FFT频率确定的问题

twel2e

求fft高手为我解答！
& [8 T* v9 g. ]8 r我现在遇到的问题是：
1 ~& D% z2 ]; \2 C! n, S+ F在matlab中我的输入信号为：
S=Adc+A1*cos(2*pi*50*t+pi*P1/180)+A2*cos(2*pi*75*t+pi*P2/180);
! f: ^: R" d% u  J取256点，通过FFT后，结果如图一所示：通常看图形的左半面，可以容易的看出50赫兹和75赫兹的信号
而我把输入信号改为：
* Z; l0 g) N- N1 K, ]S=Adc+A1*cos(2*pi*50*t+pi*P1/180)+A2*cos(2*pi*150*t+pi*P2/180);
: i( T: o  O! H8 H此时看图二，还是看左半面，我看到的是50赫兹和130赫兹的信号。。。真正的信号跑到右半面去了。。。。
如果我们的信号是未知的，那到底如何通过频谱图分析信号呢？
我就没有多少财富。。。。就50啦。。。希望高手看下。。。
我matlab的程序为：
close all; %先关闭所有图片
9 ~  U8 _. z' ]1 E/ v0 ZAdc=2; %直流分量幅度
A1=3; %频率F1信号的幅度
A2=1.5; %频率F2信号的幅度
7 A4 Z' O3 _7 L, y8 hF1=50; %信号1频率(Hz)
F2=75; %信号2频率(Hz)
Fs=256; %采样频率(Hz)
P1=-30; %信号1相位(度)
$ ?, a- m& M1 [) Y$ W0 kP2=90; %信号相位(度)
N=256; %采样点数
2 }" p  @8 Q8 b( P7 bt=[0:1/Fs:N/Fs]; %采样时刻
" b0 F! w" O; l5 E4 ?! N( r2 ]/ T%信号
7 G$ F' G9 @3 S. K3 nS=Adc+A1*cos(2*pi*F1*t+pi*P1/180)+A2*cos(2*pi*F2*t+pi*P2/180);
( A* F8 g, W" t# o! N; w, ~$ l%显示原始信号
plot(S);
title('原始信号');
figure;
Y = fft(S,N); %做FFT变换
Ayy = (abs(Y)); %取模
plot(Ayy(1:N)); %显示原始的FFT模值结果
8 W8 h3 B+ M. e$ u# G3 {" @title('图1');
figure;
! y4 ]1 F. z5 h( H9 x' y4 K2 D2 N# l; I; KAyy=Ayy/(N/2); %换算成实际的幅度
Ayy(1)=Ayy(1)/2;
F=([1:N]-1)*Fs/N; %换算成实际的频率值
plot(F(1:N/2),Ayy(1:N/2)); %显示换算后的FFT模值结果
0 B6 R9 P2 Z) e3 s; _title('幅度-频率曲线图');
figure;
! _7 D2 e2 D3 p6 w# UPyy=[1:N/2];
6 x$ \( u- f/ P, J% Kfor i=1:N/2
Pyy(i)=phase(Y(i)); %计算相位
# u7 W# W) M5 i- a! N9 o) nPyy(i)=Pyy(i)*180/pi; %换算为角度
( ~2 X2 S* Z2 Oend;
plot(F(1:N/2),Pyy(1:N/2)); %显示相位图
title('相位-频率曲线图');
! N/ z4 Y) Z: {! i7 z3 f4 V

芦根苏木

首先ADC采样要满足Nyquist采样定义，最简单就是FS>=2F0。你的FS是256Hz，而第二个仿真中出现了一个150Hz的信号分量，也就是不能满足采样定理。在仿真中的表现就是150Hz的信号混叠到采样带宽内了。
& \- A6 ?& s2 v( }

Memory00

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

ESCAPE

采样率不对，会导致混叠

land

再看看别人是怎么说的