CN113675883A - 一种并网逆变器控制延时参数辨识的方法 - Google Patents

Abstract

本发明公开了一种并网逆变器控制延时参数辨识的方法，包括以下步骤：S1，采集并网逆变器频率响应数据，利用Loewner矩阵计算初始阶数和极点；S2，利用矢量匹配修正极点，依次计算留数、一次项、常数项；S3，计算拟合均方根误差百分比，并与给定误差阈值对比，若不小于给定误差阈值，则返回步骤S2，修正初始参数，增大拟合阶数；重复步骤S2形成迭代，直至小于给定误差阈值时结束计算，输出拟合结果；S4，基于步骤S3的拟合结果，根据拟合的实际并网逆变器输出阻抗传递函数，利用数值分析方法求解控制延时传递函数，进而辨识控制延时。本发明能够实现延时参数等内部信息未知的并网逆变器的控制延时参数辨识，快速精准，为谐波谐振抑制提供一定的参考。

Description

一种并网逆变器控制延时参数辨识的方法

技术领域

本发明涉及并网逆变器参数辨识技术领域，特别是涉及一种并网逆变器控制延时参数辨识的方法。

背景技术

随着电力电子化新能源电力系统的发展，作为新能源和大电网的重要接口，多种类、多工况运行的并网逆变器应用愈加广泛。并网逆变器普遍采用数字控制器，以计算延时和脉宽调制延时为代表的控制延时环节直接改变并网逆变器的dq输出阻抗特性，是诱发并网逆变器谐波不稳定问题的关键参数之一。受控制延时的影响，并网逆变器dq输出阻抗特性曲线会出现反向谐振峰。控制延时的大小也会对并网系统的谐振频率及阻尼系数产生不同程度的影响。工程中一般选择1.5个采样周期作为控制延时。但实际系统控制延时并非如此，其数值大小会受占空比更新模式、采样方式、硬件滤波等的不同而改变。考虑实际工程现场存在大量已运行内部信息保密的商用并网逆变器，针对于这类并网逆变器怎样辨识其控制延时参数，尚没有很有效的方法。

发明内容

本发明的目的是提供一种并网逆变器控制延时参数辨识的方法，能够实现延时参数等内部信息未知的并网逆变器的控制延时参数辨识，快速精准，对于工程现场并网系统的谐波谐振抑制设计具有一定的指导意义。

为实现上述目的，本发明提供了如下方案：

一种并网逆变器控制延时参数辨识的方法，该方法包括以下步骤：

S1，采集并网逆变器频率响应数据，利用Loewner矩阵计算初始阶数和极点；

S2，利用矢量匹配修正极点，依次计算留数、一次项、常数项；

S3，计算拟合均方根误差百分比，并与给定误差阈值对比，若不小于给定误差阈值，则返回步骤S2，修正初始参数，增大拟合阶数；重复步骤S2形成迭代，直至小于给定误差阈值时结束计算，输出拟合结果；

S4，基于步骤S3的拟合结果，根据拟合的实际并网逆变器输出阻抗传递函数，利用数值分析方法求解控制延时传递函数，进而辨识控制延时。

进一步的，所述步骤S1，采集并网逆变器频率响应数据，利用Loewner 矩阵计算初始阶数和极点，具体包括：

S101，在采集的并网逆变器频率响应数据的负频率处附加复共轭频域响应数据，产生两倍的数据点，采用向量格式切线插值策略拆分数据，即：

式中：s_k，γ_i，β_j为频率采样点的复频率，S(s_k)、S(γ_i)、S(β_j)为采样点的复共轭频域响应数据，

表示复共轭函数，k＝1,2,……，2n；i＝1,2,……，n； j＝1,2,……，n；x＝1,2,……，n/2；s_x和s_n+2i/2分别为频率采样点x和(n+2i)/2的复频率；

和

分别为复频率s_x和s_n+2i/2的复共轭频率；γ_2x和β_2x分别为数据集γ和β的2x处的复频率；γ_2x-1和β_2x-1分别为数据集γ和β的2x-1处的复频率；

S102，构造矩阵L、σL、F和W如下：

式中：[L_j,i]和[σL_j,i]是n×n矩阵L和σL的第j行第i列元素，j＝1,2,……， n；i＝1,2,……，n；式(2)中构造的矩阵是复数的，用相似变换来求解得到，对应实矩阵{L_r，σL_r，F_r，W_r}如式(3)所示：

式中：kron.()表示Kronecker积函数；eye(n/2)表示(n/2)×(n/2)的单位矩阵； v是块对角矩阵V的子块矩阵；

S103，对矩阵束yL_r-σL_r进行奇异值分解：

yL_r-σL_r＝ΛΣΨ^* (4)

式中：y为除(L_r,σL_r)特征值外的样本值，Σ是包含奇异值的对角矩阵，Λ、Ψ是正交矩阵，()^*表示复共轭转置，阶数

由归一化奇异值图确定，

是图中最大下降指标；

S104，通过式(4)计算奇异值，绘制归一化奇异值变化趋势图，最大下降点即为矢量匹配计算初始阶数m，由矩阵Λ和Ψ的前

个奇异值组成矩阵Λ_R和Ψ_R，计算矩阵束

的特征值即为VF算法的初始极点 q_k。

进一步的，所述步骤S2中，利用矢量匹配修正极点，依次计算留数、一次项、常数项，具体包括：

S201，基于Loewner矩阵计算得到的初始极点q_k和初始阶数

建立目标函数，如式(5)所示：

式中：c_k是留数；q_k是初始极点；d和e为实系数；f(s)为待拟合传递函数；s＝jw,w为角频率。

S202，通过最小二乘法解式(5)所示的线性方程组求留数

计算σ(s) 的零点{p_k}即为f(s)的新的极点，实现极点修正：

{p_k}＝eig(A-b·c^T) (6)

式中：A是包含初始极点{q_k}的对角矩阵，b是单位列向量，c^T是留数

的行向量，eig.()表示特征值函数；

S203，基于修正后的极点，替代先前极点，重复最小二乘求解线性方程组，实现留数、一次项、常数项的相继计算。

进一步的，所述步骤S3，计算拟合均方根误差百分比，并与给定误差阈值对比，若不小于给定误差阈值，则返回步骤S2，修正初始参数，增大拟合阶数；重复步骤S2形成迭代，直至小于给定误差阈值时结束计算，输出拟合结果，具体包括：

S301，建立相对误差百分比r_ms作为拟合结果误差计算标准，通过式(7) 进行计算，r_ms越接近于0，表明拟合误差越小，拟合精确度越好，

式中：N为采集样本数，s_k＝j*2πf_k；f_k表示第k个频率采样点的频率；f_fit(s_k) 表示第k个频率采样点的频率拟合值；f_mea(s_k)表示第k个频率采样点的频率测量值；

S302，计算一次迭代后的r_ms，给定误差阈值r_ms0设置为1*10^-13，不满足小于误差阈值的条件时，增大方法阶数，同时将本次迭代得到的函数参数作为初始值重新带入矢量匹配进行迭代，直至拟合指标小于给定误差阈值。

进一步的，所述步骤S4，基于步骤S3的拟合结果，根据拟合的实际并网逆变器输出阻抗传递函数，利用数值分析方法求解控制延时传递函数，进而辨识控制延时，具体包括：

S401，基于并网逆变器dq输出阻抗，推导并网逆变器控制延时传递函数如式(8)所示：

其中，

Z_odd0＝sL+V_dcG_c

Z_odq＝-ω₀L

