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中设置的一致。
网友评论