matlab寻找频谱峰值,如何得到明显峰值的频谱图？

%FFT变换,获得采样数据基本信息，时域图，频域图

%这里的向量都用行向量，假设被测变量是速度，单位为m/s

clear;

close all;

A=xlsread('1.xlsx'); %将测量数据赋给A，此时A为N×2的数组

x=A(:,1);                                     %将A中的第一列赋值给x，形成时间序列

x=x';                                           %将列向量变成行向量

y=A(:,2);                                     %将B中的第二列赋值给y，形成被测量序列

y=y';                                           %将列向量变成行向量

%显示原始数据曲线图(时域)

subplot(2,1,1);

plot(x,y)                                                                                %显示原始数据曲线图

axis([min(x) max(x) 1.1*floor(min(y)) 1.1*ceil(max(y))])              %优化坐标，可有可无

xlabel('时间 (s)');

ylabel('加速度(m/s^2)');

title('原始信号(时域)');

grid on;

%傅立叶变换

y=y-mean(y);