式中：Z_odd0为不考虑控制延时逆变器的dd输出阻抗；Z_odd为实际系统存在控制延时逆变器的dd输出阻抗，即为拟合阻抗传递函数；Z_odq为不考虑控制延时逆变器的dq输出阻抗；V_dc为直流侧电压；L为滤波器电感；G_c为PI 控制器传递函数矩阵，s＝jw，w为角频率；ω₀为基波角频率；

S402，代入步骤3的并网逆变器阻抗拟合结果，求得并网逆变器控制延时传递函数G_del，将G_del用Pade近似表示为多项式形式，如式(9)所示

其中：

式中：T为实际控制延时；s＝jw，w为角频率；l和k分别为传递函数G_del分子和分母阶数，通常取l＝k；

整理化简式(8)，根据式(8)阶数确定Pade近似阶数，即：l和k的取值，计算式(9)即可进一步求解实际控制延时T具体数值。

根据本发明提供的具体实施例，本发明公开了以下技术效果：本发明提供的并网逆变器控制延时参数辨识的方法，适用并网逆变器内部信息保密的实际工程场景，以控制延时参数等内部信息未知的并网逆变器为对象，以控制延时参数辨识为目的，从阻抗拟合角度实现并网逆变器实际控制延时参数的快速精确辨识；针对传统矢量匹配的技术缺陷，本发明提出利用Loewner矩阵改进矢量匹配用于并网逆变器阻抗拟合，实现了矢量匹配初始参数的自动获取，缩短了计算耗时，可以满足实际工程现场对快速精确阻抗拟合的需求；基于本发明公开方法的计算结果可以进一步辨识系统谐波谐振点，为后续的谐波谐振抑制设计提供一定的参考，对于工程现场并网系统的谐波谐振抑制设计具有一定的指导意义。

附图说明

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

图1为本发明并网逆变器控制延时参数辨识的方法的流程图；

图2为本发明实施例三相并网逆变器拓扑图；

图3为本发明实施例归一化奇异值图；

图4为本发明实施例辨识控制延时G_del零极点图；

图5(a)为本发明实施例G_dd主导谐振零极点图；

图5(b)为本发明实施例G_dq主导谐振零极点图；

图6(a)为本发明实施例无背景谐波时PCC点处并网电流波形图；

图6(b)为本发明实施例无背景谐波时频谱图；

图6(c)为本发明实施例有背景谐波时PCC点处并网电流波形；

图6(d)为本发明实施例有背景谐波时频谱图。

具体实施方式

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

本发明的目的是提供一种并网逆变器控制延时参数辨识方法，能够实现延时参数等内部信息未知的并网逆变器的控制延时参数辨识，快速精准，对于工程现场并网系统的谐波谐振抑制设计具有一定的指导意义。

为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本发明作进一步详细的说明。

图2为一种LCL型并网逆变器典型拓扑，V_dc为直流侧电压，L为滤波器电感，L_g为电网等效电感；i_abc为dq旋转坐标系下逆变器侧输出电流，θ为公共耦合点PCC(point ofcommon coupling)处对交流电压v_pcc锁相产生的相位信号，电流环采用比例积分PI(proportion integration)控制器；i_{d_ref}和i_{q_ref}分别为d轴与q轴电流指令参考值。

如图1所示，本发明提供的并网逆变器控制延时参数辨识的方法，基于图2所示拓扑结构，该方法包括以下步骤：

S1，采集并网逆变器频率响应数据，利用Loewner矩阵计算初始阶数和极点；

S2，利用矢量匹配修正极点，依次计算留数、一次项、常数项；

S3，计算拟合均方根误差百分比，并与给定误差阈值对比，若不小于给定误差阈值，则返回步骤S2，修正初始参数，增大拟合阶数；重复步骤S2形成迭代，直至小于给定误差阈值时结束计算，输出拟合结果；

S4，基于步骤S3的拟合结果，根据拟合的实际并网逆变器输出阻抗传递函数，利用数值分析方法求解控制延时传递函数，进而辨识控制延时。

其中，所述步骤S1，采集并网逆变器频率响应数据，利用Loewner矩阵计算初始阶数和极点，具体包括：

S101，在采集的并网逆变器频率响应数据的负频率处附加复共轭频域响应数据，产生两倍的数据点，采用向量格式切线插值(Vector Format Tangential Interpolation，VFTI)策略拆分数据，即：

式中：s_k，γ_i，β_j为频率采样点的复频率，S(s_k)、S(γ_i)、S(β_j)为采样点的复共轭频域响应数据，

表示复共轭函数，k＝1,2,……，2n；i＝1,2,……，n； j＝1,2,……，n；x＝1,2,……，n/2；s_x和s_n+2i/2分别为频率采样点x和(n+2i)/2的复频率；

和

分别为复频率s_x和s_n+2i/2的复共轭频率；γ_2x和β_2x分别为数据集γ和β的2x处的复频率；γ_2x-1和β_2x-1分别为数据集γ和β的2x-1处的复频率。

S102，构造矩阵L、σL、F和W如下：

式中：[L_j,i]和[σL_j,i]是n×n矩阵L和σL的第j行第i列元素，j＝1,2,……， n；i＝1,2,……，n；式(2)中构造的矩阵是复数的，用相似变换来求解得到，对应实矩阵{L_r，σL_r，F_r，W_r}如式(3)所示：

式中：kron.()表示Kronecker积函数；eye(n/2)表示(n/2)×(n/2)的单位矩阵； v是块对角矩阵V的子块矩阵。

S103，对矩阵束yL_r-σL_r进行奇异值分解：

yL_r-σL_r＝ΛΣΨ^* (4)

式中：y为除(L_r,σL_r)特征值外的样本值，Σ是包含奇异值的对角矩阵，Λ、Ψ是正交矩阵，()^*表示复共轭转置，阶数

由归一化奇异值图确定，

是图中最大下降指标；

S104，通过式(4)计算奇异值，绘制归一化奇异值变化趋势图，如图3 所示，最大下降点即为矢量匹配计算初始阶数m，由矩阵Λ和Ψ的前

个奇异值组成矩阵Λ_R和Ψ_R，计算矩阵束

的特征值即为VF算法的初始极点q_k。

其中，所述步骤S2中，利用矢量匹配修正极点，依次计算留数、一次项、常数项，具体包括：

S201，基于Loewner矩阵计算得到的初始极点q_k和初始阶数

建立目标函数，如式(5)所示：

式中：c_k是留数；q_k是初始极点；d和e为实系数；f(s)为待拟合传递函数；s＝jw,w为角频率。

S202，通过最小二乘法解式(5)所示的线性方程组求留数

计算σ(s) 的零点{p_k}即为f(s)的新的极点，实现极点修正：

{p_k}＝eig(A-b·c^T) (6)

式中：A是包含初始极点{q_k}的对角矩阵，b是单位列向量，c^T是留数

的行向量，eig.()表示特征值函数；

S203，基于修正后的极点，替代先前极点，重复最小二乘求解线性方程组，实现留数、一次项、常数项的相继计算。

其中，所述步骤S3，计算拟合均方根误差百分比，并与给定误差阈值对比，若不小于给定误差阈值，则返回步骤S2，修正初始参数，增大拟合阶数；重复步骤S2形成迭代，直至小于给定误差阈值时结束计算，输出拟合结果，具体包括：

S301，建立相对误差百分比r_ms作为拟合结果误差计算标准，通过式(7) 进行计算，r_ms越接近于0，表明拟合误差越小，拟合精确度越好，

式中：N为采集样本数，s_k＝j*2πf_k；f_k表示第k个频率采样点的频率；f_fit(s_k) 表示第k个频率采样点的频率拟合值；f_mea(s_k)表示第k个频率采样点的频率测量值；

S302，计算一次迭代后的r_ms，给定误差阈值r_ms0设置为1*10^-13，不满足小于误差阈值的条件时，增大方法阶数，同时将本次迭代得到的函数参数作为初始值重新带入矢量匹配进行迭代，直至拟合指标小于给定误差阈值。

其中，所述步骤S4，基于步骤S3的拟合结果，根据拟合的实际并网逆变器输出阻抗传递函数，利用数值分析方法求解控制延时传递函数，进而辨识控制延时，具体包括：

S401，基于并网逆变器dq输出阻抗，推导并网逆变器控制延时传递函数如式(8)所示：

其中，

Z_odd0＝sL+V_dcG_c

Z_odq＝-ω₀L

式中：Z_odd0为不考虑控制延时逆变器的dd输出阻抗；Z_odd为实际系统存在控制延时逆变器的dd输出阻抗，即为拟合阻抗传递函数；Z_odq为不考虑控制延时逆变器的dq输出阻抗；V_dc为直流侧电压；L为滤波器电感；G_c为PI 控制器传递函数矩阵，s＝jw，w为角频率；ω₀为基波角频率。

S402，代入步骤3的并网逆变器阻抗拟合结果，求得并网逆变器控制延时传递函数G_del，将G_del用Pade近似表示为多项式形式，如式(9)所示

其中：

式中：T为实际控制延时；s＝jw，w为角频率；l和k分别为传递函数G_del分子和分母阶数，通常取l＝k；

整理化简式(8)，根据式(8)阶数确定Pade近似阶数，即：l和k的取值，计算式(9)即可进一步求解实际控制延时T具体数值。

计算式(8)得到实际并网逆变器的控制延时G_del传递函数，绘制其零极点图如图4所示。由图4分析可知，被测并网逆变器实际控制延时传递函数 G_del所对应的极点位于-1.18×10⁴-4.45×10^-14i(其他零极点均近似对消)，可以近似为-1.18×10⁴，故式(9)选择一阶Pade近似，即：l＝k＝1。计算式(9)可求解并网逆变器实际控制延时T＝169.5μs＝1.695T_s。基于阻抗判据，绘制并网系统主导谐振零极点图如图5(a)(b)所示，由图5(a)可知，对于系统阻抗G_dd，主导谐振极点位于-304±1.005×10⁴i和-147±1.030×10⁴i，谐振频率为1600Hz 和1639Hz。由图5(b)可知，对于系统阻抗G_dq，主导谐振极点位于-941±1.069×10⁴i和-783±1.090×10⁴i，谐振频率为1702Hz和1734Hz。经理论分析，受以上谐振频率处的电网背景谐波电压的激励，并网电流可能出现放大效应，并网系统产生谐波失稳现象。

为验证上述方法的有效性，基于RT-LAB搭建半实物实时仿真平台，开展如下实验：工况1：电网处于无背景谐波状态，工况1：电网处于无背景谐波状态，观察PCC点处逆变器并网电流波形，如图6(a)所示；对其进行FFT 分解，绘制频谱图如图6(b)所示。工况2：电网电压中出现1600Hz和1800Hz 背景谐波，谐波幅值为额定电压的1％，观察PCC点处逆变器并网电流波形，如图6(c)所示；对其进行FFT分解，绘制频谱图如图6(d)所示。

电网中无背景谐波时，从图6(a)中可以看出：并网逆变器并网电流波形未出现明显畸变；从图6(b)中可以看出：并网逆变器会向电网注入宽频域谐波电流，但谐波电流含量普遍较低。当电网中出现1600Hz和1800Hz背景谐波电压时，从图6(c)中可以看出：并网逆变器并网电流波形出现明显畸变；从图6(d)中可以看出：并网逆变器并网电流中的1600Hz谐波被显著放大，谐波失稳现象产生。经计算，无背景谐波电压时，并网逆变器并网电流谐波总畸变率(Total Harmonic Distortion，THD)为0.42％，存在14次背景谐波时，THD变为3.49％。对比图6(c)、图6(d)可以发现：并网逆变器并网电流THD的显著增加主要由并网电流中被放大的1600Hz谐波导致。上述实验结果恰恰证明了控制延时参数辨识的有效性，且与图5(a)分析结果几乎吻合。

