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

摘要

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

权利要求

1.一种自适应频率估计的锁相环，其特征在于，包括：

dq变换模块，用于对电网的三相电压进行dq变换，得到q轴电压分量；

级联自适应陷波器，用于对所述q轴电压分量上的扰动分量进行滤除；

锁相环节，用于根据滤除扰动分量后的q轴电压分量，得到所述三相电压的相位估计；

所述级联自适应陷波器，包括多个串联连接的用于滤除特定次扰动分量的自适应陷波器；

所述自适应陷波器为二阶自适应陷波器；所述二阶自适应陷波器的自适应算法函数表达为：；

其中， 表示状态变量， 表示需要滤除的扰动分量，表示自适应陷波器的迭代步长， 为自适应陷波器的回归因子， 表示输入信号， 表示输出信号的估计， 表示n次调节变量，a和r表示参数变量， n表示迭代次数。

2.如权利要求1所述的自适应频率估计的锁相环，其特征在于，所述锁相环节包括依次串联连接的比例积分调节器、比较器、积分环节以及比例环节；

其中，所述比例积分调节器，用于对滤除扰动分量后的q轴电压分量进PI调节，得到q轴的角频率分量；所述q轴的角频率分量经所述比较器、积分环节、比例环节后，得到所述相位估计。

3.如权利要求1所述的自适应频率估计的锁相环，其特征在于，所述自适应陷波器包括用于滤除2次扰动分量的第一自适应陷波器、用于滤除6次扰动分量的第二自适应陷波器以及用于滤除12次扰动分量的第三自适应陷波器；其中，所述第一自适应陷波器、第二自适应陷波器、第三自适应陷波器依次串联连接。

4.一种自适应频率估计的锁相方法，其特征在于，包括：对电网的三相电压进行dq变换，得到q轴电压分量；

通过预先构造的级联自适应陷波器，对所述q轴电压分量上的扰动分量进行滤除；

根据滤除扰动分量后的q轴电压分量，得到所述三相电压的相位估计；

其中，所述自适应陷波器为二阶自适应陷波器；所述二阶自适应陷波器的自适应算法函数表达为：；

其中， 表示状态变量， 表示需要滤除的扰动分量，表示自适应陷波器的迭代步长， 为自适应陷波器的回归因子， 表示输入信号， 表示输出信号的估计，表示n次调节变量，a和r表示参数变量， n表示迭代次数。

5.如权利要求4所述的自适应频率估计的锁相方法，其特征在于，所述根据滤除扰动分量后的q轴电压分量，得到所述三相电压的相位估计，包括：将滤除扰动分量后的q轴电压分量送入比例积分调节器，得到q轴的角频率分量；

所述q轴的角频率分量经比较器、积分环节、比例环节后，得到所述相位估计。

6.如权利要求4所述的自适应频率估计的锁相方法，其特征在于，所述对电网的三相电压进行dq变换，得到q轴电压分量，包括：对所述三相电压进行Clark变换，得到 坐标系下的电压分量；

对 坐标系下的电压分量进行Park变换，得到所述q轴电压分量。

7.如权利要求4所述的自适应频率估计的锁相方法，其特征在于，所述自适应陷波器包括用于滤除2次扰动分量的第一自适应陷波器、用于滤除6次扰动分量的第二自适应陷波器以及用于滤除12次扰动分量的第三自适应陷波器；其中，所述第一自适应陷波器、第二自适应陷波器、第三自适应陷波器依次串联连接。

说明书

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

技术领域

[0001] 本发明涉及电网电压相位检测技术领域，尤其涉及一种自适应频率估计的锁相环以及锁相方法。

背景技术

[0002] 随着电力系统的快速发展，大量电力电子装备得到了广泛的应用，包括交流电弧炉、连轧机组、大功率焊机、港口起重机、煤炭提升机及大型变流机组等。这些非线性、冲击性不对称负荷产生不对称无功和谐波电流注入电网、使公用电网电压产生波动、三相不平衡和谐波污染。同时，电网电压的相位作为电力电子装备控制环节和触发环节重要的参考基准，直接影响了电力电子装备的控制精度和产品性能。因此，可以快速准确跟踪被污染的电网电压相位的高性能锁相环是电力电子装备的重要组成部分。

[0003] 三相同步锁相环（SRF‑PLL）是构成高性能数字锁相环的基础，它可以快速跟踪理想电网电压的相位。但是，当电网电压被干扰后，SRF‑PLL无法准确跟踪电压相位。

发明内容

[0004] 本发明实施例提供一种自适应频率估计的锁相环以及锁相方法，通过在锁相环节引入级联自适应陷波器滤除电网电压的扰动分量，可以准确跟踪电压相位。

