基于MATLAB/Simulink的三相并网逆变器dq阻抗建模及扫频仿真

本文主要介绍三相并网逆变器dq阻抗建模的相关知识，和大家分享一下怎么使用MATLAB/Simulink来进行理论模型的搭建以及如何通过扫频获取阻抗模型，一方面是给自己留点学习资料，也希望对大家能有所帮助。如果有什么问题或者错误，欢迎留言。

本文涉及的代码和模型就不提供啦，不过后面会给出一组参数和对应阻抗的Bode图，有兴趣的可以自己做一下然后对照看看有没有问题。

整体系统介绍

图1是本文所使用的系统：三相并网逆变器；输入认为是恒定的直流电压源；滤波采用单电感，同时考虑寄生电阻；控制在dq坐标系下实现，为简单起见，使用单电流环。

图1 三相并网变换器系统框图

理论模型

理论模型不是本文叙述的重点，大家可以去找相关的论文，有比较详细的介绍（特别是一些学位论文），这里只简单提一下。

并网变换器的dq阻抗模型实际上就是一个2*2的传递函数矩阵，输入和输出分别为并网点的电压和电流

(1)

建模中，个人认为最重要的就是锁相环。由于锁相环（PLL）的存在，在有小信号扰动时，控制系统的dq坐标系和实际系统的dq坐标系之间存在着相位差Δθ。Δθ的存在引起了dq轴的耦合，使得系统产生了复杂的稳定性问题。

为此，我们需要找出主电路（上标为 s）和控制系统（上标为 c ）两个dq坐标系之间的关系。也就是两个坐标系下的电气量之间的转换关系式。根据相关的文献，有

(2)

其中 x 表示电压/电流/占空比等。

MATLAB实现

获取dq阻抗模型的控制框图或表达式后，我们可以考虑在MATLAB中用代码实现，并画出对应的Bode图。

要注意的几个点是：

1. 不需要写出各个阻抗的具体表达式，直接以矩阵的形式来进行运算；
2. 求逆矩阵时，MATLAB里面不推荐使用 inv() 函数，而是用 / 或 \ 予以替代。具体来说

(3)

3. 似乎没办法直接作出一个2*2传递函数矩阵的Bode图（也可能是我不会），所以需要分别运用 bode() 函数。另外，推荐大伙使用这么几行命令，让每次画出的图的范围都是固定的。

tt = bodeoptions;
tt.XLim={[1 10e3]}; % 横坐标范围
tt.FreqUnits='Hz'; % 设置值模式
tt.PhaseWrapping='on'; % 相频特性 限制范围
tt.PhaseWrappingBranch=-180; % ±180之间

bode(Zdd,tt);

基于Simulink的阻抗扫频仿真

整体思路

这里 阻抗扫频仿真 的意思是，我们先在Simulink中搭建电路模型，然后往PCC点注入一定的扰动信号（这里为串联电压扰动信号源），并采集对应的PCC电流，而后对电压电流信号进行处理运算，求解出系统的阻抗模型。

由于待求的阻抗有四个，但是一组扰动信号只有d轴和q轴两组方程，所以需要设计两组线性无关的扰动信号，才能建立起四个方程。

(4)

而这只是一个频率下的阻抗，还需要根据情况设计多个扰动频率点，以勾勒出一段频率范围内的阻抗特性。

整体的思路是：

1. 搭建电路模型，并且保证其可以正常工作（具有稳态工作点）；
2. 设计电压扰动信号1和电压扰动信号2；
3. 把某一频率下的信号1注入系统，采集PCC处的电压和电流，并进行dq变换；
4. 注入扰动信号2，并且采集PCC处的电压和电流；
5. 利用FFT，求出相应扰动频率下的两组信号的电压分量和电流分量的幅值和相位；
6. 根据前面的公式，得到阻抗矩阵的具体数值；
7. 注入另一频率下的扰动信号，并重复上述步骤，直至获得所有频率下的阻抗解；
8. 画图。

注意事项

这里针对上面的各个步骤介绍一些可能要加以注意的事项

1. 搭建电路模型，这一点没啥好说的，你首先需要确保自己的系统有稳态工作点，再考虑注入小信号扰动的事情；
2. 设计扰动信号，对于逆变器，一般采用串联电压源的方式，而两组线性无关的信号可以这样设计：

(5)

其中 w_g 是电网电压角频率。

3. 注入信号，采集信号，为了确保准确性，可以每次注入一个频率的一组扰动信号。采集的时候因为是dq阻抗，所以要进行dq变换。这里要注意θ的选取，需要变换到实际系统dq坐标系下（上标为 s ）；
4. 利用FFT求出幅值和相位，MATLAB中有对应的函数 fft() ，不过还需要自己补充一些内容才能符合需求；
5. 求解阻抗方程，如果需要导纳模型，同样也可以求解；
6. 注入其他频率的扰动信号，一般要避开工频（50Hz）以及其他谐波频率，免得计算结果有误。

