CN112217509B - 基于简化梳状滤波器的主动降噪锁频环装置及方法 - Google Patents
基于简化梳状滤波器的主动降噪锁频环装置及方法 Download PDFInfo
- Publication number
- CN112217509B CN112217509B CN202011050746.5A CN202011050746A CN112217509B CN 112217509 B CN112217509 B CN 112217509B CN 202011050746 A CN202011050746 A CN 202011050746A CN 112217509 B CN112217509 B CN 112217509B
- Authority
- CN
- China
- Prior art keywords
- frequency
- filter
- simplified
- frequency information
- comb filter
- 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
- 230000009467 reduction Effects 0.000 title claims abstract description 50
- 238000000034 method Methods 0.000 title claims abstract description 33
- 238000004364 calculation method Methods 0.000 claims abstract description 13
- 238000001914 filtration Methods 0.000 claims description 58
- 230000003044 adaptive effect Effects 0.000 claims description 56
- 238000010606 normalization Methods 0.000 claims description 18
- 230000009466 transformation Effects 0.000 claims description 12
- 230000003111 delayed effect Effects 0.000 claims description 10
- 238000006243 chemical reaction Methods 0.000 claims description 6
- 230000006872 improvement Effects 0.000 claims description 5
- 238000013016 damping Methods 0.000 claims description 4
- 230000007704 transition Effects 0.000 claims description 3
- 230000008030 elimination Effects 0.000 abstract description 8
- 238000003379 elimination reaction Methods 0.000 abstract description 8
- 238000005516 engineering process Methods 0.000 abstract description 8
- 238000005070 sampling Methods 0.000 abstract description 3
- 230000000694 effects Effects 0.000 description 5
- 238000001514 detection method Methods 0.000 description 4
- 230000004044 response Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 2
- 238000007781 pre-processing Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000001629 suppression Effects 0.000 description 2
- 230000001360 synchronised effect Effects 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000014509 gene expression Effects 0.000 description 1
- 230000010358 mechanical oscillation Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03L—AUTOMATIC CONTROL, STARTING, SYNCHRONISATION OR STABILISATION OF GENERATORS OF ELECTRONIC OSCILLATIONS OR PULSES
- H03L7/00—Automatic control of frequency or phase; Synchronisation
- H03L7/06—Automatic control of frequency or phase; Synchronisation using a reference signal applied to a frequency- or phase-locked loop
- H03L7/08—Details of the phase-locked loop
- H03L7/085—Details of the phase-locked loop concerning mainly the frequency- or phase-detection arrangement including the filtering or amplification of its output signal
- H03L7/093—Details of the phase-locked loop concerning mainly the frequency- or phase-detection arrangement including the filtering or amplification of its output signal using special filtering or amplification characteristics in the loop
-
- 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
- Filters That Use Time-Delay Elements (AREA)
Abstract
本发明涉及一种基于简化梳状滤波器的主动降噪锁频环装置及方法,通过减少陷波点的数量以及利用梳状滤波器中零极点和陷波点的对应关系,通过s域的二次最优平方逼近来简化梳状滤波器,使其复杂度降低并且解决存在近似小延迟和采样频率高的问题;以及通过改进型主动降噪技术实现对剩余噪声的二次消除,以这样的方式在保持原本锁频环动态特性的同时减少原本梳状滤波器所需的计算力并且具有出色的锁频结果。本发明可以被广泛应用于三相双二阶广义积分器锁频环中用来准确锁定电网频率。
Description
技术领域
本发明涉及信号处理领域,特别是涉及一种基于简化梳状滤波器的主动降噪锁频环装置及方法。
背景技术
同步技术在现代能源技术和电力应用中对信号的处理、监测、控制非常重要。电网电压相位和频率信息的准确检测是电力电子变换器并网稳定运行的关键。锁相环和锁频环用于检测电网信号的相位和频率,进而实现电网信号的同步,精确且快速地估算出电网电压的相位和频率信息是锁相环和锁频环的目标。
同步旋转坐标系下的锁相环作为最常用的锁相技术可在电网正常工作情况下获得准确的锁相结果,然而,在电网电压畸变时,传统锁相环动态响应慢,锁相误差较大,严重时三相同步锁相环(Synchronous Reference Frame Phase-Locked Loop,SRF-PLL)甚至无法工作。因此,本领域技术人员提出了各种各样的改进型锁相环和锁频环。这些改进方法本质是通过在锁定频率之前使用滤波环节对输入给锁相环或锁频环的信号进行“预处理”,通过给锁相环或锁频环更加“纯净”的信号来得到准确的电网信息。
虽然这些方案对于谐波的抑制能力出色,但这些“预处理”网络的设置,也带来了巨大的运算力需求以及较慢的动态响应问题,而这些问题是阻碍这些方案实施的重要因素。
电网发生暂升、暂降、三相不平衡和谐波污染严重时,受负序分量影响,传统锁相环的相位误差增大,不能满足并网逆变器的控制需求。为提取基波分量,需在传统锁相环之前加装低通滤波器或对结构进行改进,但仍易引起相角偏移、响应速度慢等缺点。
由于近些年在小信号建模和锁频环滤波能力方面取得的进展,使得锁频环在电力系统相关应用方面越来越受欢迎。现在,除了在功率变换器同步信号方面的应用外,锁频环还用于例如频率自适应谐振控制器、电力系统中的孤岛检测、机械振荡估计、估计两个或多个频率未知的正弦信号、检测并分离功率信号中的基波和谐波成分、抑制心电图信号中电源线干扰等等中。
锁频环用于锁定电网电压的频率而不是相位,相位计算在回路外完成。锁频环和锁相环之间的固有差异,使锁频环在电网基波检测中更具优势。锁频环通常需要利用带通滤波器的特性来提取基波成分,并衰减、消除锁频环系统中的其他谐波,因此在锁相环中广泛使用的环路内低通滤波器不在适用于锁频环。
而现有的用于锁频环的预滤波技术只能对特定的频率或者特定的某些频率进行噪声消除,因此,本领域亟需一种统一将所有噪声进行消除的技术方案。
发明内容
本发明的目的是提供一种基于简化梳状滤波器的主动降噪锁频环装置及方法,解决当前只能对特定的噪声进行消除,而且会带来比较大的计算负担的问题。
为实现上述目的,本发明提供了如下方案:
一种基于简化梳状滤波器的主动降噪锁频环装置,其特征在于,所述装置包括:简化梳状滤波器电路、双二阶广义积分器锁频环电路和主动降噪电路;
所述简化梳状滤波器电路通过减少陷波点的数量以及利用梳状滤波器中零极点和陷波点的对应关系,通过s域的二次最优平方逼近来简化;
所述双二阶广义积分器锁频环电路包括第一二阶广义积分器、第二二阶广义积分器、归一化单元和频率误差信号的运算环节;
所述主动降噪电路包括延时器、自适应滤波器和改进型基于sigmoid函数的最小均方算法;
采集到的电压信号经dq变换后输入至所述简化梳状滤波器,得到梳状滤波信号,再将所述梳状滤波信号转换至αβ坐标系下,之后输入双二阶广义积分器锁频环进行运算得到所述电压信号的频率信息,再将所述频率信息延时处理得到延时频率信息,再将所述延时频率信息分别输入自适应滤波器和改进型基于sigmoid函数的最小均方算法得到自适应滤波器预测值和改进型基于sigmoid函数的最小均方算法的输出值,从而确定自适应滤波器系数,再将所述延时频率信息输入自适应滤波器,得到最终频率信息。
可选的,所述简化梳状滤波器电路通过调整带阻滤波器的阻带减小陷波点的数量,同时利用梳状滤波器中零极点和陷波点的对应关系,通过s域的二次最优平方逼近来简化;
λ1和λ2是由梳状滤波器零极点分布的函数进行二次最优平方逼近后得到的参数,是原梳状滤波器多个零极点进行简化后的数值;其中λ1为阻带中心点,不同的λ1会导致不同的阻带中心;λ2为阻尼系数,不同的值会造成不同的阻带深度和过渡带的陡度。
可选的,所述主动降噪电路通过修改误差与滤波器之间的函数关系,使得滤波器系数更能反映误差的变换;
其中μ(n)是自适应滤波器系数,误差e(n)是实际观测值x(n)和预测值y(n)之间的偏差,自适应滤波器的系数由e(n)和x(n-d)更新,α是控制μ(n)宽窄形状的常数,β是控制μ(n)范围的常数,n为对应第n次的计算。
一种基于简化梳状滤波器的主动降噪锁频环方法,所述方法包括:
获取电压信号;
利用简化梳状滤波器对所述电压信号进行滤波,得到梳状滤波信号;所述简化梳状滤波器利用梳状滤波器中零极点和陷波点的对应关系,通过s域的二次最优平方逼近来简化;
将所述梳状滤波信号转换至αβ坐标系下,然后由双二阶广义积分器锁频环进行运算得到所述电压信号的频率信息;
将所述频率信息延时处理得到延时频率信息,并保留延时处理之前所述频率信息的观测值;
将所述延时频率信息分别输入自适应滤波器和改进型基于sigmoid函数的最小均方算法,得到自适应滤波器预测值和改进型基于sigmoid函数的最小均方算法的输出值;
根据所述改进型基于sigmoid函数的最小均方算法的输出值和所述自适应滤波器预测值与所述频率信息的观测值之间的误差确定自适应滤波器系数;
利用确定自适应滤波器系数后的自适应滤波器对所述延时频率信息进行滤波,得到最终频率信息。
可选的,所述利用简化梳状滤波器对所述电压信号进行滤波,得到梳状滤波信号之前还包括:
将所述电压信号由abc坐标经进行dq变换后得到dq坐标系下的dq电压信号。
可选的,所述由双二阶广义积分器锁频环进行运算得到所述电压信号的频率信息具体包括:
将转换至αβ坐标系下的信号第一电压信号和第二电压信号分别输入第一二阶广义积分器和第二二阶广义积分器;
将所述第一二阶广义积分器和第二二阶广义积分器的输出连接至归一化单元的输入得到归一化单元的输出值;
将所述第一二阶广义积分器和第二二阶广义积分器的输出连接至频率误差信号计算的输入端进行运算得到频率误差信号;
将所述归一化单元的输出值和所述频率误差信号相乘得出乘积值;
当所述乘积值为0时,输出所述双二阶广义积分器锁频环锁定的频率信息,所述锁定的频率信息与所述电压信号的频率信息相等。
其中μ(n)是自适应滤波器系数,误差e(n)是实际观测值x(n)和预测值y(n)之间的偏差,自适应滤波器的系数由e(n)和x(n-d)更新。
一种基于简化梳状滤波器的主动降噪锁频环系统,所述系统包括:
电压信号获取单元,用于获取电压信号;
梳状滤波单元,用于利用简化梳状滤波器对所述电压信号进行滤波,得到梳状滤波信号;所述简化梳状滤波器利用梳状滤波器中零极点和陷波点的对应关系,通过s域的二次最优平方逼近来简化;
锁频单元,用于将所述梳状滤波信号转换至αβ坐标系下,然后由双二阶广义积分器锁频环进行运算得到所述电压信号的频率信息;
延时单元,用于将所述频率信息延时处理得到延时频率信息,并保留延时处理之前所述频率信息的观测值;
延时频率信息输入单元,用于将所述延时频率信息分别输入自适应滤波器和改进型基于sigmoid函数的最小均方算法,得到自适应滤波器预测值和改进型基于sigmoid函数的最小均方算法的输出值;
自适应滤波器系数确定单元,用于根据所述改进型基于sigmoid函数的最小均方算法的输出值和所述自适应滤波器预测值与所述频率信息的观测值之间的误差确定自适应滤波器系数;
自适应滤波单元,用于利用确定自适应滤波器系数后的自适应滤波器对所述延时频率信息进行滤波,得到最终频率信息。
可选的,还包括dq变换单元,用于将所述电压信号由abc坐标经进行dq变换后得到dq坐标系下的dq电压信号。
可选的,还包括αβ坐标转换单元,用于将所述梳状滤波信号转换至αβ坐标系下。
根据本发明提供的具体实施例,本发明公开了以下技术效果:
本发明用一种新的思路简化了梳状滤波器,减少了原本梳状滤波器在实际实现中所存在近似小延迟、实现复杂度高、采样频率高的问题,同时完成了对主要谐波的消除。
本发明针对简化梳状滤波器所导致的滤波效果不足,以及可能存在的意料之外的噪声,增加了主动降噪。相较于之前的梳状滤波器或者其他的用于锁频环的预滤波技术只能对特定的频率或者特定的某些频率进行噪声消除,其不仅只能对特定的噪声进行消除,而且会带来比较大的计算负担。本发明通过针对双二阶广义积分器锁频环得到的频率信息使用主动降噪,统一的将所有噪声进行消除,从而实现了频率的确定。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例一提供的基于简化梳状滤波器的主动降噪锁频环装置的结构示意图。
符号说明:31为简化梳状滤波器,11为第一二阶广义积分器,12为第二二阶广义积分器,21为自适应滤波器,22为改进型基于sigmoid函数的最小均方算法,1为双二阶广义积分器锁频环,2为主动降噪部分,eα和eβ为稳态时α轴、β轴上的频率误差信号,να,νβ为α,β轴的输入信号,ef为频率误差的信号,ω′为二阶广义积分器的谐振频率,ω0=100π,Г为锁频环环路增益;q为90°相位滞后因子,qv′即表示该信号的相位滞后信号v'的角度为90°,να',qνα'为να通过二阶广义积分器后得到的幅值相等、相位相差90°的正交输出信号,νβ',qνβ'为νβ通过二阶广义积分器后得到的幅值相等、相位相差90°的正交输出信号,ν+和ν-分别为电压正、负序分量幅值。k为阻尼系数,
图2为本发明实施例一提供的基于简化梳状滤波器的主动降噪锁频环装置的双二阶广义积分器锁频环装置中的二阶广义积分器的放大图。
图3为本发明实施例二提供的基于简化梳状滤波器的主动降噪锁频环方法的控制流程图。
图4为本发明实施例二提供的基于简化梳状滤波器的主动降噪锁频环系统的组成示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明的目的是提供一种基于简化梳状滤波器的主动降噪锁频环装置及方法解决当前只能对特定的噪声进行消除,而且会带来比较大的计算负担的问题。
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。
实施例一:
如图1所示,一种基于简化梳状滤波器的主动降噪锁频环装置,包括:简化梳状滤波器电路、双二阶广义积分器锁频环电路和主动降噪电路;
所述简化梳状滤波器电路通过减少陷波点的数量以及利用梳状滤波器中零极点和陷波点的对应关系,通过s域的二次最优平方逼近来简化;
传统的数字梳状滤波器在Z域中具有以下传递函数:
式中,M是滤波器的阶数,z是在复数平面上定义的一个复变量,通常称为z变换算子。为了使梳状滤波器更有选择性或更接近期望的行为,将极点从原点移向单位圆:
其中参数r为梳状滤波器在陷波频率处的衰减系数,且0≤r<1,Pk是指第k个极点,Zk指第k个零点,极坐标系下(1,0)为k=1。通过极点位置的改变,获得了在通频带上更为平坦的增益和陷波处更好的衰减。通过在Z域上的传递函数,从而得到了改进后的梳状滤波器在s域的传递函
其中Tw表示梳状滤波器的窗口长度,fn表示梳状滤波器的陷波频率fn=n/Tw(n∈{1,2,...,M})。
为了在选定的陷波频率上提供较大的衰减,同时在通频带上获得平坦的增益而不影响信号的相位,本发明选择r为接近但小于1。基波经dq变换后可成为直流分量,所以梳状滤波器只用于取消一组位于1∠0°线上的零极点对就可完成滤波功能。
下面简要说明简化梳状滤波器的原理:
梳状滤波器与串联的多个陷波滤波器的幅值响应相当。每一个陷波滤波器是一个阻带很小的带阻滤波器(BSF)。通过在主要谐波含量区域扩大阻带,在其他区域减小或者取消阻带,从而通过调整带阻滤波器的阻带减少陷波点的数量,进而实现梳状滤波器的简化。简化梳状滤波器基于这样的思想,同时利用梳状滤波器中零极点和陷波点的对应关系,通过s域的二次最优平方逼近来简化。
λ1和λ2是由梳状滤波器零极点分布的函数进行二次最优平方逼近后得到的参数,是原梳状滤波器多个零极点进行简化后的数值;其中λ1为阻带中心点,不同的λ1会导致不同的阻带中心;λ2为阻尼系数,不同的值会造成不同的阻带深度和过渡带的陡度。
简化梳状滤波器是为了尽可能的保留原本梳状滤波器的滤波能力的同时,尽可能的降低原本梳状滤波器所需要的运算力和硬件需求。
所述双二阶广义积分器锁频环电路包括第一二阶广义积分器、第二二阶广义积分器、归一化单元和频率误差信号的运算环节;
所述第一二阶广义积分器和第二二阶广义积分器的具体结构见附图2;
双二阶广义积分器锁频环根据给定的三相正弦输入信号,可以得到输入三相正弦信号的频率。其锁定频率的基本原理是:首先将abc坐标系下三相正弦输入信号转换至αβ坐标系下得到να,νβ,然后分别经过二阶广义积分器和归一化单元的运算后,将频率误差信号ef和归一化单元的数值相乘。当正弦输入信号的频率ω小于双二阶广义积分器锁频环的输出频率时,频率误差信号ef和归一化单元的数值乘积小于0,当正弦输入信号的频率ω大于双二阶广义积分器锁频环的输出频率时,频率误差信号ef和归一化单元的数值乘积大于0,当正弦输入信号的频率ω等于双二阶广义积分器锁频环的输出频率时,频率误差信号ef和归一化单元的数值乘积等于0,双二阶广义积分器锁频环中的积分器根据频率误差信号ef和归一化单元的数值乘积的变化,逐渐调整锁频环输出频率,最终使得输入频率和双二阶广义积分器锁频环的输出频率相等,实现频率的锁定。
所述主动降噪电路包括延时器、自适应滤波器和改进型基于sigmoid函数的最小均方算法;
经由简化梳状滤波器滤波后的信号的频率信息,引入改进型的主动降噪对得到的频率信号进行处理,
为改进处,即为改进型基于sigmoid函数的最小均方算法(M-SVSLMS),通过修改误差与滤波器系数之间的函数关系,使得滤波器系数更能反映误差的变换。通过主动降噪电路将所得到的频率信号中的谐波影响消除掉,从而得到准确的电网频率信息。
公式(4)。下面是主动降噪电路的数学表达式(1-4)。
y(n)=wT(n)x(n) (1)
e(n)=x(n)-y(n) (2)
w(n+1)=w(n)+2μ(n)e(n)x(n) (3)
其中,公式(4)为上述提及到的改进处,主动降噪系统中,动态调整自适应滤波器系数和误差之间的关系从而保证快速性和最佳滤波效果。如何动态调整μ(n)和e(n)之间的关系则是由公式(4)指明的。x(n)是观察到的频率信号,ω(n)为权向量,ωT(n)为权向量的转置,α是控制(4)形状的常数,其控制(4)的宽窄形状,α越大(4)越宽大,μ(n)在不同误差作用下变化的更缓慢,β是控制(4)范围的常数,即(4)的高低,β越大(4)越高,最大值越大,n为对应第n次的计算,d为延时器的延时常数,μ(n)是自适应滤波器系数,输入频率信号通过延迟器将得到x(n-d)作为自适应滤波器的输入,输出y(n)。y(n)是自适应滤波器从过去的信号预测的当前值,误差e(n)是实际观测值x(n)和预测值y(n)之间的偏差,自适应滤波器的系数由e(n)和x(n-d)更新。通过主动降噪部分的作用,从而得到最终的频率信息。其中μ(n)是自适应滤波器系数,误差e(n)是实际观测值x(n)和预测值y(n)之间的偏差,自适应滤波器的系数由e(n)和x(n-d)更新。
采集到的电压信号经dq变换后输入至所述简化梳状滤波器,得到梳状滤波信号,再将所述梳状滤波信号转换至αβ坐标系下,之后输入双二阶广义积分器锁频环进行运算得到所述电压信号的频率信息,再将所述频率信息延时处理得到延时频率信息,再将所述延时频率信息分别输入自适应滤波器和改进型基于sigmoid函数的最小均方算法得到自适应滤波器预测值和改进型基于sigmoid函数的最小均方算法的输出值,从而确定自适应滤波器系数,再将所述延时频率信息输入自适应滤波器,得到最终频率信息。
本发明实施例提供的基于简化梳状滤波器的主动降噪锁频环装置用一种新的思路简化了梳状滤波器,减少了原本梳状滤波器在实际实现中所存在近似小延迟、实现复杂度高、采样频率高的问题,同时完成了对主要谐波的消除。并且针对简化梳状滤波器所导致的滤波效果不足,以及可能存在的意料之外的噪声,增加了主动降噪。相较于之前的梳状滤波器或者其他的用于锁频环的预滤波技术只能对特定的频率或者特定的某些频率进行噪声消除,其不仅只能对特定的噪声进行消除,而且会带来比较大的计算负担。本发明实施例通过针对双二阶广义积分器锁频环得到的频率信息使用主动降噪,统一的将所有噪声进行消除,从而实现了频率的确定。
实施例二:
如图2所示,一种基于简化梳状滤波器的主动降噪锁频环方法,其特征在于,所述方法包括:
S1、获取电压信号;Vabc为获取到的电压信号。
将所述电压信号由abc坐标经进行dq变换后得到dq坐标系下的dq电压信号。
S2、利用简化梳状滤波器对所述电压信号进行滤波,得到梳状滤波信号,从而提升了基波在信号中的占比成分;所述简化梳状滤波器利用梳状滤波器中零极点和陷波点的对应关系,通过s域的二次最优平方逼近来简化;
S3、将所述梳状滤波信号转换至αβ坐标系下,然后由双二阶广义积分器锁频环进行运算得到所述电压信号的频率信息;
其具体过程为:
将所述梳状滤波信号转换至αβ坐标系下,得到Vα、Vβ,将Vα、Vβ分别输入第一二阶广义积分器和第二二阶广义积分器;
分别采集两个二阶广义积分器的输出να',qνα',νβ'和qνβ';
同时,将采集到的να与να'相加后与qνα'相乘得到α轴上的频率误差信号eα,νβ与νβ'相加后与qνβ'相乘得到β轴上的频率误差信号eβ,eα和eβ相加后得到误差频率信号ef,误差频率信号依次与锁频环的环路增益、归一化单元的输出值相乘后,将得到的值送给积分器,积分器的输出加上ω0得到的值,得到所述电压信号的频率信息。
S4、将所述频率信息延时处理得到延时频率信息,并保留延时处理之前所述频率信息的观测值;
S5、将所述延时频率信息分别输入自适应滤波器和改进型基于sigmoid函数的最小均方算法,得到自适应滤波器预测值和改进型基于sigmoid函数的最小均方算法的输出值;
S6、根据所述改进型基于sigmoid函数的最小均方算法的输出值和所述自适应滤波器预测值与所述频率信息的观测值之间的误差确定自适应滤波器系数;
S7、利用确定自适应滤波器系数后的自适应滤波器对所述延时频率信息进行滤波,得到最终频率信息。通过对梳状滤波器的简化,对基于sigmoid函数的最小均方算法的改进以及加入改进型的主动降噪部分实现对噪声的双次消除,以这样的方式在保持原本锁频环动态特性的同时减少原本梳状滤波器所需的计算力并且具有出色的锁频结果。
针对上述方法,本发明还提供了一种基于简化梳状滤波器的主动降噪锁频环系统,如图4所示,包括:
电压信号获取单元,用于获取电压信号;
梳状滤波单元,用于利用简化梳状滤波器对所述电压信号进行滤波,得到梳状滤波信号;所述简化梳状滤波器利用梳状滤波器中零极点和陷波点的对应关系,通过s域的二次最优平方逼近来简化;
锁频单元,用于将所述梳状滤波信号转换至αβ坐标系下,然后由双二阶广义积分器锁频环进行运算得到所述电压信号的频率信息;
延时单元,用于将所述频率信息延时处理得到延时频率信息,并保留延时处理之前所述频率信息的观测值;
延时频率信息输入单元,用于将所述延时频率信息分别输入自适应滤波器和改进型基于sigmoid函数的最小均方算法,得到自适应滤波器预测值和改进型基于sigmoid函数的最小均方算法的输出值;
自适应滤波器系数确定单元,用于根据所述改进型基于sigmoid函数的最小均方算法的输出值和所述自适应滤波器预测值与所述频率信息的观测值之间的误差确定自适应滤波器系数;
自适应滤波单元,用于利用确定自适应滤波器系数后的自适应滤波器对所述延时频率信息进行滤波,得到最终频率信息。
作为一种可选的实施方式,所述系统还包括dq变换单元,用于将所述电压信号由abc坐标经进行dq变换后得到dq坐标系下的dq电压信号;
还包括αβ坐标转换单元,用于将所述梳状滤波信号转换至αβ坐标系下。
本发明实施例提供的一种基于简化梳状滤波器的主动降噪锁频环方法和系统涉及一种简化梳状滤波器与主动降噪技术被应用于三相双二阶广义积分器锁频环中用来准确锁定电网频率。首先通过简化梳状滤波器对信号噪声的初级滤除,然后再通过改进型主动降噪技术实现对剩余噪声的二次消除,以这样的方式在保持原本锁频环动态特性的同时减少原本梳状滤波器所需的计算力并且具有出色的锁频结果。
相较于之前的梳状滤波器或者其他的用于锁频环的预滤波技术只能对特定的频率或者特定的某些频率进行噪声消除,本发明实施例通过针对双二阶广义积分器锁频环得到的频率信息使用主动降噪,统一的将所有噪声进行消除,从而实现了频率的确定。
本说明书中各个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。对于实施例公开的系统而言,由于其与实施例公开的方法相对应,所以描述的比较简单,相关之处参见方法部分说明即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处。综上所述,本说明书内容不应理解为对本发明的限制。
Claims (10)
1.一种基于简化梳状滤波器的主动降噪锁频环装置,其特征在于,所述装置包括:简化梳状滤波器电路、双二阶广义积分器锁频环电路和主动降噪电路;
所述简化梳状滤波器电路通过减少陷波点的数量以及利用梳状滤波器中零极点和陷波点的对应关系,通过s域的二次最优平方逼近来简化;
所述双二阶广义积分器锁频环电路包括第一二阶广义积分器、第二二阶广义积分器、归一化单元和频率误差信号的运算环节;
所述主动降噪电路包括延时器、自适应滤波器和改进型基于sigmoid函数的最小均方算法;
采集到的电压信号经dq变换后输入至所述简化梳状滤波器,得到梳状滤波信号,再将所述梳状滤波信号转换至αβ坐标系下,之后输入双二阶广义积分器锁频环进行运算得到所述电压信号的频率信息,再将所述频率信息延时处理得到延时频率信息,再将所述延时频率信息分别输入自适应滤波器和改进型基于sigmoid函数的最小均方算法得到自适应滤波器预测值和改进型基于sigmoid函数的最小均方算法的输出值,从而确定自适应滤波器系数,再将所述延时频率信息输入自适应滤波器,得到最终频率信息。
4.一种基于简化梳状滤波器的主动降噪锁频环方法,其特征在于,所述方法包括:
获取电压信号;
利用简化梳状滤波器对所述电压信号进行滤波,得到梳状滤波信号;所述简化梳状滤波器利用梳状滤波器中零极点和陷波点的对应关系,通过s域的二次最优平方逼近来简化;
将所述梳状滤波信号转换至αβ坐标系下,然后由双二阶广义积分器锁频环进行运算得到所述电压信号的频率信息;
将所述频率信息延时处理得到延时频率信息,并保留延时处理之前所述频率信息的观测值;
将所述延时频率信息分别输入自适应滤波器和改进型基于sigmoid函数的最小均方算法,得到自适应滤波器预测值和改进型基于sigmoid函数的最小均方算法的输出值;
根据所述改进型基于sigmoid函数的最小均方算法的输出值和所述自适应滤波器预测值与所述频率信息的观测值之间的误差确定自适应滤波器系数;
利用确定自适应滤波器系数后的自适应滤波器对所述延时频率信息进行滤波,得到最终频率信息。
5.根据权利要求4所述的一种基于简化梳状滤波器的主动降噪锁频环方法,其特征在于,所述利用简化梳状滤波器对所述电压信号进行滤波,得到梳状滤波信号之前还包括:
将所述电压信号由abc坐标经进行dq变换后得到dq坐标系下的dq电压信号。
6.根据权利要求4所述的一种基于简化梳状滤波器的主动降噪锁频环方法,其特征在于,所述由双二阶广义积分器锁频环进行运算得到所述电压信号的频率信息具体包括:
将转换至αβ坐标系下的信号第一电压信号和第二电压信号分别输入第一二阶广义积分器和第二二阶广义积分器;
将所述第一二阶广义积分器和第二二阶广义积分器的输出连接至归一化单元的输入得到归一化单元的输出值;
将所述第一二阶广义积分器和第二二阶广义积分器的输出连接至频率误差信号计算的输入端进行运算得到频率误差信号;
将所述频率误差信号依次与锁频环环路增益、所述归一化单元的输出值相乘得出乘积值;
将所述乘积值送给积分器;
将所述积分器的输出加上ω0得到的频率信息与所述电压信号的频率信息相等,其中ω0=100π。
8.一种基于简化梳状滤波器的主动降噪锁频环系统,其特征在于,所述系统包括:
电压信号获取单元,用于获取电压信号;
梳状滤波单元,用于利用简化梳状滤波器对所述电压信号进行滤波,得到梳状滤波信号;所述简化梳状滤波器利用梳状滤波器中零极点和陷波点的对应关系,通过s域的二次最优平方逼近来简化;
锁频单元,用于将所述梳状滤波信号转换至αβ坐标系下,然后由双二阶广义积分器锁频环进行运算得到所述电压信号的频率信息;
延时单元,用于将所述频率信息延时处理得到延时频率信息,并保留延时处理之前所述频率信息的观测值;
延时频率信息输入单元,用于将所述延时频率信息分别输入自适应滤波器和改进型基于sigmoid函数的最小均方算法,得到自适应滤波器预测值和改进型基于sigmoid函数的最小均方算法的输出值;
自适应滤波器系数确定单元,用于根据所述改进型基于sigmoid函数的最小均方算法的输出值和所述自适应滤波器预测值与所述频率信息的观测值之间的误差确定自适应滤波器系数;
自适应滤波单元,用于利用确定自适应滤波器系数后的自适应滤波器对所述延时频率信息进行滤波,得到最终频率信息。
9.根据权利要求8所述的一种基于简化梳状滤波器的主动降噪锁频环系统,其特征在于,
还包括dq变换单元,用于将所述电压信号由abc坐标经进行dq变换后得到dq坐标系下的dq电压信号。
10.根据权利要求8所述的一种基于简化梳状滤波器的主动降噪锁频环系统,其特征在于,
还包括αβ坐标转换单元,用于将所述梳状滤波信号转换至αβ坐标系下。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011050746.5A CN112217509B (zh) | 2020-09-29 | 2020-09-29 | 基于简化梳状滤波器的主动降噪锁频环装置及方法 |
AU2021100247A AU2021100247A4 (en) | 2020-09-29 | 2021-01-15 | Frequency-locked loop (fll) device and method for active noise control (anc) based on simplified comb filter |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011050746.5A CN112217509B (zh) | 2020-09-29 | 2020-09-29 | 基于简化梳状滤波器的主动降噪锁频环装置及方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112217509A CN112217509A (zh) | 2021-01-12 |
CN112217509B true CN112217509B (zh) | 2022-11-22 |
Family
ID=74052433
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011050746.5A Active CN112217509B (zh) | 2020-09-29 | 2020-09-29 | 基于简化梳状滤波器的主动降噪锁频环装置及方法 |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN112217509B (zh) |
AU (1) | AU2021100247A4 (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114235143B (zh) * | 2021-12-16 | 2023-09-22 | 常州大学 | 一种磁悬浮多跨转子不对中振动的自适应检测系统 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103199779A (zh) * | 2013-04-22 | 2013-07-10 | 哈尔滨工业大学 | 基于自适应滤波的内置式永磁同步电机转子位置观测装置及观测方法 |
CN106953634A (zh) * | 2017-03-03 | 2017-07-14 | 燕山大学 | 一种基于双自调谐二阶广义积分器的锁频环方法 |
CN107196329A (zh) * | 2017-05-12 | 2017-09-22 | 上海电力学院 | 一种电气化铁路电能治理调节装置的并网锁相方法 |
CN107706929A (zh) * | 2016-07-26 | 2018-02-16 | 南京工程学院 | 基于最小方差滤波的自适应锁相环方法及系统 |
CN107786201A (zh) * | 2017-09-30 | 2018-03-09 | 中国农业大学 | 一种基于锁频环的二阶广义积分器结构及锁相环同步方法 |
CN107846218A (zh) * | 2017-12-08 | 2018-03-27 | 广东电网有限责任公司电力科学研究院 | 一种基于梳状滤波器的锁相环 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6977827B2 (en) * | 2004-03-22 | 2005-12-20 | American Superconductor Corporation | Power system having a phase locked loop with a notch filter |
US7425874B2 (en) * | 2006-06-30 | 2008-09-16 | Texas Instruments Incorporated | All-digital phase-locked loop for a digital pulse-width modulator |
US10386396B2 (en) * | 2015-10-06 | 2019-08-20 | Sparq Systems Inc. | Hybrid PLL for grid synchronization in distributed generation |
-
2020
- 2020-09-29 CN CN202011050746.5A patent/CN112217509B/zh active Active
-
2021
- 2021-01-15 AU AU2021100247A patent/AU2021100247A4/en not_active Ceased
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103199779A (zh) * | 2013-04-22 | 2013-07-10 | 哈尔滨工业大学 | 基于自适应滤波的内置式永磁同步电机转子位置观测装置及观测方法 |
CN107706929A (zh) * | 2016-07-26 | 2018-02-16 | 南京工程学院 | 基于最小方差滤波的自适应锁相环方法及系统 |
CN106953634A (zh) * | 2017-03-03 | 2017-07-14 | 燕山大学 | 一种基于双自调谐二阶广义积分器的锁频环方法 |
CN107196329A (zh) * | 2017-05-12 | 2017-09-22 | 上海电力学院 | 一种电气化铁路电能治理调节装置的并网锁相方法 |
CN107786201A (zh) * | 2017-09-30 | 2018-03-09 | 中国农业大学 | 一种基于锁频环的二阶广义积分器结构及锁相环同步方法 |
CN107846218A (zh) * | 2017-12-08 | 2018-03-27 | 广东电网有限责任公司电力科学研究院 | 一种基于梳状滤波器的锁相环 |
Non-Patent Citations (3)
Title |
---|
Single-Phase Frequency-Locked Loops: A Comprehensive Review;Saeed Golestan等;《IEEE》;20190411;全文 * |
一种基于自适应滤波器的新型单相锁相环技术;冬雷等;《电工技术学报》;20161226;全文 * |
抑制直流偏移的改进型二阶广义积分器的锁频环设计;胡挜喆等;《电测与仪表》;20200510;全文 * |
Also Published As
Publication number | Publication date |
---|---|
AU2021100247A4 (en) | 2021-04-15 |
CN112217509A (zh) | 2021-01-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104578172B (zh) | 一种带线性锁频环的光伏逆变调节器控制方法 | |
CN109193707B (zh) | 基于虚拟同步发电机的负序电压补偿双环控制方法及系统 | |
CN109831184B (zh) | 一种带通滤波方法、系统、存储介质及终端设备 | |
CN109245103B (zh) | 基于改进滑动平均值滤波器的软件锁相环实现方法及装置 | |
CN105449718B (zh) | 基于改进型串联信号延迟对消算法的并网同步锁相方法 | |
CN110943632B (zh) | 基于级联广义积分器的储能变换器虚拟惯量控制方法 | |
CN112217509B (zh) | 基于简化梳状滤波器的主动降噪锁频环装置及方法 | |
CN109406854B (zh) | 单相有源滤波器的谐波指令电流获取方法、装置及设备 | |
CN109358228B (zh) | 基于双增强型锁相环的电网电压正负序分量实时估计方法 | |
CN111693813A (zh) | 一种基于复合广义积分的fbd谐波检测方法 | |
CN109193793B (zh) | 一种变流器免电压检测的并网控制系统和方法 | |
CN110661263B (zh) | 含有自适应延时滤波器的锁频环及基于该锁频环的并网逆变器控制方法 | |
CN113014250A (zh) | 一种可消除直流偏移电压的锁相环及其锁相控制方法 | |
CN110365018B (zh) | 一种自适应宽频带正负序分离方法 | |
CN115912489A (zh) | 一种适用于非理想电网的lms-sogi三相锁相环设计方法及系统 | |
CN116112008A (zh) | 一种改进型dsogi-pll锁相环 | |
CN109599867A (zh) | 基于双环结构的三相并网变流器电流环干扰抑制控制方法 | |
CN114629112A (zh) | 基于二阶广义积分器的锁频环及其控制方法 | |
CN110702987B (zh) | 电网电压信号的正负序基波分量提取系统 | |
Sozański | Three phase active power filter with selective harmonics elimination | |
CN113346899A (zh) | 一种基于级联型滤波器的三相并网软件锁相环 | |
CN107798162B (zh) | 一种电网正序电压分量与负序电压分量估计方法及装置 | |
CN112332809A (zh) | 一种幅度非衰减均衡相位的分数阶椭圆滤波器设计方法 | |
CN109617077B (zh) | 一种全数字化电网同步锁相方法 | |
CN113472346B (zh) | 一种基于复合型滤波器的电网同步软件锁相环 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |