原创 matlab中针对txt数据文件的FFT变换

最近研究了matlab中的FFT变换，通过收集资料

整理了自己的频谱提取程序

a=textread('jb4.txt'); 

%读取幅值数据  幅值数据存在txt文本中

fs=5000;    采样频率

L=10000;  数据长度

NFFT =5000; 

% FFT点数，暂取5000，取一半，有资料说取越多越好

y = fft(a,NFFT)/L;

%进行fft变换（除以总采样点数，是为了后面精确看出原始信号幅值）

f = fs/2*linspace(0,1,NFFT/2+1);

%频率轴(只画到Fs/2即可，由于y为实数，后面一半是对称的)

plot(f,2*abs(y(1:NFFT/2+1)));