CN104979829A - 一种基于改进的FastICA算法的谐波电流估计方法 - Google Patents
一种基于改进的FastICA算法的谐波电流估计方法 Download PDFInfo
- Publication number
- CN104979829A CN104979829A CN201510332421.9A CN201510332421A CN104979829A CN 104979829 A CN104979829 A CN 104979829A CN 201510332421 A CN201510332421 A CN 201510332421A CN 104979829 A CN104979829 A CN 104979829A
- Authority
- CN
- China
- Prior art keywords
- harmonic
- newton
- matrix
- harmonic current
- fastica
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
- Y02E40/40—Arrangements for reducing harmonics
Landscapes
- Measurement Of Resistance Or Impedance (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了一种基于改进的FastICA算法的谐波电流估计方法,其特征在于,采用拟牛顿法代替牛顿法作为FastICA的优化方法。本发明的有益效果:(1)本发明能在谐波阻抗未知情况下准确估计电网中的谐波电流;(2)克服了系统实际谐波阻抗难以获得和谐波测量装置成本昂贵的难题;(3)原有方法相比,提高了算法的精度,更好的解决了谐波阻抗难以确定的难题。
Description
技术领域
本发明涉及一种基于改进的FastICA算法的谐波电流估计方法,属于电力系统电能质量技术领域。
背景技术
在谐波研究领域中,谐波源的识别是一个非常重要的问题,它对谐波责任的划分和谐波的治理有着重要的意义。现有的方法大多是基于谐波阻抗已知情况下的谐波源识别,然而在实际系统中,要得到电网的准确谐波参数是比较困难的。因此,需要加强研究如何在谐波阻抗未知情况下,利用已知的谐波量测信息估计谐波源的状态。牛顿算法的优点是具有二阶局部收敛特性,缺点是对初始点目标函数的要求高,使计算量和储存量增大。
发明内容
本发明的目的克服现有技术存在的不足,提出一种基于改进的FastICA算法的谐波电流估计方法,解决了谐波测量设备成本昂贵和系统实际谐波阻抗难以获得的技术问题。
本发明采用如下技术方案:一种基于改进的FastICA算法的谐波电流估计方法,其特征在于,采用拟牛顿法代替牛顿法作为FastICA的优化方法。
优选地,所述拟牛顿法的迭代公式:
式(7)中:αk为线性搜索步长因子,一般情况下取αk=1,Bk为拟牛顿矩阵。
所述拟牛顿法的限定条件:
式(8)中:sk=Xk+1-Xk,yk=f(Xk+1)-f(Xk)。
优选地,具体包括如下步骤:
SS1测量母线的谐波电压;
SS2利用线性均值滤波器把谐波电压分解成为快速变化分量和缓慢变化分量;
SS3对谐波电压的快速变化分量,利用Quasi-Newton FastICA算法估计出系统的谐波阻抗矩阵;
SS4根据缓慢变化分量、原信号的先验知识、谐波阻抗矩阵W,利用Quasi-Newton FastICA算法估计出系统的谐波电流。
优选地,所述步骤SS1具体包括:选择谐波注入电流作为谐波潮流方程中的状态变量,而以节点的谐波电压作为量测量,忽略量测误差,谐波电流估计的模型为:
Uh=ZhIh (2)
式中,Uh——n维节点谐波电压矩阵;Zh——n×n阶谐波阻抗矩阵;Ih——n维节点谐波注入电流矩阵;h是谐波次数;
如果谐波电压矩阵Uh和谐波阻抗矩阵Zh已知,则谐波电流矩阵的提取如式(3)或式(4)所示;
Ih=YhUh (4)
Yh——系统的谐波导纳矩阵;
如果谐波阻抗矩阵Zh未知,则可以利用式(5)估计谐波电流;
Iest-h=WhUh (5)
Iest-h为估计的谐波电流,Wh为待估计的h次谐波节点导纳矩阵。
优选地,所述Quasi-Newton FastICA算法包括:输入量为预处理后的混合信号X=(x1,x2,…,xm);输出量为源信号的估计值Y=(y1,y2,…,yn);具体步骤如下:
步骤(1)令k=0(k=0,1,…,n),初始化权向量矩阵W0;
步骤(2)k=k+1;
步骤(3)对W进行调整
步骤(4)归一化处理Wk+1=Wk+1/||Wk+1||;
步骤(5)若算法没有收敛,则回到步骤(3),否则转入步骤(6);
步骤(6)若算法收敛,则求出一个独立分量y1=wx。
本发明所达到的有益效果:(1)本发明能在谐波阻抗未知情况下准确估计电网中的谐波电流;(2)克服了系统实际谐波阻抗难以获得和谐波测量装置成本昂贵的难题;(3)原有方法相比,提高了算法的精度,更好的解决了谐波阻抗难以确定的难题。
名词解释:Newton法—牛顿法;Quasi-Newton法—拟牛顿法。
附图说明
图1是本发明的一种基于改进的FastICA算法的谐波电流估计方法的流程框图。
图2是本发明的一个实施例的结构示意图。
图3是节点7的电压波形图。
图4是节点16的电压波形图。
图5是节点28的电压波形图。
图6是节点7采用Newton FastICA算法估计结果的谐波电流波形图。
图7是节点7采用Quasi-Newton FastICA算法估计结果的谐波电流波形图。
图8是节点16采用Newton FastICA算法估计结果的谐波电流波形图。
图9是节点16采用Quasi-Newton FastICA算法估计结果的谐波电流波形图。
图10是节点28采用Newton FastICA算法估计结果的谐波电流波形图。
图11是节点28采用Quasi-Newton FastICA算法估计结果的谐波电流波形图。
图12是电压采集数量分别对Newton FastICA算法和Quasi-Newton FastICA算法的误差影响图。
具体实施方式
下面结合附图对本发明作进一步描述。以下实施例仅用于更加清楚地说明本发明的技术方案,而不能以此来限制本发明的保护范围。
图1是本发明的一种基于改进的FastICA算法的谐波电流估计方法的流程框图,本发明提出一种基于改进的FastICA算法的谐波电流估计方法,其特征在于,采用拟牛顿法代替牛顿法作为FastICA的优化方法。
具体包括如下步骤:
SS1测量母线的谐波电压:选择谐波注入电流作为谐波潮流方程中的状态变量,而以节点的谐波电压作为量测量,忽略量测误差,谐波电流估计的模型为:
Uh=ZhIh (2)
式中,Uh——n维节点谐波电压矩阵;Zh——n×n阶谐波阻抗矩阵;Ih——n维节点谐波注入电流矩阵;h是谐波次数;
如果谐波电压矩阵Uh和谐波阻抗矩阵Zh已知,则谐波电流矩阵的提取如式(3)或式(4)所示;
Ih=YhUh (4)
Yh——系统的谐波导纳矩阵;
如果谐波阻抗矩阵Zh未知,则可以利用式(5)估计谐波电流;
Iest-h=WhUh (5)
Iest-h为估计的谐波电流,Wh为待估计的h次谐波节点导纳矩阵。
SS2利用线性均值滤波器把谐波电压分解成为快速变化分量和缓慢变化分量;
SS3对谐波电压的快速变化分量,利用Quasi-Newton FastICA算法估计出系统的谐波阻抗矩阵;
SS4根据缓慢变化分量、原信号的先验知识、谐波阻抗矩阵W,利用Quasi-Newton FastICA算法估计出系统的谐波电流。
最优化问题常用的一些方法都是基于牛顿算法的优化方法,FastICA也通常以牛顿法作为优化算法。牛顿法的迭代公式:
Xk+1=Xk-[f'(Xk)]-1f(Xk)] (6)
但是牛顿算法的缺点是对初始点目标函数的要求高,使计算量和储存量增大。而拟牛顿法克服了以上缺点,其原理是用Hesse近似Bk代替实际目标函数真实二阶导数。拟牛顿法的迭代公式:
式(7)中:αk为线性搜索步长因子,一般情况下取αk=1,Bk为拟牛顿矩阵。
拟牛顿法的限定条件:
式(8)中:sk=Xk+1-Xk,yk=f(Xk+1)-f(Xk)
拟牛顿法最大优点是计算准确,具有很强的收敛性,是一种有效的优化算法。以拟牛顿法作为FastICA的优化算法,可以改进原有FastICA的性能。
FastICA算法是独立分量分析的常用算法,通常用牛顿法作为它的迭代算法,目标函数采用最大化负熵的形式,然后对测量量X进行批量处理,算法每次计算都从观测量信号中分离出一个独立成分分量。
随机向量y的负熵定义为:
J(y)=H(yG)-H(y) (9)
式中,yG为与y具有相同的均值与方差的高斯变量;H(y)为概率密度为p(y)的向量y的信息熵:
H(y)=-∫p(y)logp(y)dy (10)
根据信息论的原理,在相同方差的随机信号中,高斯信号的信息熵总是最大的,由式(9)和(10)可知负熵总是具有非负的值,随机变量的非高斯性越强则负熵J(y)越大,因此可以用负熵J(y)来反映随机变量的非高斯性。因为信息熵p(y)难以求出,于是经常采用如式(11)所示的近似公式代替式(10)来计算负熵J(y):
J(y)={E[g(y)]-E[g(yG)]}2 (11)
式中,E[·]表示均值运算,g(·)代表了非线性函数。
FastICA算法的目的在于找到一个最优方向,使输出的随机变量y=WTx达到最大的非高斯性。而非高斯性通常用式子(11)中的负熵J(WTx)的近似值来衡量。算法的目的就是要使J(y)=J(WTx)达到最大值,即E[g(y)]=E[g(WTx)]取最大值。根据拟牛顿法得到Quasi-NewtonFastICA的迭代公式:
所述Quasi-Newton FastICA算法包括:输入量为预处理后的混合信号X=(x1,x2,…,xm);输出量为源信号的估计值Y=(y1,y2,…,yn);具体步骤如下:
步骤(1)令k=0(k=0,1,…,n),初始化权向量矩阵W0;
步骤(2)k=k+1;
步骤(3)对W进行调整
步骤(4)归一化处理Wk+1=Wk+1/||Wk+1||;
步骤(5)若算法没有收敛,则回到步骤(3),否则转入步骤(6);
步骤(6)若算法收敛,则求出一个独立分量y1=wx。
图2是本发明的一个实施例的结构示意图,采用39节点系统模型,在节点7,节点16和节点28处接入了谐波负荷(以5次谐波为例),用MATLAB软件进行仿真;由于系统含有三个独立的谐波源,可以选择三个母线作为谐波电压的测量点。以节点7、16、28为测量点,从这三个节点得到的电压数据作为量测量。
图3是节点7的电压波形图,图4是节点16的电压波形图,图5是节点28的电压波形图,图中横坐标表示采样点数n;本发明的所做的仿真,负荷的电压曲线满足随机性的特点,所以按快速变化分量来处理,谐波信号中不包含缓慢变化分量。所以不需要用线性均值滤波器分离出缓慢变化分量和快速变化分量,对所测的谐波信号可以直接运用FastICA算法进行处理。
根据图2中的谐波电压数据,分别运用Newton FastICA算法和Quasi-Newton FastICA算法对节点7、节点16和节点28的谐波电流进行估计。估计结果如图6~11所示,图中实线表示实测到的谐波电流,虚线表示用算法估计出的电流。
图6是节点7采用Newton FastICA算法估计结果的谐波电流波形图,图7是节点7采用Quasi-Newton FastICA算法估计结果的谐波电流波形图;图8是节点16采用Newton FastICA算法估计结果的谐波电流波形图,图9是节点16采用Quasi-Newton FastICA算法估计结果的谐波电流波形图;图10是节点28采用Newton FastICA算法估计结果的谐波电流波形图,图11是节点28采用Quasi-NewtonFastICA算法估计结果的谐波电流波形图。
FastICA算法估计得到的电流与实际谐波电流之间的误差用相关系数来表示。相关系数越接近1,表示算法的误差越小。表1计算出了取500个电压采集点时的相关系数。表2计算出了取1500个电压采集点时的相关系数。
表1电压采集样本为500时的相关系数
表2电压采集样本为1500时的相关系数
由表1和表2可知,与Newton FastICA算法相比,Quasi-NewtonFastICA算法估计结果精度更高,算法性能更优。
由表1和表2可知,电压采集点数会对算法的精度产生影响。图12是电压采集数量分别对Newton FastICA算法和Quasi-NewtonFastICA算法的误差影响图,由结果可知,与Newton FastICA算法相比,Quasi-Newton FastICA算法受电压采集点数的影响较小。因而在相同的电压采集数量条件下,Quasi-Newton FastICA算法估计谐波电流的结果更加准确。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明技术原理的前提下,还可以做出若干改进和变形,这些改进和变形也应视为本发明的保护范围。
Claims (5)
1.一种基于改进的FastICA算法的谐波电流估计方法,其特征在于,采用拟牛顿法代替牛顿法作为FastICA的优化方法。
2.根据权利要求1所述的一种基于改进的FastICA算法的谐波电流估计方法,其特征在于,所述拟牛顿法的迭代公式:
式(7)中:αk为线性搜索步长因子,一般情况下取αk=1,Bk为拟牛顿矩阵。
所述拟牛顿法的限定条件:
式(8)中:sk=Xk+1-Xk,yk=f(Xk+1)-f(Xk)。
3.根据权利要求1所述的一种基于改进的FastICA算法的谐波电流估计方法,其特征在于,具体包括如下步骤:
SS1测量母线的谐波电压;
SS2利用线性均值滤波器把谐波电压分解成为快速变化分量和缓慢变化分量;
SS3对谐波电压的快速变化分量,利用Quasi-Newton FastICA算法估计出系统的谐波阻抗矩阵W;
SS4根据缓慢变化分量、原信号的先验知识、谐波阻抗矩阵W,利用Quasi-Newton FastICA算法估计出系统的谐波电流。
4.根据权利要求3所述的一种基于改进的FastICA算法的谐波电流估计方法,其特征在于,所述步骤SS1具体包括:选择谐波注入电流作为谐波潮流方程中的状态变量,而以节点的谐波电压作为量测量,忽略量测误差,谐波电流估计的模型为:
Uh=ZhIh (2)
式中,Uh——n维节点谐波电压矩阵;Zh——n×n阶谐波阻抗矩阵;Ih——n维节点谐波注入电流矩阵;h是谐波次数;
如果谐波电压矩阵Uh和谐波阻抗矩阵Zh已知,则谐波电流矩阵的提取如式(3)或式(4)所示;
Ih=YhUh (4)
Yh——系统的谐波导纳矩阵;
如果谐波阻抗矩阵Zh未知,则可以利用式(5)估计谐波电流;
Iest-h=WhUh (5)
Iest-h为估计的谐波电流,Wh为待估计的h次谐波节点导纳矩阵。
5.根据权利要求3所述的一种基于改进的FastICA算法的谐波电流估计方法,其特征在于,所述Quasi-Newton FastICA算法包括:输入量为预处理后的混合信号X=(x1,x2,…,xm);输出量为源信号的估计值Y=(y1,y2,…,yn);具体步骤如下:
步骤(1)令k=0(k=0,1,…,n),初始化权向量矩阵W0;
步骤(2)k=k+1;
步骤(3)对W进行调整
步骤(4)归一化处理Wk+1=Wk+1/||Wk+1||;
步骤(5)若算法没有收敛,则回到步骤(3),否则转入步骤(6);
步骤(6)若算法收敛,则求出一个独立分量y1=wx。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510332421.9A CN104979829B (zh) | 2015-06-16 | 2015-06-16 | 一种基于改进的FastICA算法的谐波电流估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510332421.9A CN104979829B (zh) | 2015-06-16 | 2015-06-16 | 一种基于改进的FastICA算法的谐波电流估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104979829A true CN104979829A (zh) | 2015-10-14 |
CN104979829B CN104979829B (zh) | 2018-07-24 |
Family
ID=54276067
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510332421.9A Expired - Fee Related CN104979829B (zh) | 2015-06-16 | 2015-06-16 | 一种基于改进的FastICA算法的谐波电流估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104979829B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104280612A (zh) * | 2014-10-28 | 2015-01-14 | 清华大学 | 一种基于单频电流传输特性的分布式谐波源辨识方法 |
CN105911334A (zh) * | 2016-04-19 | 2016-08-31 | 国网上海市电力公司 | 一种谐波阻抗未知条件下的谐波电流估计方法 |
CN107688135A (zh) * | 2017-07-04 | 2018-02-13 | 广西大学 | 一种基于FastICA的小电流接地系统的单相接地故障的选线方法 |
CN109193787A (zh) * | 2018-10-28 | 2019-01-11 | 国网吉林省电力有限公司四平供电公司 | 一种含新能源配电网的谐波路径识别方法 |
CN110161373A (zh) * | 2019-06-06 | 2019-08-23 | 广东电网有限责任公司 | 一种不平衡源定位方法、系统及装置 |
CN114336795A (zh) * | 2022-01-18 | 2022-04-12 | 华南理工大学 | 一种适用于交直流混联电网谐波状态的综合优化配置方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103743949A (zh) * | 2014-01-06 | 2014-04-23 | 国家电网公司 | 一种基于单通道FastICA的谐波与间谐波的检测方法 |
US20140341236A1 (en) * | 2013-05-03 | 2014-11-20 | Zte (Usa) Inc. | Method and apparatuses for algorithm on qam coherent optical detection |
-
2015
- 2015-06-16 CN CN201510332421.9A patent/CN104979829B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20140341236A1 (en) * | 2013-05-03 | 2014-11-20 | Zte (Usa) Inc. | Method and apparatuses for algorithm on qam coherent optical detection |
CN103743949A (zh) * | 2014-01-06 | 2014-04-23 | 国家电网公司 | 一种基于单通道FastICA的谐波与间谐波的检测方法 |
Non-Patent Citations (1)
Title |
---|
任立志等: ""基于ICA和MI的谐波源识别研究"", 《电测与仪表》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104280612A (zh) * | 2014-10-28 | 2015-01-14 | 清华大学 | 一种基于单频电流传输特性的分布式谐波源辨识方法 |
CN105911334A (zh) * | 2016-04-19 | 2016-08-31 | 国网上海市电力公司 | 一种谐波阻抗未知条件下的谐波电流估计方法 |
CN107688135A (zh) * | 2017-07-04 | 2018-02-13 | 广西大学 | 一种基于FastICA的小电流接地系统的单相接地故障的选线方法 |
CN107688135B (zh) * | 2017-07-04 | 2020-04-17 | 广西大学 | 一种基于FastICA的小电流接地系统的单相接地故障的选线方法 |
CN109193787A (zh) * | 2018-10-28 | 2019-01-11 | 国网吉林省电力有限公司四平供电公司 | 一种含新能源配电网的谐波路径识别方法 |
CN109193787B (zh) * | 2018-10-28 | 2021-07-06 | 国网吉林省电力有限公司四平供电公司 | 一种含新能源配电网的谐波路径识别方法 |
CN110161373A (zh) * | 2019-06-06 | 2019-08-23 | 广东电网有限责任公司 | 一种不平衡源定位方法、系统及装置 |
CN110161373B (zh) * | 2019-06-06 | 2022-06-21 | 广东电网有限责任公司 | 一种不平衡源定位方法、系统及装置 |
CN114336795A (zh) * | 2022-01-18 | 2022-04-12 | 华南理工大学 | 一种适用于交直流混联电网谐波状态的综合优化配置方法 |
CN114336795B (zh) * | 2022-01-18 | 2024-03-22 | 华南理工大学 | 一种适用于交直流混联电网谐波状态的综合优化配置方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104979829B (zh) | 2018-07-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104979829A (zh) | 一种基于改进的FastICA算法的谐波电流估计方法 | |
Evrenosoglu et al. | Travelling wave based fault location for teed circuits | |
CN103412242B (zh) | 一种基于快速独立分量分析和互信息的谐波源定位方法 | |
CN110082640B (zh) | 一种基于长短时记忆网络的配网单相接地故障辨识方法 | |
CN102135569B (zh) | 基于波动量法的用户侧谐波发射水平实用化的估计方法 | |
CN110618316B (zh) | 一种谐波源辨识方法、装置、终端设备和存储介质 | |
CN108155648A (zh) | 基于自适应h无穷扩展卡尔曼滤波的状态估计方法 | |
CN109444505B (zh) | 一种电动汽车充电站谐波电流检测算法 | |
CN107658881A (zh) | 基于戴维南等值方法的电压稳定临界点判断方法 | |
CN105548718A (zh) | 一种基于混合整体最小二乘法的系统谐波阻抗计算方法 | |
CN107370150A (zh) | 基于同步相量量测的电力系统状态估计不良数据处理方法 | |
Yang et al. | Non‐linear characteristic quantity extraction of ferroresonance overvoltage time series | |
CN107064629A (zh) | 一种基于频率相对偏差预估的分段综合单频信号频率估计方法 | |
CN101576586A (zh) | 一种基于自适应滤波的相位检测方法 | |
CN103823158B (zh) | 采用不变矩的谐振接地系统故障选线方法 | |
CN110161375B (zh) | 一种基于分布电阻参数的高压直流输电线路计算模型 | |
CN111537776A (zh) | 一种行波波头标定方法、装置、终端以及介质 | |
CN115201563A (zh) | 一种基于联合熵的多谐波源定位方法及系统 | |
Sangeeth et al. | High impedance fault detection using multi-domain feature with artificial neural network | |
Davis et al. | Estimation of transmission line parameters from historical data | |
CN109638811B (zh) | 基于模型等值的配电网电压功率灵敏度鲁棒估计方法 | |
Li et al. | A new method to estimate the state of charge of the green battery | |
CN104167732B (zh) | 一种基于联络线相角差的电网等值方法 | |
Carta et al. | Harmonics detector in distribution systems based on compressive sensing | |
Amin et al. | Kalman filter estimation of impedance parameters for medium transmission line |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180724 Termination date: 20210616 |