本发明提供的并网逆变器控制延时参数辨识的方法，适用并网逆变器内部信息保密的实际工程场景，以控制延时参数等内部信息未知的并网逆变器为对象，以控制延时参数辨识为目的，从阻抗拟合角度实现并网逆变器实际控制延时参数的快速精确辨识；针对传统矢量匹配的技术缺陷，本发明提出利用Loewner矩阵改进矢量匹配用于并网逆变器阻抗拟合，实现了矢量匹配初始参数的自动获取，缩短了计算耗时，可以满足实际工程现场对快速精确阻抗拟合的需求；基于本发明公开方法的计算结果可以进一步辨识系统谐波谐振点，为后续的谐波谐振抑制设计提供一定的参考，对于工程现场并网系统的谐波谐振抑制设计具有一定的指导意义。

本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。对于实施例公开的系统而言，由于其与实施例公开的方法相对应，所以描述的比较简单，相关之处参见方法部分说明即可。

本申请中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处。综上所述，本说明书内容不应理解为对本发明的限制。

Claims (5)

1.一种并网逆变器控制延时参数辨识的方法，其特征在于，包括以下步骤：

S1，采集并网逆变器频率响应数据，利用Loewner矩阵计算初始阶数和极点；

S2，利用矢量匹配修正极点，依次计算留数、一次项、常数项；

S3，计算拟合均方根误差百分比，并与给定误差阈值对比，若不小于给定误差阈值，则返回步骤S2，修正初始参数，增大拟合阶数；重复步骤S2形成迭代，直至小于给定误差阈值时结束计算，输出拟合结果；

S4，基于步骤S3的拟合结果，根据拟合的实际并网逆变器输出阻抗传递函数，利用数值分析方法求解控制延时传递函数，进而辨识控制延时。

2.根据权利要求1所述的并网逆变器控制延时参数辨识的方法，其特征在于，所述步骤S1，采集并网逆变器频率响应数据，利用Loewner矩阵计算初始阶数和极点，具体包括：

S101，在采集的并网逆变器频率响应数据的负频率处附加复共轭频域响应数据，产生两倍的数据点，采用向量格式切线插值策略拆分数据，即：

式中：s_k，γ_i，β_j为频率采样点的复频率，S(s_k)、S(γ_i)、S(β_j)为采样点的复共轭频域响应数据，

表示复共轭函数，k＝1,2,……，2n；i＝1,2,……，n；j＝1,2,……，n；x＝1,2,……，n/2；s_x和s_n+2i/2分别为频率采样点x和(n+2i)/2的复频率；

和

分别为复频率s_x和s_n+2i/2的复共轭频率；γ_2x和β_2x分别为数据集γ和β的2x处的复频率；γ_2x-1和β_2x-1分别为数据集γ和β的2x-1处的复频率。

S102，构造矩阵L、σL、F和W如下：

式中：[L_j,i]和[σL_j,i]是n×n矩阵L和σL的第j行第i列元素，j＝1,2,……，n；i＝1,2,……，n；式(2)中构造的矩阵是复数的，用相似变换来求解得到对应实矩阵{L_r，σL_r，F_r，W_r}如式(3)所示：

式中：kron.()表示Kronecker积函数；eye(n/2)表示(n/2)×(n/2)的单位矩阵；v是块对角矩阵V的子块矩阵；

S103，对矩阵束yL_r-σL_r进行奇异值分解：

yL_r-σL_r＝ΛΣΨ^* (4)

式中：y为除(L_r,σL_r)特征值外的样本值，Σ是包含奇异值的对角矩阵，Λ、Ψ是正交矩阵，()^*表示复共轭转置，阶数

