Question

matlab做快速傅里叶变换，用FFT函数，程序看不懂，有高手能解释一下吗？

这是论文里的一段程序
Y = fft(ft,n);% 返回n点的离散傅里叶变换，如果向量ft的长度小于n，函数要将向量ft补零到长度n；如果向量ft的长度大于n，则函数ft使之长度为n。若ft 是矩阵，按相同方法对ft 进行处理。
halfLength = floor(n/2);
Pyy = Y(1:halfLength+1);%只选取前半截部分
Pyy = abs(Pyy1);%用于计算向量y的振幅
f =((0:halfLength)+1)*Fs1/n;
plot(f1,Pyy1),title('频谱'),grid;
我看不懂f为什么这样算，以及Pyy取前半截部分为什么Y(1:halfLength+1)这里要加1，不加1正好是一半吧。有人能给我解答一下吗？非常感谢！

• 你的回答被采纳后将获得：
• 系统奖励15（财富值+成长值）+难题奖励10（财富值+成长值）+提问者悬赏10（财富值+成长值）

Answer 1

Y(1:halfLength+1)中1:halfLength+1是索引，而索引必须是正整数，所以，不能从0开始，要都加1，但是指却是从0值开始到最后一个值。
至于f的算法是不一样的，在f =((0:halfLength)+1)*Fs1/n这里，0:halfLength)+1是数值，不是索引，结果是一个矢量。