流程框图

如果每次注入信号都手动操作的话，未免过于繁琐和费时（简直非人哉），所以强烈建议大家把这一个过程自动化。

得益于强大的MATLAB，我们最终可以实现这样的效果：在设计好系统和扰动频率点后，待程序运行结束后便可以自动生成理论模型和扫频仿真的Bode图。

个人的思路如下，供大家参考。

在Simulink中搭建模型，扰动信号的频率、幅值等将从MATLAB中获得（同理，也可以设置其他的模块参数从MATLAB中获得）；这一功能利用如下函数实现

% 载入仿真模型文件
load_system('xxx.slx');

% 设定参数
set_param('xxx/Constant','Value','3');

% 运行仿真
sim('xxx', [t_start, t_end]);

具体可以参考这篇文章：MATLAB脚本调用simulink仿真文件及simulink模块参数修改-load_system函数-sim函数-set_param函数

另一方面，采集到的电压电流信号也要传到MATLAB中，以便于进一步的处理。这里可以使用 To Workspace 或 To File 模块；

点击 运行 后，首先启动Simulink模型，并选择注入扰动信号1（这里要设置好合适的仿真时间）。仿真结束后，把数据先存放到一个数组里，然后再注入扰动信号2，并且也给它存起来。获得两组数据后，利用FFT求得幅值和相位，利用阻抗表达式求得这一扰动频率下的阻抗解。为了获得其他频率下的阻抗，可以设计一个扰动频率数组，利用循环语句每次注入不同的信号。最后，把得到的阻抗解利用 plot() 给它画出来，并利用 hold on ， 对比理论模型和扫频模型。

图2 流程框图

其他

这里提供了一组参数和阻抗Bode图，具体系统框图如前面的图1所示，供大家参考。

描述	参数	数值
直流电压	Udc	1000V
电网电压有效值	Ug	220V
有功功率	P	100kW
无功功率	Q	20kVar
滤波电感	L	1mH
滤波电感寄生电阻	R	1mΩ
锁相环比例系数	Kp_PLL	1
锁相环积分系数	Ki_PLL	100
电流环比例系数	Kp_i	0.01
电流环积分系数	Ki_i	0.1

图3 dq阻抗模型Bode图

其中 蓝色线 表示理论模型，红色线 表示扫频仿真模型。
2023.08.23

就个人的理解，扰动信号的注入时间不能太短，同时还要满足一定的要求，否则可能影响FFT的计算结果。举例来说，如果扰动信号的频率为2Hz，那么用来进行分析的这段信号至少要为(1/2)s，也就是0.5s。但是如果我们选0.6s，似乎也不太行(忘记有没有试过），因为在0.6s不是2Hz(0.5s)的整数倍，波形是不完整的。

如果我们选取0.5s为用于分析的信号时长，那么这段信号的分辨率就是2Hz，所以在选取扰动信号的频率时，最好都是2Hz的整数倍，也就是偶数。

此外，不能直接把仿真一开始的一段信号拿去分析，因为此时系统还未稳定。

主要参考文献：

1.曾志杰，并网逆变器DQ阻抗模型的建立方法研究，东南大学，2021.
2.李杨，多变换器系统小信号稳定性分析与提升方法研究，湖南大学，2021.
3.黄耀，新能源并网逆变器序阻抗建模及其阻抗测量研究，合肥工业大学，2020.
4.刘家源，新能源并网逆变器的阻抗测量方法研究与装备研制，湖南大学，2019.
5.伍文华，新能源发电接入弱电网的宽频带振荡机理及抑制方法研究，湖南大学，2019.
2023.10.08

在建模时，一开始没必要做逆变器，可以先做无源负载，比如电网接一个RL负载，测量该电路的阻抗特性，这样会简单很多；

如果测量逆变器的阻抗特性，也没必要一上来就是完全体，比如很复杂的滤波电路和控制策略，可以先做本文所用的这样一个简单的电路和控制；

当采集到PCC处的三相电压/电流并打算将其转换到dq坐标系时，Park变换的相角不能和控制环一样，用锁相环出来的相角，因为那就跑到控制环下的dq坐标系了。可以利用电网电压建一个新的锁相环，这样得到的相角没有动态误差。
2023.12.19

这里统一回复一下，本文涉及的代码和模型也许以后会公开，但是目前不太方便。有偿也不考虑哈，谢谢。
2024.01.06