1.一种基于APF直流电压信息反演的工业园区谐波监测方法，其特征在于：其步骤是：S1、三相之间采用星型连接方式，其中us为三相电网电压，is为网侧电流，iL为负载电流，ic为APF输出补偿电流，udc为直流侧电容电压，idc为流过电容的电流，LC为连接电抗；根据电力系统基本原理得到APF系统输出到电网的三相瞬时有功功率为：其中Us为ABC三相相电压有效值，I(n)和 为n次输出电流对应的有效值和相位，n＝±(6k±1)以及n＝‑1时对应三相瞬时功率的交流部分，其中k＝1,2,3…，n＝1时对应三相瞬时功率的直流部分；

APF稳态运行时，直流侧电容电压由直流恒定分量与交流分量组成，即：根据APF交直流两侧瞬时功率相等的原则得：

令m＝n‑1，由式(3)进一步得出直流侧电容电压中角频率为mω的交流纹波分量的具体表达式：通过以上对直流侧电容电压纹波的理论分析，得出以下结论：(1)当APF补偿m+1次正序谐波电流时，其中m+1＞0，或者|‑m+1|次负序谐波电流时，其中‑m+1＜0，直流侧电容电压将会产生m倍频的波动，直流侧电容电压纹波频率与补偿电流频率的对应关系(2)在APF输出补偿电流相序及频率不变的前提下，APF输出补偿电流幅值I(m+1)或I(‑m+1)越大，其对应的直流侧电容电压纹波幅值也越大；

S2、获得直流侧电容电压信号的奇异本征模函数

(1)利用CEEMD分解直流侧电容电压信号，获得n个本征模函数，并按照式(5)、式(6)计算每个IMF的其中：T为傅里叶频谱的周期，SlnT,i为以lnT为自变量的第i个IMF的傅里叶频谱， 为从第i个IMF的频谱中计算出的平均周期，Ei为第i个IMF的能量密度，N表示第i个IMF的数据长度；

(2)判断行程点

令满足式(7)所示条件的第一个IMF为行程点p，并根据式(8)计算该行程点之前所有IMF的 的平均值其中， 为第i个IMF的能量，i＝1,2,…,n‑1；

其中， 为行程点p之前所有IMF的平均能量；

(3)确定奇异本征模函数并进行阈值优化

令满足式(9)所示条件的本征模函数为奇异本征模函数，在确定奇异本征模函数后按照式(10)对其进行阈值优化其中i＝p+1,p+2,…,n‑1，p为行程点；

其中：IMFi′(t)为优化阈值后的奇异本征模函数，τi为阈值，σi为标准差；

S3、对优化后的IMF应用希尔伯特变换

直流侧电容电压信号经过CEEMD处理后表示为：其中，rn为残差，m为奇异本征模函数的个数；

对不同频率的IMF进行希尔伯特变换可以得到：

构建解析信号：

根据式(14)、式(15)得到IMFi′的瞬时幅值、瞬时相位以及瞬时频率：S4、通过查表法确定配网内谐波次数

根据式(16)获得直流侧电容电压纹波的频次及幅值，通过与步骤S1的直流侧电容电压纹波频率与补偿电流频率的对应关系对照比对，根据直流侧电容电压纹波频次反演出配网内谐波次数，并且根据直流侧电容电压各频次纹波的幅值可以反演出内谐波负载的变化情况，继而掌握系统受谐波的污染程度。