我知道规范会给我该分量的频率和相应的幅度，而 fft 将计算信号的 DFT 并抛出我必须根据物理频率绘制的每个分量的复数。

y 是时间序列信号 Fs 是采样频率。

N = length(y);
f1=fft(y)
fa = seq(0,N-1,by=1)*(Fs/N)
plot(fa,abs(f1),type='l',xlim = c(0,3),ylim = c(0,50))
s1=spec(y,f=1000,plot = FALSE)
s1$x<-s1$x*1000   (KHz to Hz)
plot(s1$x,s1$y,type='l',xlim = c(0,3))

红色一个来自 spec 函数，蓝色来自 fft 函数。