CN103439567B - 一种moa阻性电流在线监测系统的工作方法 - Google Patents

Abstract

本发明公开了一种MOA阻性电流在线监测系统的工作方法，其阻性泄漏电流基波参数通常按照投影法来获得，具体过程为：同时采样电网电压和MOA泄漏电流信号；应用可变栅栏准同步DFT谐波分析技术获取电网电压的基波幅值和初相角应用可变栅栏准同步DFT谐波分析技术获取MOA泄漏电流的基波幅值和初相角按照投影法获取泄漏电流基波在电网电压基波上的投影角获取MOA阻性泄漏电流基波

Description

一种MOA阻性电流在线监测系统的工作方法

本申请是分案申请，原申请的申请号：201110322641.5，发明创造名称：一种MOA阻性电流在线监测系统的工作方法，申请日：2011-10-21。

技术领域

本发明涉及MOA阻性电流测试的技术领域，具体是一种基于高精度的谐波分析方法的MOA阻性电流在线监测系统的工作方法。

背景技术

金属氧化物避雷器(以下简称MOA)因其优越的过电压保护特性在电力系统中得到广泛应用，但MOA电阻片老化以及经受热和冲击破坏会引起故障，严重可能会导致其爆炸，避雷器击穿还会导致变电站母线短路，影响系统安全运行，因此必须对运行中的MOA进行严格有效的检测和定期预防性试验。在氧化锌避雷器的检测和试验中，交流运行电压下的泄漏电流测量是一个重要的项目，泄漏电流中的阻性电流基波成分的大小能较准确的反映氧化锌避雷器受潮、阀片的老化和内部绝缘受损等缺陷。因此，对阻性泄漏电流基波等参量进行准确测量就可以较为准确地对MOA性能进行判别。

谐波分析技术在电能质量监控、电子产品生产检验、电器设备监控等众多领域应用广泛，是进行电网监控、质量检验、设备监控的重要技术手段。目前谐波分析应用最广泛的技术是离散傅里叶变换(DFT)和快速傅里叶变换(FFT)。准同步采样技术和DFT技术相结合的谐波分析技术(即：准同步DFT谐波分析技术)能够提高谐波分析的精度，其算式为：

