虚拟同步发电机VSG仿真：光储并网模型搭建与参数整定全攻略

这两年做光储并网仿真，我越来越强烈地感觉到一个矛盾：电网里的同步发电机占比越来越低，电力电子变流器占比越来越高，系统频率特性正在发生根本性变化。同步发电机自带一个巨大的机械转子，本质上就是一个大飞轮，在功率波动时能天然提供惯量支撑，给一次调频争取时间和缓冲空间。而光伏和储能接口的逆变器是静止设备，没有任何旋转部件，传统PQ控制策略下它们对电网频率几乎没有任何支撑能力，甚至会因为响应太快而加剧频率恶化。这个问题的解决方案之一就是虚拟同步发电机——用控制算法在软件层面"虚拟"出同步发电机的转子和励磁绕组特性，让逆变器从外特性上看起来像一台真实的同步发电机。

这篇文章要聊的就是基于Matlab/Simulink搭建一套完整的"光伏储能+虚拟同步发电机+并网"仿真模型。我在实际项目里用这套模型做过并网暂态分析、调频调压性能验证和参数整定测试，整个过程踩了不少坑，也总结了一套可以复用的搭建思路。如果你正在做新能源并网控制研究、微电网/储能变流器控制，或者写论文需要一套能跑的仿真模型做验证，这篇文章应该能帮你省下大量对着文档摸索的时间。

1. 光储并网的核心矛盾：逆变器为什么要有"机械惯性"

1.1 高比例新能源并网后的频率支撑困境

先聊一个基本现象。传统电力系统里，同步发电机的转子具有真实的转动惯量，当系统出现有功缺额时，发电机转子会释放或吸收动能，减缓频率变化速度。这个"减缓变化速度"的能力叫惯量支撑，它和一次调频、二次调频是三个不同层级的概念——惯量支撑是毫秒到秒级瞬时响应，一次调频是秒级到十秒级的功率调节，二次调频则是分钟级的自动发电控制恢复。

到了高比例新能源系统里，问题就出现了。光伏并网通过逆变器，储能系统通过变流器，两者都没有旋转部件，系统里"看得见的飞轮"越来越少了。同样是100MW的负荷突变，过去同步机组占比高的时候频率变化率可能只有0.1Hz/s，现在可能直接飙到0.5Hz/s以上。频率变化率太快，继电保护里的低频减载装置可能来不及正确动作，光伏逆变器自身的频率保护也可能误触发，导致连锁脱网。

这个问题的根源不在于"新能源发不出电"，而在于逆变器和电网之间的接口特性。传统PQ控制的逆变器只有电流源特性，它锁相跟踪电网电压，按照给定有功无功指令输出电流，电网频率变化对它来说只是一个"观测信号"，并不会主动改变输出功率来支撑频率。换句话说，它不是一个参与者，只是一个跟随者。

1.2 VSG的解决思路：给逆变器装一个虚拟飞轮

虚拟同步发电机（Virtual Synchronous Generator，VSG）的思路很直接：既然逆变器没有物理转子，那我就用控制算法模拟一个转子出来。核心是一个转子运动方程：

J * dω/dt = Pm - Pe - D * (ω - ωn)

这个方程和同步发电机的摇摆方程在数学形式上是完全一致的。J是虚拟惯量，D是虚拟阻尼，Pm是虚拟机械功率（由有功指令和调频环节产生），Pe是逆变器实际输出的电磁功率，ω是虚拟角速度，ωn是额定角速度。把ω积分就得到虚拟功角θ，也就是逆变器输出电压相位的参考值。

这个"虚拟飞轮"的效果是：当电网频率瞬间跌落时，VSG感受到ω和ωn的偏差，通过阻尼项D*(ω-ωn)自动调整输出功率，同时惯量项J*dω/dt让相位和功率变化过程呈现出与同步发电机类似的惯性特征。电网侧看到的不是一个"瞬间变化功率"的电流源，而是一个"缓缓变化、有惯性、有阻尼"的电压源。

