matlab fft变换求出频谱后，求具体点的幅值的问题

程序如下
fs=250;
N=1024;
n=0:N-1;
t=n/fs;
y=fft(x(:,1),N);
mag=abs(y);
f=n*fs/N;
plot(f,mag);
set(gca,'xtick',0:10:250)
由原始常微分方程的结果进行fft变换，求得频谱如下图：

我想求得频谱图中20,40,60,80,100,120hz处的幅值，然后求他们的比值，我该如何编写程序呢？望大神指教。

你的mag中保存的就是幅值信息，你可以算一下频率对应的点，然后取该带点的值
比如 mag(20)应该就是，如果不是，看看是不是20-1 或20+1的点

fft结果是关于采样频率一半共轭对称的 入如考虑幅值，就是关于采样频率一半对称的，只能给出0频到采样频率一半的频谱信息 因此fft结果通常只取前一半，而真实信号的幅度被平均分到前后各一半 所以要乘以2 而fft是积分变换，随着采样点数的增加，信号的总量会增加 例如同一个信号，如果由n个采样点变成2n个采样点 整个信号的积分就会成倍的增加，所以结果要除以n以抵销采样数对结果的影响 所以最后要乘以2，除以n，也就是除以 (n/2)