CN103245830B - 一种结合ar谱估计与非线性优化的间谐波检测方法 - Google Patents
一种结合ar谱估计与非线性优化的间谐波检测方法 Download PDFInfo
- Publication number
- CN103245830B CN103245830B CN201310115657.8A CN201310115657A CN103245830B CN 103245830 B CN103245830 B CN 103245830B CN 201310115657 A CN201310115657 A CN 201310115657A CN 103245830 B CN103245830 B CN 103245830B
- Authority
- CN
- China
- Prior art keywords
- parameter
- harmonic
- formula
- nonlinear optimization
- iteration
- 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.)
- Active
Links
Landscapes
- Measuring Frequencies, Analyzing Spectra (AREA)
- Complex Calculations (AREA)
Abstract
一种结合AR谱估计与非线性优化的间谐波检测方法,其特征应用AR模型谱估计计算谐波初始参数来建立信号模型:选取待分析的采样信号序列,对其进行AR模型谱估计,给出精度较低的谐波参数(频率、幅值和相位)信息;再取模型,提出求解信号模型参数的非线性优化算法,对相位值进行计算。本发明的有益效果是:解决了非线性优化算法的参数建模问题,克服了AR谱估计不能计算谐波幅值和相位的缺点,进一步提高了频率的计算精度;降低迭代算法对初值的敏感度,提高了迭代稳定性和计算效率;且有一定的抗噪能力。
Description
技术领域
本发明涉及一种电网信号间谐波分析领域,具体的说,涉及一种结合AR谱估计与非线性优化的间谐波检测方法。
背景技术
电力电子等非线性设备在电力系统中的广泛应用,谐波和间谐波日益增多,因而分析谐波和间谐波对电力系统的经济安全运行有着重要意义。目前的谐波分析主要通过对电网信号采样和数字化处理实现的,如快速傅里叶变换(FFT)和加窗插值类FFT方法。但这些方法在检测谐波时不可避免会产生频谱泄漏和栅栏效应,影响检测精度,用于间谐波(非整数次谐波)检测时效果更差,原因在于间谐波常出现在整数次谐波附近的频点且幅值较小,易被整数次谐波的旁瓣淹没,导致其频率分辨率大大降低,而低分辨率造成的谱线识别错误又会进一步地降低幅值与相位的检测精度。
现有技术的缺点是:难以获得电力线路中各次谐波的频率、幅值和相位的准确值,无法描述出含谐波干扰的电路波形。
发明内容
本发明针对现有技术的不足,提供了一种结合AR谱估计与非线性优化的间谐波检测方法,用于计算谐波参数的非线性优化算法,实现电网信号的频率、幅值和相位的准确计算。
为了实现上述发明目的,本发明提供一种结合AR谱估计与非线性优化的间谐波检测方法。
一种结合AR谱估计与非线性优化的间谐波检测方法,其特征在于,按以下步骤进行:
步骤一:AR谱初步估计
1.1):根据谐波参数的初步信息建立该谐波的参数模型,取参数模型如下:
式中y(t)为待分析信号,Am为各次谐波和间谐波的幅值,frm为对应的频率,θm为对应的相位,m为对应的谐波个数,t为时间;
1.2):对求解信号模型参数进行非线性优化算法:将所述模型(1)作如下处理,将幅值、相位(Am,θm)等价转化为求解两个幅值(Cm,Dm):
Amsin(2πfrmt+θm)=Amsin(2πfrmt)cos(θm)+Amcos(2πfrmt)sin(θm)
=Cmsin(2πfrmt)+Dmcos(2πfrmt) (2)
其中,Cm=Amcos(θm),Dm=Amsin(θm);
1.3):针对步骤一中的参数模型(1),将式(2)代入式(1),建立最小二乘方程如下:
式中n为采样点数;目标函数表征了信号计算结果和实际采样结果的误差,令y(tn)为信号第n点的计算值、y1(tn)为信号第n点实际采样值,将参数fm、Cm和Dm统一用向量x表示,对应于各采样时刻的误差向量用f(x)表示:
f(x)=[f1,f2,...,fn]
=[y(t1)-y1(t1),(t2)-y1(t2),...,(tn)-y1(tn)]
代入(3)式可得:
式中[f(x)]T表示f(x)的转置,则参数的求解转化为求如下最值的问题:
步骤二:共轭梯度法一次迭代,也就是参数初值进行预先处理,使参数迭代值接近全局最优解。
2.1)对参数初值进行迭代,迭代方程为:
xk+1=xk+λkpk (5)
其中,k为迭代次数,λk为最优步长,pk为迭代方向,
2.2)pk满足下式(7)
置初始近似值x0∈Rn(Rn为实数域),一次迭代误差限ε1>0
2.3):计算,令0→k,p0是迭代方向矩阵pk的初值;
2.4):求最优步长λk满足,计算迭代值xk+1=xk+λkpk以及
2.5):计算,当迭代到时,转步骤三,否则,转2.6;
2.6):计算k+1→k,转2.4;
步骤三:当后,再用阻尼最小二乘法二次迭代:即采用阻尼最小二乘法快速寻找到全局最优解,即最终的间谐波参数计算值。
3.1):为使沿着其梯度方向减小,必有的法方程等于零:
可得:
式中:G(xb)=Df(xb)TDf(xb),g(xb)=Df(xb)Tf(xb),p(xb)=-G(xb)-1g(xb);
3.2):考虑阻尼时有:
其中μb为第b次迭代的阻尼因子,I是与G(xb)同维数的单位对角阵;则最终迭代公式为:
当μb>0总可以保证是收敛的,原则上μb=10-4~10-2;
3.3):计算得出f(xb),Df(xb),令0→j(j为3.7步骤所需的一个判断变量);
3.4):解方程组求得
3.5):计算xb+1=xb+p(xb)及
3.6):若则取(缩放因子v>1,可取2,5,10),转3.3),否则转3.7);
3.7):若j=0,则取μb=vμb,1→j,转3.3),否则转3.8);
3.8):若||p(xb)||≤ε2,则xb+1为极小值点x*的近似,算法停止,输出xb+1;否则取xb+1→xb,b+1→b,转3.3)。
本发明所述步骤一中的参数模型式(1)中的m次数可根据功率谱中局部最大值个数确定,然后依次确定功率谱中每一个局部最大值处对应的频率和局部最大值的开方,把它们分别作为参数fm、Am的初值,然后用0~2π内的随机值作为相位参数θm的初值。
本发明所述步骤三中的3.2)阻尼因子μb,当μb>0总可以保证因而是收敛的,当μb太大时会收敛速度下降,若μb太小则收敛域过小,原则上μb=10-4~10-2,通常取μb=10-2。
本发明AR模型谱估计有较好的频率分辨率,能检测出间谐波频率,但其幅值参数计算结果与实际值相距甚远且不能给出相位参数,所以需结合非线性优化算法来搜寻准确的谐波参数;应用AR谱估计初步估算信号的谐波参数,给出频率和幅值估值,相位初始参数在0~2π内随机选取;然后基于前述参数建立信号模型,再应用非线性优化算法求取最终的模型参数,给出间谐波检参数检测结果。在非线性优化算法中,为克服阻尼最小二乘法对初值的敏感而选择共轭梯度法对参数初值进行预先处理,使参数迭代值接近全局最优解,再用阻尼最小二乘法快速寻找到全局最优解,即最终的间谐波参数计算值。
式(1)中各次谐波的幅值Am和频率frm是可以准确获得,但相位θm为假设值,后续步骤就是为使相位θm与真实值相吻合。
步骤三:当后,再用阻尼最小二乘法二次迭代:即采用阻尼最小二乘法快速寻找到全局最优解,即最终的间谐波参数计算值。
3.1):由于Df(x)为f(x)的偏导矩阵:
则可得到:
为了使沿着其梯度方向减小,必有的法方程等于零:
将f(x)在xb(b为迭代次数)处作泰勒公式线性展开得:
又[f″(xb)]T(x-xb)2以后各项的值都远小于前两项,故近似有:
f(x)≈f(xb)+Df(xb)T(x-xb) (9)
令lb(x)=Df(xb)T(x-xb)+f(xb),代入法方程(8)得到:
0=Df(x)Tf(x)≈Df(xb)Tlb(x)
=Df(xb)T[Df(xb)(x-xb)+f(xb)]
=Df(xb)TDf(xb)(x-xb)+Df(xb)Tf(xb)
则迭代公式为:
为表示方便,上式可进一步写为:
式中:G(xb)=Df(xb)TDf(xb),g(xb)=Df(xb)Tf(xb),p(xb)=-G(xb)-1g(xb);
3.2):考虑阻尼时有:
其中μb为第b次迭代的阻尼因子,I是与G(xb)同维数的单位对角阵;则最终迭代公式为:
当μb>0总可以保证是收敛的;
3.3)计算得出f(xk),Df(xk),令0→j(j为3.7步骤所需的一个判断变量);
3.4)解方程组求得
3.5)计算xb+1=xb+p(xb)及
3.6)若则取(缩放因子v>1,可取2,5,10),转3.3),否则转3.7);
3.7)若j=0,则取μb=vμb,1→j,转3.3),否则转3.8);
3.8)若||p(xb)||≤ε2,则xb+1为极小值点x*的近似,算法停止,输出xb+1;否则取xb+1→xb,b+1→b,转3.3)。
本发明是通过应用AR模型谱估计计算谐波初始参数来建立信号模型的方法:将待分析的采样信号序列进行AR模型谱估计,给出精度较低的谐波参数(频率、幅值和相位)信息;然后基于初步的谐波信息建立信号的参数模型,再取模型进行式(1)计算。模型(1)为非线性模型,需采用非线性优化算法;同时,为减小模型的非线性,将参数模型进行处理,将幅值、相位(Am,θm)等价转化为求解两个幅值(Cm,Dm)进行计算。并且将式中参数frm、Cm和Dm统一用向量x表示,则对应于各采样时刻的误差向量用f(x)表示,进而将非线性优化问题转化最小值问题进行求解。
其中目标函数表征了信号计算结果和实际采样结果的误差,差别越小表明参数计算精度越高。
所述参数模型式(1)中的m次数可根据功率谱中局部最大值个数确定,然后依次确定功率谱中每一个局部最大值处对应的频率和局部最大值的开方,分别作为参数frm、Am的初值,然后用0~2π内的随机值作为相位参数θm的初值。
由于谐波的出现使待分析信号波形产生了畸变,所以通过式(1)中参数的确定可通过选取AR功率谱中的局部最大值实现,然后用0~2π内的随机值作为相位参数θm的迭代初值;
所述参数模型式(1)中的m次数可根据功率谱中局部最大值个数确定,然后依次确定功率谱中每一个局部最大值处对应的频率和局部最大值的开方,把它们分别作为参数frm、Am的初值,然后用0~2π内的随机值作为相位参数θm的初值。
所述步骤三中的3.2的阻尼因子μk,当μb>0总可以保证因而是收敛的,当μb太大时会收敛速度下降,若μb太小则收敛域过小,原则上μb=10-4~10-2,通常取μb=10-2;
本发明的显著效果是:
利用AR谱估计给出模型的结构和初始参数,解决了非线性优化算法的参数建模问题,也给出了参数的初值;结合非线性优化算法,克服了AR谱估计不能计算谐波幅值和相位的缺点,且进一步提高了频率的计算精度。
非线性优化算法的设计中,首先选用共轭梯度法对参数初值进行预处理,使之接近全局最优解;随后再使用阻尼最小二乘法能使参数迭代的谐波参数计算值),降低迭代算法对初值的敏感度,提高了迭代稳定性和计算效率;
AR谱估计与非线性优化结合方法的参数计算精度高出加Hanning窗插值法1~4个数量级,且有一定的抗噪能力。
附图说明
图1为本发明流程图。
具体实施方式
下面结合附图和实施例对本发明作进一步说明:
如图1所示,本发明实现的流程图按以下步骤进行:
(1)输入采样数据y(n)并作AR谱估计运算给出精度较低的谐波参数信息;
(2)基于初步谐波信息,通过选取功率谱局部最大值建立参数模型;
(3)应用非线性优化算法确定准确的谐波参数,步骤包括:
1)置初始近似值x0∈Rn,误差限ε1>0;
2)计算令0→k;
3)求最优步长tk满足计算迭代值xk+1=xk+λkpk以及
4)计算若转步骤6),否则,转步骤5);
5)计算k+1→k,转步骤3;
6)置ε2>0,阻尼因子μ0=10-2,缩放因子v>1(可取2,5,10),0→k;
7)计算f(xk),Df(xk),0→j;
8)解方程组求得p(xk)=p(μk);
9)计算xk+1=xk+p(xk)及
10)若则取转步骤7),否则转步骤11);
11)若j=0,则取μk=vμk,1→j,转步骤7),否则转步骤12);
12)若||p(xk)||≤ε2,则xk+1为极小值点x*的近似,算法停止,输出xk+1;否则取xk+1→xk,k+1→k,转步骤7);
算例1:取含基波与相近频率间谐波的电网信号模型如下:
其中,v(t)是信噪比为60dB的高斯白噪声;取采样频率fs=1kHz,采样点数N=512个,AR模型的阶数p=176;在非线性优化法中,取误差限ε1=1e-2,ε2=1e-6,缩放因子v=2,阻尼因子μ0=10-2。
表1谐波与间谐波参数计算精度相对误差表
注:加窗指加Hanning窗插值FFT算法,AR指AR谱估计与非线性优化结合方法。
算例2:取含多个谐波与间谐波的电网信号模型如下:
其中,v(t)是信噪比为60dB的高斯白噪声;取采样频率fs=1kHz,采样点数N=512个,AR模型的阶数p=176;在非线性优化法中,取误差限ε1=1e-2,ε2=1e-6,缩放因子v=2,阻尼因子μ0=10-2。
表2谐波与间谐波参数计算精度相对误差表
注:加窗法指加Hanning窗插值FFT算法,结合法指AR谱估计与非线性优化结合方法。
从表1和2的结果可见,噪声下的谐波和间谐波的检测中,AR谱估计和非线性优化的结合方法能准确地检测出信号的谐波与间谐波分量,且参数计算精度较加Hanning窗插值法提高了约1~4个数量级。
本发明的技术方案可应用于电网谐波分析、电能计量和电能质量监测。
Claims (3)
1.一种结合AR谱估计与非线性优化的间谐波检测方法,其特征在于,按以下步骤进行:
步骤一:AR谱初步估计
1.1):根据谐波参数的初步信息建立该谐波的参数模型,取参数模型如下:
式中y(t)为待分析信号,Am为各次谐波和间谐波的幅值,frm为对应的频率,θm为对应的相位,m为对应的谐波和间谐波个数,t为时间;
1.2):对求解信号模型参数进行非线性优化算法:将所述模型(1)作如下处理,将幅值、相位(Am,θm)等价转化为求解两个幅值(Cm,Dm):
其中,Cm=Amcos(θm),Dm=Amsin(θm);
1.3):针对步骤一中的参数模型(1),将式(2)代入式(1),建立最小二乘方程如下:
式中n为采样点数;目标函数表征了信号计算结果和实际采样结果的误差,令y(tn)为信号第n点的计算值、y1(tn)为信号第n点实际采样值,将参数frm、Cm和Dm统一用向量x表示,对应于各采样时刻的误差向量用f(x)表示:
代入(3)式可得:
式中[f(x)]T表示f(x)的转置,则参数的求解转化为求如下最值的问题:
步骤二:共轭梯度法一次迭代,也就是参数初值进行预先处理,使参数迭代值接近全局最优解:
2.1)对参数初值进行迭代,迭代方程为:
xk+1=xk+λkpk (5)
其中,k为迭代次数,λk为最优步长,pk为迭代方向,
2.2)pk满足下式(7)
置初始近似值x0∈Rn,Rn为实数域,一次迭代误差限ε1>0,
2.3):计算令0→k,p0是迭代方向矩阵pk的初值;
2.4):求最优步长λk满足计算迭代值
2.5):计算当迭代到时,转步骤三,否则,转2.6);
2.6):计算k+1→k,转2.4);
步骤三:当后,再用阻尼最小二乘法二次迭代:即采用阻尼最小二乘法快速寻找到全局最优解,即最终的间谐波参数计算值:
3.1):为使沿着其梯度方向减小,必有的法方程等于零:
可得:
式中:G(xb)=Df(xb)TDf(xb),g(xb)=Df(xb)Tf(xb),p(xb)=-G(xb)-1g(xb);
3.2):考虑阻尼时有:
其中μb为第b次迭代的阻尼因子,I是与G(xb)同维数的单位对角阵;则最终迭代公式为:
当μb>0总可以保证是收敛的,原则上μb=10-4~10-2;
3.3):计算得出令0→j,j为3.7)步骤所需的一个判断变量;
3.4):解方程组求得
3.5):计算xb+1=xb+p(xb)及
3.6):若则取缩放因子v>1,可取2,5,10,转3.3),否则转3.7);
3.7):若j=0,则取μb=vμb,1→j,转3.3),否则转3.8);
3.8):若||p(xb)||≤ε2,则xb+1为极小值点x*的近似,算法停止,输出xb+1;否则取xb+1→xb,b+1→b,转3.3)。
2.根据权利要求1所述一种结合AR谱估计与非线性优化的间谐波检测方法,其特征在于:所述步骤一中的参数模型式(1)中的m次数可根据功率谱中局部最大值个数确定,然后依次确定功率谱中每一个局部最大值处对应的频率和局部最大值的开方,把它们分别作为参数frm、Am的初值,然后用0~2π内的随机值作为相位参数θm的初值。
3.根据权利要求1所述一种结合AR谱估计与非线性优化的间谐波检测方法,其特征在于:所述步骤三中的3.2)阻尼因子μb,当μb>0总可以保证 因而是收敛的,当μb太大时会收敛速度下降,若μb太小则收敛域过小,原则上μb=10-4~10-2,通常取μb=10-2。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310115657.8A CN103245830B (zh) | 2013-04-03 | 2013-04-03 | 一种结合ar谱估计与非线性优化的间谐波检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310115657.8A CN103245830B (zh) | 2013-04-03 | 2013-04-03 | 一种结合ar谱估计与非线性优化的间谐波检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103245830A CN103245830A (zh) | 2013-08-14 |
CN103245830B true CN103245830B (zh) | 2015-08-19 |
Family
ID=48925475
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310115657.8A Active CN103245830B (zh) | 2013-04-03 | 2013-04-03 | 一种结合ar谱估计与非线性优化的间谐波检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103245830B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104655928A (zh) * | 2013-11-21 | 2015-05-27 | 国家电网公司 | 一种电动汽车充电机的输入电压间谐波的检测方法 |
CN105259546A (zh) * | 2015-09-30 | 2016-01-20 | 西安电子科技大学 | 一种基于ar谱分形的海面微弱雷达目标检测方法 |
CN106909736B (zh) * | 2017-02-24 | 2020-07-14 | 北京航天控制仪器研究所 | 一种扰动运动统计特性的解析计算方法 |
CN110161311B (zh) * | 2019-05-17 | 2020-09-18 | 华中科技大学 | 一种谐波与间谐波的检测方法 |
CN111443241A (zh) * | 2020-04-27 | 2020-07-24 | 南京工程学院 | 一种高压岸电电力系统谐波检测方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102222911A (zh) * | 2011-04-19 | 2011-10-19 | 哈尔滨工业大学 | 基于ar模型和卡尔曼滤波的电力系统间谐波估计方法 |
CN102323477A (zh) * | 2011-06-09 | 2012-01-18 | 南京信息工程大学 | 一种电网间谐波测量方法及测量仪 |
CN202119835U (zh) * | 2011-05-20 | 2012-01-18 | 四川电力科学研究院 | 非稳定谐波及间谐波测量仪 |
CN102636693A (zh) * | 2012-05-04 | 2012-08-15 | 重庆大学 | 一种结合fft与非线性最小二乘的谐波分析算法 |
-
2013
- 2013-04-03 CN CN201310115657.8A patent/CN103245830B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102222911A (zh) * | 2011-04-19 | 2011-10-19 | 哈尔滨工业大学 | 基于ar模型和卡尔曼滤波的电力系统间谐波估计方法 |
CN202119835U (zh) * | 2011-05-20 | 2012-01-18 | 四川电力科学研究院 | 非稳定谐波及间谐波测量仪 |
CN102323477A (zh) * | 2011-06-09 | 2012-01-18 | 南京信息工程大学 | 一种电网间谐波测量方法及测量仪 |
CN102636693A (zh) * | 2012-05-04 | 2012-08-15 | 重庆大学 | 一种结合fft与非线性最小二乘的谐波分析算法 |
Also Published As
Publication number | Publication date |
---|---|
CN103245830A (zh) | 2013-08-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102435844B (zh) | 一种频率无关的正弦信号相量计算方法 | |
Chen et al. | Adaptive clutter suppression and detection algorithm for radar maneuvering target with high-order motions via sparse fractional ambiguity function | |
Su et al. | Power harmonic and interharmonic detection method in renewable power based on Nuttall double‐window all‐phase FFT algorithm | |
CN103245830B (zh) | 一种结合ar谱估计与非线性优化的间谐波检测方法 | |
Banerjee et al. | A subspace-based dynamic phasor estimator for synchrophasor application | |
CN105137180B (zh) | 基于六项余弦窗四谱线插值的高精度谐波分析方法 | |
CN110007141B (zh) | 一种基于电压电流谐波相似度的谐振点检测方法 | |
CN103308766A (zh) | 一种基于凯撒自卷积窗双谱线插值fft谐波分析方法及其装置 | |
CN106501602B (zh) | 一种基于滑窗频谱分离的基波参数测量方法 | |
CN103308804A (zh) | 基于快速k-s变换电能质量扰动信号时频参数提取方法 | |
CN107102255A (zh) | 单一adc采集通道动态特性测试方法 | |
CN102636693A (zh) | 一种结合fft与非线性最小二乘的谐波分析算法 | |
CN106324340A (zh) | 一种同步相量和频率测量动态性能的方法 | |
CN103941090A (zh) | 基于谱线能量插值的谐波测量方法 | |
CN107392123A (zh) | 一种基于相参积累消噪的射频指纹特征提取和识别方法 | |
CN105445541A (zh) | 一种任意频率下自适应功率计算方法 | |
CN103543331B (zh) | 一种计算电信号谐波和间谐波的方法 | |
Ma et al. | Harmonic and interharmonic analysis of mixed dense frequency signals | |
Li et al. | Frequency estimation based on modulation FFT and MUSIC algorithm | |
Zhan et al. | Improved WLS-TF algorithm for dynamic synchronized angle and frequency estimation | |
CN105372492B (zh) | 基于三条dft复数谱线的信号频率测量方法 | |
CN105738698B (zh) | 一种基于中心频移的谐波参数估计算法 | |
CN103592514A (zh) | 一种新型谐波高精度检测方法 | |
CN105606893B (zh) | 基于空间平滑修正music的电力间谐波检测方法 | |
CN104950215A (zh) | 一种微机保护方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |