離散樣本序列求諧波的前提是離散樣本序列為周期序列，最好是已知周期。

　　取一個周期的樣本，假設樣本為一維矩陣X。

　　在matlab的命令行輸入下述命令即可：

　　y=fft(x,N);    %對信號進行快速Fourier變換

　　mag=abs(y)     %求得Fourier變換后的振幅

　　矩陣mag的第一個數(shù)是0次諧波（直流分量），第二個數(shù)為基波（一次諧波幅值，第n+1個數(shù)就是n次諧波的幅值。

　　上述運算方法存在一定的缺陷，就是一個周期（或整數(shù)個周期）樣本的點數(shù)應該是2的N次冪，否則，F(xiàn)FT運算時，會在后面自動補零，影響頻譜的真實性，改善的方法是增加樣本的周期數(shù)，當樣本點數(shù)足夠多時，補零造成的影響可以忽略不計。

　　WP4000變頻功率分析儀采用DFT算法，犧牲了運算速度，但是，對運算點數(shù)沒有2的N次冪的限制，提高了基波測量及頻譜運算精度。

---