同样的道理，无功-电压控制也可以用类似方法实现。同步发电机的励磁调节器根据机端电压偏差调整励磁电流，从而改变无功输出。VSG里就用一个无功-电压下垂方程来模拟这个行为：E = E_ref + n*(Q_ref - Q)。E是虚拟内电势幅值，n是无功下垂系数。这样逆变器在外特性上同时具备了有功-频率和无功-电压两个维度的自我调节能力。

1.3 这套模型能做什么、适合谁来用

这套光储VSG并网仿真模型，核心覆盖了从直流侧光伏发电、储能稳压，到交流侧虚拟同步机控制、并网运行的全链路。具体能做的实验包括：

• 并网稳态运行：验证VSG控制下的有功/无功输出与电网交换功率是否跟踪指令。
• 电网频率扰动响应：在电网侧注入频率阶跃或缓慢波动，观察VSG输出的有功功率变化是否具备惯量支撑特征。
• 负载突变/有功指令阶跃：验证虚拟惯量和阻尼系数对暂态振荡的抑制效果。
• 无功电压调节：改变无功指令，观察并网点电压和无功功率的跟随情况。
• 光伏功率波动场景：光照突变导致光伏输出波动，储能通过直流母线稳压和VSG并网控制平滑功率送入电网。
• 并网与离网切换/预同步：验证预同步控制在并网瞬间对冲击电流的抑制。

适合的读者很明确：做并网逆变器/储能变流器控制的研究生和工程师，需要做VSG控制算法验证的科研人员，以及想快速搭一套可复现仿真模型用于课程设计或论文佐证的学生。对Simulink基础要求不需要太高，但至少要知道Simscape Electrical的基本模块操作和S-function的基本概念——这些在后面的章节里我会逐个讲清楚。

2. VSG控制核心：转子运动方程、调频调压与双闭环

2.1 有功-频率控制：虚拟调速器与虚拟转子

VSG控制器内部可以拆成三个层级：功率外环（模拟同步发电机的调速器和励磁调节器）、虚拟转子运动（模拟机械惯量）、电压电流双闭环（模拟定子绕组与输出电压控制）。

功率外环的第一部分是虚拟调速器，它的作用是建立有功-频率下垂关系。同步发电机的原动机带调速器，当系统频率降低时调速器会增加进汽/进水量来增加机械功率。VSG里用这样一个方程模拟：

Pm = P_ref + Kf * (fn - f)

其中P_ref是有功功率指令，f是当前频率（通常由PLL锁相得到，或直接使用VSG内部的虚拟频率ω），fn是额定频率，Kf是频率响应系数，也叫有功调频系数。这个Kf与后面要讲的下垂系数m互为倒数关系。实际工程里常用百分数表示下垂：2%~5%的频率偏差对应100%的功率变化，即Kf = Pn / (2%~5%的fn)。

虚拟转子运动方程我已经在1.2节写过了，这里重点说它和调速器指令之间的关系。机械功率Pm由调速器产生，电磁功率Pe由主电路实际输出的有功功率反馈回来。两者之差驱动dω/dt变化，进而改变虚拟角速度ω，再积分得到相位θ。注意，实际的同步发电机转子方程里是转矩平衡 Jdω/dt = Tm - Te - D(ω - ωg)，其中Tm和Te是机械转矩和电磁转矩，功率与转矩相差一个ω。在VSG的一些实现里为了简化，直接把功率当作转矩用力矩平衡的方式处理，这在ω接近额定值的情况下误差很小，工程上完全能接受。但在推导小信号模型时要注意区分，不然参数整定公式会对不上。

2.2 无功-电压控制：模拟励磁调节器

同步发电机励磁系统的核心功能是维持机端电压并参与无功分配。VSG的无功-电压外环模拟了这个特性：

E = E_ref + n * (Q_ref - Q)

