（matlab）地震数据频谱分析-频谱图代码

figure
plot(DATA)
title(‘原始地震道’)
DATA_fft_temp = abs(fft(DATA));
fs = 1/dt;
sample_point = length(DATA_fft_temp);
f_x = 0:fs/(sample_point-1):fs; % 频率数组
DATA_fft = DATA_fft_temp(1:end/2);
figure
plot(f_x(1:end/2),DATA_fft,‘k’) % 展示频谱图
grid on
title(‘频谱图’)