[0005] 第一方面，本发明实施例提供了一种自适应频率估计的锁相环，包括：

[0006] dq变换模块，用于对电网的三相电压进行dq变换，得到q轴电压分量；

[0007] 级联自适应陷波器，用于对所述q轴电压分量上的扰动分量进行滤除；

[0008] 锁相环节，用于根据滤除扰动分量后的q轴电压分量，得到所述三相电压的相位估计。

[0009] 作为上述方案的改进，所述锁相环节包括依次串联连接的比例积分调节器、比较器、积分环节以及比例环节；

[0010] 其中，所述比例积分调节器，用于对滤除扰动分量后的q轴电压分量进PI调节，得到q轴的角频率分量；所述q轴的角频率分量经所述比较器、积分环节、比例环节后，得到所述相位估计。

[0011] 作为上述方案的改进，所述级联自适应陷波器，包括多个串联连接的用于滤除特定次扰动分量的自适应陷波器。

[0012] 作为上述方案的改进，所述自适应陷波器包括用于滤除2次扰动分量的第一自适应陷波器、用于滤除6次扰动分量的第二自适应陷波器以及用于滤除12次扰动分量的第三自适应陷波器；其中，所述第一自适应陷波器、第二自适应陷波器、第三自适应陷波器依次串联连接。

[0013] 作为上述方案的改进，所述自适应陷波器为二阶自适应陷波器；

[0014] 所述二阶自适应陷波器的自适应算法函数表达为：

[0015] ；

[0016] 其中， 表示状态变量， 表示需要滤除的扰动分量，表示自适应陷波器的迭代步长， 为自适应陷波器的回归因子， 表示输入信号， 表示输出信号的估计，表示n次调节变量，a和r表示参数变量，n表示迭代次数。

[0017] 第二方面，本发明实施例提供了一种自适应频率估计的锁相方法，包括：

[0018] 对电网的三相电压进行dq变换，得到q轴电压分量；

[0019] 通过预先构造的级联自适应陷波器，对所述q轴电压分量上的扰动分量进行滤除；

[0020] 根据滤除扰动分量后的q轴电压分量，得到所述三相电压的相位估计。

[0021] 作为上述方案的改进，所述根据滤除扰动分量后的q轴电压分量，得到所述三相电压的相位估计，包括：

[0022] 将滤除扰动分量后的q轴电压分量送入比例积分调节器，得到q轴的角频率分量；

[0023] 所述q轴的角频率分量经比较器、积分环节、比例环节后，得到所述相位估计。

[0024] 作为上述方案的改进，所述对电网的三相电压进行dq变换，得到q轴电压分量，包括：

[0025] 对所述三相电压进行Clark变换，得到 坐标系下的电压分量；

[0026] 对 坐标系下的电压分量进行Park变换，得到所述q轴电压分量。

[0027] 作为上述方案的改进，所述级联自适应陷波器，包括多个串联连接的用于滤除特定次扰动分量的自适应陷波器。

[0028] 作为上述方案的改进，所述自适应陷波器包括用于滤除2次扰动分量的第一自适应陷波器、用于滤除6次扰动分量的第二自适应陷波器以及用于滤除12次扰动分量的第三自适应陷波器；其中，所述第一自适应陷波器、第二自适应陷波器、第三自适应陷波器依次串联连接。

[0029] 相对于现有技术，本发明实施例的有益效果在于：通过对电网的三相电压进行dq变换，得到q轴电压分量；然后使用预先构造的级联自适应陷波器，对所述q轴电压分量上的扰动分量进行滤除；再根据滤除扰动分量后的q轴电压分量，得到所述三相电压的相位估计；通过在锁相环引入级联自适应陷波器滤除电网电压的扰动分量，可以准确跟踪电压相位。

附图说明

[0030] 为了更清楚地说明本发明的技术方案，下面将对实施方式中所占据要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施方式，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。

[0031] 图1是本发明实施例提供的自适应频率估计的锁相环的示意图；

[0032] 图2是一阶陷波器的迭代结构示意图；

[0033] 图3是二阶陷波器的迭代结构示意图；

[0034] 图4是本发明实施例提供的自适应频率估计的锁相环的整体结构示意图；

[0035] 图5是本发明实施例提供的一种自适应频率估计的锁相方法的流程图。

具体实施方式

[0036] 下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

[0037] 实施例一

[0038] 请参阅图1，其是本发明实施例提供一种自适应频率估计的锁相环的结构示意图，所述自适应频率估计的锁相环包括：

[0039] dq变换模块1，用于对电网的三相电压进行dq变换，得到q轴电压分量；