E_ref是空载电势参考值，Q_ref是无功功率指令，Q是实际无功输出，n是无功下垂系数。当实际无功大于指令时，内电势幅值会适当降低，自动减少无功输出；当系统电压跌落需要更多无功支撑时，电压偏差项会通过无功调节改变E，提高无功输出。

还有一种常见做法是引入积分环节实现无差调节：E = E_ref + n*(Q_ref - Q) + Ki*∫(Q_ref - Q)dt。这种方式在孤岛或微网模式下可以精确跟踪无功指令，但并网模式下可能会和电网电压调节相互作用，注意权衡使用。我的模型里保留的是带小比例积分的版本，既兼顾稳态精度又避免积分饱和导致电压越限。

2.3 电压电流双闭环与调制信号生成

VSG功率外环输出的是同步发电机的虚拟内电势幅值E和相位θ，但是逆变器不能直接输出电压源，它只能通过开关管占空比来逼近期望的电压波形。所以需要把E∠θ转换为三相电压参考信号，然后用电压电流双闭环来控制逆变器的输出电压跟踪这个参考。

电压外环的输出是电流内环的参考值，电流内环的输出被限制在逆变器允许的最大电流以内以避免过流损坏。电压外环通常采用PR（比例谐振）调节器或者dq坐标系下的PI调节器。PR调节器在静止坐标系下对正弦参考可以做到无静差跟踪，实现简单，无需坐标变换；PI调节器在dq旋转坐标系下效果也很好，但需要PLL或VSG内部参考角做Park变换。

电流内环输出的是电压调制波指令，经过SVPWM调制生成IGBT的开关脉冲。SVPWM相比SPWM的优点是直流电压利用率更高、谐波含量更低，对于三相三线系统这是最常用的选择。VSG中电压电流双闭环的带宽要求：电流内环要快到能跟随开关频率变化，一般是开关频率的1/10到1/5；电压外环比电流内环慢一个数量级，保证两者解耦，不互相干扰。

2.4 用S-function实现VSG核心状态方程

VSG控制核心——包括虚拟调速器、转子运动方程、无功-电压方程——是一个强耦合的动态系统，在Simulink里直接画框图当然可以做，但参数修改和复用性很差。更专业的做法是用S-function把这部分封装成MEX文件或Level-2 MATLAB S-function，把状态方程写在代码里。

以Level-1 C MEX S-function为例，核心部分就是四个状态量：Δω（相对于额定角速度的偏差）、θ（相位）、积分项（无功无差调节用）、以及滤波后的功率值（如果功率滤波状态放在这里做）。微分方程就是转子运动方程加相位微分方程：

d(Δω)/dt = (Pm - Pe - D*Δω) / J
dθ/dt = ωn + Δω

这里单独解释一下为什么相位微分里要加回额定角速度ωn。因为ωn=314.159rad/s，如果没有这项，相位就不会正常旋转，整个电压参考信号就是静止的，完全错乱。很多初学VSG的人仿真一跑就发疯，最后发现就是这个ωn漏加了。我当年也在这上面浪费过半天。

如果是在S-function Builder里写，可以在mdlDerivatives里填连续状态方程，在mdlOutputs里输出状态变量。但是注意，如果powergui选择的是离散仿真模式，S-function的连续状态应该改成离散差分方程，否则求解器匹配会出问题。这一段在5.3节再展开讲。

3. Simulink整体模型架构：四个子系统逐一拆解

3.1 系统拓扑总览与信号流向

整个模型的顶层拓扑可以这样设计：光伏阵列经过Boost变换器接入直流母线，储能电池通过双向Buck-Boost变换器也接入直流母线，直流母线再经过三相两电平逆变器变成交流，通过LCL滤波器、并网开关和公共连接点接到三相电网。

