在MATLAB中用FFT把时域信号转化为频域信号

MATLAB实例

fe=12800;%采样频率
N=1280;%采样点数
dt=1/fe;%时间间隔
t=0:dt:(N-1)dt;
a=sin(2pi960t)+sin(2pi1200t)+sin(2pi1040t);%信号
PH2=(fft(a)); %a为时域信，作fft变换

%下面为幅值修正
P2 = (PH2/N);
P1 = P2(1:N/2+1);%此时选取前半部分，因为fft之后为对称的双边谱
P1(2:end-1) = 2*P1(2:end-1);

%画图
f = fe*(0:(N/2))/N;
plot(f,abs(P1))
title(‘时域信号’)
xlabel(‘f (Hz)’)
ylabel(’|P1(f)|’)

对比时域信号图和频域信号图：

可以实现由时域信号转换为频域信号，频域图中峰值所在频率与程序中信号a中设置的一致。