基于小波变换与正弦曲线拟合的非稳态畸变噪声检测方法技术

本发明专利技术属于电能计量领域，具体涉及一种基于小波变换与曲线拟合的非稳态畸变信号噪声检测方法；该方法首先进行非稳态畸变信号小波系数求解，然后利用小波变换方法重构基波信号，在进行基波信号正弦曲线拟合，最后计算非稳态畸变噪声；本发明专利技术基于小波变换与曲线拟合的非稳态畸变信号噪声检测方法不仅能够提高非稳态噪声信号的检测精度，而且能够减小不同层数之间分析结果的差异，具有更广泛的层数适用范围。

【技术实现步骤摘要】

本专利技术属于电能计量领域，具体涉及一种基于小波变换与曲线拟合的非稳态畸变信号噪声检测方法。
技术介绍
随着电网中非线性负载的增加，电网基波信号上会叠加非稳态噪声形成电网非稳态畸变信号。非稳态噪声信号由于具有时变、频率复杂特性，增加了电能计量难度，降低了电能计量的准确性，不仅使国家电力部门利益受到损失，以及电网系统安全受到影响，而且容易烧毁对电网信号变化敏感的负荷设备。解决这些问题需要提出一种能够准确检测非稳态噪声信号的新方法。噪声信号检测方法中,早期采用模拟滤波器原理进行分析，这种方法电路原理简单，但电路元件参数发生变化时，检测结果不理想。20世纪80年代中期出现了基于传统傅里叶变换的数字滤波器，该方法有计算量小、实现简单等优势，但由于不具有频率局部化能力，可能会漏采非稳态噪声。小波变换方法是近些年快速发展起来的一种变换分析新方法，随着小波变换方法的完善和发展，弥补了奈奎斯特采样定理条件限定下传统方法的不足，可以从非稳态畸变信号中分离基波信号与噪声信号，以便检测。而现有大多数小波变换方法研究，都是建立在已知非稳态噪声频率能量分布基础上，选取合适的小波变换方法分解重构层数，才能比较准确的检测非稳态噪声，分解层数过多或过少都可能会降低检测非稳态噪声的准确性。而实际电网非稳态噪声频率复杂并且未知，传统方法可能降低非稳态噪声检测精度，所以就需要改进传统方法准确检测噪声信号。
技术实现思路
针对上述问题，本专利技术公开了一种基于小波变换与正弦曲线拟合的非稳态畸变噪声检测方法，该方法不仅能够提高非稳态噪声信号的检测精度，而且能够减小不同层数之间分析结果的差异，具有更广泛的层数适用范围。本专利技术的目的是这样实现的：基于小波变换与正弦曲线拟合的非稳态畸变噪声检测方法，包括以下步骤：S1、非稳态畸变信号小波系数求解采用小波变换Mallat算法对非稳态畸变信号进行分解，分解算法如下所示，A0[i(t)]=i(t)Aj[i(t)]=ΣkH(2t-k)Aj-1[i(t)]Dj[i(t)]=ΣkG(2t-k)Aj-1[i(t)]]]>式中，t为离散时间序列，i(t)为非稳态畸变信号，j为分解层数，k为位移因子，H和G均为时域小波分析滤波器，Aj为信号i(t)在第j层的低频系数，Dj为信号第j层的高频系数；S2、小波变换方法重构基波信号采用小波变换Mallat算法对非稳态畸变信号中的基波信号进行重构，重构后基波信号ia(t)如下所示，ia(t)＝Aj[i(t)]H(2k-t)式中，t为离散时间序列，i(t)为非稳态畸变信号，j为分解层数，k为位移因子，H和G均为时域小波分析滤波器，Aj为信号i(t)在第j层的低频系数；S3、基波信号正弦曲线拟合按照如下公式拟合基波信号i'a(t)i′a(t)=ω2π∫02π/ωi0(t)dt+ωπ∫t0t0+2π/ωia(t)cos(ωt)dtcos(ωt)+ωπ∫t0t0+2π/ωia(t)sin(ωt)dtsin(ωt)]]>式中，ω为基波信号的角频率，t0表示一个周期的起始时间；S4、计算非稳态畸变噪声按照公式id(t)＝i(t)-i'a(t)计算非稳态畸变噪声id(t)。有益效果：第一、本专利技术基于小波变换与正弦曲线拟合的非稳态畸变噪声检测方法，首先通过小波变换方法计算非稳态畸变信号的高频系数和低频系数，再根据低频系数重构出基波信号，重要的是，在重构得到基波信号的基础上，进行正弦曲线拟合，因此能够提高基波信号拟合结果与真值之间的一致性，进而通过非稳态畸变信号与基波信号拟合结果做差，得到更准确的非稳态噪声信号；第二、由于在重构得到基波信号的基础上，进行了正弦曲线拟合，使得拟合后的基波信号与真值之间的一致性进一步提高，因此降低了不同层数分解下得到结果的差异，使得本专利技术基于小波变换与正弦曲线拟合的非稳态畸变噪声检测方法适用的层数更广泛。附图说明图1是本专利技术基于小波变换与正弦曲线拟合的非稳态畸变噪声检测方法的流程图。图2是标准电网基波信号。图3是电网冲击噪声信号。图4是电网噪声信号与电网基波信号叠加后的非稳态畸变信号。图5是6层分解重构条件下单独小波分析方法得到的基波信号。图6是6层分解重构条件下单独小波分析方法得到的噪声信号。图7是本专利技术方法重构得到的基波信号。图8是本专利技术方法重构得到的噪声信号。具体实施例下面结合附图对本专利技术具体实施例作进一步详细描述。具体实施例一本实施例的基于小波变换与正弦曲线拟合的非稳态畸变噪声检测方法，流程图如图1所示，该方法包括以下步骤：S1、非稳态畸变信号小波系数求解采用小波变换Mallat算法对非稳态畸变信号进行分解，分解算法如下所示，A0[i(t)]=i(t)Aj[i(t)]=ΣkH(2t-k)Aj-1[i(t)]Dj[i(t)]=ΣkG(2t-k)Aj-1[i(t)]]]>式中，t为离散时间序列，i(t)为非稳态畸变信号，j为分解层数，k为位移因子，H和G均为时域小波分析滤波器，Aj为信号i(t)在第j层的低频系数，Dj为信号第j层的高频系数；S2、小波变换方法重构基波信号采用小波变换Mallat算法对非稳态畸变信号中的基波信号进行重构，重构后基波信号ia(t)如下所示，ia(t)＝Aj[i(t)]H(2k-t)式中，t为离散时间序列，i(t)为非稳态畸变信号，j为分解层数，k为位移因子，H和G均为时域小波分析滤波器，Aj为信号i(t)在第j层的低频系数；S3、基波信号正弦曲线拟合按照如下公式拟合基波信号i'a(t)i′a(t)=ω2π∫02π/ωi0(t)dt+ωπ∫t0t0+2π/ωia(t)cos(ωt)dtcos(ωt)+ωπ∫t0t0+2π/ωia(t)sin(ωt)dtsin(ωt)]]>式中，ω为基波信号的角频率，t0表示一个周期的起始时间；S4、计算非稳态畸变噪声按照公式id(t)＝i(t)-i'a(t)计算非稳态畸变噪声id(t)。为了进一步验证本专利技术的有益效果，进行了如下仿真实验：电网冲击信号模型的仿真结果中的电网基波信号如图2所示。电网冲击噪声信号如图3所示，冲击噪声信号发生在0.05s-0.15s之间。噪声信号与基波信号叠加后的非稳态畸变信号如图4所示。以6层分解为例，分别通过本专利技术方法和单独小波变换方法对图中的电网冲击信号模型进行分解重构，单独小波分析方法得到的基波信号与噪声信号分别如图5和图6所示，本专利方法重构得到的基波信号与噪声信号如图本文档来自技高网...