信号流方向是：光伏Boost用MPPT控制，目标是让光伏工作在最大功率点，它只负责向直流母线"推"功率。储能双向DC-DC的目的是稳定直流母线电压——光伏功率供大于求时储能吸收多余功率，光伏功率不足时储能释放功率补缺口。直流母线电压稳定后，逆变器作为前级"电源"把直流功率送往交流侧，VSG控制根据并网点的电压电流算出有功无功，再根据转子运动方程决定输出电压的相位和幅值，最后通过PWM驱动开关管。

在Simulink里，我把整个模型拆成四个子系统来管理，避免顶层布线乱成一团。

3.2 子系统一：光伏阵列与MPPT

光伏部分我直接用Simscape Electrical自带的PV Array模块，选型号时可以选厂商参数，也可以自己设定开路电压、短路电流、最大功率点电压电流。关键是MPPT算法要自己搭：最常用的是扰动观察法，每一步给工作电压一个小扰动，比较扰动前后功率变化方向，决定下一步扰动的方向。Simulink里可以用MATLAB Function写，也可以用S-function，或者直接用基本数学模块搭逻辑。

Boost变换器把光伏电压升到直流母线电压水平。MPPT的输出是参考电压Vpv_ref，电压外环PI调节器比较Vpv和Vpv_ref，输出电流内环的参考值，电流内环再输出PWM占空比。这里有个细节：Boost电感要按电流纹波要求设计，一般取连续工作模式，电感值L = Vpv * (Vdc - Vpv) / (fsw * ΔI * Vdc)，其中fsw是开关频率，ΔI是允许的电流纹波。

我建议在光伏子系统里把MPPT的使能开关预留出来。做实验时经常要验证"光伏功率突变"的场景，直接给一个阶跃扰动修改光照辐射度，比临时改控制逻辑方便得多。

3.3 子系统二：储能电池与双向DC-DC稳压

储能电池用Simscape Electrical的Battery模块，可以选标准预置模型，设置额定容量Ah、额定电压、SOC初始值。双向Buck-Boost变换器是两个IGBT/二极管对组成的半桥结构，通过控制两个开关管的互补导通实现功率双向流动。

控制目标是维持直流母线电压。控制方法用直流母线电压外环，比较Vdc和Vdc_ref，输出电流内环参考值，内环控制储能电感电流。当Vdc高于参考时，外环输出负的电流参考，储能进入充电模式，把母线多余功率吸收掉；当Vdc低于参考时，外环输出正的电流参考，储能放电补足缺额。

这个环节有两个容易出问题的地方。第一是PI参数的积分限幅要设置好，否则直流母线电压会出现低频振荡，并网功率也会跟着波动。第二是Battery模块本身带SOC上下限，SOC接近边界时内阻变化大，母线电压可能撑不住，建议仿真场景里SOC初始值设在50%~80%之间，留出充分的工作区间。我在5.5节会再细说SOC和母线电压耦合的问题。

3.4 子系统三：VSG控制器与PWM生成

这是整个模型的核心子系统。输入信号包括：并网点三相电压Va/Vb/Vc、逆变器输出三相电流Ia/Ib/Ic、有功功率指令Pref、无功功率指令Qref，以及一些控制参数（J、D、Kf、n等）。输出是逆变器的PWM脉冲。

VSG控制器的数据处理流程：

第一步，用瞬时功率法或者dq变换计算有功和无功。瞬时功率法最简单：P = VaIa + VbIb + VcIc，Q = (VbcIa + VcaIb + VabIc)/√3，注意这里用的是相电流和线电压的瞬时值乘积组合，不需要坐标变换。算出来的P和Q含有高次谐波，必须经过低通滤波再送进VSG功率外环，滤波时间常数一般取1~10ms，太大会让惯量响应变慢，太小则滤波不干净会导致调制波畸变。

第二步，按2.1和2.2节的方程计算虚拟机械功率Pm和虚拟内电势幅值E。转子运动方程求出θ，电压参考信号就是：

Va_ref = E * sin(θ)
Vb_ref = E * sin(θ - 2π/3)
Vc_ref = E * sin(θ + 2π/3)

