3. 阻抗建模的数学基础

新能源并网系统的阻抗建模，本质上是在稳态工作点附近建立“电压小扰动—电流小扰动”之间的线性映射。要理解后续的阻抗扫频、奈奎斯特判据和振荡分析，必须先掌握四块数学基础：小信号线性化、坐标系变换、谐波线性化，以及阻抗/导纳矩阵的定义。本节按“从非线性到线性、从三相到两相、从时域到频域、从标量到矩阵”的顺序展开。

一、小信号线性化：在稳态工作点上“切一刀”

并网逆变器的控制环路、PLL、坐标变换都含乘法、三角函数和饱和环节，整体是非线性系统。设系统状态方程为 ẋ = f(x, u)，在稳态工作点 (x₀, u₀) 处做泰勒展开并忽略高阶项，得到小信号模型：

Δẋ = A·Δx + B·Δu
Δy = C·Δx + D·Δu
其中 A = ∂f/∂x |₀ ，B = ∂f/∂u |₀

这一步的关键是：线性化只在工作点附近有效。工作点由并网电压幅值、有功/无功指令、直流母线电压共同决定。例如某 100 kW 逆变器在 V_dc = 800 V、V_g = 311 V(相电压峰值)、P = 100 kW、Q = 0 时，其 d 轴电流指令约为 214 A（按三相功率换算），q 轴约为 0。若把工作点改到 Q = 30 kvar，线性化后的 A、B 矩阵会变化，阻抗特性也随之改变。因此工程上做阻抗扫频时，必须记录并固定工作点，否则不同工况下的阻抗曲线不可直接比较。

二、dq 与 αβ 坐标系变换：把交流量变成直流量

三相静止坐标系 abc 下的电压电流是 50 Hz 正弦量，直接做小信号分析很麻烦。Clarke 变换把 abc 投影到两相静止坐标系 αβ：

[xα]  2 [1   -1/2   -1/2 ] [xa]
[xβ] = - [0  √3/2  -√3/2 ] [xb]
       3                     [xc]

Park 变换再把 αβ 旋转到与电网电压同步的 dq 坐标系，旋转角 θ 由 PLL 输出：

[xd]   [ cosθ  sinθ] [xα]
[xq] = [-sinθ  cosθ] [xβ]

变换后，基波正序分量变成 dq 下的直流量，控制器用 PI 就能实现无静差跟踪。但要注意：dq 变换是时变的，θ 本身受 PLL 动态影响。当电网电压出现小扰动时，θ 也会抖动，这个抖动会耦合进 dq 电流，形成 dq 阻抗矩阵中的交叉项（如 Z_dq、Z_qd）。这正是 dq 域阻抗建模比 αβ 域更复杂、但也更贴近控制器实际结构的原因。

αβ 域则保持了两相正交的静止特性，适合分析不平衡和諧波问题；dq 域适合分析对称基波附近的稳定性。两者通过旋转矩阵相互转换，建模时需明确参照系，全文保持一致。

三、谐波线性化：把扰动放到频域里看

小信号线性化在时域完成，谐波线性化则把扰动写成不同频率的正弦分量。设并网电压在基波上叠加一个小幅值、频率为 f_p 的正弦扰动：

v(t) = V₁·cos(2πf₁t) + V_p·cos(2πf_p t + φ_p)

由于系统非线性，响应电流中不仅含 f_p，还会出现 f₁±f_p 的边带分量。谐波线性化就是分别提取这些频率分量的幅值和相位，定义各频率点的阻抗：

Z(f_p) = V_p / I_p ， 同时记录 I(f₁±f_p)

工程上常用“电压扰动注入法”：在并网点串联一个宽频扰动源，逐点扫描 10 Hz～2 kHz，测量每个频点的电压电流相量，得到 Z(f) 曲线。若某频点阻抗幅值出现尖峰或相位穿越 -90°，往往对应系统潜在谐振点。例如在 f_p = 650 Hz 附近测得 |Z| = 2.1 Ω、∠Z = -88°，而逆变器输出阻抗在该频点呈感性，两者可能形成并联谐振，需要重点关注。

四、阻抗与导纳矩阵：从标量到多输入多输出

三相系统在 dq 域是两输入两输出，阻抗不再是标量，而是 2×2 矩阵。以电压为输入、电流为输出，导纳矩阵定义为：

[Δid]   [Ydd  Ydq] [Δvd]
[Δiq] = [Yqd  Yqq] [Δvq]

其逆矩阵即为阻抗矩阵 Z = Y⁻¹。对角元素 Ydd、Yqq 表示同轴电压扰动引起的同轴电流响应；非对角元素 Ydq、Yqd 表示交叉耦合，主要来自 PLL、坐标变换和控制环路的交叉项。若忽略交叉项，矩阵退化为两个独立标量，会丢失 dq 耦合引起的振荡模式，因此实际建模中一般不能随意忽略。

稳定性判据通常基于阻抗比矩阵：把逆变器输出阻抗 Z_inv 与电网阻抗 Z_grid 组成回比矩阵 L = Z_grid · Z_inv⁻¹，通过广义奈奎斯特判据判断系统稳定性。矩阵形式下，需要同时检查两个特征值的轨迹，而不是只看单条曲线。

五、建模时的几个实操要点

第一，工作点必须记录。同一台逆变器在 P=50 kW 和 P=100 kW 下的 dq 阻抗曲线不同，比较前先确认 V_dc、V_g、P、Q 一致。第二，坐标系约定要统一。d 轴对齐电网电压矢量、q 轴滞后 90°，还是反过来，会影响交叉项符号；本文采用 d 轴对齐电压、q 轴滞后 90° 的约定。第三，扫频扰动幅值要足够小，一般取基波幅值的 1%～5%，否则会激发非线性边带，测得的阻抗不再代表小信号特性。第四，矩阵求逆要检查条件数，若 Ydd 与 Yqq 接近相等且 Ydq、Yqd 很小，矩阵接近奇异，数值求逆会放大误差，此时可改用解析表达式或对矩阵做正则化处理。