由归一化奇异值图确定，

是图中最大下降指标；

S104，通过式(4)计算奇异值，绘制归一化奇异值变化趋势图，最大下降点即为矢量匹配计算初始阶数m，由矩阵Λ和Ψ的前

个奇异值组成矩阵Λ_R和Ψ_R，计算矩阵束

的特征值即为VF算法的初始极点q_k。

3.根据权利要求2所述的并网逆变器控制延时参数辨识的方法，其特征在于，所述步骤S2中，利用矢量匹配修正极点，依次计算留数、一次项、常数项，具体包括：

S201，基于Loewner矩阵计算得到的初始极点q_k和初始阶数

建立目标函数，如式(5)所示：

式中：c_k是留数；q_k是初始极点；d和e为实系数；f(s)为待拟合传递函数；s＝jw,w为角频率；

S202，通过最小二乘法解式(5)所示的线性方程组求留数

计算σ(s)的零点{p_k}，即为f(s)的新的极点，实现极点修正：

{p_k}＝eig(A-b·c^T) (6)

式中：A是包含初始极点{q_k}的对角矩阵，b是单位列向量，c^T是留数

的行向量，eig.()表示特征值函数；

S203，基于修正后的极点，替代先前极点，重复最小二乘求解线性方程组，实现留数、一次项、常数项的相继计算。

4.根据权利要求3所述的并网逆变器控制延时参数辨识的方法，其特征在于，所述步骤S3，计算拟合均方根误差百分比，并与给定误差阈值对比，若不小于给定误差阈值，则返回步骤S2，修正初始参数，增大拟合阶数；重复步骤S2形成迭代，直至小于给定误差阈值时结束计算，输出拟合结果，具体包括：

S301，建立相对误差百分比r_ms作为拟合结果误差计算标准，通过式(7)进行计算，r_ms越接近于0，表明拟合误差越小，拟合精确度越好，

式中：N为采集样本数，s_k＝j*2πf_k；f_k表示第k个频率采样点的频率；f_fit(s_k)表示第k个频率采样点的频率拟合值；f_mea(s_k)表示第k个频率采样点的频率测量值；

S302，计算一次迭代后的r_ms，给定误差阈值r_ms0设置为1*10^-13，不满足小于误差阈值的条件时，增大方法阶数，同时将本次迭代得到的函数参数作为初始值重新带入矢量匹配进行迭代，直至拟合指标小于给定误差阈值。

5.根据权利要求4所述的并网逆变器控制延时参数辨识的方法，其特征在于，所述步骤S4，基于步骤S3的拟合结果，根据拟合的实际并网逆变器输出阻抗传递函数，利用数值分析方法求解控制延时传递函数，进而辨识控制延时，具体包括：

S401，基于并网逆变器dq输出阻抗，推导并网逆变器控制延时传递函数如式(8)所示：

其中，

Z_odd0＝sL+V_dcG_c

Z_odq＝-ω₀L

式中：Z_odd0为不考虑控制延时逆变器的dd输出阻抗；Z_odd为实际系统存在控制延时逆变器的dd输出阻抗，即为拟合阻抗传递函数；Z_odq为不考虑控制延时逆变器的dq输出阻抗；V_dc为直流侧电压；L为滤波器电感；G_c为PI控制器传递函数矩阵，s＝jw，w为角频率；ω₀为基波角频率。

S402，代入步骤3的并网逆变器阻抗拟合结果，求得并网逆变器控制延时传递函数G_del，将G_del用Pade近似表示为多项式形式，如式(9)所示

其中：

式中：T为实际控制延时；s＝jw，w为角频率；l和k分别为传递函数G_del分子和分母阶数，通常取l＝k；

整理化简式(8)，根据式(8)阶数确定Pade近似阶数，即：l和k的取值，计算式(9)即可进一步求解实际控制延时T具体数值。

Images