第三步，电压参考送进电压电流双闭环。在dq坐标系下，电压外环的PI调节器输出电流参考，电流内环PI调节器输出电压指令，再经过Park反变换得到三相调制波。SVPWM模块可以用Simulink自带的，也可以用S-function自己写。

PWM生成部分，逆变器开关频率通常选10kHz~20kHz，SVPWM载波频率和采样时间要匹配，离散仿真步长建议在载波周期的1/20以下，即5e-5~1e-4秒量级。

3.5 子系统四：主电路、LCL滤波与并网开关

主电路用Simscape Electrical的三相两电平桥（Three-Level Bridge / Two-Level Bridge）实现，开关管选IGBT/Diode模型。直流侧母线电容要有足够的容量来抑制母线电压纹波，一般按支撑持续时间来选：C = P * Δt / (Vdc * ΔVdc)，比如10kW系统允许母线电压纹波±10V，支撑时间5ms，C = 100000.005/(80020) ≈ 3125μF，取3300μF。

LCL滤波器是并网逆变器与电网之间的接口，它的参数直接影响谐波抑制效果和控制稳定性。LCL滤波的谐振频率要设计在开关频率的1/10到1/5之间，同时避开电网基波频率。典型设计：逆变器侧电感L1取0.1~0.15pu，电网侧电感L2取0.05~0.1pu，滤波电容Cf取额定功率条件下的5%以下容性无功功率对应的电容值。谐振峰需要加阻尼，最常用的是在滤波电容支路串联无源阻尼电阻，或者用控制算法做有源阻尼。

并网开关用Three-Phase Breaker模块控制，断开状态下模拟离网/孤岛运行，闭合后模拟并网运行。实际仿真中并网瞬间的冲击电流是重点考察指标，所以并网开关的闭合时刻需要和预同步完成信号联动，不是随便给个时间触发的。

3.6 直流母线电压稳定的控制逻辑

直流母线电压是整个系统的"枢纽变量"。光伏波动时储能必须立即响应，但储能响应前的一小段时间里母线电压会出现波动，这个波动通过VSG输出功率波动传导到交流侧。为了保证并网电能质量，我在模型里做了一个前馈补偿：把逆变器输出功率（交流侧）折算成直流电流前馈量，叠加到储能电流内环的参考值上。

前馈的逻辑不复杂：直流母线电容上的瞬时功率差等于光伏输出功率减去储能吸收功率再减去逆变器送入电网的功率。如果逆变器功率突然变化，前馈环节提前改变储能电流指令，而不是等母线电压偏差出来后再由PI调节器慢慢纠正。这样做母线电压的动态跌落可以减小一半以上。对于并网电能质量的改善非常明显，建议在做光储模型时一定把这个环节加上。

4. 关键参数整定：J、D与下垂系数不是拍脑袋定的

4.1 虚拟惯量和阻尼系数的物理意义

VSG的参数整定是整个仿真模型最重要的环节，也是最容易"调了半天还是振荡"的部分。先说虚拟惯量J。J的单位是kg·m²，但在仿真模型里它的数值不是来自物理转动部件，而是人为设定的"等效转动惯量"。J的作用是控制频率变化率：J越大，同样的有功不平衡引起的dω/dt越小，频率变化越平缓，但动态响应也越慢——频率要经过更长时间才能回到新的稳态。

阻尼系数D的作用是衰减振荡和提供稳态功率偏差。D越大，频率偏差引起的功率修正越强，系统越"硬"，振荡衰减越快。但D过大会让VSG在频率偏差时输出功率过大，影响调频机构的稳定性。

J和D必须一起匹配，不是单独调的。从二阶系统看，转子运动方程对应的特征方程是 Js² + Ds + K = 0（K是与调频系数和功率耦合相关的增益），阻尼比ζ = D / (2√(JK))。如果J设得很大但D不动，阻尼比会变小，系统反而更容易振荡；反过来J很小D很大，系统虽然稳定但惯量支撑能力就没有意义了。

4.2 有功下垂系数和无功下垂系数的选取

