VSG（虚拟同步发电机）做并网稳定性分析时，最怕两件事：一是拿不到准确的行为模型，二是拿到了也不知道怎么判断稳不稳。前两年我搞一块VSG样机的时候，就栽在低频振荡上，当时折腾了半天才意识到，只要把正负序阻抗扫出来，再加上时域下的阻抗稳定性判断，整个问题就清楚一大半。这篇文章就把我在Simulink里做VSG序阻抗扫频、电压电流双闭环建模，以及用序阻抗判据做时域稳定性分析的全套过程搬出来，适合正在做逆变器并网仿真、尤其被弱电网振荡问题折磨的人。
1. 项目背景与分析思路
1.1 为什么选用序阻抗法做VSG稳定性分析
VSG通常用在并网逆变器里，通过模拟同步发电机的转子运动方程，给电力电子变换器提供惯性和阻尼。很多人觉得只要把转子运动方程写对、参数整定好，系统就稳定。实际接入弱电网后，电网阻抗和控制器带宽耦合在一起，经常在10~300Hz这个频段出现振荡。前几年做样机时遇到过大约23Hz的功率振荡，当时第一反应是调惯量J、调阻尼D，调来调去就是压不住，后来换成阻抗视角才定位到问题：逆变器输出阻抗在中频段的相位和电网感性阻抗形成了负相位裕度，本质是控制环节和LCL滤波器共同造成的谐振风险。
阻抗法的优势在于它是一种“模块化”的分析思路。电网、逆变器、控制器可以分开处理，不需要建立超大的状态空间方程，也不需要对内部所有参数做小信号线性化。对VSG这种功率外环加电压电流双闭环的控制来说，只要把正序、负序输出阻抗的频率特性扫出来或建模出来，稳定性结论就落在几条波特图和奈奎斯特曲线上。
这个项目最终目标拆成三块：
- 在Simulink中搭建VSG主电路和电压电流双闭环控制；
- 通过扫频扰动获得正序和负序输出阻抗；
- 用阻抗比判据判断不同电网阻抗下的稳定性，再用时域仿真验证。
1.2 双闭环控制与阻抗之间的耦合关系
VSG完整控制链路是：功率计算 → VSG功率外环 → 参考电压 → 电压外环 → 电流内环 → PWM。其中电压外环和电流内环构成双闭环，直接决定逆变器输出阻抗的中高频段形状，而VSG功率外环决定低频段形状。
在分析阻抗时，功率外环的虚拟惯量J和阻尼D都会映射到输出阻抗的实部和虚部上，尤其在有功-频率、无功-电压两条通道上，低频阻抗往往表现出“负电阻”或“非线性”特征。双闭环则主导中高频段，电压外环的PI参数、电流内环的带宽、控制延时、滤波器参数，都会在波特图上留下可观察的特征。因此做“电压电流双闭环”的序阻抗扫频，本质上是在测整套控制策略在频域上的“对外表现”。
2. VSG控制原理与双闭环参数设计
2.1 VSG功率外环的原理
VSG的功率外环有两个核心方程。有功频率部分，虚拟转子运动方程一般写成：
J * dω/dt = Pm - Pe - D * (ω - ωref)
其中J是虚拟惯量，D是虚拟阻尼，Pe是实测有功功率，Pm是机械功率参考。有功调频部分可以再加一个有功下垂系数Dp：Pm = Pref + Dp * (ωref - ω)。无功电压部分，参考电压幅值可写成：
E = Eref + n * (Qref - Q)
n是无功下垂系数，Q是实测无功功率。两条通道给出的结果一个是相角θ，一个是电压幅值E，合成三相参考电压e_abc。
这里有一个容易被忽略的点：VSG功率外环本身是一个非线性、有一定带宽的环节，它的带宽通常只有几赫兹到十几赫兹。对比之下，电流内环带宽是1~2kHz，电压外环带宽是200~500Hz。所以不同频段的阻抗特性由不同控制环节主导，如果只拿传统的下垂控制模型去等效VSG，低频段阻抗必然对不上。
2.2 电压电流双闭环设计与参数范围
电压电流双闭环的共同作用是把VSG给出的三相参考电压跟踪成实际输出电压。电压外环的误差量是参考电压与实际输出电压的比较，PI调节器输出电流内环的参考值；电流内环再把电流参考和实际电感电流做比较，输出调制波。
电流内环一般按带宽1~2kHz整定。以L=3mH、R=0.1Ω的滤波电感为例，电流内环开环传递函数近似为G_i(s) = (Kp_is + Ki_i) / (sL + R)，再叠加上数字控制约1.5拍延时。按简化零极点对消法整定：
ω_bw ≈ 2π * 1500 rad/s Kp_i ≈ L * ω_bw ≈ 28 Ki_i ≈ R * ω_bw ≈ 940
实际工程中因为离散化、PWM调制增益、延时等因素，Kp_i一般不会取这么大，我项目中最后调试下来用的Kp_i=15、Ki_i=300，系统依然有足够带宽且阻尼更充分。电压外环带宽取电流环的1/5~1/10，约200~400Hz，我常用Kp_v=0.6、Ki_v=40。
参数整定确实没有唯一答案，关键是要理解“内环快、外环慢”的原则。内环不够快，外环可能引入振荡；外环不够慢，又会在某段频率上形成不利的相位滞后。这也是双闭环VSG阻抗扫频结果里中频段经常出现谐振峰的原因之一。
2.3 为什么正序和负序阻抗特性会不同
正序和负序阻抗不一样，直接原因就是控制的非对称性。VSG功率外环只建立在正序旋转坐标系中，负序分量基本不参与功率外环的调节，所以负序扰动进入后，响应路径主要是电压外环、电流内环以及LC滤波器的自然特性。正序扰动则会被功率外环先调一遍，多了一个低频段的“搬运”环节。
这意味着做稳定性分析时不能只测一个正序阻抗就完事。尤其是电网不对称、单相接地、负载不平衡等场景下，负序阻抗直接决定负序电流分配，影响保护和控制效果。这也是标题里强调“序阻抗扫频”而不是单纯“阻抗扫频”的原因。
3. 序阻抗扫频的实现方法
3.1 扫频的基本思想
阻抗扫频的本质是叠加信号和系统辨识。在50Hz基波工作点附近，加入一个小幅值正弦扰动，扰动频率从一个较低值逐步扫到上千赫兹。在每个频点，测并网点电压和并网电流的同频分量，做复数除法，就得到该频点下逆变器从端口看进去的等效阻抗。
输出阻抗定义成Zout = ΔV_pcc / ΔI_out，其中ΔI_out是逆变器流向电网方向的电流扰动。这里的方向符号特别重要，我前面提过，定义反了会导致后面所有相位结果偏差180°，奈奎斯特判据的结论完全相反。
扫频频率范围通常覆盖1Hz到1kHz，或者根据研究对象细化到10~500Hz。对于VSG，重点关注5~200Hz的次同步/中频段，所以扫频点在这个区间要加密，高频段可以稀疏一些。
3.2 扰动注入方式对比
扫频扰动有两种常见注入方式，我在项目里都试过。
串联电压源扰动：在PCC点和理想电网之间串入一个受控电压源，受控源的控制信号中包含扫频分量。这个方式的好处是物理直观，可以直接测逆变器输出阻抗；缺点是受控源串在功率回路里，如果处理不好，两个理想电压源串联可能造成代数奇异，需要在受控源两端并一个极小电阻或加缓冲阻抗。
并联电流源扰动：在PCC点并联一个受控电流源，注入扫频电流，测量PCC电压响应。这个方式不用断开主回路，但测得的结果是逆变器阻抗和电网阻抗的并联特征，除非把电网支路断开或做额外处理，否则不能直接得到Zinv。
综合来看，我推荐串联电压源方式。在Simulink建模时，受控源输出等于扫频扰动电压，基波部分由理想电网提供，这样受控源在无扰动时等同一个零电压源，不改变系统工作点。为了防止求解器报奇异的麻烦，我会在受控源模型里串入一个0.01Ω的阻尼电阻，对阻抗测量结果的影响可以忽略。
3.3 正负序扰动注入与对称分量处理
正序扰动和负序扰动差别只在于三相相位顺序。以注入80Hz扰动为例：
正序扰动：A相相位0°，B相-120°，C相+120°。 负序扰动：A相相位0°，B相+120°，C相-120°。
在Simulink里，可以用两个并联的受控电压源组分别产生正序和负序扰动，然后通过手动选择接入哪个扰动源。我习惯做成一个外部开关信号控制选择。
数据处理时使用对称分量变换。设某个频点处三相电压测得的复幅值为Va、Vb、Vc，则正序电压为Vpos = (Va + aVb + a^2Vc) / 3，负序电压为Vneg = (Va + a^2Vb + aVc) / 3，其中a = e^(j*2π/3)。电流也做同样的变换，正序阻抗就是Vpos / Ipos，负序阻抗是Vneg / Ineg。
这里有一个细节：FFT分析时通常推荐用矩形窗并保证截取的是扰动周期的整数倍，如果仿真时间波动、截取不整齐，就换Hann窗。由于阻抗是电压除以电流，同一个窗函数对两者产生相同加权，比值不受影响，所以用Hann窗也能获得相对准确的幅相结果。
4. 时域下阻抗扫频稳定性分析方法
4.1 阻抗比判据的推导思路
单机并网系统可以简化为逆变器等效阻抗Zinv与电网阻抗Zg串联的结构。如果Zg和Zinv形成一个闭合回路，系统稳定性就取决于这个回路的开环增益特性。常用的开环传递函数是L(s) = Zg(s) / Zinv(s)，然后用奈奎斯特判据判断L(jω)是否环绕(-1, j0)点。
对于三相平衡系统，正序和负序解耦，可以分别对正序回路和负序回路做独立的单输入单输出奈奎斯特判断。这正是先做序阻抗扫频，再做稳定性分析的原因。如果直接在dq坐标系下分析，就要面对一个2×2的阻抗矩阵，需要做广义奈奎斯特判据，复杂度高不少。
4.2 波特图交点与相位裕度
在不画奈奎斯特曲线的情况下，快速判断法是在幅频特性中找|Zg| = |Zinv|的交点，然后看交点处的相位关系。电网阻抗Zg通常是感性，相位大致在90°附近，如果逆变器阻抗在这一频率也呈感性，相位差不大，系统稳定裕度足够；如果逆变器阻抗呈容性，相位差接近180°，就可能在交点频率附近形成负阻尼振荡。
相位裕度的定义：
PM = 180° + ∠Zg(jωc) - ∠Zinv(jωc)
其中ωc是幅值交点频率。PM大于0且有一定裕量，系统一般认为稳定；PM接近0或为负，就需要警惕。
实际操作中发现，扫频得到的相位曲线在控制带宽附近往往有一个快速下降的过程，这也是VSG中频振荡的高发区。比如电流内环带宽约1.5kHz，那么在几百赫兹频率处，相位就可能从接近0°快速掉到-90°甚至更低，如果在这个频段恰好和电网阻抗的幅值相交，相位裕度很容易吃紧。
4.3 时域验证设计
扫频得到的是频域数据，属于“预测”。要说服自己，还得回到时域做验证。我的验证套路很固定：
先按阻抗判据找到临界状态，比如电网电感Lg增大到某个值时相位裕度接近0。然后在Simulink模型里把电网电感设成这个临界值，在0.2s时投入一个负荷阶跃或并网断路器动作，观察功率、电流、电压波形。
如果理论判据说稳定，时域波形应该快速衰减；如果判据说失稳，时域波形应该出现增幅振荡或等幅振荡。振荡频率应该和阻抗幅值交点频率接近，这一步可以对上号。
5. 基于Simulink的建模仿真全流程
5.1 模型整体架构
Simulink模型我习惯按五块搭：
- 主电路部分：直流电压源、三相桥、LC滤波器、串联扰动源、理想电网。
- VSG功率外环：有功功率计算、无功功率计算、转子运动方程、参考电压生成。
- 电压外环和电流内环：各轴PI调节、坐标变换、电流内环受控源输出。
- 测量与信号记录：电压电流测量、To Workspace记录。
- 扫频控制接口：接收基座工作空间变量fp_scan、Vscan等。
为了扫频效率和数值稳定性，我强烈建议先搭平均模型，也就是用受控电压源代替IGBT桥臂，调制信号乘以Vdc/2就是输出电压。平均模型没有开关纹波，定步长可以放到几十微秒，扫30个频点很快就能跑完。需要看高频谐波或验证EMC时再用开关模型，但阻抗扫频的核心频段（10~1000Hz）平均模型和开关模型的幅相特性基本一致。
5.2 主要参数参考表
我用一组典型的10kW/380V VSG参数作为示例，方便直接复现：
| 参数名称 | 符号 | 数值 | 说明 |
|---|---|---|---|
| 直流母线电压 | Vdc | 750V | 三相桥直流侧 |
| 滤波电感 | Lf | 3mH | 含内阻Rf=0.1Ω |
| 滤波电容 | Cf | 20μF | 输出滤波 |
| 电网相电压有效值 | Vg | 220V | 单相 |
| 电网频率 | fg | 50Hz | |
| 开关频率 | fsw | 10kHz | 平均模型可放宽 |
| 电流内环PI | Kp_i / Ki_i | 15 / 300 | |
| 电压外环PI | Kp_v / Ki_v | 0.6 / 40 | |
| VSG惯量 | J | 0.8 | kg·m²量级 |
| VSG阻尼 | D | 25 | |
| 扫频扰动幅值 | Vscan | 6V | 约为额定相电压2% |
| 扫频范围 | fp | 10~1000Hz | 对数分布 |
电网电感Lg是稳定性分析中要变化的量，建议用1mH、3mH、5mH、8mH几个档位分别扫。
5.3 扰动源与数据记录
串联扰动源的控制信号生成很简单。以正序扫频为例，用一个正弦信号发生器，频率设为fp_scan（从工作空间读取），A、B、C三相的相位偏移分别为0、-2π/3、+2π/3。控制信号乘以Vscan后送给受控电压源。
在Simulink中，受控电压源模块（Controlled Voltage Source）的输入就是扰动电压瞬时值。由于受控源串联在PCC和电网之间，基波情况下不输出电压，所以它不影响稳态工作点。注意我在该受控源支路串入0.01Ω电阻，可以避免理想电压源串联的数值问题。
数据记录用To Workspace模块，输出变量名设为v_pcc和i_out，采样时间设为固定步长，比如1e-6秒或2e-5秒。所有扫频点的数据都保存后，脚本统一做FFT处理。
5.4 自动化扫频脚本
扫频全过程用MATLAB脚本批量跑。模型里的扰动频率通过assignin注入到工作空间，Simulink里的正弦源直接引用这个变量。脚本核心代码如下：
%% VSG序阻抗扫频主脚本 clear; clc; close all; % 扫频点：10Hz 到 1000Hz，对数分布取30个点 freqs = logspace(log10(10), log10(1000), 30); Vscan = 6; % 扰动幅值，约为相电压幅值311V的2% simStop = 0.4; % 每个频点仿真时长 % 预分配结果存储 Z_p = zeros(1, length(freqs)); Z_n = zeros(1, length(freqs)); for k = 1:length(freqs) fp = freqs(k); % 把扰动频率和幅值注入到模型工作空间 assignin('base', 'fp_scan', fp); assignin('base', 'Vscan_', Vscan); % 运行仿真 simOut = sim('VSG_ImpedanceScan', ... 'StopTime', num2str(simStop), ... 'Solver', 'ode3', ... 'SolverType', 'Fixed-step', ... 'FixedStep', '1e-6'); % 取数据 t = simOut.tout; v_pcc = simOut.v_pcc; % N x 3 i_out = simOut.i_out; % N x 3 % 跳过暂态，从0.2s之后取数据 idx = t > 0.2; t = t(idx); v_pcc = v_pcc(idx,:); i_out = i_out(idx,:); % 按扰动周期截取整数倍 Tp = 1 / fp; tend = t(1) + floor((t(end)-t(1)) / Tp) * Tp; idx2 = t <= tend; t = t(idx2); v_pcc = v_pcc(idx2,:); i_out = i_out(idx2,:); % 加窗做FFT N = length(t); win = hann(N); V = fft(v_pcc .* repmat(win,1,3), N, 1); I = fft(i_out .* repmat(win,1,3), N, 1); % 频率轴 dt = t(2) - t(1); freq_axis = (0:N-1) / (N*dt); [~, kf] = min(abs(freq_axis - fp)); % 取扰动频率点的复数幅值 Va = V(kf,1); Vb = V(kf,2); Vc = V(kf,3); Ia = I(kf,1); Ib = I(kf,2); Ic = I(kf,3); % 对称分量变换 a = exp(1i * 2 * pi / 3); Vpos = (Va + a*Vb + a^2*Vc) / 3; Vneg = (Va + a^2*Vb + a*Vc) / 3; Ipos = (Ia + a*Ib + a^2*Ic) / 3; Ineg = (Ia + a^2*Ib + a*Ic) / 3; % 记录阻抗 Z_p(k) = Vpos / Ipos; Z_n(k) = Vneg / Ineg; fprintf('f=%.2fHz Zp=%.3f∠%.1f°, Zn=%.3f∠%.1f°\n', ... fp, abs(Z_p(k)), angle(Z_p(k))*180/pi, ... abs(Z_n(k)), angle(Z_n(k))*180/pi); end % 绘制正序阻抗波特图 figure; subplot(2,1,1); semilogx(freqs, 20*log10(abs(Z_p)), 'b-o'); grid on; ylabel('幅值 (dBΩ)'); title('VSG正序输出阻抗'); subplot(2,1,2); semilogx(freqs, angle(Z_p)*180/pi, 'b-o'); grid on; ylabel('相角 (deg)'); xlabel('频率 (Hz)'); % 绘制负序阻抗波特图 figure; subplot(2,1,1); semilogx(freqs, 20*log10(abs(Z_n)), 'r-o'); grid on; ylabel('幅值 (dBΩ)'); title('VSG负序输出阻抗'); subplot(2,1,2); semilogx(freqs, angle(Z_n)*180/pi, 'r-o'); grid on; ylabel('相角 (deg)'); xlabel('频率 (Hz)');扫频时要注意，负序扰动需要把B、C两相相位对调。自动化脚本里可以在模型内部加一个“正负序选择”开关，通过assignin('base','scan_mode', mode)来切换。负序模式下的对称分量变换逻辑完全一样，不需要改脚本。
6. 仿真结果分析与案例展示
6.1 正序阻抗扫频结果
以我项目中的参数为例，正序输出阻抗波特图通常长这样：
低频段（1~30Hz），由于VSG功率外环的积分特性，阻抗幅值整体偏低，相位接近0°，甚至可能出现小幅的“负电阻”区域，反映的是虚拟惯量和阻尼的参数映射。中频段（50~300Hz），电压外环作用开始显现，阻抗曲线会出现一个明显的谐振峰，相位从感性向容性过渡。高频段（500Hz以上），电流内环带宽进入衰减区，滤波电感主导，阻抗幅值按照+20dB/十倍频程上升，相位回到90°附近。
这些特征和理论小信号模型基本吻合。如果扫出来曲线形状不一样，先检查参数整定是否合理，再检查电压外环和电流内环的带宽是否满足“内环快、外环慢”。
6.2 负序阻抗扫频结果
负序阻抗曲线和正序有一个明显区别：低频段没有VSG功率环带来的“负电阻”区域，而是直接受电压外环控制。一般来说，负序阻抗在低频段表现为容性或电感性的组合特征，相位更接近-90°到+90°之间平滑变化，高频段与正序阻抗逐渐重合。
这也说明，如果只在正序域做稳定性判断，会漏掉负序通路带来的隐患。例如电网发生不对称故障或单相负载较重时，负序阻抗决定了负序电压如何分配，尽管功率外环不感知负序，但电压外环和电流内环仍然会响应并影响负序阻抗，进而影响故障电流和并网点电压质量。
6.3 稳定性判定与时域验证
我取了三组电网电感来演示稳定性判断：
| 电网电感Lg | 幅值交点频率 | 相位裕度PM | 时域结果 |
|---|---|---|---|
| 1mH | 约210Hz | 约48° | 稳定 |
| 4mH | 约130Hz | 约6° | 临界稳定 |
| 8mH | 约80Hz | 约-15° | 失稳 |
判断方法就是画L(s)=Zg(s)/Zinv(s)的奈奎斯特曲线，或者干脆用|Zg|=|Zinv|的交点算相位裕度。在时域验证时，把电网电感改成8mH，模型在0.2s投入一个小负荷扰动，很快就看到功率波形出现约80Hz的增幅振荡，和阻抗交点频率对应。
这个案例是典型的VSG弱电网中频振荡场景。很多人在现场看到几十赫兹的功率振荡，第一反应是逆变器参数不对，但通过阻抗扫频能快速定位到哪个频段、哪个控制环节在贡献负阻尼。这就是序阻抗扫频的价值。
7. 常见问题与排查技巧
7.1 阻抗曲线毛刺多怎么办
扫频结果毛刺多，最常见原因是FFT截取的不是扰动周期的整数倍，频谱泄漏严重。解决方式是保持每个频点仿真时间足够长，截取数据时精确按Tp的整数倍截取。如果还是抖，就把Hann窗加上，虽然牺牲了一点频谱分辨率，但幅相提取更稳定。
另一个原因是扰动频率和电网基波或谐波太靠近。比如扫到49Hz、51Hz这类频率时，响应分量和50Hz基波混叠，很难分离。我会在扫频点里主动排除基频附近的频点，或者设置一个“死区”，例如10~40Hz、60~100Hz这样分段扫。
7.2 扰动幅值怎么选
扰动幅值太大会使系统进入非线性区，扫出来的“阻抗”其实是大信号等效阻抗，不是小信号阻抗；太小时信噪比低，FFT提取不准确。我的经验是取额定相电压幅值的2%~5%，比如311V的系统用6~15V扰动。
电网电感大的时候可以适当减小扰动幅值，因为弱电网对大信号扰动更敏感；电网强时可以用偏大一点的幅值提高信噪比。
7.3 与理论模型对不上
扫频结果和解析小信号模型对不上，先怀疑三个地方：
第一，电流方向定义。前面反复强调过，输出阻抗定义中电流方向是“逆变器流向电网”，如果模型里把电流测量模块的极性接反，相位会整体差180°。
第二，控制延时模型。数字控制器的延迟对高频段阻抗影响很大，1拍延时和1.5拍延时在高频段的相位差异可能达到几十度，和理论模型对比时要把控制延时计入。
第三，正负序分离逻辑。对称分量变换公式中a的幂次写错，或者负序扰动相位没反转，都会导致结果错乱。建议先对一个纯阻感负载模型做扫频，验证算法正确性，再上完整VSG控制。
7.4 仿真速度太慢
平均模型是首选。如果坚持用开关模型，在Simulink里用离散定步长，步长至少小于开关周期十分之一，10kHz开关频率下建议1e-7到1e-6秒，但每个频点仿真0.4s会比较慢，30个频点可能要跑很久。
优化办法是减少每个频点仿真时长。低频频点需要更多时间建立稳态，可以跑0.6s；高频频点暂态建立快，跑0.15s就够。另外可以用parfor并行跑扫频频点，只要每个频点之间没有状态依赖，并行效率很高。
8. 个人经验与总结
阻抗扫频这条路，我踩过最大的坑就是方向符号和对称分量提取的细节，看似小问题，却能让后续所有稳定性结论反转。后来我养成了一个习惯：每次搭好模型，先拿一个已知RL负载做扫频验证，扫出来的阻抗曲线和解析值对得上，再往VSG上换。
还有一点体会，扫频只是工具，判断稳定性最终要看奈奎斯特曲线或相位裕度，但同时一定要配一组时域验证。频域预测说稳定，时域却振荡了，那一定是有序阻抗法没覆盖的因素，比如负序通路、控制非线性、锁相环带宽等。VSG相对传统下垂控制多了一个功率外环，低频段阻抗特性更复杂，扫频范围一定要覆盖到几赫兹，否则会漏掉低频振荡风险。
另外，采样数据做FFT时，建议把每个频点的仿真时长和扰动幅值都自动记录到结果文件里，方便回溯。阻抗扫频和稳定性分析往往要迭代好几轮参数，有完整的记录能省去大量重复劳动。
这套方法不只适用于VSG，任何并网变流器、甚至在电机驱动并网场景里，序阻抗扫频加时域稳定性验证都是一套通用打法。如果你正在被弱电网振荡或模型不透明的问题困扰，建议按这个链路试一遍：搭模型、扫正负序阻抗、算相位裕度、时域验证，四步走下来，问题基本能定位到具体频段和控制环节。
