做新能源仿真这块，风光储联合发电系统的Matlab模型我前前后后搭过好几个版本，最让我觉得值得拿出来分享的，是带超级电容储能的那一版。风力发电和光伏发电天生看天吃饭，出力曲线忽高忽低，直流母线电压跟着晃；想把这种晃动压住，又不让储能系统成为整个系统的短板，超级电容是个非常有意思的选择。这篇笔记把整个风光储Matlab仿真模型的思路、关键模块搭建过程、控制器设计和调试经验都捋一遍，适合正在做微电网、可再生能源发电系统仿真，或者对超级电容储能建模感兴趣的同学参考。既讲清楚为什么这么设计，也会把实际操作里的配置细节、参数计算过程放出来，照着能搭出来，少走我的弯路。
1. 风光储系统在做什么，超级电容为什么必须在场
1.1 风光出力的波动到底有多频繁
很多刚开始接触风光储模型的人，容易把风电和光伏当成一个"稳定的电源"去处理。实际仿真里把风速模型一加，光伏辐照度曲线一导入，你就会发现根本不是这么回事。风功率和风速是三次方关系，风速稍微动一动，出力就是大幅度的起伏；光伏更不用说，云层只要从组件上方飘过，辐照度可以在几秒到几十秒内快速跌落，输出功率直接掉一半。
这种波动在直流微电网里会直接反映为母线电压的波动和负载供电的不稳定。如果系统是并网的，情况还会更复杂——功率振荡会通过逆变器传导到电网侧，严重的时候会引起保护动作。我做模型时最直观的感受是：不给储能时，母线电压波形就像心电图一样乱跳；给上储能后，波形立刻被"熨平"，这就是储能在这个系统里的价值，不是简单的存电放点，而是在几十毫秒到几十秒的时间尺度上吸收和释放功率。
1.2 超级电容和电池，谁更适合干这活
储能设备选型这个问题，做仿真的人必须想明白，否则后面搭模型会陷入"为什么我要用超容"的自我怀疑。先说结论：在平抑高频波动这件事上，超级电容比电池合适得多。
电池的能量密度高，意味同样体积能存的电多，但它的功率密度有限，大电流充放电会发热、折寿。你用锂电池去吸收风光系统几秒钟一次的高频功率抖动，循环寿命会掉得很快。超级电容正好反过来，功率密度可以做到几千瓦每千克，充放电响应在毫秒级，循环寿命动辄几十万甚至上百万次。我在Simulink里做过一个对比：把同样的高频波动功率指令给电池模型和超容模型，电池的端电压和SOC曲线明显跟不上参考，而超容几乎能贴着指令走。
所以工程上常见的做法是"电池+超容"混合储能，电池管长时间尺度的能量转移，超容管短时间尺度的功率平抑。我这套模型先把焦点放在超级电容上，把功率型储能建模和控制做扎实，之后再扩展混合储能就顺理成章了。
| 特性 | 超级电容 | 锂电池 |
|---|---|---|
| 功率密度 | 高，可达数千W/kg | 较低，数百W/kg级别 |
| 能量密度 | 低，不适合长时间储能 | 高 |
| 响应速度 | 毫秒级 | 秒级 |
| 循环寿命 | 数十万到百万次 | 数千次 |
| 典型适用场景 | 高频功率波动平抑、瞬时支撑 | 削峰填谷、长时间备用 |
1.3 这个模型要回答什么问题
我在设计这套仿真模型时，给自己定了三个必须回答的问题。第一，风光出力波动在直流母线上怎么传递，储能介入前后波形有什么变化。第二，超级电容通过双向DC-DC变换器接入母线后，控制策略怎么设计，才能既平抑高频波动又稳住母线电压。第三，控制器里的参数——低通滤波时间常数、PI参数、容量配置——具体怎么计算和整定。
这三个问题解决掉，模型也就搭成了。后面所有模块，都是围绕这三个问题展开的。
2. 超级电容储能系统的模型架构，跟它背后的原理
2.1 整体拓扑和能量流向
这套模型采用的拓扑是经典的直流微网结构：风力发电单元（风机+永磁同步发电机+整流器）、光伏发电单元（光伏阵列+Boost变换器）、超级电容储能单元（超容模组+双向DC-DC变换器）全部挂在一条共同的直流母线上，母线再通过逆变器连接交流负载或电网。
从仿真角度看，每个发电单元都在向母线注入电流，储能单元是唯一既能从母线吸收功率、又能向母线注入功率的环节。我用一句话总结这个拓扑的本质：直流母线是一个功率汇集与平衡节点，谁出力多谁出力少，最终都要靠储能把差值吃掉。仿真时只需要关注母线电压是否稳定，就能间接判断功率平衡是否成立。
风力和光伏单元在模型中都被设计成可控电流源的形式，这样便于和储能控制器配合。超级电容单元则通过双向DC-DC变换器接入，变换器的作用是转换电压等级，同时也是控制功率流向的执行机构。
2.2 超级电容的等效电路：怎么从物理变成仿真
超级电容建模是整个模型里最需要仔细琢磨的部分。它不像电阻电容那么简单，但也不像电池模型那么复杂。我做模型时用的是经典RC等效电路：一个理想电容C，一个串联等效电阻R_esr，一个并联等效电阻R_ep。
理想电容C对应超级电容存储电荷的能力，串联电阻R_esr代表极板、端子和电解质带来的欧姆压降，并联电阻R_ep代表漏电流，主要影响长时间静置时的自放电。在秒级到百秒级的储能动态仿真里，R_ep的影响非常小，我通常直接忽略，只保留C和R_esr，这样模型更简洁，状态变量也少。
对应的数学模型很简单：
% 超级电容RC模型状态方程 dUc_dt = I_uc / C; % 忽略漏电阻时的电容电压变化率 V_uc = Uc + I_uc * R_esr; % 端口电压 = 电容电压 + 内阻压降注意符号约定：I_uc的正方向定义为超容放电时电流流出端口。放电时I_uc为正，R_esr上的压降会拉低端口电压；充电时I_uc为负，端口电压会高于电容电压。这个细节决定后续DC-DC变换器模型的输入电压是否正确，很容易搞反。
储能能量的计算是另一个关键点。超级电容存储能量与电压的平方成正比，E = 0.5CU²。假设模组最高工作电压是200V，如果放电到100V，释放的能量不是一半，而是75%，因为0.5C(200²-100²) = 0.5C30000，占满充能量40000的75%。这直接影响容量配置计算，我在下一节实操里会给出具体算例。
2.3 双向DC-DC变换器：能量的"双向门"
超级电容的端电压很低，一个几十伏的模组不可能直接挂在400V的直流母线上，中间一定要有双向DC-DC变换器。这里我用的是经典的两相交错Buck-Boost拓扑，超容侧接低压端，母线侧接高压端。
Boost模式用于超容放电：开关管导通时，超容给电感充磁；开关管关断时，电感通过续流二极管向母线侧释放能量，母线电压被抬升。Buck模式用于超容充电：母线侧功率通过开关管斩波，降压后给超容充电。仿真模型里如果直接用理想开关器件建模，开关频率会迫使求解器采用微秒级步长，整套风光储模型跑起来非常慢。
所以我把变换器做成平均模型，只保留能量平均关系，用受控电压源和受控电流源替代开关器件。平均模型的核心关系在Boost模式下是：
- 母线侧电压 = 超容侧电压 / (1 - d)，其中d是占空比
- 超容侧电流 = 母线侧电流 / (1 - d)
这样我只需要控制占空比d，就能控制功率流向和大小，仿真速度能提高一个数量级，而系统级动态特性基本不受影响。这个取舍在工程仿真里非常常见——研究储能对整个系统的影响时，没人在意每个开关周期内电流怎么纹波。
2.4 控制策略：低通滤波功率分配加双闭环稳压
储能控制器是整套模型的灵魂。我采用的策略是低通滤波功率分配 + 母线电压外环、电流内环双闭环。
先说低通滤波功率分配。风光发电总功率P_total会先经过一阶低通滤波器，得到平滑参考功率P_smooth。高频波动分量就是P_total减去P_smooth，这部分指令交给超级电容。也就是说，超容吸收或释放的功率参考值P_uc_ref = P_total - P_smooth。这个策略的逻辑本质是把功率按频率分开：低频分量由原动机和电网承担，高频分量由超容承担。低通滤波器的时间常数T决定了"高频"的边界，T越大，超容承担波动的频带越宽，容量消耗也越快。我做仿真时T取10秒到20秒之间，能覆盖大部分风光波动周期。
再说母线电压控制。在独立微网场景下，母线电压需要储能来维持。我使用电压外环加电流内环：电压外环检测母线电压误差，输出超容电流参考值；电流内环跟踪参考电流，输出占空比给DC-DC变换器。外环带宽要比内环低一个数量级，避免两个环路打架。母线电压控制在400V，电压外环带宽设20Hz，电流内环带宽设1kHz，这个主从分离的设计能让控制稳定性和动态响应兼得。
3. 在Matlab/Simulink里一步步把模型搭出来
3.1 准备工作：环境与基础模块选型
在动手之前，先把环境准备好。Simulink模型用Matlab R2023b或更新的版本都能跑，Simscape Electrical工具箱不是必须的——因为我自己搭的等效电路模型更轻，也更容易调试。如果大家想省事，用Simscape里的现成模块也行，但注意模型的运行速度会有差异。
新建模型时我建议先按子系统把骨架搭好：wind_system、pv_system、supercap_storage、inverter、controller五个主要模块。先放空壳，再逐个填内容，这样调试时能单独验证每个子系统，不会一上来就面对一团乱线。
所有仿真步长先统一设置为1毫秒采样，有助于后续控制器离散化设计。
3.2 光伏单元与MPPT：让光伏先"听话"
光伏阵列模块我推荐直接用Simscape Electrical中的PV Array模块，参数设置简单：光伏组件开路电压Voc、短路电流Isc、最大功率点电压Vmp、最大功率点电流Imp，然后把串并联数量配好。例如，要构成一套额定功率2kW的光伏阵列，组件Voc=45V，Isc=8A，Vmp=36V，Imp=7A，用8块串联组成360V阵列，再并联1串，就能接近这个功率等级。
MPPT控制我用经典的扰动观察法。它实现简单、鲁棒性好，光伏仿真中非常成熟。核心逻辑是周期性扰动工作点电压，比较功率变化方向，决定下一步扰动方向。在Matlab Function里写这个逻辑非常方便。
function D_ref = mppt_pao(V_pv, I_pv, D_old, dD, V_old, P_old) P_now = V_pv * I_pv; if P_now > P_old if V_pv > V_old D_ref = D_old - dD; else D_ref = D_old + dD; end else if V_pv > V_old D_ref = D_old + dD; else D_ref = D_old - dD; end end endBoost变换器的占空比由MPPT输出直接控制，光伏阵列运行在最大功率点附近。这样光伏单元在模型中就是一个"能跟随环境变化的受控电流源"。
3.3 风电单元：风速模型到直流侧的简化
风电单元如果做全细节机电模型（风速、桨距角控制、永磁同步发电机、背靠背变流器），模型复杂度会成倍增加。我做系统级储能仿真时，用的是简化模型：风速模型 + 空气动力学转矩公式 + PMSG整流后等效到直流侧的受控功率源。
风速模型用四分量组合法，基本风加阵风加渐变风加随机风。基本风设为10m/s，阵风幅值3m/s，渐变风持续5秒，随机风用白噪声滤波生成。这样合成的风速曲线接近真实情况，又不会复杂到难以调试。风机输出机械功率按经典公式P = 0.5ρACpv³计算，Cp取0.4，风轮半径2米，功率就能和2kW的光伏系统匹配起来。
在仿真模型里，我直接把风电机组的输出接到直流母线，用一个受控电流源注入功率。整流器和变流器的动态响应比储能控制的响应快很多，在系统级时间尺度上可以忽略，所以这个简化既能体现风速波动，又不拖慢仿真速度。
3.4 超级电容模组建模：三个关键参数
超级电容模组是我亲手搭的最重要的模块。在Simulink里用积分器加增益模块实现RC等效电路：输入为超容电流I_uc，输出为电容电压Uc和端口电压V_uc。
容量配置的算例我放在这里，非常实用。假设直流母线400V，超容模组最高工作电压250V，最低工作电压125V（半压），需要支撑1000W的功率差额持续10秒。所需能量E = 1000W × 10s = 10000J。根据能量公式反推容量：
C = 2*E / (Umax^2 - Umin^2) C = 2*10000 / (250^2 - 125^2) C = 20000 / 46875 ≈ 0.43F数值不大对吧？这就是功率型储能的特性，平抑秒级波动所需的电容能量并不算大，但瞬时功率能力要足够强。模型里我取1F作为设计值，留出裕量，模组串联等效电阻ESR按0.05Ω设置。如果读者做时间更长的支撑，比如30秒，容量就要按比例提升到1.3F左右。
建模时一定记得在积分器的初始条件里设置电容初始电压，否则仿真第一步就是NaN。初始电压建议设在0.7倍的额定电压，比如250V的0.7倍就是175V，这样可以兼顾充放电裕量。
3.5 双向DC-DC与控制器参数整定
双向DC-DC平均模型的搭建方法：用Controlled Voltage Source表达母线侧电压，用Controlled Current Source表达输入电流。占空比d由控制器输出。在Boost平均模型中，母线侧电压受控源电压 = V_uc / (1-d)，超容侧电流受控源电流 = I_bus / (1-d)。
控制器采用离散PI实现，在Simulink中直接用Discrete PID Controller模块。电流内环的PI参数我按频域带宽法整定。假设电感L=2mH，期望内环带宽1kHz，比例增益Kp可以按L*带宽估算，大约0.2左右，积分增益Ki取Kp的百分之一到十分之一，先设20。电压外环带宽20Hz，比例增益Kp=2，积分Ki=10。这套参数只是起始值，跑阶跃响应后还要微调。
调试时我有个经验：先断开MPPT和功率分配，只让储能做母线电压控制，给母线接一个固定负载，观察电压阶跃响应。如果母线电压振荡，说明外环Kp偏大；如果稳态有静差，加大Ki。电流内环同理，电流响应慢就加Kp，但注意别让振荡。
3.6 仿真结果怎么看：波形与分析要点
模型搭好后，仿真时长设120秒，求解器选ode23t，相对误差1e-3。先在光照和风速都变化的条件下运行，重点看四组波形：
第一组是风力和光伏的输出功率曲线，可以看到明显的波动。第二组是超级电容的功率指令与实际功率，观察超容是否跟踪了高频波动分量。第三组是直流母线电压，这个波形最能说明储能效果——不加储能时母线电压可能摆动在380V到420V之间，加上储能控制后应稳定在400V左右，误差不超过2%。第四组是超级电容的SOC变化曲线，用电压代替SOC观察，确认超容电压始终在最低工作电压之上。
这四个波形全部正常，模型就算基本跑通了。
4. 调试实录：踩坑记录与排查方法
4.1 模型一运行就报NaN
我最初把整个系统连起来时，仿真第一步就报"Input data contains NaN"。排查发现是超级电容模块的问题——积分器初始条件没设置，或者双向DC-DC模块里出现了除零操作。占空比初始值设为0.5，所有状态变量的初始条件都显式赋好值，问题立刻消失。
4.2 功率分配振荡，母线电压还是稳不住
现象是母线电压周期性波动，超容电流来回充放，但平抑效果很差。我看波形后发现低通滤波器时间常数T设成了2秒，和电压外环带宽20Hz对应的控制频率靠得太近，功率分配环路与电压控制环路形成了正反馈。把T调到10秒以上，振荡就消失了。
4.3 超级电容电压掉到下限以下
有几次跑长时段仿真，超容电压一路掉到125V以下，这时候超容基本失去了储能能力，母线电压开始失控。原因是功率指令限幅没做，低通滤波后的高频分量过大，导致超容过放。解决办法两招：一是给超容功率参考加限幅器；二是在控制器里加过压欠压保护逻辑，用Switch判断超容电压范围，越界时将功率参考强制归零。
| 常见问题 | 可能原因 | 解决办法 |
|---|---|---|
| 仿真NaN | 积分器初始条件缺失、除零 | 显式设置初始条件，占空比初始0.5 |
| 功率分配振荡 | 滤波器时间常数太小、环路带宽重叠 | T调到10s~20s，添加限幅 |
| 超容电压越限 | 容量配置不足、指令无限制 | 加大容量，加功率限幅和保护逻辑 |
| 仿真非常慢 | 开关级DC-DC模型、步长过小 | 改用平均模型，用ode23t |
| 波形毛刺明显 | PI输出变化过快、占空比跳变 | 给占空比加速率限制器 |
4.4 仿真跑得太慢怎么办
如果大家用理想开关器件搭DC-DC，仿真会非常卡。我一开始也犯过这个错误，开关频率设20kHz，固定步长必须压到1e-6秒才能收敛，跑120秒仿真要等很久。后来把DC-DC改成2.3节描述的平均模型，仿真速度提升了将近一个数量级。
4.5 控制器参数的一个通病
电流内环PI参数如果Kp设太大，占空比会饱和，波形出现严重毛刺。这时候不要盲目调参数，先检查限幅器是否生效，再给占空比输出加一个速率限制器，比如每个控制周期变化不超过0.01。这个小技巧能解决大部分波形质量问题。
5. 我的一些实操体会
做这套超级电容储能模型，我最深的感受是：仿真模型的难点不在模块本身，而在参数之间的匹配。低通滤波器时间常数、PI带宽、电容容量、功率指令限幅，这些参数单独看都合理，但合在一起需要反复联调。我建议大家一开始先跑一个小规模的直流微网——一个DC源加一个超容储能加一个负载，把母线电压控制调到满意，再接风电和光伏，不然问题混在一起很难定位。
参数设置的先后顺序也很重要。我的习惯是先定母线电压和功率等级，再定超容容量和电压范围，然后定DC-DC电感参数，最后整定控制器。这个顺序保证每个参数都有依据，不会出现"调Kp调Ki调了一个小时却发现是容量不够"这种尴尬。
还有一点值得提醒：超容容量不是越大越好。容量大确实让平抑效果更好，但成本高、体积大，而且大容量会让低通滤波器时间常数牵一发而动全身。根据实际波动的能量频谱去配置容量，才是专业做法。这套模型后续如果要扩展，很自然的路线是加入锂电池做混合储能，让电池管长时间能量转移、超容管短时间功率波动，再用一个能量管理策略把两条储能支路协调起来。有兴趣的话，可以沿着这个方向继续玩下去。