下垂系数的物理含义是"每偏差1Hz或每偏差100%功率对应的调节量"。有功下垂系数Kf（单位W/Hz）通常按下垂百分比来取：Kf = Pn / (下垂百分比 × fn)。以10kW系统为例，下垂取2%，50Hz额定频率，Kf = 10000 / (0.02×50) = 10000 W/Hz。意思是电网频率每下降1Hz，VSG输出有功增加10kW。这个值在工程上是合理的，2%~5%是电网调频的常规整定范围。

无功下垂系数n（单位V/Var）类似：n = ΔU / ΔQ。一般取额定电压的3%~5%作为无功从零到满发对应的电压偏差。在10kW/220V系统中，如果电压额定相电压有效值220V，无功下垂3%，即电压偏差6.6V时无功功率从0变到10kVar，n = 6.6/10000 ≈ 6.6e-4 V/Var。

4.3 一个10kW系统的完整参数算例

拿我验证过的一个10kW光储VSG系统举例，各环节参数如下：

系统额定功率10kW，交流线电压有效值380V，频率50Hz，直流母线电压800V，开关频率10kHz。

虚拟同步机参数：

• J = 2.5 kg·m²。计算依据：假设极端情况下有功阶跃5kW，要求频率变化率不超过1Hz/s（即6.28 rad/s²），由J = ΔP / (ωn × dω/dt) = 5000 / (314×6.28) ≈ 2.53 kg·m²。实际取2.5留一点裕度。
• D = 30。按阻尼比ζ≈0.7来反推：先根据系统增益估算K，然后由D = 2ζ√(JK)得出。如果仿真中发现振荡衰减偏慢，可以把D逐步加到40~50。
• Kf = 10000 W/Hz（下垂2%）。
• n = 6.6e-4 V/Var（下垂3%）。

功率滤波时间常数：5ms。

LCL滤波器参数：

• L1 = 2.5mH（逆变器侧）
• L2 = 1mH（电网侧）
• Cf = 10μF
• 阻尼电阻Rd = 2Ω（无源阻尼）

MPPT扰动步长：0.5%的参考电压，速度适中，不会在稳态时引起过大振荡。

储能电池参数：100Ah，额定电压540V（通过双向DC-DC升压到800V母线），SOC初始65%。

4.4 不同参数下的动态响应对比

下面这张表是我在参数扫描时记录的实际现象，可以直接作为你调参时的经验参考：

具体到仿真波形：J=1kg·m²时，5kW有功阶跃下频率能冲到51.2Hz，且有明显的一段时间振荡；J=3kg·m²时，峰值降到50.8Hz，但恢复时间反而变长。这说明J不是越大越好，它改变的是惯量频率动态的"时间尺度"，要和D配合起来才能同时获得小的频率偏移和快的恢复速度。

我在实际调试时习惯先固定下垂系数，然后调J和D：先给一个有功阶跃，只调D让系统振荡衰减到临界阻尼，再调J优化频率变化率指标。反复迭代三四轮就能找出一组不错的参数。如果你不想一步步试，可以用Matlab的Parameter Estimation工具箱把J和D作为待辨识参数，以实测/期望的功率阶跃波形为目标做自动整定，但前提是模型的初始波形已经基本稳定，否则辨识结果没有意义。

5. 仿真调试中的五个典型坑与排查路径

5.1 代数环：瞬时功率反馈与控制环路形成死锁

VSG模型最常见的仿真崩溃原因是代数环。代数环的产生路径是：VSG输出决定电压参考 → 逆变器输出电压电流变化 → 功率计算变化 → 功率反馈回VSG改变电压参考。这个过程在现实中是有物理延时的（PWM载波周期、线路传播、滤波电容充电都需要时间），但在Simulink连续仿真里，如果我们把功率计算模块的输出直接接入VSG功率外环的输入，会在同一仿真步内形成一个无延时的闭环，求解器必须迭代求解这个环，出现代数环，轻则仿真速度极慢，重则直接报错"Algebraic loop detected"。

