1．前言

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

　　2．傅里葉算法原理

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

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

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

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

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

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

　　3．軟件實(shí)現(xiàn)

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

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

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

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

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

基礎(chǔ)軟件開(kāi)發(fā)部  羊陽(yáng)