CN114784810A - 一种自适应频率估计的锁相环以及锁相方法 - Google Patents

一种自适应频率估计的锁相环以及锁相方法 Download PDF

Info

Publication number
CN114784810A
CN114784810A CN202210683832.2A CN202210683832A CN114784810A CN 114784810 A CN114784810 A CN 114784810A CN 202210683832 A CN202210683832 A CN 202210683832A CN 114784810 A CN114784810 A CN 114784810A
Authority
CN
China
Prior art keywords
phase
adaptive
component
notch filter
disturbance
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
Application number
CN202210683832.2A
Other languages
English (en)
Other versions
CN114784810B (zh
Inventor
吴亚楠
谭亮
傅鹏
黄懿赟
卢晶
茆智伟
李俊
茆华风
田贇祥
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hefei Institutes of Physical Science of CAS
Original Assignee
Hefei Institutes of Physical Science of CAS
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Hefei Institutes of Physical Science of CAS filed Critical Hefei Institutes of Physical Science of CAS
Priority to CN202210683832.2A priority Critical patent/CN114784810B/zh
Publication of CN114784810A publication Critical patent/CN114784810A/zh
Application granted granted Critical
Publication of CN114784810B publication Critical patent/CN114784810B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/12—Circuit arrangements for ac mains or ac distribution networks for adjusting voltage in ac networks by changing a characteristic of the network load
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R25/00—Arrangements for measuring phase angle between a voltage and a current or between voltages or currents
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/30—Circuit design
    • G06F30/32—Circuit design at the digital level
    • G06F30/337—Design optimisation
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/01—Arrangements for reducing harmonics or ripples
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/26—Arrangements for eliminating or reducing asymmetry in polyphase networks
    • H—ELECTRICITY
    • H03—ELECTRONIC CIRCUITRY
    • H03H—IMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
    • H03H21/00—Adaptive networks
    • H03H21/0012—Digital adaptive filters
    • H03H21/002—Filters with a particular frequency response
    • H03H21/0021—Notch filters

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Theoretical Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Supply And Distribution Of Alternating Current (AREA)
  • Stabilization Of Oscillater, Synchronisation, Frequency Synthesizers (AREA)

Abstract

本发明公开了一种自适应频率估计的锁相环以及锁相方法，所述锁相环包括：dq变换模块，用于对电网的三相电压进行dq变换，得到q轴电压分量；级联自适应陷波器，用于对所述q轴电压分量上的扰动分量进行滤除；锁相环节，用于根据滤除扰动分量后的q轴电压分量，得到所述三相电压的相位估计；本发明通过在锁相环引入级联自适应陷波器滤除电网电压的扰动分量，可以准确跟踪电压相位。

Description

一种自适应频率估计的锁相环以及锁相方法
技术领域
本发明涉及电网电压相位检测技术领域，尤其涉及一种自适应频率估计的锁相环以及锁相方法。
背景技术
随着电力系统的快速发展，大量电力电子装备得到了广泛的应用，包括交流电弧炉、连轧机组、大功率焊机、港口起重机、煤炭提升机及大型变流机组等。这些非线性、冲击性不对称负荷产生不对称无功和谐波电流注入电网、使公用电网电压产生波动、三相不平衡和谐波污染。同时，电网电压的相位作为电力电子装备控制环节和触发环节重要的参考基准，直接影响了电力电子装备的控制精度和产品性能。因此，可以快速准确跟踪被污染的电网电压相位的高性能锁相环是电力电子装备的重要组成部分。
三相同步锁相环（SRF-PLL）是构成高性能数字锁相环的基础，它可以快速跟踪理想电网电压的相位。但是，当电网电压被干扰后，SRF-PLL无法准确跟踪电压相位。
发明内容
本发明实施例提供一种自适应频率估计的锁相环以及锁相方法，通过在锁相环节引入级联自适应陷波器滤除电网电压的扰动分量，可以准确跟踪电压相位。
第一方面，本发明实施例提供了一种自适应频率估计的锁相环，包括：
dq变换模块，用于对电网的三相电压进行dq变换，得到q轴电压分量；
级联自适应陷波器，用于对所述q轴电压分量上的扰动分量进行滤除；
锁相环节，用于根据滤除扰动分量后的q轴电压分量，得到所述三相电压的相位估计。
作为上述方案的改进，所述锁相环节包括依次串联连接的比例积分调节器、比较器、积分环节以及比例环节；
其中，所述比例积分调节器，用于对滤除扰动分量后的q轴电压分量进PI调节，得到q轴的角频率分量；所述q轴的角频率分量经所述比较器、积分环节、比例环节后，得到所述相位估计。
作为上述方案的改进，所述级联自适应陷波器，包括多个串联连接的用于滤除特定次扰动分量的自适应陷波器。
作为上述方案的改进，所述自适应陷波器包括用于滤除2次扰动分量的第一陷波器、用于滤除6次扰动分量的第二陷波器以及用于滤除12次扰动分量的第三陷波器；其中，所述第一陷波器、第二陷波器、第三陷波器依次串联连接。
作为上述方案的改进，所述自适应陷波器为二阶自适应陷波器；
所述二阶自适应陷波器的自适应算法函数表达为：
Figure 195604DEST_PATH_IMAGE001
；
其中，
Figure 194784DEST_PATH_IMAGE002
表示状态变量，
Figure DEST_PATH_IMAGE003
表示需要滤除的扰动分量，
Figure 354501DEST_PATH_IMAGE004
表示自适应陷波器的迭代步长，
Figure DEST_PATH_IMAGE005
为自适应陷波器的回归因子，
Figure 203246DEST_PATH_IMAGE006
表示输入信号，
Figure DEST_PATH_IMAGE007
表示输出信号的估计，
Figure 106611DEST_PATH_IMAGE008
表示n次调节变量，a和r表示参数变量，n表示迭代次数。
第二方面，本发明实施例提供了一种自适应频率估计的锁相方法，包括：
对电网的三相电压进行dq变换，得到q轴电压分量；
通过预先构造的级联自适应陷波器，对所述q轴电压分量上的扰动分量进行滤除；
根据滤除扰动分量后的q轴电压分量，得到所述三相电压的相位估计。
作为上述方案的改进，所述根据滤除扰动分量后的q轴电压分量，得到所述三相电压的相位估计，包括：
将滤除扰动分量后的q轴电压分量送入比例积分调节器，得到q轴的角频率分量；
所述q轴的角频率分量经比较器、积分环节、比例环节后，得到所述相位估计。
作为上述方案的改进，所述对电网的三相电压进行dq变换，得到q轴电压分量，包括：
对所述三相电压进行Clark变换，得到
Figure DEST_PATH_IMAGE009
坐标系下的电压分量；
对
Figure 225877DEST_PATH_IMAGE009
坐标系下的电压分量进行Park变换，得到所述q轴电压分量。
作为上述方案的改进，所述级联自适应陷波器，包括多个串联连接的用于滤除特定次扰动分量的自适应陷波器。
作为上述方案的改进，所述自适应陷波器包括用于滤除2次扰动分量的第一陷波器、用于滤除6次扰动分量的第二陷波器以及用于滤除12次扰动分量的第三陷波器；其中，所述第一陷波器、第二陷波器、第三陷波器依次串联连接。
相对于现有技术，本发明实施例的有益效果在于：通过对电网的三相电压进行dq变换，得到q轴电压分量；然后使用预先构造的级联自适应陷波器，对所述q轴电压分量上的扰动分量进行滤除；再根据滤除扰动分量后的q轴电压分量，得到所述三相电压的相位估计；通过在锁相环引入级联自适应陷波器滤除电网电压的扰动分量，可以准确跟踪电压相位。
附图说明
为了更清楚地说明本发明的技术方案，下面将对实施方式中所占据要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施方式，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1是本发明实施例提供的自适应频率估计的锁相环的示意图；
图2是一阶陷波器的迭代结构示意图；
图3是二阶陷波器的迭代结构示意图；
图4是本发明实施例提供的自适应频率估计的锁相环的整体结构示意图；
图5是本发明实施例提供的一种自适应频率估计的锁相方法的流程图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
实施例一
请参阅图1，其是本发明实施例提供一种自适应频率估计的锁相环的结构示意图，所述自适应频率估计的锁相环包括：
dq变换模块1，用于对电网的三相电压进行dq变换，得到q轴电压分量；
级联自适应陷波器2，用于对所述q轴电压分量上的扰动分量进行滤除；
锁相环节3，用于根据滤除扰动分量后的q轴电压分量，得到所述三相电压的相位估计。
进一步，所述锁相环节包括依次串联连接的比例积分调节器、比较器、积分环节以及比例环节；
其中，所述比例积分调节器，用于对滤除扰动分量后的q轴电压分量进PI调节，得到q轴的角频率分量；所述q轴的角频率分量经所述比较器、积分环节、比例环节后，得到所述相位估计。
进一步，所述级联自适应陷波器，包括多个串联连接的用于滤除特定次扰动分量的自适应陷波器。
其中，所述自适应陷波器包括用于滤除2次扰动分量的第一陷波器、用于滤除6次扰动分量的第二陷波器以及用于滤除12次扰动分量的第三陷波器；其中，所述第一陷波器、第二陷波器、第三陷波器依次串联连接。
在一种可选的实施例中，所述自适应陷波器为二阶自适应陷波器。
所述二阶自适应陷波器的函数表达为：
Figure 795311DEST_PATH_IMAGE001
；
其中，
Figure 695134DEST_PATH_IMAGE002
表示状态变量，
Figure 198927DEST_PATH_IMAGE003
表示需要滤除的扰动分量，
Figure 500596DEST_PATH_IMAGE004
表示自适应陷波器的迭代步长，
Figure 798853DEST_PATH_IMAGE005
为自适应陷波器的回归因子，
Figure 622190DEST_PATH_IMAGE006
表示输入信号，
Figure 867358DEST_PATH_IMAGE007
表示输出信号的估计，
Figure 757954DEST_PATH_IMAGE008
表示n次调节变量，a和r表示参数变量，n表示迭代次数。
其中，级联自适应陷波器构造原理如下：
对于连续系统或离散系统，自适应陷波器可以采用梯度下降算法在误差性能面上寻找最优解。最小均方（LMS）算法被广泛用作梯度下降算法，以寻找二次误差最小为目标的函数。因此可以选择合适的代价函数来搜索最小点，使得自适应陷波器利用LMS算法识别一个输入信号频率。
示例性，将电网的三相电压看作是单一正弦波和背景白噪声的组合，对自适应陷波器的构造原理进行说明。三相电压表达为单一正弦波和背景白噪声的组合的函数为：
Figure 492691DEST_PATH_IMAGE010
；（1）
其中，p1表示正弦信号的幅值，
Figure DEST_PATH_IMAGE011
表示角频率，
Figure 475429DEST_PATH_IMAGE012
表示相位，
Figure 648921DEST_PATH_IMAGE013
表示白噪声（等同于下文所述的扰动分量），n表示迭代次数。白噪声
Figure 659602DEST_PATH_IMAGE013
与
Figure 565241DEST_PATH_IMAGE014
在统计学上互相独立，则自适应陷波器可以表达为：
Figure 599057DEST_PATH_IMAGE015
；（2）
其中，H(z)表示传递函数，z表示变换因子。
通过LMS算法可以得到代价函数，所述代价函数的表达如下：
Figure 310661DEST_PATH_IMAGE016
；（3）
其中，
Figure 612067DEST_PATH_IMAGE017
表示背景白噪声输出的方差，
Figure 688607DEST_PATH_IMAGE018
表示
Figure 272035DEST_PATH_IMAGE019
的范数，
Figure 459434DEST_PATH_IMAGE020
表示幅值响应的平方；
Figure 444708DEST_PATH_IMAGE019
表示传递函数。
理想自适应陷波器的频率响应为：
Figure 754466DEST_PATH_IMAGE021
；（4）
由公示（4）可以得知，除了
Figure 762873DEST_PATH_IMAGE022
处，
Figure 816280DEST_PATH_IMAGE023
为陷波频率，沿单位圆的幅值响应值始终等于1。
将公式（4）代入公式（3）可以得到：
Figure 656060DEST_PATH_IMAGE024
；（5）
代价函数
Figure 567079DEST_PATH_IMAGE025
在
Figure 859520DEST_PATH_IMAGE022
处取得最小值。
自适应陷波器的直接形式如图1所示，用于实现传递函数
Figure 716618DEST_PATH_IMAGE019
。
Figure 348587DEST_PATH_IMAGE026
；（6）
Figure 148DEST_PATH_IMAGE027
；（7）
Figure 983148DEST_PATH_IMAGE028
；（8）
自适应陷波器是带宽比较窄的带阻陷波器，可以快速滤波输入信号中的特定次频率分量，比如电网三相电压中的不平衡分量、谐波分量等扰动分量。在公式（6）中，分子多项式的零点在单位圆上，分母多项式的极点在单位圆上，通过调整略小于1来保持系统稳定。在本发明实施例中可通过参数a调整陷波频率。
其中，a和r表示参数变量，z表示变换因子；通过选择合适的参数a和参数r调节幅值响应
Figure 378357DEST_PATH_IMAGE029
。比如当r，陷波带宽减小，因此，在不考虑噪声的情况下，
Figure 927150DEST_PATH_IMAGE030
近似为理想响应，并且代价函数的梯度下降具有最小点。
将LMS算法引用于最小化代价函数，则有：
Figure 952875DEST_PATH_IMAGE031
；（9）
对于-2<a<2，则有：
Figure 219908DEST_PATH_IMAGE032
；（10）
此时，代价函数
Figure 418808DEST_PATH_IMAGE025
在a=-2cos
Figure 258326DEST_PATH_IMAGE033
处可以获得全局最小值。通过公式（9）和（10）可以得知，求取
Figure 517269DEST_PATH_IMAGE034
的解相当找到代价函数
Figure 271599DEST_PATH_IMAGE025
对于参数
Figure 946294DEST_PATH_IMAGE035
的全局最小值。
依据LMS直接型结构形式，一阶自适应陷波器的自适应迭代算法的表达为：
Figure 469679DEST_PATH_IMAGE036
；（11）
Figure 899523DEST_PATH_IMAGE037
；（12）
其中，
Figure 875569DEST_PATH_IMAGE038
表示自适应陷波器的迭代步长，
Figure 353955DEST_PATH_IMAGE005
为自适应陷波器的回归因子。
因此，可以得到二阶自适应陷波滤波器的迭代结构示意图2，二阶自适应陷波器的自适应算法的表达为：
Figure 997426DEST_PATH_IMAGE001
；（13）
其中，
Figure 270276DEST_PATH_IMAGE002
表示状态变量，
Figure 999197DEST_PATH_IMAGE003
表示需要滤除的扰动分量，
Figure 514230DEST_PATH_IMAGE004
表示自适应陷波器的迭代步长，
Figure 12207DEST_PATH_IMAGE005
为自适应陷波器的回归因子，
Figure 518275DEST_PATH_IMAGE006
表示输入信号，
Figure 672176DEST_PATH_IMAGE007
表示输出信号的估计，
Figure 492364DEST_PATH_IMAGE008
表示n次调节变量，a和r表示参数变量，n表示迭代次数。
在发明实施例中，以非理想电网的三相电压作为输入信号，由于三相电压含有不平衡分量、谐波分量等干扰，不平衡分量经过dq变换模块在q轴上面变成2次扰动分量，5和7次谐波分量经过dq变换模块在q轴上面变成6次扰动分量，同理，11和13次谐波分量经过dq变换模块在q轴上面变成12次扰动分量。为了滤除三相电压q轴上面的干扰，保证相位跟踪的准确性，需要在q轴上面加入滤波器消除干扰项。单独自适应陷波器只对单一频率的干扰具有很好的滤除效果。在本发明实施例中为了消除多次干扰分量，采用自适应陷波器级联的结构方式，根据q轴信号上面的不同干扰量，设置不同自适应陷波器的组合，每个自适应陷波滤波器负责滤除特定次扰动分量。比如基于上述自适应滤波器的原理，可以构造出三个二阶自适应陷波器串联的结构，分别用于滤除2次扰动分量V2、6次扰动分量V6、12次扰动分量V12。
自适应频率估计锁相环：
将输入信号
Figure 844848DEST_PATH_IMAGE039
的电压相角记为
Figure 787397DEST_PATH_IMAGE040
，角频率记为
Figure 163014DEST_PATH_IMAGE033
。当dq坐标系与
Figure 583631DEST_PATH_IMAGE041
按照同步转速
Figure 492419DEST_PATH_IMAGE042
旋转时，
Figure 340289DEST_PATH_IMAGE039
与d轴之前的相位是固定的，记为
Figure 265520DEST_PATH_IMAGE043
。
Figure 427511DEST_PATH_IMAGE039
在dq坐标系上的投影
Figure 754587DEST_PATH_IMAGE044
，
Figure 38938DEST_PATH_IMAGE045
。若
Figure 389148DEST_PATH_IMAGE043
=0，
Figure 151568DEST_PATH_IMAGE046
，此时
Figure 333150DEST_PATH_IMAGE047
=0，d轴与
Figure 726086DEST_PATH_IMAGE039
重合。
因此，
Figure 625908DEST_PATH_IMAGE047
作为被控变量，相位估计
Figure 634097DEST_PATH_IMAGE048
对三相电压
Figure 670186DEST_PATH_IMAGE039
进行dq变换，通过反馈信号对
Figure 296339DEST_PATH_IMAGE048
进行修正，保证
Figure 621141DEST_PATH_IMAGE047
=0。最终
Figure 725363DEST_PATH_IMAGE048
成为输入信号
Figure 881538DEST_PATH_IMAGE039
的相位估计。
自适应频率估计锁相环的结构如图4所示，包括dq变换模块，级联自适应陷波器，比例积分调节器(PI)、比较器、积分环节、比例环节。通过将级联自适应陷波器（G2（z）、G6（z）、G12（z））嵌入三相同步锁相环中，滤除三相电压
Figure 616276DEST_PATH_IMAGE039
中的干扰量。
为了更清楚地说明本发明的工作原理，下面以电网接入变流器为例，分析电网公共连接点的电网（PCC）电压的扰动分量，其中，电网公共连接点的电网（PCC）电压可以表示为：
Figure 225112DEST_PATH_IMAGE049
；（14）
其中，V1表示电压幅值，
Figure 398604DEST_PATH_IMAGE040
表示电压相角，
Figure 346969DEST_PATH_IMAGE043
表示B相不平衡因子，（k表示数值1到无穷递增），
Figure 314925DEST_PATH_IMAGE050
表示C相不平衡因子；V6k+1表示谐波电压幅值；
VPCC为三相电压，VPCC首先经过Clark变换，abc坐标系变为
Figure 847275DEST_PATH_IMAGE009
坐标系的表达，然后以相位估计
Figure 293300DEST_PATH_IMAGE048
进行Park变换，得到VPCC在q轴的表达如公式（15）所示：
Figure 424067DEST_PATH_IMAGE051
；（15）
其中，
Figure 500607DEST_PATH_IMAGE052
，
Figure 84035DEST_PATH_IMAGE053
，E1表示
Figure 333751DEST_PATH_IMAGE009
坐标系下电压幅值，
Figure 256708DEST_PATH_IMAGE054
表示相角误差；
根据工程经验，PCC点电网电压在abc坐标系的参数设定如表1所示。
表1-输入信号参数设置
Figure 566466DEST_PATH_IMAGE055
上述V1表示电压理想状态（即无干扰）的值；(1+
Figure 574874DEST_PATH_IMAGE043
)V1、(1+
Figure 628280DEST_PATH_IMAGE050
)V1表示电压的不平衡分量，V5、V7、V11、V13对应表示5、7、11、13次谐波分量。
则通过上述公式（15），可以得到dq坐标的参数如表2所示。
表2-Vq轴扰动分量
Figure 468060DEST_PATH_IMAGE056
根据上述Vq轴扰动分量V2、V6、V12，级联自适应陷波器由3个自适应陷波滤波器串联构成，分别滤除2次谐波，6次谐波和12次谐波扰动分量。具体结构如图4所示。为了滤除对应的扰动分量，需要调节自适应陷波滤波器的中心频率，分别对应2次，6次和12次谐波、频带宽度（BW2、BW6、BW12）和迭代因子（
Figure 384939DEST_PATH_IMAGE057
、
Figure 677380DEST_PATH_IMAGE058
、
Figure 534477DEST_PATH_IMAGE059
）。
在本发明实例中，自适应陷波器使用的自适应递归算法可以在没有参考估计频率的情况下跟踪频率变化，因此，可以提高自适应陷波器的品质因数，有效地提高滤波能力。对于固定陷波器，必须调整较小的Q值以增加滤波器带宽并减少频率变化的干扰。系统相位延迟与自适应陷波器的参数和Q值有关。因此，必须设计自适应陷波器的参和Q以获得合适的相位延迟。在本发明实施例中调谐参数如表3所示。
表3-自适应陷波器的参数
Figure 166447DEST_PATH_IMAGE060
为了验证上述自适应陷波器的参数是否合适，需要对系统参数进行仿真验证，确认系统参数是否满足锁相环的性能要求。
其中，输入信号中叠加不平衡和谐波干扰分量见上述表1。频率的波动变化设定为47Hz-53Hz，即输入信号频率分别设定为47Hz和53Hz，通过Vq轴上扰动分量的消除比例来确认级联自适应滤波器性能。根据上述表2所示的Vq轴中的扰动分量，47Hz和53Hz扰动滤除比例如下表4所示。
表4
Figure 818008DEST_PATH_IMAGE061
如上表所示，本发明实施例涉及的锁相环具有较好的扰动滤除性能，能够准确的反应输入信号的相位信息。
相对于现有技术，本发明实施例提供一种包括dq变换模块，级联自适应陷波器，比例积分调节器(PI)、比较器、积分环节、比例环节的自适应频率估计的锁相环，通过将级联自适应陷波器嵌入传统的三相同步锁相环（SRF-PLL）中，构成自适应频率估计的锁相环，通过级联自适应陷波滤波器可以自动跟踪电网电压频率的变化，同时消除非理想电网电压的扰动，包括三相不平衡产生的负序分量，谐波分量，能够快速准确估计电网电压频率波动。
实施例二
请参阅图5，本发明实施例提供了一种自适应频率估计的锁相方法，包括：
S1：对电网的三相电压进行dq变换，得到q轴电压分量；
S2：通过预先构造的级联自适应陷波器，对所述q轴电压分量上的扰动分量进行滤除；
S3：根据滤除扰动分量后的q轴电压分量，得到所述三相电压的相位估计。
作为上述方案的改进，所述根据滤除扰动分量后的q轴电压分量，得到所述三相电压的相位估计，包括：
将滤除扰动分量后的q轴电压分量送入比例积分调节器，得到q轴的角频率分量；
所述q轴的角频率分量经比较器、积分环节、比例环节后，得到所述相位估计。
作为上述方案的改进，所述对电网的三相电压进行dq变换，得到q轴电压分量，包括：
对所述三相电压进行Clark变换，得到坐标系下的电压分量；
对
Figure 801008DEST_PATH_IMAGE009
坐标系下的电压分量进行Park变换，得到所述q轴电压分量。
作为上述方案的改进，所述三相电压的函数表达为：
Figure 196217DEST_PATH_IMAGE049
；
其中，V1表示电压幅值，
Figure 745010DEST_PATH_IMAGE040
表示电压相角，
Figure 770735DEST_PATH_IMAGE043
表示B相不平衡因子，（k表示数值1到无穷递增），
Figure 37768DEST_PATH_IMAGE050
表示C相不平衡因子；V6k+1表示谐波电压幅值；
所述q轴电压分量的函数表达为：
Figure 971089DEST_PATH_IMAGE051
；
其中，
Figure 76186DEST_PATH_IMAGE052
，
Figure 335129DEST_PATH_IMAGE053
；E1表示
Figure 823879DEST_PATH_IMAGE009
坐标系下电压幅值，
Figure 764153DEST_PATH_IMAGE054
表示相角误差。
在一种可选的实施例中，所述级联自适应陷波器，包括多个串联连接的用于滤除特定次扰动分量的自适应陷波器。
进一步，所述自适应陷波器包括用于滤除2次扰动分量的第一陷波器、用于滤除6次扰动分量的第二陷波器以及用于滤除12次扰动分量的第三陷波器；其中，所述第一陷波器、第二陷波器、第三陷波器依次串联连接。
进一步，所述自适应陷波器为二阶自适应陷波器；
所述二阶自适应陷波器的自适应算法函数表达为：
Figure 287538DEST_PATH_IMAGE001
；
其中，
Figure 389487DEST_PATH_IMAGE002
表示状态变量，
Figure 631112DEST_PATH_IMAGE003
表示需要滤除的扰动分量，
Figure 171815DEST_PATH_IMAGE004
表示自适应陷波器的迭代步长，
Figure 487390DEST_PATH_IMAGE005
为自适应陷波器的回归因子，
Figure 88135DEST_PATH_IMAGE006
表示输入信号，
Figure 551478DEST_PATH_IMAGE007
表示输出信号的估计，
Figure 349668DEST_PATH_IMAGE008
表示n次调节变量，a和r表示参数变量，n表示迭代次数。
其中，自适应陷波器的参数设置以及工作原理请参见实施例一，在这里不再重复说明。
相对于现有技术，本发明实施例的有益效果在于：通过对电网的三相电压进行dq变换，得到q轴电压分量；然后使用预先构造的级联自适应陷波器，对所述q轴电压分量上的扰动分量进行滤除；再根据滤除扰动分量后的q轴电压分量，得到所述三相电压的相位估计；通过在锁相环节引入级联自适应陷波器滤除电网电压的扰动分量，可以准确跟踪电压相位。
需说明的是，以上所描述的装置实施例仅仅是示意性的，其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。另外，本发明提供的装置实施例附图中，模块之间的连接关系表示它们之间具有通信连接，具体可以实现为一条或多条通信总线或信号线。本领域普通技术人员在不付出创造性劳动的情况下，即可以理解并实施。
以上所述是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出多台改进和润饰，这些改进和润饰也视为本发明的保护范围。

Claims (10)

1.一种自适应频率估计的锁相环，其特征在于，包括：
dq变换模块，用于对电网的三相电压进行dq变换，得到q轴电压分量；
级联自适应陷波器，用于对所述q轴电压分量上的扰动分量进行滤除；
锁相环节，用于根据滤除扰动分量后的q轴电压分量，得到所述三相电压的相位估计。
2.如权利要求1所述的自适应频率估计的锁相环，其特征在于，所述锁相环节包括依次串联连接的比例积分调节器、比较器、积分环节以及比例环节；
其中，所述比例积分调节器，用于对滤除扰动分量后的q轴电压分量进PI调节，得到q轴的角频率分量；所述q轴的角频率分量经所述比较器、积分环节、比例环节后，得到所述相位估计。
3.如权利要求1所述的自适应频率估计的锁相环，其特征在于，所述级联自适应陷波器，包括多个串联连接的用于滤除特定次扰动分量的自适应陷波器。
4.如权利要求3所述的自适应频率估计的锁相环，其特征在于，所述自适应陷波器包括用于滤除2次扰动分量的第一陷波器、用于滤除6次扰动分量的第二陷波器以及用于滤除12次扰动分量的第三陷波器；其中，所述第一陷波器、第二陷波器、第三陷波器依次串联连接。
5.如权利要求3所述的自适应频率估计的锁相环，其特征在于，所述自适应陷波器为二阶自适应陷波器；
所述二阶自适应陷波器的自适应算法函数表达为：
Figure 877616DEST_PATH_IMAGE001
；
其中，
Figure 765937DEST_PATH_IMAGE002
表示状态变量，
Figure 586126DEST_PATH_IMAGE003
表示需要滤除的扰动分量，
Figure 876293DEST_PATH_IMAGE005
表示自适应陷波器的迭代步长，
Figure 192742DEST_PATH_IMAGE006
为自适应陷波器的回归因子，
Figure 630677DEST_PATH_IMAGE007
表示输入信号，
Figure 51294DEST_PATH_IMAGE008
表示输出信号的估计，
Figure 461547DEST_PATH_IMAGE009
表示n次调节变量，a和r表示参数变量，n表示迭代次数。
6.一种自适应频率估计的锁相方法，其特征在于，包括：
对电网的三相电压进行dq变换，得到q轴电压分量；
通过预先构造的级联自适应陷波器，对所述q轴电压分量上的扰动分量进行滤除；
根据滤除扰动分量后的q轴电压分量，得到所述三相电压的相位估计。
7.如权利要求6所述的自适应频率估计的锁相方法，其特征在于，所述根据滤除扰动分量后的q轴电压分量，得到所述三相电压的相位估计，包括：
将滤除扰动分量后的q轴电压分量送入比例积分调节器，得到q轴的角频率分量；
所述q轴的角频率分量经比较器、积分环节、比例环节后，得到所述相位估计。
8.如权利要求6所述的自适应频率估计的锁相方法，其特征在于，所述对电网的三相电压进行dq变换，得到q轴电压分量，包括：
对所述三相电压进行Clark变换，得到
Figure 574996DEST_PATH_IMAGE010
坐标系下的电压分量；
对
Figure 437910DEST_PATH_IMAGE010
坐标系下的电压分量进行Park变换，得到所述q轴电压分量。
9.如权利要求6所述的自适应频率估计的锁相方法，其特征在于，所述级联自适应陷波器，包括多个串联连接的用于滤除特定次扰动分量的自适应陷波器。
10.如权利要求9所述的自适应频率估计的锁相方法，其特征在于，所述自适应陷波器包括用于滤除2次扰动分量的第一陷波器、用于滤除6次扰动分量的第二陷波器以及用于滤除12次扰动分量的第三陷波器；其中，所述第一陷波器、第二陷波器、第三陷波器依次串联连接。
CN202210683832.2A 2022-06-17 2022-06-17 一种自适应频率估计的锁相环以及锁相方法 Active CN114784810B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210683832.2A CN114784810B (zh) 2022-06-17 2022-06-17 一种自适应频率估计的锁相环以及锁相方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210683832.2A CN114784810B (zh) 2022-06-17 2022-06-17 一种自适应频率估计的锁相环以及锁相方法

Publications (2)

Publication Number Publication Date
CN114784810A true CN114784810A (zh) 2022-07-22
CN114784810B CN114784810B (zh) 2022-09-16

Family

ID=82420231

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210683832.2A Active CN114784810B (zh) 2022-06-17 2022-06-17 一种自适应频率估计的锁相环以及锁相方法

Country Status (1)

Country Link
CN (1) CN114784810B (zh)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100312365A1 (en) * 2009-06-05 2010-12-09 University Of Southern California Adaptive notch filter
CN102590829A (zh) * 2012-03-14 2012-07-18 西安电子科技大学 用于卫星导航系统的完全自适应陷波器及其陷波方法
CN103227640A (zh) * 2013-04-16 2013-07-31 苏州汇川技术有限公司 谐波电网相位检测方法及锁相环
CN103929151A (zh) * 2014-04-21 2014-07-16 北京航空航天大学 自适应最优相角陷波滤波器设计方法
CN104901685A (zh) * 2015-05-20 2015-09-09 南京南瑞太阳能科技有限公司 一种基于自适应陷波器的锁相环

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100312365A1 (en) * 2009-06-05 2010-12-09 University Of Southern California Adaptive notch filter
CN102590829A (zh) * 2012-03-14 2012-07-18 西安电子科技大学 用于卫星导航系统的完全自适应陷波器及其陷波方法
CN103227640A (zh) * 2013-04-16 2013-07-31 苏州汇川技术有限公司 谐波电网相位检测方法及锁相环
CN103929151A (zh) * 2014-04-21 2014-07-16 北京航空航天大学 自适应最优相角陷波滤波器设计方法
CN104901685A (zh) * 2015-05-20 2015-09-09 南京南瑞太阳能科技有限公司 一种基于自适应陷波器的锁相环

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
JING LU: "Research on Optimization Method of Passive Resonance Damping of Hybrid Active Power Filter", 《2022 INTERNATIONAL CONFERENCE ON POWER ENERGY SYSTEMS AND APPLICATIONS》 *
田贇祥: "托卡马克电源窄带宽近基频间谐波检测方法", 《核 聚 变 与 等 离 子 体 物 理》 *
解朦朦等: "基于复数陷波器的窄带干扰抑制研究", 《电子科技》 *

Also Published As

Publication number Publication date
CN114784810B (zh) 2022-09-16

Similar Documents

Publication Publication Date Title
CN106655276B (zh) 一种适用于三相电网电压的锁相方法
Wang et al. A method to improve the dynamic performance of moving average filter-based PLL
CN106953634B (zh) 一种基于双自调谐二阶广义积分器的锁频环方法
CN108599261B (zh) 基于非线性pi和解耦双同步坐标系锁相环的锁相方法
CN105680854B (zh) 一种锁相环及锁相方法
CN110557118B (zh) 一种锁相装置及锁相方法
Dinh et al. Adaptive notch filter solution under unbalanced and/or distorted point of common coupling voltage for three‐phase four‐wire shunt active power filter with sinusoidal utility current strategy
CN106936427B (zh) 一种具有谐波抑制功能的单相锁相环优化设计方法
CN107591809A (zh) 多周期并联重复控制谐波和间谐波指定次补偿方法
Sridharan et al. Accurate phase detection system using modified SGDFT-based PLL for three-phase grid-interactive power converter during interharmonic conditions
CN107863774B (zh) 谐波指令电流的获取方法、装置、系统及可读存储介质
CN116094025A (zh) 一种含跟网型和构网型变流器功率系统暂态稳定分析方法
CN109358228B (zh) 基于双增强型锁相环的电网电压正负序分量实时估计方法
Gautam et al. Enhanced single-phase phase locked loop based on complex-coefficient filter
CN105182072B (zh) 一种基于虚拟滤波因子的谐波检测方法与系统
CN110749769A (zh) 一种三相电压跌落的快速检测方法
US10215784B1 (en) Measuring apparatus including phase locked loop and measuring method thereof
CN114784810B (zh) 一种自适应频率估计的锁相环以及锁相方法
CN108110802A (zh) 一种并网功率控制方法
CN116760070A (zh) 直驱风机并网系统宽频振荡抑制方法、装置、介质及设备
CN110048442B (zh) 模块化多电平变换器微分平滑非线性控制方法及装置
CN106451573A (zh) 一种多变量反馈控制的三相lcl型联网变换器及方法
Kang et al. Improving harmonic rejection capability of OSG based on n-th order bandpass filter for single-phase system
da Fonseca et al. Three-phase phase-locked loop algorithms for the detection of positive-sequence and negative-sequence components
CN112702058A (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