基于FREEDM架构的微电网稳定性与参数灵敏度分析系统

资 源 简 介

本项目针对基于FREEDM（未来可再生电能传输与管理）概念的先进微电网架构，开发了一套完整的数学建模与稳定性分析平台。项目首先基于固态变压器（SST）、分布式可再生能源（DRER）及储能装置（DESD）构建了微电网的高阶非线性动态模型，并利用小信号平均法将其线性化为状态空间模型。核心功能在于对系统进行多维度的稳定性分析，包括通过计算系统雅可比矩阵的特征值来评估不同工况（如并网、孤岛、负载突变）下的系统暂态稳定性。在此基础上，项目深度集成了灵敏度分析算法，能够定量计算控制器参数（如电压/电流环PI增益）、电路硬件参数（如滤波电感、直流母线电容）以及线路阻抗对系统主导极点（特征值）实部和虚部的影响程度。通过计算参与因子和灵敏度系数，项目能够自动识别出对系统稳定性影响最大的关键参数，并绘制参数变化的根轨迹图，从而为微电网控制器的鲁棒设计和参数优化提供直接的理论依据和量化指标，确保系统在宽参数变化范围内的稳定运行。

详 情 说 明

基于FREEDM架构的微电网稳定性与参数灵敏度分析系统

项目简介

本项目针对基于FREEDM（未来可再生电能传输与管理）概念的先进微电网架构，开发了一套数学建模与稳定性分析平台。系统核心聚焦于固态变压器（SST）逆变级的动态特性，构建了基于dq旋转坐标系的高阶小信号状态空间模型。

该工具旨在解决微电网中电力电子变换器（SST）在弱电网或参数摄动下的稳定性评估问题。通过集成特征值分析、参与因子计算、参数灵敏度量化分析以及根轨迹扫描功能，项目能够定量评估控制器参数（如PI增益）和电路硬件参数（如滤波电感）对系统主导极点的影响，为微电网控制器的鲁棒设计和参数优化提供理论依据。

功能特性

本项目实现了以下核心功能：

• 高阶系统建模：基于电路拓扑和双闭环控制策略，自动构建包含LC滤波器动态与PI控制器状态的8阶小信号状态空间模型。
• 稳定性量化评估：计算系统雅可比矩阵的特征值，并输出系统的最小阻尼比和自然频率，判断系统在当前工况下的暂态稳定性。
• 模态参与分析：通过计算参与因子矩阵，识别主导极点（距离虚轴最近的不稳定/低阻尼极点）与系统状态变量（如电感电流、电容电压、控制器积分态）之间的关联度。
• 参数灵敏度分析：定量计算关键控制参数（$K_{pv}, K_{iv}, K_{pi}$）及硬件参数（$L_f$）对主导极点实部的灵敏度系数，直观展示参数变化对稳定性的影响方向与程度。
• 根轨迹扫描：支持对指定参数（如电压环比例增益 $K_{pv}$）进行大范围扫描，绘制系统特征值的根轨迹图，揭示参数变化过程中的稳定性边界。
• 时域仿真验证：基于线性化模型进行时域仿真（Lsim），模拟电压参考阶跃及负载电流突变工况下的系统动态响应。
• 多维可视化：自动生成包含特征值分布、参与因子柱状图、灵敏度系数对比、根轨迹图及电压时域响应的综合分析图表。

系统要求

• MATLAB R2018b 或更高版本
• Control System Toolbox（控制系统工具箱）

使用方法

1. 配置环境：确保MATLAB的工作路径包含项目脚本所在文件夹。
2. 运行程序：直接运行主脚本（main），程序将自动执行从参数定义到绘图的全流程。
3. 查看结果：控制台将输出系统最小阻尼比等关键指标，随后弹出图形窗口展示5个子图的分析结果。

代码实现逻辑与算法分析

系统主体逻辑主要包含在主执行脚本和核心建模函数中，具体实现细节如下：

1. 系统参数初始化