【技术保护点】
基于小波变换与正弦曲线拟合的非稳态畸变噪声检测方法，其特征在于，包括以下步骤：S1、非稳态畸变信号小波系数求解采用小波变换Mallat算法对非稳态畸变信号进行分解，分解算法如下所示，A0[i(t)]=i(t)Aj[i(t)]=ΣkH(2t-k)Aj-1[i(t)]Dj[i(t)]=ΣkG(2t-k)Aj-1[i(t)]]]>式中，t为离散时间序列，i(t)为非稳态畸变信号，j为分解层数，k为位移因子，H和G均为时域小波分析滤波器，Aj为信号i(t)在第j层的低频系数，Dj为信号第j层的高频系数；S2、小波变换方法重构基波信号采用小波变换Mallat算法对非稳态畸变信号中的基波信号进行重构，重构后基波信号ia(t)如下所示，ia(t)＝Aj[i(t)]H(2k‑t)式中，t为离散时间序列，i(t)为非稳态畸变信号，j为分解层数，k为位移因子，H和G均为时域小波分析滤波器，Aj为信号i(t)在第j层的低频系数；S3、基波信号正弦曲线拟合按照如下公式拟合基波信号i'a(t)i′a(t)=ω2π∫02π/ωia(t)dt+ωπ∫t0t0+2π/ωia(t)cos(ωt)dt cos(ωt)+ωπ∫t0t0+2π/ωia(t)sin(ωt)dt sin(ωt)]]>式中，ω为基波信号的角频率，t0表示一个周期的起始时间；S4、计算非稳态畸变噪声按照公式id(t)＝i(t)‑i'a(t)计算非稳态畸变噪声id(t)。...

【技术特征摘要】
1.基于小波变换与正弦曲线拟合的非稳态畸变噪声检测方法，其特征在于，包括以下步骤：S1、非稳态畸变信号小波系数求解采用小波变换Mallat算法对非稳态畸变信号进行分解，分解算法如下所示，A0[i(t)]=i(t)Aj[i(t)]=ΣkH(2t-k)Aj-1[i(t)]Dj[i(t)]=ΣkG(2t-k)Aj-1[i(t)]]]>式中，t为离散时间序列，i(t)为非稳态畸变信号，j为分解层数，k为位移因子，H和G均为时域小波分析滤波器，Aj为信号i(t)在第j层的低频系数，Dj为信号第j层的高频系数；S2、小波变换方法重构基波信号采用小波变换Mallat算法对非稳态畸变信号中的基波信号进行重构，重构后基波信号ia(t)如下所示，ia(t)＝Aj[i(t)]H...

【专利技术属性】
技术研发人员：于晓洋，朱凯鹏，赵烟桥，
申请(专利权)人：哈尔滨理工大学，
类型：发明
国别省市：黑龙江;23