[0040] 级联自适应陷波器2，用于对所述q轴电压分量上的扰动分量进行滤除；

[0041] 锁相环节3，用于根据滤除扰动分量后的q轴电压分量，得到所述三相电压的相位估计。

[0042] 进一步，所述锁相环节包括依次串联连接的比例积分调节器、比较器、积分环节以及比例环节；

[0043] 其中，所述比例积分调节器，用于对滤除扰动分量后的q轴电压分量进PI调节，得到q轴的角频率分量；所述q轴的角频率分量经所述比较器、积分环节、比例环节后，得到所述相位估计。

[0044] 进一步，所述级联自适应陷波器，包括多个串联连接的用于滤除特定次扰动分量的自适应陷波器。

[0045] 其中，所述自适应陷波器包括用于滤除2次扰动分量的第一自适应陷波器、用于滤除6次扰动分量的第二自适应陷波器以及用于滤除12次扰动分量的第三自适应陷波器；其中，所述第一自适应陷波器、第二自适应陷波器、第三自适应陷波器依次串联连接。

[0046] 在一种可选的实施例中，所述自适应陷波器为二阶自适应陷波器。

[0047] 所述二阶自适应陷波器的函数表达为：

[0048] ；

[0049] 其中， 表示状态变量， 表示需要滤除的扰动分量，表示自适应陷波器的迭代步长， 为自适应陷波器的回归因子， 表示输入信号， 表示输出信号的估计， 表示n次调节变量，a和r表示参数变量，n表示迭代次数。

[0050] 其中，级联自适应陷波器构造原理如下：

[0051] 对于连续系统或离散系统，自适应陷波器可以采用梯度下降算法在误差性能面上寻找最优解。最小均方（LMS）算法被广泛用作梯度下降算法，以寻找二次误差最小为目标的函数。因此可以选择合适的代价函数来搜索最小点，使得自适应陷波器利用LMS算法识别一个输入信号频率。

[0052] 示例性，将电网的三相电压看作是单一正弦波和背景白噪声的组合，对自适应陷波器的构造原理进行说明。三相电压表达为单一正弦波和背景白噪声的组合的函数为：

[0053] ；（1）

[0054] 其中，p1表示正弦信号的幅值， 表示角频率， 表示相位， 表示白噪声（等同于下文所述的扰动分量），n表示迭代次数。白噪声 与 在统计学上互相独立，则自适应陷波器可以表达为：

[0055] ；（2）

[0056] 其中，H(z)表示传递函数，z表示变换因子。

[0057] 通过LMS算法可以得到代价函数，所述代价函数的表达如下：

[0058] ；（3）

[0059] 其中， 表示背景白噪声输出的方差， 表示 的范数， 表示幅值响应的平方； 表示传递函数。

[0060] 理想自适应陷波器的频率响应为：

[0061] ；（4）

[0062] 由公示（4）可以得知，除了 处， 为陷波频率，沿单位圆的幅值响应值始终等于1。

[0063] 将公式（4）代入公式（3）可以得到：

[0064] ；（5）

[0065] 代价函数 在 处取得最小值。

[0066] 自适应陷波器的直接形式如图2所示，用于实现传递函数 。

[0067] ；（6）

[0068] ；（7）

[0069] ；（8）

[0070] 自适应陷波器是带宽比较窄的带阻陷波器，可以快速滤波输入信号中的特定次频率分量，比如电网三相电压中的不平衡分量、谐波分量等扰动分量。在公式（6）中，分子多项式的零点在单位圆上，分母多项式的极点在单位圆上，通过调整略小于1来保持系统稳定。在本发明实施例中可通过参数a调整陷波频率。

[0071] 其中，a和r表示参数变量，z表示变换因子；通过选择合适的参数a和参数r调节幅值响应 。比如当r，陷波带宽减小，因此，在不考虑噪声的情况下， 近似为理想响应，并且代价函数的梯度下降具有最小点。

[0072] 将LMS算法引用于最小化代价函数，则有：

[0073] ；（9）

[0074] 对于‑2

[0075] ；（10）

[0076] 此时，代价函数 在a=‑2cos 处可以获得全局最小值。通过公式（9）和（10）可以得知，求取 的解相当找到代价函数 对于参数 的全局最小值。

[0077] 依据LMS直接型结构形式，一阶自适应陷波器的自适应迭代算法的表达为：

[0078] ；（11）

[0079] ；（12）

[0080] 其中，表示自适应陷波器的迭代步长， 为自适应陷波器的回归因子。

[0081] 因此，可以得到二阶自适应陷波滤波器的迭代结构示意图3，二阶自适应陷波器的自适应算法的表达为：