解决办法有几个：一是用瞬时功率法的输出加一阶低通滤波器，滤波器本身提供了延迟，破坏了代数环。二是在功率计算输出端直接加一个Memory模块或Unit Delay模块，打一拍延迟。三是在离散仿真模式下，Simulink天然自带一拍延迟，代数环的出现概率会大幅降低。我强烈建议功率信号务必加滤波，这不仅解决代数环，还让功率反馈更平滑，否则反馈信号里的PWM谐波会导致VSG输出电压参考畸变，并网电流THD会明显超标。

5.2 并网瞬间的相位同步与预同步控制

如果你直接在仿真开始就让并网开关闭合，大概率会看到逆变器输出电流瞬间飙到几十上百安培——这是并网冲击电流。原因很简单：VSG输出的虚拟电势相位θ和电网相位的差值在并网瞬间不是零，相当于把一个"不匹配的电压源"硬接到电网上，压差全部砸在LCL滤波器上。

解决办法是做预同步控制。并网开关闭合前，VSG先空载运行，控制θ跟踪电网相位。预同步的实现方法：用PLL锁相得到电网相位θg，和VSG的θ做差，如果偏差不为零，通过一个PI调节器给虚拟角速度ω叠加一个修正量，让θ逐渐逼近θg，直到两者误差进入死区（比如0.5°以内），再闭合并网开关。

在Simulink里实现预同步控制时，注意需要在并网信号里加一个逻辑判断：只有预同步完成标志为true时，Breaker闭合信号才有效。否则会出现开关动作时刻相位还没对准的尴尬情况。我自己就遇到过"预同步逻辑写了但没加使能条件"的问题，仿真结果看起来并网开关在0.3s闭合，但相位对齐过程还没结束，冲击电流照样存在。

5.3 S-function在连续/离散仿真模式下的配置差异

VSG控制如果用了S-function，需要注意powergui的仿真模式选择。Simulink默认为连续仿真，S-function里的连续状态（mdlDerivatives）会用变步长积分器求解，这在电力电子电路里会导致一个严重的效率问题——开关管的非线性动态迫使求解器把步长缩到极小，仿真跑几秒钟往往要几个小时。

更合理的做法是：powergui选择Discrete模式，仿真步长设为5e-5秒（20kHz采样，与开关频率匹配）。此时S-function不能再写连续微分方程，要把转子运动方程离散化：

Δω[k+1] = Δω[k] + (Δt/J) * (Pm[k] - Pe[k] - D*Δω[k])
θ[k+1] = θ[k] + Δt * (ωn + Δω[k+1])

在S-function里就是把mdlDerivatives换成mdlUpdate，状态用离散状态定义。这个改动虽然代码差异不大，但很多人会忽略。我见过不少同学在连续模式下用S-function搭VSG模型，仿真卡到怀疑人生，最后全部重写成离散模式，问题立刻解决。

另一个S-function相关的坑是采样时间的设置。如果S-function的采样时间设为连续（-1），在离散模式下会出现采样率不匹配的警告；如果设为固定步长，要和powergui的步长保持一致。建议在S-function的mdlInitializeSizes里明确设置继承式采样时间，或者干脆用固定采样时间等于仿真步长。

5.4 powergui仿真步长与求解器选择

电力电子仿真里powergui的设置直接决定仿真能否顺利跑通。我常用的配置是：Simulation type选择Discrete，Solver选择ode45（变步长）或固定步长ode4，步长5e-5~1e-4秒。注意Simulink对电力电子的离散仿真要求整个模型统一步长，所以MCU里的PWM频率、S-function采样时间、滤波器数字控制采样时间全部要和powergui的步长对齐，不能出现"PWM用1e-5步长但powergui步长设成1e-4"这种错配，否则波形失真、谐波完全不能反映真实情况。

