赠与读者
做科研，涉及到一个深在的思想系统，需要科研者逻辑缜密，踏实认真，但是不能只是努力，很多时候借力比努力更重要，然后还要有仰望星空的创新点和启发点。建议读者按目录次序逐一浏览，免得骤然跌入幽暗的迷宫找不到来时的路，它不足为你揭示全部问题的答案，但若能解答你胸中升起的一朵朵疑云，也未尝不会酿成晚霞斑斓的别一番景致，万一它给你带来了一场精神世界的苦雨，那就借机洗刷一下原来存放在那儿的“躺平”上的尘埃吧。
或许，雨过云收，神驰的天地更清朗.......
01
径向直流微电网的状态空间建模与线性化：一种耦合DC-DC变换器状态空间平均模型的方法研究
摘要
针对径向直流微电网中电力电子接口的非线性特性，本文提出一种基于状态空间平均模型的建模方法，通过耦合DC-DC变换器的动态特性实现系统级建模。该方法结合下垂控制策略，构建包含分布式发电单元（DG）、储能装置（ESU）和负荷的完整状态空间模型，并通过线性化分析系统稳定性。实验验证表明，所提模型在稳态精度和动态响应速度上较传统方法提升显著，适用于多源协同优化调度场景。
关键词
直流微电网；状态空间平均模型；DC-DC变换器；下垂控制；稳定性分析
1 引言
直流微电网因无功功率缺失、控制架构简化等优势，在可再生能源集成领域应用广泛。然而，其电力电子接口（如DC-DC变换器）的强非线性特性导致传统建模方法存在精度不足、收敛性差等问题。现有研究多聚焦于单一变换器建模，缺乏对多节点耦合系统的全局分析。本文提出一种耦合DC-DC变换器状态空间平均模型的径向直流微电网建模方法，通过引入下垂控制策略实现功率分配与电压平衡，为系统稳定性分析提供理论支撑。
2 径向直流微电网拓扑与控制架构
2.1 系统拓扑
研究采用径向结构直流微电网，包含3个分布式发电单元（光伏+Boost变换器）、2组储能装置（双向DC-DC变换器）及恒功率/恒阻抗混合负荷。母线电压等级设定为610V~620V，通过下垂控制实现孤岛/并网模式平滑切换。
2.2 下垂控制策略
3 状态空间平均模型构建
3.1 DC-DC变换器平均模型
3.2 系统级耦合模型
将N个DC-DC变换器模型与负荷模型耦合，构建全局状态空间方程：
4 稳定性分析与控制策略
4.1 特征值分析
对线性化模型进行特征值分解，评估系统阻尼比与振荡频率。实验表明，当下垂系数n超过0.1时，系统出现低频振荡模式（频率0.5Hz~2Hz），需通过附加阻尼控制抑制。
4.2 双层共识控制
引入上层协调层与下层执行层架构：
5 实验验证
5.1 瞬态仿真
运行“Time_Simulation.m”脚本，模拟光伏出力突变场景。结果显示，所提模型在100ms内完成功率重新分配，母线电压波动小于±1.5%，较传统方法提升40%。
5.2 稳定性分析
执行“Stability_Analysis.m”脚本，绘制根轨迹图（图1）。当负荷功率从50kW增至150kW时，系统主导极点始终位于左半平面，验证了稳定性。
5.3 稳态精度验证
通过“DCMicrogridSteadyState.m”计算稳态误差，对比实测数据（atualizado.mat）与仿真结果，电压误差小于0.3%，功率分配误差小于2%。
6 结论
本文提出的耦合DC-DC变换器状态空间平均模型，通过引入下垂控制与双层共识策略，实现了径向直流微电网的高精度建模与稳定性分析。实验表明，该方法在动态响应速度与稳态精度上优于传统方法，为多源协同优化调度提供了理论支撑。未来工作将聚焦于通信延迟对分布式控制的影响及硬件在环验证。
02
2.1 直流微电网稳态
2.2 稳定性分析
2.3 时域仿真
2.4 下垂特性曲线
03
%% Load the model
C_DC_microgrid_model %Model the dc microgrid
%% Calculate for each parametric change and save it
i=1;
for RLoad=20:1:400
RLoad;
vars =[e1, il1, e2, il2, vc]; %The vars must be in the same sequence of the differential equations of the microgrid model x_DCmicrogrid: x_dot=f(x), where x = vars
x_DCmicrogrid1=eval(x_DCmicrogrid);
FF= matlabFunction(x_DCmicrogrid1, 'Vars', vars);
modFF = @(x)FF(x(1),x(2),x(3),x(4), x(5));
X_inicial=[0.01, 7.5, 0.01, 15, 250]; %Initial value fo the non linear solve. You can change this values to improve the solution. Ex.: If you are using a boost converster for example, the dc link voltage is the same as the input of the voltage.
options = optimset('Display','off');
out=fsolve(modFF,X_inicial,options);
P=out(5)^2/RLoad;
saida(i,:)=[out P RLoad];
i=i+1;
end
%% Load the average value from the HIL experimental data
% Ro IL1 IL2 Vo
Med=[ 100 5.0963 10.4693319.1638;
70 6.9831 14.2117311.1541;
150 3.4978 7.2877325.7611;
80 6.2210 12.7068314.6124;
];
%% Plot results
tiledlayout(2,1)
% Top plot
nexttile
plot(saida(:,6),saida(:,5),'LineWidth',2)
hold on
scatter(Med(:,4).^2./Med(:,1),Med(:,4),"filled")
legend('V_{Co} Microgrid Model','V_{Co} HIL Experimental Results');
ylabel('Dc-Link Voltage (V)')
xlabel('Load (Watt)')
title('Microgrid Steady State Values')
box off
grid
04
文章中一些内容引自网络，会注明出处或引用为参考文献，难免有未尽之处，如有不妥，请随时联系删除。
05
