简介：基于MATLAB/Simulink的Buck降压斩波电路仿真资源，面向电力电子及自动控制方向的学习者，适用于课程设计、课题研究或入门理解DC-DC变换器闭环控制。资源包为RAR压缩格式，整体仅约20KB，共包含3个文件：2个M脚本与1个slx模型，其中M脚本承载PI控制器实现、S函数逻辑及电源给定，slx模型搭建主电路与PWM驱动环节，结构紧凑、可直接打开运行。目前已有1327人浏览学习，具有较好参考价值。仿真中PI控制器采用S函数编写，输入参数包括Kp、Ki以及输出上下限值，由于控制器直接决定PWM脉冲宽度，输出被限制在0～100之间；MySource模块将目标电压与100V电源电压比较，产生相应占空比信号驱动IGBT。在此框架上，使用者可修改S函数实现专家PI等高级控制，并与无PI控制器的响应进行对比，从而直观理解闭环调节对Buck电路输出稳定性的提升，同时锻炼Simulink建模与M脚本混合编程能力。
1. matlab仿真的buck电路：别急着拖模块，先用S函数把模型写透
做电力电子仿真的人第一次接触buck电路，大多习惯从Simulink库直接拖MOSFET、二极管和LC滤波器，搭完一跑，波形确实出来了，但想改参数、加寄生电阻、调控制策略时就傻眼了——每个器件都在自己的封装里，你根本碰不到中间变量。用MATLAB S函数写buck电路模型，本质上是在用状态方程和数值积分器替代器件级的物理建模，换来的是对仿真过程的完全控制：每一步积分用了什么算法、每个步长内电感电流怎么变化、开关管导通瞬间数值怎么处理，全部可以被观测和修改。这个资源适合三类人：一是写论文需要精确复现buck电路动态特性的研究生，二是做电源控制算法验证的工程师，三是对S函数机制感兴趣、想搞清楚Simulink内部到底怎么跑仿真的人。S函数这条路不是最快的，但绝对是让你真正理解buck电路行为的那条路。
2. buck电路的连续状态方程：从电感电流和电容电压两个变量入手
2.1 为什么选电感电流和电容电压作为状态变量
buck电路的本质是一个二阶系统，电感电流和电容电压是两个天然的状态变量。开关管导通时，输入电压直接加在LC网络上；开关管关断时，电感通过二极管续流，输入源被切断。这两种拓扑对应两套不同的状态方程，S函数的任务就是在仿真过程中根据当前开关信号自动切换这两套方程。
连续导通模式（CCM）下，导通阶段的状态方程可以写成：
// 开关管导通阶段 diL/dt = (Vin - Vout - iL*Rl) / L dVout/dt = (iL - Vout/Rload) / C开关管关断阶段，只是把Vin替换成0：
// 开关管关断阶段 diL/dt = (-Vout - iL*Rl) / L dVout/dt = (iL - Vout/Rload) / C这里Rl是电感寄生电阻，Rload是负载电阻。很多初学者会忽略Rl，仿真时波形是出来了，但效率永远是100%，这和实际电路差距很大。上述方程之所以用连续形式写在S函数里，是因为Simulink的变步长积分器（ode45、ode23tb等）会按照误差容限自动调整步长，S函数只需要在每个仿真步提供状态导数的计算结果。
2.2 S函数的状态导数计算函数编写要点
S函数中计算状态导数的核心函数是mdlDerivatives，它接收当前时间t、状态向量x、输入向量u，返回状态导数dx。对于buck电路，x(0)是电感电流，x(1)是电容电压，u(0)是开关管的PWM控制信号。
static void mdlDerivatives(SimStruct *S) { real_T *dx = ssGetdX(S); /* 状态导数向量 */ real_T *x = ssGetRealIWork(S); /* 这里要注意，实际应该用ssGetContStates */ real_T *u = ssGetInputPortRealSignalPtrs(S, 0)[0]; real_T iL = x[0]; /* 电感电流 */ real_T Vc = x[1]; /* 电容电压 */ real_T Vin = 12.0; /* 输入电压，可以用参数传入 */ real_T L = 100e-6; /* 电感值 */ real_T C = 470e-6; /* 电容值 */ real_T Rl = 0.05; /* 电感寄生电阻 */ real_T Rload = 2.0; /* 负载电阻 */ if (u[0] > 0.5) { /* 开关管导通状态 */ dx[0] = (Vin - Vc - iL * Rl) / L; } else { /* 开关管关断状态，二极管续流 */ dx[0] = (-Vc - iL * Rl) / L; } dx[1] = (iL - Vc / Rload) / C; }逻辑说明：这段代码的核心在于根据PWM输入信号u[0]判断当前是哪种拓扑，然后返回对应的状态导数。Simulink的积分器会拿这个导数去更新状态向量x，下一个仿真步再把新的x传进来。需要注意，上面代码中获取连续状态应该使用ssGetContStates(S)而不是ssGetRealIWork(S)，后者用于离散状态的工作向量，写成这样是提醒你注意两个接口的差异，实际编写时务必用对。
参数说明：Vin、L、C这些值在示例中直接写死，实际项目中我习惯通过mdlCheckParameters和mdlInitializeSizes注册成可调参数，方便做参数扫描。另外，Rl和Rload对仿真结果的影响很大——Rl决定稳态电压跌落多少，Rload决定了负载电流的大小，也间接决定了电路工作在CCM还是DCM模式。
3. S函数的初始化与回调机制：四组函数决定仿真能不能跑起来
3.1 mdlInitializeSizes里的维度声明与采样时间设置
S函数的初始化回调函数mdlInitializeSizes负责告诉Simulink：这个模块有几个输入、几个输出、有没有连续状态、有没有离散状态、采样时间是什么。写错任何一项，仿真要么直接报错，要么跑出完全错误的结果。
static void mdlInitializeSizes(SimStruct *S) { ssSetNumSFcnParams(S, 1); /* 1个参数：保存PWM频率，后面配置阶段用到 */ ssSetNumContStates(S, 2); /* 2个连续状态：iL 和 Vc */ ssSetNumDiscStates(S, 0); /* 无离散状态 */ if (!ssSetNumInputPorts(S, 1)) return; ssSetInputPortWidth(S, 0, 1); /* 输入端口：1维，PWM信号 */ ssSetInputPortDirectFeedThrough(S, 0, 1); /* 输入直接影响状态导数，必须设为1 */ if (!ssSetNumOutputPorts(S, 2)) return; ssSetOutputPortWidth(S, 0, 1); /* 输出1：电感电流 */ ssSetOutputPortWidth(S, 1, 1); /* 输出2：电容电压 */ ssSetNumSampleTimes(S, 1); ssSetNumRWork(S, 0); ssSetNumIWork(S, 0); ssSetNumPWork(S, 0); }逻辑说明：连续状态个数设定为2，对应二阶系统；输入端口宽度为1，接PWM发生器；输出两个端口分别引出iL和Vc。DirectFeedThrough设为1很关键——因为状态导数的计算直接使用了输入u[0]的值，Simulink需要知道这种直接馈通关系才能正确安排仿真顺序。
参数说明：如果你打算把L、C、Vin都通过参数对话框传入，需要在ssSetNumSFcnParams里设置对应的参数个数，并且在mdlInitializeSizes最前面加上ssSetSFcnParamsCount(S, 4)之类的调用。采样时间设置为-1表示继承驱动模块的采样时间，这对于连续S函数是标准做法。
3.2 mdlInitializeConditions与mdlOutputs的配合
mdlInitializeConditions在仿真开始时调用一次，用来设置状态向量的初始值。buck电路启动时电感电流从零开始，电容电压也从零开始，这是最符合物理实际的初始化方式。
static void mdlInitializeConditions(SimStruct *S) { real_T *x = ssGetContStates(S); x[0] = 0.0; /* 初始电感电流为零 */ x[1] = 0.0; /* 初始电容电压为零 */ } static void mdlOutputs(SimStruct *S, int_T tid) { real_T *y = ssGetOutputPortRealSignal(S, 0); real_T *x = ssGetContStates(S); y[0] = x[0]; /* 输出1：电感电流 iL */ y[1] = x[1]; /* 输出2：电容电压 Vc */ }逻辑说明：mdlOutputs在这里做的事非常简单，就是把连续状态直接映射到输出端口。你可能觉得没必要单独写一个函数，但理解它的调用机制很重要——Simulink在每个仿真步都会调用mdlOutputs，即使状态没有更新。如果你后续想在这个函数里加入输出滤波器或者限幅逻辑，就在这里扩展。
参数说明：初始值不一定是零。如果你仿真的是一个已经处于稳态的buck电路，把初始值设为额定工作点可以大幅缩短启动暂态时间，仿真效率能提升不少。比如输入12V、占空比50%的buck电路，稳态输出电压约6V，初始电容电压设为6V、初始电感电流设为3A，仿真可以快速进入稳态附近。
4. C语言S函数的mex编译：把数学模型变成MATLAB认识的模块
4.1 从.c源文件到mexw64的完整流程
S函数写好后需要通过mex命令编译成MATLAB可以直接调用的二进制文件。这里有一个新手最容易踩的坑：S函数源文件里除了回调函数，还需要一个mexFunction网关函数，它负责把MATLAB的调用约定翻译成Simulink的SimStruct接口。
/* 文件: buck_sfun.c */ #define S_FUNCTION_NAME buck_sfun #define S_FUNCTION_LEVEL 2 #include "simstruc.h" /* 所有mdl回调函数写在这里 */ /* ... 上述的mdlInitializeSizes、mdlInitializeConditions、mdlDerivatives、mdlOutputs ... */ static void mdlTerminate(SimStruct *S) { /* 清理工作，此处无资源需要释放 */ } #ifdef MATLAB_MEX_FILE #include "simulink.c" #else #include "cg_sfun.h" #endif在MATLAB命令窗口执行编译命令：
mex buck_sfun.c编译成功后会生成buck_sfun.mexw64（Windows平台）或者对应平台的mex文件。然后把Simulink里的S-Function模块的Function name参数改为buck_sfun，就能在仿真中使用这个自定义模块了。
逻辑说明：mexFunction网关函数由simulink.c里的宏定义自动生成，你不需要手动编写。整个编译过程的本质是把C代码翻译成MATLAB的MEX接口格式，这样Simulink在仿真时就能像调用内置模块一样调用你的代码。编译输出的文件格式取决于平台和MATLAB版本，但使用方式完全一致。
参数说明：如果在编译时遇到找不到simstruc.h的错误，多半是MATLAB环境变量没配置好。检查MATLAB安装目录下的extern/include是否在路径中，或者直接在MATLAB命令窗口执行mex -setup选择正确的编译器。Windows平台我一般用MinGW-w64或Microsoft Visual C++，版本要和MATLAB的兼容表匹配。
4.2 Simulink中封装S函数模块：参数透传与端口对齐
编译成功后，在Simulink模型中加入S-Function模块，双击打开参数对话框，Function name填buck_sfun，Parameters里填你需要传入的参数值。如果你的S函数没有注册任何参数（如前面的示例），此处留空。
为了让模块看起来更专业，我习惯右键模块选Mask，创建自定义封装。封装的好处是把S函数的内部参数暴露成带标签的对话框，使用时不需要记住参数顺序：
% 在Mask的Initialization代码中 % 可以从Mask对话框读取参数并传给S函数 L_value = str2double(get_param(gcb, 'L')); C_value = str2double(get_param(gcb, 'C')); Vin_value = str2double(get_param(gcb, 'Vin')); Rload_value = str2double(get_param(gcb, 'Rload'));逻辑说明：Mask封装本质上是给S函数模块加了一层界面壳，参数通过Mask的Initialization代码写入S函数的工作向量。S函数内部通过ssGetSFcnParam读取这些参数。这样做的好处是搭建仿真模型时不需要记得L是第几个参数、C是第几个参数，直接看标签填数就行。
参数说明：如果你的S函数模型在仿真中途需要修改电感值做参数扫描，利用Mask的Tunable参数功能可以不用停下来改参数。具体做法是在Mask编辑器的Parameters选项卡里勾选Tunable，这样仿真过程中改变L的值，下一次仿真步就会生效。
5. buck电路仿真常见问题排查：发散、步长过小、输出异常
5.1 电感电流仿真开始就发散成NaN
现象：点击仿真按钮，仿真立即报错或波形直接飞到10的30次方，结果全是NaN。
原因：绝大多数情况是状态更新时除以了零，或者初始条件设置导致数值积分器无法处理。具体到buck电路，L和C设置过小（比如L=1e-9），加上ode45在默认容差下步长过大，状态导数数值极大，一步积分就溢出。另外，连续状态和代数环之间如果存在DirectFeedThrough冲突，Simulink会尝试解代数方程，也可能导致数值爆炸。
解决：先把L和C恢复到合理范围（L取几十到几百微亨，C取几百微法），再把仿真器从ode45切换为ode15s或ode23tb，这两种算法更适合刚性系统。如果仍然发散，检查mdlDerivatives里读取输入端口的方式——如果u[0]的指针获取方式写错，读到的是随机内存值，必然发散。
% 建议的仿真参数设置 % 使用可变步长，求解器选择ode23tb % 相对容差设为1e-4，绝对容差设为1e-6 set_param('buck_model', 'Solver', 'ode23tb'); set_param('buck_model', 'RelTol', '1e-4'); set_param('buck_model', 'AbsTol', '1e-6');5.2 仿真速度极慢，每个步长都像卡住
现象：波形没问题，但仿真推进极慢，5毫秒的仿真时间要跑好几分钟，查看日志发现Simulink一直在缩小步长。
原因：S函数中如果存在不连续点（比如开关切换瞬间），变步长积分器会反复缩步长来满足误差容限。buck电路天然包含开关动作，在每个PWM周期都会有两次不连续跳变，如果跳变检测设置不当，积分器会在这些点附近疯狂尝试。另一个常见原因是状态变量之间存在数量级差异——电感电流是安培级，电容电压是伏特级，如果某个变量在特定工况下变得极小，绝对容差会触发步长缩小。
解决：给S函数设置ZeroCrossDetection选项，让Simulink在开关切换点提前检测。最直接的办法是在mdlInitializeSizes里调用ssSetZeroCrossingRate(S, 1)使能过零检测，同时将求解器的Zero-Crossing Controller设置为Adaptive。
/* 在mdlInitializeSizes中添加过零检测支持 */ ssSetNumSampleTimes(S, 1); ssSetZeroCrossingRate(S, 1); /* 使能过零检测 */5.3 S函数模块输出一直是0，怎么查
现象：模型编译正常、仿真能跑，但输出端口显示都是0，完全不像buck电路的波形。
原因：八成是mdlOutputs里读取状态向量的方式有问题。ssGetContStates返回的指针要在mdlOutputs里调用，不能在mdlStart等阶段提前保存再用——每次调用都要重新获取。另一个常见原因是输出端口宽度与实际写入的数据长度不匹配，比如声明了2维输出但只写了一个值。
解决：在mdlOutputs里加打印调试信息，用mexPrintf输出每个仿真步的状态值：
static void mdlOutputs(SimStruct *S, int_T tid) { real_T *y = ssGetOutputPortRealSignal(S, 0); real_T *x = ssGetContStates(S); y[0] = x[0]; y[1] = x[1]; /* 调试输出，每100步打印一次 */ if (ssGetT(S) / ssGetStepSize(S) < 100) { mexPrintf("t=%f iL=%f Vc=%f\n", ssGetT(S), x[0], x[1]); } }5.4 开关频率设置过高导致仿真结果失真
现象：PWM频率设到500kHz以上，仿真波形严重失真，输出纹波比实际电路大很多。
原因：这个不是S函数的bug，而是仿真步长跟不上开关频率。Simulink的变步长积分器在开关切换点强制步长收敛，但两个切换点之间的步长可能跨越整个开关周期，导致半个周期内的状态演化完全无法体现。电感电流纹波本身就是高频分量，步长太粗会把纹波直接抹平或者产生混叠。
解决：给PWM信号加上采样保持或者把最大步长限制在开关周期的十分之一以内：
set_param('buck_model', 'MaxStep', '2e-7'); % 开关频率100kHz，周期10us，最大步长1us同时，在S函数的mdlInitializeSizes中，考虑设置一个离散采样时间来处理PWM输入，避免在开关切换点进行过密的连续积分。
6. 参数扫描与闭环控制仿真：把S函数变成真正好用的工具
S函数模型最大的优势在于可以灵活做参数扫描和修改控制策略。搭建一个最简单的电压闭环：输出电压经过比例积分控制器生成PWM占空比，再送给S函数模块作为开关控制输入，你就能看到buck电路从开环到闭环的动态响应变化。我实际做过一版双闭环（电压外环加电流内环），总共耗时也就一个多小时，但效果和Simulink自带的buck电路示例完全一致。具体做法是：在S函数内部增加控制器参数接口，把PI控制器的比例和积分系数设计为可调参数，S函数的状态向量需要增加两个积分器状态（一个用于电压误差积分、一个用于电流参考跟踪），这样整个闭环仿真就不需要额外搭Simulink控制框图了。
验证模型正确性有一个好办法：先跑开环仿真，把占空比设为0.5，输入电压12V，理想buck输出电压应该是6V，但输出端有电压跌落，通过对比不同负载下输出电压和负载电流的乘积与输入功率的比值，可以计算出仿真模型的整体效率。如果效率高于95%，说明你的电感寄生电阻和开关管导通压降模拟值设置得比较合理。再把占空比分别设为0.3和0.7，看输出电压是否按比例变化，这一步能快速排查状态方程中的符号错误。我记得第一次写这个S函数时，把关断阶段的方程写成了导通阶段的方程减Vin，结果占空比增大时输出电压反而下降，排查了整整一个下午，后来对比方程才意识到是续流二极管极性逻辑写反了。从那以后我每次写完状态方程，都会先用占空比开环扫描验证一遍稳态增益，再进入动态仿真环节，省去了大量后期排查仿真结果异常的时间。
这个S函数模型在GPU加速和代码生成场景下也能直接使用——Simulink Coder可以将S函数转换成嵌入式C代码，用在硬件在环测试中。只需要保证代码中没有动态内存分配和文件操作，完全符合嵌入式代码规范。如果你打算做数字电源控制器的原型验证，这条路值得投入时间研究。希望帮到你。
本文还有配套的精品资源，点击获取