[0082] ；（13）

[0083] 其中， 表示状态变量， 表示需要滤除的扰动分量，表示自适应陷波器的迭代步长， 为自适应陷波器的回归因子， 表示输入信号， 表示输出信号的估计， 表示n次调节变量，a和r表示参数变量，n表示迭代次数。

[0084] 在发明实施例中，以非理想电网的三相电压作为输入信号，由于三相电压含有不平衡分量、谐波分量等干扰，不平衡分量经过dq变换模块在q轴上面变成2次扰动分量，5和7次谐波分量经过dq变换模块在q轴上面变成6次扰动分量，同理，11和13次谐波分量经过dq变换模块在q轴上面变成12次扰动分量。为了滤除三相电压q轴上面的干扰，保证相位跟踪的准确性，需要在q轴上面加入滤波器消除干扰项。单独自适应陷波器只对单一频率的干扰具有很好的滤除效果。在本发明实施例中为了消除多次干扰分量，采用自适应陷波器级联的结构方式，根据q轴信号上面的不同干扰量，设置不同自适应陷波器的组合，每个自适应陷波滤波器负责滤除特定次扰动分量。比如基于上述自适应滤波器的原理，可以构造出三个二阶自适应陷波器串联的结构，分别用于滤除2次扰动分量V2、6次扰动分量V6、12次扰动分量V12。

[0085] 自适应频率估计锁相环：

[0086] 将输入信号 的电压相角记为 ，角频率记为 。当dq坐标系与 按照同步转速 旋转时， 与d轴之前的相位是固定的，记为 。 在dq坐标系上的投影， 。若 =0， ，此时 =0，d轴与 重合。

[0087] 因此， 作为被控变量，相位估计 对三相电压 进行dq变换，通过反馈信号对进行修正，保证 =0。最终 成为输入信号 的相位估计。

[0088] 自适应频率估计锁相环的结构如图4所示，包括dq变换模块，级联自适应陷波器，比例积分调节器(PI)、比较器、积分环节、比例环节。通过将级联自适应陷波器（G（2 z）、G6（z）、G1（2 z））嵌入三相同步锁相环中，滤除三相电压 中的干扰量。

[0089] 为了更清楚地说明本发明的工作原理，下面以电网接入变流器为例，分析电网公共连接点的电网（PCC）电压的扰动分量，其中，电网公共连接点的电网（PCC）电压可以表示为：

[0090] ；（14）

[0091] 其中，V1表示电压幅值，表示电压相角，表示B相不平衡因子，（k表示数值1到无穷递增），表示C相不平衡因子；V6k+1表示谐波电压幅值；

[0092] VPCC为三相电压，VPCC首先经过Clark变换，abc坐标系变为 坐标系的表达，然后以相位估计 进行Park变换，得到VPCC在q轴的表达如公式（15）所示：

[0093] ；（15）

[0094] 其中， ， ，E1表示 坐标系下电压幅值，表示相角误差；

[0095] 根据工程经验，PCC点电网电压在abc坐标系的参数设定如表1所示。

[0096] 表1‑输入信号参数设置

[0097]

[0098] 上述V1表示电压理想状态（即无干扰）的值；(1+ )V1、(1+ )V1表示电压的不平衡分量，V5、V7、V11、V13对应表示5、7、11、13次谐波分量。

[0099] 则通过上述公式（15），可以得到dq坐标的参数如表2所示。

[0100] 表2‑Vq轴扰动分量

[0101]

[0102] 根据上述Vq轴扰动分量V2、V6、V12，级联自适应陷波器由3个自适应陷波滤波器串联构成，分别滤除2次谐波，6次谐波和12次谐波扰动分量。具体结构如图4所示。为了滤除对应的扰动分量，需要调节自适应陷波滤波器的中心频率，分别对应2次，6次和12次谐波、频带宽度（BW2、BW6、BW12）和迭代因子（ 、 、 ）。

[0103] 在本发明实例中，自适应陷波器使用的自适应递归算法可以在没有参考估计频率的情况下跟踪频率变化，因此，可以提高自适应陷波器的品质因数，有效地提高滤波能力。对于固定陷波器，必须调整较小的Q值以增加滤波器带宽并减少频率变化的干扰。系统相位延迟与自适应陷波器的参数和Q值有关。因此，必须设计自适应陷波器的参和Q以获得合适的相位延迟。在本发明实施例中调谐参数如表3所示。

[0104] 表3‑自适应陷波器的参数

[0105]

[0106] 为了验证上述自适应陷波器的参数是否合适，需要对系统参数进行仿真验证，确认系统参数是否满足锁相环的性能要求。