程序首先定义了微电网的物理参数，包括电网电压（380V）、频率（50Hz）、SST的额定功率以及LC滤波器的无源器件参数（$L_f, R_f, C_f, R_c$）。控制部分定义了基于dq坐标系的电压外环和电流内环的双闭环PI参数，以及数字控制延迟（虽然定义了延迟参数，但在当前状态空间模型中主要体现为无延迟的连续域近似或特定处理）。

2. S-SST 状态空间建模 (build_sst_model)

这是系统的核心算法函数。它构建了SST逆变级在dq同步旋转坐标系下的小信号模型 $dot{x} = Ax + Bu$。

• 状态变量：选取了8个状态变量 $x = [i_{Ld}, i_{Lq}, v_{Cd}, v_{Cq}, gamma_d, gamma_q, phi_d, phi_q]$。其中 $i_{L}, v_{C}$ 为LC滤波器的物理状态，$gamma$ 为电流环积分器状态，$phi$ 为电压环积分器状态。
• 系统矩阵 (A)：推导了包含电路微分方程（如 $-R/L$ 阻尼项、$omega_0$ 耦合项）和闭环控制方程（PI增益、$v_{ref}$ 前馈等）的雅可比矩阵。代码中详细实现了d轴与q轴之间的交叉耦合项（如 $omega L i_q$ 和 $omega C v_q$）。
• 输入矩阵 (B)：定义了4个系统输入 $u = [v_{ref_d}, v_{ref_q}, i_{load_d}, i_{load_q}]^T$，分别对应参考电压设定值和负载电流扰动。
• 输出矩阵 (C)：提取 $[v_{Cd}, v_{Cq}, i_{Ld}, i_{Lq}]$ 作为观测输出。

3. 特征值与参与因子分析

• 利用 eig(A) 计算全系统特征值。
• 计算阻尼比 $zeta = -real(lambda) / |lambda|$，用于评估振荡衰减特性。
• 通过左特征向量 $Psi$ 和右特征向量 $Phi$ 的Hadamard积计算参与因子矩阵 $P = |Phi circ Psi^T|$，并针对主导极点（实部最大的负根）进行归一化处理，从而找出导致系统不稳定的主要状态变量。

4. 关键参数灵敏度计算

程序调用灵敏度计算模块（在代码结构中体现为 calculate_sensitivity），针对主导极点，计算其特征值实部 $sigma$ 对参数 $p$ 的偏导数 $frac{partial sigma}{partial p}$。在主流程中，提取了电压环PI参数（$K_{pv}, K_{iv}$）、滤波电感（$L_f$）和电流环比例系数（$K_{pi}$）的灵敏度结果并进行对比展示。

5. 根轨迹参数扫描

为了观察电压环比例增益 $K_{pv}$ 的影响，代码设定了一个从 0.01 到 0.5 的线性区间。在循环中，程序反复修改临时参数结构体并重新调用 build_sst_model 生成模型，记录这一过程中所有特征值的移动轨迹。

6. 负载阶跃时域仿真

利用 lsim 函数对构建的线性系统 (sys_ss) 进行时域仿真。

• 工况设置：在 $t=0$ 时刻施加 d轴电压参考值（311V）。
• 扰动设置：在 $t=0.05s$ 时刻，模拟负载电流 d轴分量突增 10A（负载突变）。
• 仿真结果用于验证频域分析（特征值）与时域表现（超调、调节时间）的一致性。

7. 结果可视化

程序最终生成一个包含5个子图的综合窗口：

1. 特征值分布图：展示所有极点在复平面上的位置，标示稳定性边界。
2. 参与因子分析图：柱状图显示各状态变量对主导极点的贡献度。
3. 灵敏度系数图：对比 $K_{pv}, K_{iv}, L_f, K_{pi}$ 对稳定性的敏感程度。
4. 根轨迹图：展示随着 $K_{pv}$ 增大，系统极点的移动趋势。
5. 时域响应图：绘制 d轴和 q轴输出电压随时间变化的曲线，直观展示负载突变时的电压跌落与恢复过程。