如果你真的需要跑连续仿真，建议先把S-function写成连续形式，然后powergui选择Continuous，Solver选ode23tb或ode15s（刚性问题求解器）。但实际经验告诉我：对于带IGBT开关的并网逆变器模型，连续模式基本没法用，一次并网暂态仿真能跑一个小时；离散模式下5秒仿真时间可能只需要几分钟。所以在做光伏储能VSG并网模型时，直接选离散模式是明智的。

5.5 储能SOC与直流母线电压的耦合问题

储能电池的SOC初始值如果设置不合理，会在仿真中途触发一个很隐蔽的问题：SOC接近下边界（比如20%以下）时电池内阻显著增大，双向DC-DC输出能力下降，直流母线电压开始缓慢跌落，VSG的并网有功功率跟着下降。从波形上看很像VSG控制参数有问题——输出功率慢慢偏离指令值——但根因其实在直流侧。

排查方法很简单：在仿真模型里把SOC波形、直流母线电压波形、VSG输出有功功率波形同时拉出来对比。如果SOC在下降、直流母线电压同步下降、有功功率也跟着跌，问题基本锁定在储能侧。解决办法要么提高SOC初始值，要么修改电池容量/放电倍率参数，要么优化DC-DC的电压外环PI增益，让母线电压在SOC变化时仍能稳定在允许范围内。

5.6 一次完整排查示例：VSG输出功率低频振荡

最后分享一个真实调试案例，完整呈现我从发现问题到定位根因的全过程。

现象：10kW光储VSG模型在0.5s时给Pref加了3kW阶跃，并网有功功率在阶跃后产生约2Hz的低频振荡，持续1秒多才衰减，峰值还超过了设定值30%。

排查第一步，先排除S-function或控制算法的问题。我单独给VSG控制子系统加了阶跃输入，断开与主电路的交互，用纯控制仿真验证转子运动方程的输出——如果控制回路自身就振荡，那问题在主电路；如果控制回路稳定，问题在功率耦合或直流侧。实测控制回路本身输出平滑，排除了算法bug。

排查第二步，检查直流母线电压波形。发现有功功率振荡时，直流母线电压也出现同频振荡。再检查储能电池电流，同样存在振荡分量。这说明问题不在交流侧VSG参数，也不在算法，而是直流母线的功率平衡被打破了——储能电流调节跟不上逆变器输出功率的变化，母线电容在充放电过程中产生电压波动，这个波动又通过VSG的有功-母线电压耦合回路反馈回交流侧，形成闭环振荡。

排查第三步，定位到储能DC-DC电压外环PI参数偏慢。把电压外环的带宽从5Hz提高到15Hz，储能电流的响应速度加快后，直流母线电压振荡幅度明显下降；再结合前馈补偿进一步改善，有功功率振荡问题完全消除。整个过程下来，问题根源是储能电流环的响应速度无法匹配VSG惯量响应的时间尺度。

这个案例给的教训是：VSG不是孤立的控制算法，它和直流侧储能、母线电容构成一个完整动态系统。调参时不能只盯着J和D，也要检查直流侧控制回路的带宽是否能覆盖VSG的动态范围。如果你在仿真里看到莫名其妙的低频振荡，先去看看储能电流环和母线电压波形，往往比你调半天VSG参数有效得多。

这套模型我从最开始的纯VSG控制，陆陆续续加进光伏MPPT、储能双向DC-DC、预同步和直流母线前馈补偿，前前后后用了大半年才稳定。现在它的可扩展性很好：想加故障穿越、低电压穿越逻辑，可以在VSG控制里加附加功角稳定器；想验证微网孤岛运行，把电网换成可编程三相电源做负荷突变就行；想研究多机并联，复制一套VSG子系统通过线路阻抗并联起来，正好可以研究无功功率分配问题。我用它复现过好几篇论文里的实验结果，波形趋势完全对得上，省去了大量重复造轮子的时间。如果你也正在搭类似的模型，希望这篇内容能帮你少踩几个我踩过的坑。