1．前言

　　在微機保護裝置中，首先要對反映被保護設(shè)備的電氣量（電壓、電流）進行采集，然后對這些采集的數(shù)據(jù)進行數(shù)字濾波，再對這些經(jīng)過數(shù)字濾波的數(shù)字信號進行數(shù)學(xué)運算、邏輯運算，并進行分析判斷，最終輸出跳閘命令、信號命令或計算結(jié)果，以實現(xiàn)各種繼電保護功能。目前，實現(xiàn)模擬量計算幅值的算法有很多，而傅立葉算法在微機保護中作為計算信號幅值的算法被廣泛采用，不但可以求出各次諧波分量的幅值和相角，而且具有一定的濾波作用，可有效濾除恒定直流分量和各正次諧波分量。

　　2．傅里葉算法原理

　　那么，到底什么是傅里葉變換算法？傅里葉變換（Fourier transform）是一種線性的積分變換。因其基本思想首先由法國學(xué)者傅里葉系統(tǒng)提出，所以以其名字來命名以示紀(jì)念。

　　傅里葉變換一般分為連續(xù)傅里葉變換和離散傅里葉變換，電力系統(tǒng)中電氣量（電壓、電流）在進入微機保護裝置后，軟件通過AD對模擬量進行固定采樣頻率的等間隔采樣，會得到一組周期性的離散信號，我們以此信號來計算出相應(yīng)幅值和相角。因此，這里應(yīng)用的是離散傅里葉變換，其變換公式如下：

　　傅氏算法的基本思路來自傅里葉級數(shù)，一個周期函數(shù)滿足狄里赫利條件，就可以將這個周期函數(shù)分解為一個級數(shù)，即一個周期性函數(shù)可以分解為直流分量、基波分量及各次諧波的無窮級數(shù)，如下表示：

　　據(jù)傅氏級數(shù)的原理，可以求出 、 分別如下：

　　上述公式中 表示基波角頻率； 和 分別是各次諧波的正弦和余弦的幅值，其中比較特殊的有： 表示直流分量， 表示基波分量正、余弦項的幅值。根據(jù)上述公式計算各次諧波時，是通過對余弦各值進行周期性的累加，所以就消除了直流分量 的影響。

　　已知獲得的各次諧波的正弦和余弦的幅值 和 ，即可得到諧波有效的幅值及應(yīng)用相角θn。

　　3．軟件實現(xiàn)

　　軟件通過AD對模擬量進行固定采樣頻率的等間隔采樣，設(shè)每周波采樣點為N，即w1=2π/N，采樣次數(shù)為k。 和 可變換為：

　　這里我們來計算信號的基波分量，首先求出基波分量（n=1）的實部和虛部 。以我公司最新繼電保護裝置為例，采用32位高性能SoC嵌入式微機處理器，大容量的RAM和Flash Memory，高精度AD，每周波采樣點數(shù)N=128，則每周波采樣值分別是i0、i1、i2……i127，同時其對應(yīng)的基波正弦和余弦系數(shù)如下表所示。

　　基波正弦和余弦系數(shù)（N=128）

　　根據(jù)上述公式，進行基波正、余弦的系數(shù)累計相加，于是可得：

　　由此即可進一步計算出相應(yīng)的幅值和相位。

基礎(chǔ)軟件開發(fā)部  羊陽