先进行傅里叶变换，将时域数据转化为频域数据，例如：

快速傅里叶变换
fftValue = abs(fftshift(fft(rawdata))); %快速傅里叶变换：将时域转化为频域

全部数据：

频谱图1：全部数据
figure % 频谱图1
set(gcf,'units','normalized','pos',[0 0 1 1]); % 最大化图形

Amp1 = (fftValue)/datalen;  % 纵坐标：实际幅值(每个幅度值除以datalen归一化)
f1 = (1:datalen)*Fs/datalen;  % 横坐标：生成对称频率坐标（数组索引必须是正整数或者逻辑值）

plot(f1,Amp1);    
xlim([0 8192]);%限制横坐标
ylim([0 0.006]);%限制横坐标
set(gca,'Fontsize',10);
title('频谱图1 全部');
xlabel('频率/HZ');ylabel('Amplitude');