[0107] 其中，输入信号中叠加不平衡和谐波干扰分量见上述表1。频率的波动变化设定为47Hz‑53Hz，即输入信号频率分别设定为47Hz和53Hz，通过Vq轴上扰动分量的消除比例来确认级联自适应滤波器性能。根据上述表2所示的Vq轴中的扰动分量，47Hz和53Hz扰动滤除比例如下表4所示。

[0108] 表4

[0109]

[0110] 如上表所示，本发明实施例涉及的锁相环具有较好的扰动滤除性能，能够准确的反应输入信号的相位信息。

[0111] 相对于现有技术，本发明实施例提供一种包括dq变换模块，级联自适应陷波器，比例积分调节器(PI)、比较器、积分环节、比例环节的自适应频率估计的锁相环，通过将级联自适应陷波器嵌入传统的三相同步锁相环（SRF‑PLL）中，构成自适应频率估计的锁相环，通过级联自适应陷波滤波器可以自动跟踪电网电压频率的变化，同时消除非理想电网电压的扰动，包括三相不平衡产生的负序分量，谐波分量，能够快速准确估计电网电压频率波动。

[0112] 实施例二

[0113] 请参阅图5，本发明实施例提供了一种自适应频率估计的锁相方法，包括：

[0114] S1：对电网的三相电压进行dq变换，得到q轴电压分量；

[0115] S2：通过预先构造的级联自适应陷波器，对所述q轴电压分量上的扰动分量进行滤除；

[0116] S3：根据滤除扰动分量后的q轴电压分量，得到所述三相电压的相位估计。

[0117] 作为上述方案的改进，所述根据滤除扰动分量后的q轴电压分量，得到所述三相电压的相位估计，包括：

[0118] 将滤除扰动分量后的q轴电压分量送入比例积分调节器，得到q轴的角频率分量；

[0119] 所述q轴的角频率分量经比较器、积分环节、比例环节后，得到所述相位估计。

[0120] 作为上述方案的改进，所述对电网的三相电压进行dq变换，得到q轴电压分量，包括：

[0121] 对所述三相电压进行Clark变换，得到坐标系下的电压分量；

[0122] 对 坐标系下的电压分量进行Park变换，得到所述q轴电压分量。

[0123] 作为上述方案的改进，所述三相电压的函数表达为：

[0124] ；

[0125] 其中，V1表示电压幅值，表示电压相角，表示B相不平衡因子，（k表示数值1到无穷递增），表示C相不平衡因子；V6k+1表示谐波电压幅值；

[0126] 所述q轴电压分量的函数表达为：

[0127] ；

[0128] 其中， ， ；E1表示 坐标系下电压幅值，表示相角误差。

[0129] 在一种可选的实施例中，所述级联自适应陷波器，包括多个串联连接的用于滤除特定次扰动分量的自适应陷波器。

[0130] 进一步，所述自适应陷波器包括用于滤除2次扰动分量的第一自适应陷波器、用于滤除6次扰动分量的第二自适应陷波器以及用于滤除12次扰动分量的第三自适应陷波器；其中，所述第一自适应陷波器、第二自适应陷波器、第三自适应陷波器依次串联连接。

[0131] 进一步，所述自适应陷波器为二阶自适应陷波器；

[0132] 所述二阶自适应陷波器的自适应算法函数表达为：

[0133] ；

[0134] 其中， 表示状态变量， 表示需要滤除的扰动分量，表示自适应陷波器的迭代步长， 为自适应陷波器的回归因子， 表示输入信号， 表示输出信号的估计， 表示n次调节变量，a和r表示参数变量，n表示迭代次数。

[0135] 其中，自适应陷波器的参数设置以及工作原理请参见实施例一，在这里不再重复说明。

[0136] 相对于现有技术，本发明实施例的有益效果在于：通过对电网的三相电压进行dq变换，得到q轴电压分量；然后使用预先构造的级联自适应陷波器，对所述q轴电压分量上的扰动分量进行滤除；再根据滤除扰动分量后的q轴电压分量，得到所述三相电压的相位估计；通过在锁相环节引入级联自适应陷波器滤除电网电压的扰动分量，可以准确跟踪电压相位。

[0137] 需说明的是，以上所描述的装置实施例仅仅是示意性的，其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。另外，本发明提供的装置实施例附图中，模块之间的连接关系表示它们之间具有通信连接，具体可以实现为一条或多条通信总线或信号线。本领域普通技术人员在不付出创造性劳动的情况下，即可以理解并实施。

[0138] 以上所述是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出多台改进和润饰，这些改进和润饰也视为本发明的保护范围。