$\left\{\begin{array}{c}{a}_k={2F}_{\mathrm{ak}}^n=\frac{2}{Q}\underset{i=0}{\overset{W}{\mathrm{\Σ}}}{\mathrm{\γ}}_{i}f\left(i\right)\cos \left(k\frac{2\mathrm{\π}}{N}i\right)\\ b_k={2F}_{\mathrm{bk}}^n=\frac{2}{Q}\underset{i=0}{\overset{W}{\mathrm{\Σ}}}{\mathrm{\γ}}_{i}f\left(i\right)\sin \left(k\frac{2\mathrm{\π}}{N}i\right)\end{array}\right.,$

式中：k为需要获得的谐波的次数(如基波k=1，3次谐波k=3)；sin和cos分别为正弦和余弦函数；而a_k和b_k分别为k次谐波的实部和虚部；n为迭代次数；W由积分方法决定，采用复化梯形积分方法时，W=nN；γ_i为一次加权系数；为所有加权系数之和；f(i)为分析波形的第i个采样值；N为周期内采样次数。

在工程应用中，谐波分析总是进行有限点的采样和难以做到严格意义的同步采样。这样，在应用准同步DFT进行谐波分析时，就会存在由于截断效应导致的长范围泄漏和由于栅栏效应导致的短范围泄漏，使得分析结果精度不高，甚至不可信。

发明内容

本发明要解决的技术问题是提供一种精度较高的MOA阻性电流在线监测系统的工作方法，以有效改进准同步DFT谐波分析技术的分析误差，获得高精度的谐波分析结果，从而提高基于谐波分析理论的MOA阻性电流测试的可靠性。

为解决上述技术问题，本发明提供的MOA阻性电流在线监测系统包括：MOA阻性电流在线测试单元、与该MOA阻性电流在线测试单元通过有线或无线网络进行数据通讯的监测服务器；所述MOA阻性电流在线监测系统的工作方法包括以下步骤：

(1)、所述MOA阻性电流在线测试单元的CPU同时采样电网电压V和MOA的泄漏电流I，并分别对电网电压V和泄漏电流I等间隔采样W+2个采样点数据：{f(i)，i=0，1，…，w+1}；

(2)、所述CPU从所述电网电压V的采样点i=0开始应用准同步DFT公式：

$\left\{\begin{array}{c}{a}_k={2F}_{\mathrm{ak}}^n=\frac{2}{Q}\underset{i=0}{\overset{W}{\mathrm{\Σ}}}{\mathrm{\γ}}_{i}f\left(i\right)\cos \left(k\frac{2\mathrm{\π}}{N}i\right)\\ b_k={2F}_{\mathrm{bk}}^n=\frac{2}{Q}\underset{i=0}{\overset{W}{\mathrm{\Σ}}}{\mathrm{\γ}}_{i}f\left(i\right)\sin \left(k\frac{2\mathrm{\π}}{N}i\right)\end{array}\right.,$ 分析W+1个数据获得基波信息和

然后，所述CPU从所述电网电压V的采样点i=1开始应用准同步DFT公式：

$\left\{\begin{array}{c}{a}_k={2F}_{\mathrm{ak}}^n=\frac{2}{Q}\underset{i=1}{\overset{W+1}{\mathrm{\Σ}}}{\mathrm{\γ}}_{i}f\left(i\right)\cos \left(k\frac{2\mathrm{\π}}{N}(i-1)\right)\\ b_k={2F}_{\mathrm{bk}}^n=\frac{2}{Q}\underset{i=1}{\overset{W+1}{\mathrm{\Σ}}}{\mathrm{\γ}}_{i}f\left(i\right)\sin \left(k\frac{2\mathrm{\π}}{N}(i-1)\right)\end{array}\right.,$ 分析W+1个数据获得基波信息和

应用公式：计算所述电网电压V的频率漂移μ；

应用公式： $\left\{\begin{array}{c}{a}_k=\frac{2}{Q}\underset{i=0}{\overset{W}{\mathrm{\Σ}}}{\mathrm{\γ}}_{i}f\left(i\right)\cos \left(k\frac{\mathrm{\μ}2\mathrm{\π}}{N}i\right)\\ b_k=\frac{2}{Q}\underset{i=0}{\overset{W}{\mathrm{\Σ}}}{\mathrm{\γ}}_{i}f\left(i\right)\sin \left(k\frac{\mathrm{\μ}2\mathrm{\π}}{N}i\right)\end{array}\right.,$ 得出所述电网电压V的基波的实部a_k和虚部b_k，然后计算出所述电网电压V的基波的幅值和初相角

(3)、所述CPU从所述泄漏电流I的采样点i=0开始应用准同步DFT公式：

$\left\{\begin{array}{c}{a}_k=2F_{\mathrm{ak}}^n=\frac{2}{Q}\underset{i=0}{\overset{W}{\mathrm{\Σ}}}{\mathrm{\γ}}_{i}f\left(i\right)\cos \left(k\frac{2\mathrm{\π}}{N}i\right)\\ b_k=2F_{\mathrm{bk}}^n=\frac{2}{Q}\underset{i=0}{\overset{W}{\mathrm{\Σ}}}{\mathrm{\γ}}_{i}f\left(i\right)\sin \left(k\frac{2\mathrm{\π}}{N}i\right)\end{array}\right.,$ 分析W+1个数据获得基波信息和

然后，所述CPU从泄漏电流I的采样点i=1应用准同步DFT公式：

$\left\{\begin{array}{c}{a}_k=2F_{\mathrm{ak}}^n=\frac{2}{Q}\underset{i=1}{\overset{W+1}{\mathrm{\Σ}}}{\mathrm{\γ}}_{i}f\left(i\right)\cos \left(k\frac{2\mathrm{\π}}{N}(i-1)\right)\\ b_k=2F_{\mathrm{bk}}^n=\frac{2}{Q}\underset{i=1}{\overset{W+1}{\mathrm{\Σ}}}{\mathrm{\γ}}_{i}f\left(i\right)\sin \left(k\frac{2\mathrm{\π}}{N}(i-1)\right)\end{array}\right.,$

分析W+1个数据获得基波信息和

应用公式：计算所述泄漏电流I的频率漂移μ；

应用公式： $\left\{\begin{array}{c}{a}_k=\frac{2}{Q}\underset{i=0}{\overset{W}{\mathrm{\Σ}}}{\mathrm{\γ}}_{i}f\left(i\right)\cos \left(k\frac{\mathrm{\μ}2\mathrm{\π}}{N}i\right)\\ b_k=\frac{2}{Q}\underset{i=0}{\overset{W}{\mathrm{\Σ}}}{\mathrm{\γ}}_{i}f\left(i\right)\sin \left(k\frac{\mathrm{\μ}2\mathrm{\π}}{N}i\right)\end{array}\right.$

，分别得出所述泄漏电流I的基波的实部a_k和虚部b_k，然后计算出所述泄漏电流I的基波的幅值和初相角

(4)、所述CPU按照投影法获取泄漏电流I的基波在电网电压V的基波上的投影角

(5)、所述CPU获取MOA阻性泄漏电流的基波并输出显示，同时发送至所述监测服务器。

所述MOA阻性电流在线监测系统的MOA阻性电流在线测试单元的工作方法中，采用了一种谐波分析方法，其应用准同步DFT进行谐波分析时频域抽样的位置根据信号频率的漂移而改变，即所述频域抽样位置为μ2π／N，其中：μ为信号频率的漂移，无漂移时μ为1。

本发明的所述谐波分析方法基于可变栅栏的思想，是通过5个分析步骤实现的。

可变栅栏的思想：准同步DFT分析误差的主要原因是信号频率的漂移导致频谱峰值出现的位置与理想位置发生偏差，如果仍然按照2π／N在频域中以进行抽样的话得到的分析结果极不正确。可变栅栏指的是：频域抽样的位置的并不是固定的2π／N，而是根据信号频率的漂移而改变，即频率抽样位置为μ2π／N(μ为信号频率的漂移)。频域抽样栅栏随着信号频率的漂移而改变可以准确估计出各次谐波峰值出现的位置，进而获取高精度的幅值和相角信息。

准同步DFT谐波分析可以有效地抑制长范围泄漏，其频谱泄漏的主要原因是信号频率漂移导致的短范围泄漏，而信号频率漂移导致的短范围泄漏的主要特征是谱峰峰值出现位置随着信号频率漂移而同步改变，所以可变栅栏频域采样能够有效根据信号漂移捕捉谱峰峰值出现的位置，从而获得高精度的谐波信息。

等间隔采样是根据进行谐波分析的理想信号的周期T和频率f(如工频信号频率f为50Hz，周期为20mS)，在一个周期内采样N点，即采样频率为f_s=Nf，且N≥64。

所述的采样W+2个采样点数据是根据所选择的积分方法而作相应选择，若采用复化梯形积分方法，则W=nN；若采用复化矩形积分方法，则W=n(N-1)；若采用复化辛普森积分方法，则W=n(N-1)／2。然后根据采样频率f_s=Nf，获得采样点数据序列{f(i)，i=0，1，…，w+1}，n≥3，最后对该数据序列进行谐波分析。

一次迭代系数γ_i由积分方法、理想周期采样点N和迭代次数n决定，具体推导过程参见文献【戴先中.准同步采样应用中的若干问题[J].电测与仪表，1988，(2)：2-7.】。

a_k和b_k为k次谐波的虚部和实部，根据a_k和b_k就可以获得谐波幅值和初相角。

信号频率的漂移μ是根据相邻采样点基波相角差与理想周期内采样点数N的固定关系而获得的，信号频率的漂移μ也可用于修正基波的频率f₁和高次谐波的频率f_k。

本发明的MOA阻性电流在线测试单元应用可变栅栏准同步DFT谐波分析技术对运行中的MOA的阻性泄漏电流基波等参量进行检测，进而判断MOA的性能。MOA阻性电流在线测试单元由MOA泄漏泄漏电流传感器、电网电压传感器、与该MOA泄漏电流传感器和电网电压传感器相连的信号调理电路、与该信号调理电路相连的数据采集电路、与该数

据采集电路相连的CPU、与该CPU相连的有线或无线网络通讯单元、与该CPU相连的LCD单元和相应的分析软件组成。所述无线网络包括GSM或CDMA或3G无线通讯网络；MOA泄漏电流传感器一般采用电流互感器，串装在MOA的下端的接地线中，获取MOA的泄漏电流；电网电压传感器一般采用电压互感器，用于获取电网的电压信号；信号调理电路将传感器送来的电压信号、泄漏电流信号变换为适合采集电路采样的电压信号；由数据采集电路将调理后的信号进行采样，再由CPU通过分析软件来完成对谐波的分析与检测。

本发明的MOA阻性电流在线测试单元的工作方法，用于改进和提高谐波分析的质量，获取高精度的电网电压和MOA泄漏电流的基波幅值、相角。

本发明的MOA阻性电流在线测试单元采用的谐波分析方法，为基于可变栅栏思想的谐波分析技术，具有以下技术优势：

(1)高精度的谐波分析结果。本发明所述的谐波分析技术获得的分析结果无论是幅值还是相角误差提高4个数量级以上。

(2)本发明所述的谐波分析技术从根本上解决了准同步DFT分析精度低的问题，而无需进行复杂的反演和修正，算法简单。

(3)相对于准同步DFT，本发明所述的谐波分析技术只需要增加一个采样点就解决了准同步DFT分析误差大的问题，易于实现。

(4)应用本发明来改进现有的仪器设备，技术上是可行，并且不需要增加任何的硬件开销就可使分析结果可以提高4个数量级以上。

(5)可变栅栏思想也同样也适用于进行多次迭代而非一次迭代的谐波分析过程，此时只需要把一次迭代分解成多次迭代实现就可以了。一次迭代和多次迭代本质上是一样的，只是在计算时多次迭代进行分步计算，而一次迭代是把多次迭代的过程合并到迭代系数γ_i中一次计算完成，所以本发明同样适用于多次迭代过程。

具体实施方式

本发明的一种MOA阻性电流在线监测系统包括：MOA阻性电流在线测试单元、与该MOA阻性电流在线测试单元通过有线或无线网络进行数据通讯的监测服务器。

所述MOA阻性电流在线监测系统的工作方法，包括以下步骤：

(1)、所述MOA阻性电流在线测试单元的CPU同时采样电网电压V和MOA的泄漏电流I，并分别对电网电压V和泄漏电流I等间隔采样W+2个采样点，以获得被分析信号的离散序列{f(k)，k＝0，1，…，w+1}。W由积分方法、迭代次数n和理想周期内采样点数N共同决定。等间隔采样指的是根据进行谐波分析的理想信号的频率f(如工频信号频率为50Hz，周期为20mS)确定采样频率f_s＝Nf，在采样频率f_s的作用下在一个周期内均匀地采样N点。一般地，周期采样点N＝64或以上就能获得较好的谐波分析结果，而迭代次数n＝3-5就能获得较理想的谐波分析结果。积分方法有复化梯形积分方法W＝nN、复化矩形积分方法W＝n(N-1)、辛普森积分方法W＝n(N-1)/2等多种，可以根据实际情况进行选择。

(2)、所述CPU从电网电压V的采样点i＝0开始应用准同步DFT公式

$\left\{\begin{array}{c}{a}_k={2F}_{\mathrm{ak}}^n=\frac{2}{Q}\underset{i=0}{\overset{W}{\mathrm{\Σ}}}{\mathrm{\γ}}_{i}f\left(i\right)\cos \left(k\frac{2\mathrm{\π}}{N}i\right)\\ b_k={2F}_{\mathrm{bk}}^n=\frac{2}{Q}\underset{i=0}{\overset{W}{\mathrm{\Σ}}}{\mathrm{\γ}}_{i}f\left(i\right)\sin \left(k\frac{2\mathrm{\π}}{N}i\right)\end{array}\right.,$ 分析W+1个数据获得基波信息和

其中，一次迭代系数γ_i由积分方法、理想周期采样点N和迭代次数n决定，而为所有加权系数之和。

所述CPU从所述电网电压V的采样点i＝1开始应用准同步DFT公式：

$\left\{\begin{array}{c}{a}_k={2F}_{\mathrm{ak}}^n=\frac{2}{Q}\underset{i=1}{\overset{W+1}{\mathrm{\Σ}}}{\mathrm{\γ}}_{i}f\left(i\right)\cos \left(k\frac{2\mathrm{\π}}{N}(i-1)\right)\\ b_k={2F}_{\mathrm{bk}}^n=\frac{2}{Q}\underset{i=1}{\overset{W+1}{\mathrm{\Σ}}}{\mathrm{\γ}}_{i}f\left(i\right)\sin \left(k\frac{2\mathrm{\π}}{N}(i-1)\right)\end{array}\right.,$ 分析W+1个数据获得基波信息和

应用公式：计算所述电网电压V的频率漂移μ：

应用公式： $\left\{\begin{array}{c}{a}_k=\frac{2}{Q}\underset{i=0}{\overset{W}{\mathrm{\Σ}}}{\mathrm{\γ}}_{i}f\left(i\right)\cos \left(k\frac{\mathrm{\μ}2\mathrm{\π}}{N}i\right)\\ b_k=\frac{2}{Q}\underset{i=0}{\overset{W}{\mathrm{\Σ}}}{\mathrm{\γ}}_{i}f\left(i\right)\sin \left(k\frac{\mathrm{\μ}2\mathrm{\π}}{N}i\right)\end{array}\right.,$ 得出所述电网电压V的基波的实部a_k和虚部b_k，然后计算出所述电网电压V的基波的幅值和初相角

(3)、所述CPU从所述泄漏电流I的采样点i＝0开始应用准同步DFT公式：

$\left\{\begin{array}{c}{a}_k={2F}_{\mathrm{ak}}^n=\frac{2}{Q}\underset{i=0}{\overset{W}{\mathrm{\Σ}}}{\mathrm{\γ}}_{i}f\left(i\right)\cos \left(k\frac{2\mathrm{\π}}{N}i\right)\\ b_k={2F}_{\mathrm{bk}}^n=\frac{2}{Q}\underset{i=0}{\overset{W}{\mathrm{\Σ}}}{\mathrm{\γ}}_{i}f\left(i\right)\sin \left(k\frac{2\mathrm{\π}}{N}i\right)\end{array}\right.,$ 分析W+1个数据获得基波信息和

所述CPU从泄漏电流I的采样点i=1应用准同步DFT公式：

$\left\{\begin{array}{c}{a}_k=2F_{\mathrm{ak}}^n=\frac{2}{Q}\underset{i=1}{\overset{W+1}{\mathrm{\Σ}}}{\mathrm{\γ}}_{i}f\left(i\right)\cos \left(k\frac{2\mathrm{\π}}{N}(i-1)\right)\\ b_k=2F_{\mathrm{bk}}^n=\frac{2}{Q}\underset{i=1}{\overset{W+1}{\mathrm{\Σ}}}{\mathrm{\γ}}_{i}f\left(i\right)\sin \left(k\frac{2\mathrm{\π}}{N}(i-1)\right)\end{array}\right.,$

分析W+1个数据获得基波信息和

应用公式：计算所述泄漏电流I的频率漂移μ；

应用公式： $\left\{\begin{array}{c}{a}_k=\frac{2}{Q}\underset{i=0}{\overset{W}{\mathrm{\Σ}}}{\mathrm{\γ}}_{i}f\left(i\right)\cos \left(k\frac{\mathrm{\μ}2\mathrm{\π}}{N}i\right)\\ b_k=\frac{2}{Q}\underset{i=0}{\overset{W}{\mathrm{\Σ}}}{\mathrm{\γ}}_{i}f\left(i\right)\sin \left(k\frac{\mathrm{\μ}2\mathrm{\π}}{N}i\right)\end{array}\right.,$

，分别得出所述泄漏电流I的基波的实部a_k和虚部b_k，然后计算出所述泄漏电流I的基波的幅值和初相角

(4)、所述CPU按照投影法获取泄漏电流I的基波在电网电压V的基波上的投影角

(5)、所述CPU获取MOA阻性泄漏电流的基波并输出显示，同时发送至所述监测服务器。

本技术领域的普通技术人员应当认识到，以上的实施例仅是用来说明本发明，而并非作为对本发明的限定，本发明还可以变化成更多的方式，只要在本发明的实质精神范围内，对以上所述实施例的变化、变型都将落在本发明的权利要求书范围内。

Claims (1)

1.一种MOA阻性电流在线监测系统的工作方法，所述MOA阻性电流在线监测系统包括：MOA阻性电流在线测试单元、与该MOA阻性电流在线测试单元通过有线或无线网络进行数据通讯的监测服务器；其特征在于包括以下步骤：

(1)、所述MOA阻性电流在线测试单元的CPU同时采样电网电压V和MOA的泄漏电流I，并分别对电网电压V和泄漏电流I等间隔采样W+2个采样点数据：{f(i)，i=0，1，…，w+1}；

(2)、所述CPU从所述电网电压V的采样点i=0开始应用准同步DFT公式：

$\left\{\begin{array}{c}{a}_k={2F}_{\mathrm{ak}}^n=\frac{2}{Q}\underset{i=0}{\overset{W}{\mathrm{\Σ}}}{\mathrm{\γ}}_{i}f\left(i\right)\cos \left(k\frac{2\mathrm{\π}}{N}i\right)\\ b_k={2F}_{\mathrm{bk}}^n=\frac{2}{Q}\underset{i=0}{\overset{W}{\mathrm{\Σ}}}{\mathrm{\γ}}_{i}f\left(i\right)\sin \left(k\frac{2\mathrm{\π}}{N}i\right)\end{array}\right.,$ 分析W+1个数据获得基波信息和

然后，所述CPU从所述电网电压V的采样点i=1开始应用准同步DFT公式：

$\left\{\begin{array}{c}{a}_k={2F}_{\mathrm{ak}}^n=\frac{2}{Q}\underset{i=1}{\overset{W+1}{\mathrm{\Σ}}}{\mathrm{\γ}}_{i}f\left(i\right)\cos \left(k\frac{2\mathrm{\π}}{N}(i-1)\right)\\ b_k={2F}_{\mathrm{bk}}^n=\frac{2}{Q}\underset{i=1}{\overset{W+1}{\mathrm{\Σ}}}{\mathrm{\γ}}_{i}f\left(i\right)\sin \left(k\frac{2\mathrm{\π}}{N}(i-1)\right)\end{array}\right.,$ 分析W+1个数据获得基波信息和

应用公式：计算所述电网电压V的频率漂移μ；

应用公式： $\left\{\begin{array}{c}{a}_k=\frac{2}{Q}\underset{i=0}{\overset{W}{\mathrm{\Σ}}}{\mathrm{\γ}}_{i}f\left(i\right)\cos \left(k\frac{\mathrm{\μ}2\mathrm{\π}}{N}i\right)\\ b_k=\frac{2}{Q}\underset{i=0}{\overset{W}{\mathrm{\Σ}}}{\mathrm{\γ}}_{i}f\left(i\right)\sin \left(k\frac{\mathrm{\μ}2\mathrm{\π}}{N}i\right)\end{array}\right.,$ 得出所述电网电压V的基波的实部a_k和虚部b_k，然后计算出所述电网电压V的基波的幅值和初相角

(3)、所述CPU从所述泄漏电流I的采样点i＝0开始应用准同步DFT公式：

$\left\{\begin{array}{c}{a}_k={2F}_{\mathrm{ak}}^n=\frac{2}{Q}\underset{i=0}{\overset{W}{\mathrm{\Σ}}}{\mathrm{\γ}}_{i}f\left(i\right)\cos \left(k\frac{2\mathrm{\π}}{N}i\right)\\ b_k={2F}_{\mathrm{bk}}^n=\frac{2}{Q}\underset{i=0}{\overset{W}{\mathrm{\Σ}}}{\mathrm{\γ}}_{i}f\left(i\right)\sin \left(k\frac{2\mathrm{\π}}{N}i\right)\end{array}\right.,$ 分析W+1个数据获得基波信息和

然后，所述CPU从泄漏电流I的采样点i=1应用准同步DFT公式：

$\left\{\begin{array}{c}{a}_k={2F}_{\mathrm{ak}}^n=\frac{2}{Q}\underset{i=1}{\overset{W+1}{\mathrm{\Σ}}}{\mathrm{\γ}}_{i}f\left(i\right)\cos \left(k\frac{2\mathrm{\π}}{N}(i-1)\right)\\ b_k={2F}_{\mathrm{bk}}^n=\frac{2}{Q}\underset{i=1}{\overset{W+1}{\mathrm{\Σ}}}{\mathrm{\γ}}_{i}f\left(i\right)\sin \left(k\frac{2\mathrm{\π}}{N}(i-1)\right)\end{array}\right.,$

分析W+1个数据获得基波信息和

应用公式：计算所述泄漏电流I的频率漂移μ；

应用公式： $\left\{\begin{array}{c}{a}_k=\frac{2}{Q}\underset{i=0}{\overset{W}{\mathrm{\Σ}}}{\mathrm{\γ}}_{i}f\left(i\right)\cos \left(k\frac{\mathrm{\μ}2\mathrm{\π}}{N}i\right)\\ b_k=\frac{2}{Q}\underset{i=0}{\overset{W}{\mathrm{\Σ}}}{\mathrm{\γ}}_{i}f\left(i\right)\sin \left(k\frac{\mathrm{\μ}2\mathrm{\π}}{N}i\right)\end{array}\right.,$

，分别得出所述泄漏电流I的基波的实部a_k和虚部b_k，然后计算出所述泄漏电流I的基波的幅值和初相角

(4)、所述CPU按照投影法获取泄漏电流I的基波在电网电压V的基波上的投影角

(5)、所述CPU获取MOA阻性泄漏电流的基波并输出显示，同时发送至所述监测服务器；

所述频率漂移μ是根据相邻采样点基波相角差与理想周期内采样点数N的固定关系而获得的，频率漂移μ用于修正基波的频率f₁。