[发明专利]基于多复合优化算法的电能质量检测分析方法有效
| 申请号: | 201910501158.X | 申请日: | 2019-06-11 | 
| 公开(公告)号: | CN110221147B | 公开(公告)日: | 2021-03-30 | 
| 发明(设计)人: | 李征;詹振辉;刘帅;孟浩 | 申请(专利权)人: | 东华大学 | 
| 主分类号: | G01R31/00 | 分类号: | G01R31/00 | 
| 代理公司: | 上海申汇专利代理有限公司 31001 | 代理人: | 翁若莹;柏子雵 | 
| 地址: | 201600 上*** | 国省代码: | 上海;31 | 
| 权利要求书: | 查看更多 | 说明书: | 查看更多 | 
| 摘要: | 本发明涉及一种多复合优化算法的电能质量检测分析方法。本发明先采用EEMD方法对采样的电信号进行去噪、分解。然后采用Prony算法来处理计算分解得到的每一个IMF的幅度、频率、衰减因子和初始相位角。为了提高计算的精度,准确度,本发明在计算的过程中加入了改进后的蜻蜓算法来与Prony算法相结合得到新的智能算法(本发明暂且称之为DP算法),来对IMF进行求解。为了得到更好的效果,本发明又在DP算法中引入了信噪比,以此来进行算法精确度的评判,使最终的结果更精确。 | ||
| 搜索关键词: | 基于 复合 优化 算法 电能 质量 检测 分析 方法 | ||
【主权项】:
                1.一种多复合优化算法的电能质量检测分析方法,其特征在于,包括以下步骤:步骤1、采用EEMD方法对采样的电力故障暂态信号去噪过后进行EEMD分解,EEMD分解后得到各个IMF分量,将第i个IMF分量记为cki(t),将每次得到的IMF分量做集成平均处理,则处理后得到的第i个IMF分量记为ci(t),![]() N为加入高斯白噪声的总体个数,N越大对应的高斯白噪声的IMF和将趋于0,分解结果也越趋近于真实值;步骤2、设置最大迭代次数Tmax,设置信噪比阈值SNR',设置迭代次数t=0;步骤3、将步骤1得到的每个IMF分量作为蜻蜓算法的初始值代入计算得到每个蜻蜓个体的分离度、对齐度、内聚度、食物吸引力及天敌排斥力,其中:第i个蜻蜓个体的分离度为Si,则有:
N为加入高斯白噪声的总体个数,N越大对应的高斯白噪声的IMF和将趋于0,分解结果也越趋近于真实值;步骤2、设置最大迭代次数Tmax,设置信噪比阈值SNR',设置迭代次数t=0;步骤3、将步骤1得到的每个IMF分量作为蜻蜓算法的初始值代入计算得到每个蜻蜓个体的分离度、对齐度、内聚度、食物吸引力及天敌排斥力,其中:第i个蜻蜓个体的分离度为Si,则有:![]() 式中,X表示第i个蜻蜓个体的位置,Xj表示与X相邻的蜻蜓个体的位置,与X相邻的蜻蜓个体共有N个;第i个蜻蜓个体的对齐度为Ai,则有:
式中,X表示第i个蜻蜓个体的位置,Xj表示与X相邻的蜻蜓个体的位置,与X相邻的蜻蜓个体共有N个;第i个蜻蜓个体的对齐度为Ai,则有:![]() 式中,Vj表示与X相邻的蜻蜓个体的速度;第i个蜻蜓个体的内聚度为Ci,则有:
式中,Vj表示与X相邻的蜻蜓个体的速度;第i个蜻蜓个体的内聚度为Ci,则有:![]() 第i个蜻蜓个体的食物吸引力为Fi,则有:Fi=X+‑X式中,X+表示食物源位置;第i个蜻蜓个体的天敌排斥力为Ei,则有:Ei=X‑+X式中,X‑表示天敌位置;步骤4、设置对齐度权重a、分离度权重s、内聚度权重c、惯性因子w,计算得到改进后的天敌因子e'及食物因子f':f'=λ·fe'=β·e式中,f表示改进前的食物因子,e表示改进前的天敌因子,
第i个蜻蜓个体的食物吸引力为Fi,则有:Fi=X+‑X式中,X+表示食物源位置;第i个蜻蜓个体的天敌排斥力为Ei,则有:Ei=X‑+X式中,X‑表示天敌位置;步骤4、设置对齐度权重a、分离度权重s、内聚度权重c、惯性因子w,计算得到改进后的天敌因子e'及食物因子f':f'=λ·fe'=β·e式中,f表示改进前的食物因子,e表示改进前的天敌因子,![]()
![]() 步骤5、根据改进后的天敌因子e'及食物因子f'得到改进后蜻蜓的更新位置Xt+1,则有:ΔXt+1=(sSi+aAi+cCi+f'Fi+e'Ei)+wΔXtXt+1=Xt+ΔXt+1式中,ΔXt表示更新步长,Xt当前第i个蜻蜓个体的位置;步骤6、并将新的位置参数Xt+1带入
步骤5、根据改进后的天敌因子e'及食物因子f'得到改进后蜻蜓的更新位置Xt+1,则有:ΔXt+1=(sSi+aAi+cCi+f'Fi+e'Ei)+wΔXtXt+1=Xt+ΔXt+1式中,ΔXt表示更新步长,Xt当前第i个蜻蜓个体的位置;步骤6、并将新的位置参数Xt+1带入![]() 来求解
来求解![]() 将结果带入
将结果带入![]() 来求取目标函数,其中:
来求取目标函数,其中:![]() 表示x(n)的拟合估计值,ak表示特征方程系数,bk表示x(n)拟合估计值对应的复数,
表示x(n)的拟合估计值,ak表示特征方程系数,bk表示x(n)拟合估计值对应的复数,![]() 表示x(n)拟合估计值对应的复数,z表示目标函数;步骤7、计算得到ak、Ak、θk及fk,Ak为第k阶振幅,fk为第k阶频率,θk为第k阶相位:Ak=|bk|
表示x(n)拟合估计值对应的复数,z表示目标函数;步骤7、计算得到ak、Ak、θk及fk,Ak为第k阶振幅,fk为第k阶频率,θk为第k阶相位:Ak=|bk|![]()
![]()
![]() 式中,zk表示多项式求解得到的根,Δt表示采样间隔;步骤8、计算得到SNR,
式中,zk表示多项式求解得到的根,Δt表示采样间隔;步骤8、计算得到SNR,![]() t=t+1;步骤9、若t≥Tmax或SNR<SNR',则结束本方法,否则,返回步骤3。
t=t+1;步骤9、若t≥Tmax或SNR<SNR',则结束本方法,否则,返回步骤3。
            
                    下载完整专利技术内容需要扣除积分,VIP会员可以免费下载。
                
                
            该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于东华大学,未经东华大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/patent/201910501158.X/,转载请声明来源钻瓜专利网。





