✅作者简介：热爱科研的Matlab仿真开发者，擅长数学建模、模型改进、算法创新、程序设计科研仿真。
🍎 往期回顾关注个人主页：完整代码获取 定制创新 论文复现私信
🍊个人信条：做科研，博学之、审问之、慎思之、明辨之、笃行之，是为：博学慎思，明辨笃行。
🔥 内容介绍
1型糖尿病患者的日常血糖管理长期以来依赖每日多次手动胰岛素注射，很难精准匹配人体动态变化的血糖代谢特性，非常容易出现血糖剧烈波动，长期血糖不达标会引发视网膜病变、肾病、神经病变等严重慢性并发症，而不当的大剂量注射又极易诱发危及生命的严重低血糖。传统纯PID闭环胰岛素输注系统在患者进餐摄入碳水化合物的场景下，很容易出现控制量叠加过度的问题：进餐前注射大剂量胰岛素后，PID反馈控制依然持续追加输注，短时间内体内胰岛素浓度过高，直接导致餐后2~3小时出现严重低血糖，临床使用安全性存在明显缺陷。
本研究基于真实1型糖尿病患者的历史连续血糖监测（CGM）数据、碳水摄入记录、胰岛素输注记录，搭建符合人体生理特性的HUPA代谢仿真模型，设计进餐事件触发的PID自动启停切换机制，构建“进餐大剂量+临时关闭PID+血糖峰值后重启PID”的混合闭环人工胰腺系统，在有效提升血糖达标时间（TIR）的同时，从机制层面大幅降低餐后低血糖的发生风险，完全符合当前临床人工胰腺产品的核心设计逻辑。
人体血糖代谢仿真建模
本研究采用经过临床验证的HUPA 1型糖尿病人体生理代谢模型，完整刻画人体的胰岛素药代动力学与碳水化合物代谢动力学特性，所有模型参数与真实1型糖尿病患者的生理特征完全对齐：
胰岛素效应FIR滤波器建模：采用高斯型因果有限脉冲响应（FIR）滤波器精准模拟皮下注射外源性胰岛素后的体内动态作用过程，完整覆盖胰岛素从皮下吸收、进入血液循环、到最终作用于肝脏与外周组织降低血糖的全时间轴效应，胰岛素的作用时间曲线完全贴合临床公认的速效胰岛素门冬氨酸的药效特性，作用峰值出现在注射后约90分钟，总作用时长约4~5小时。
碳水化合物效应FIR滤波器建模：专门构建碳水化合物吸收代谢FIR滤波器，模拟进餐摄入的碳水化合物从胃肠道逐步吸收入血、升高血糖的动态过程，碳水的血糖响应曲线与临床实测数据高度吻合，摄入后约30分钟开始出现明显升糖效应，120~180分钟左右达到升糖峰值。
随机噪声模型：仿真系统同时引入过程噪声与测量噪声，过程噪声模拟人体代谢的日常随机波动特性，测量噪声模拟真实CGM传感器的测量误差特性，两组仿真对照组共用完全相同的随机噪声序列，保证两组方案的性能对比完全公平，排除随机因素对结果的干扰。
仿真系统的时间步长设置为5分钟，与主流商用连续血糖监测设备的采样周期完全一致，仿真时长默认覆盖24小时完整昼夜周期，完整包含3次常规进餐时段与夜间空腹时段，完全复现真实患者一天的血糖动态变化过程。
进餐切换式混合PID控制器设计
本研究针对性设计带进餐事件自动切换机制的混合PID反馈控制器，从根源上解决传统闭环系统餐后低血糖的痛点，控制器的完整工作逻辑如下：
基础PID参数整定：采用位置式PID控制结构，控制器的输入为当前实时血糖值与目标血糖参考值（临床推荐设置为5.5~6.5mmol/L）的偏差，控制器的输出为基础胰岛素的每分钟输注速率。控制器内置抗积分饱和机制，当胰岛素输注量达到泵的最大基础速率上限时，自动停止积分项的累加，彻底避免传统PID的积分饱和超调问题。
进餐事件识别与PID自动关闭逻辑：系统自动识别所有与进餐相关的大剂量胰岛素注射事件，当用户手动输入进餐碳水总量、注射进餐大剂量之后，系统自动临时关闭PID反馈控制回路，将PID的输出强制锁定为预设的基础胰岛素输注基线值，同时暂停PID控制器的积分项更新，避免反馈控制在碳水吸收上升阶段额外追加胰岛素，防止胰岛素过量叠加。
PID自动重启机制：系统持续实时监测血糖的动态变化趋势，当检测到血糖经过峰值开始出现连续下降趋势，或者距离进餐大剂量注射的时间超过预设的安全超时阈值（临床推荐设置为150分钟）后，自动重新开启PID反馈控制回路，同时将重启时刻的积分项初始值设置为匹配当前基础输注率的对应数值，避免PID重启瞬间出现输出跳变，保证胰岛素输注过程全程平滑过渡。
整套切换逻辑完全不需要用户的额外干预，系统自动完成PID的关闭与重启，操作流程和当前主流商用胰岛素泵的使用习惯完全兼容，患者学习成本极低。
⛳️ 运行结果
🔗 参考文献
🍅更多免费数学建模和仿真教程关注领取
🏆团队擅长辅导定制多种科研领域MATLAB仿真，助力科研梦：
#各类智能优化算法改进及应用
#生产调度 #经济调度#装配线调度#充电优化#车间调度#发车优化#水库调度#三维装箱#物流选址#货位优化#公交排班优化#充电桩布局优化#车间布局优化#集装箱船配载优化#水泵组合优化#解医疗资源分配优化#设施布局优化#可视域基站和无人机选址优化#背包问题#风电场布局#时隙分配优化#最佳分布式发电单元分配#多阶段管道维修#工厂-中心-需求点三级选址问题 #应急生活物质配送中心选址#基站选址#道路灯柱布置#枢纽节点部署#输电线路台风监测装置#集装箱调度 #机组优化 #投资优化组合 #云服务器组合优化#天线线性阵列分布优化#CVRP问题#VRPPD问题#多中心VRP问题#多层网络的VRP问题#多中心多车型的VRP问题 # 动态VRP问题 #双层车辆路径规划（2E-VRP） #充电车辆路径规划（EVRP） #油电混合车辆路径规划#混合流水车间问题#订单拆分调度问题#公交车的调度排班优化问题#航班摆渡车辆调度问题#选址路径规划问题#港口调度#港口岸桥调度#停机位分配#机场航班调度#泄漏源定位#冷链#时间窗#多车场等#选址优化#港口岸桥调度优化#交通阻抗#重分配#停机位分配#机场航班调度#通信上传下载分配优化
#机器学习和深度学习时序#回归#分类#聚类和降维
#bp时序#回归预测和分类
#ENS声神经网络时序#回归预测和分类
#SVM#CNN-SVM#LSSVM#RVM支持向量机系列时序
#CNN#TCN#GCN卷积神经网络系列时序
#ELM#KELM#RELM#DELM极限学习机系列时序
#GRU#Bi-GRU#CNN-GRU#CNN-BiGRU门控神经网络时序
#ELMAN递归神经网络时序
#LSTM#BiLSTM#CNN-LSTM#CNN-BiLSTM/长短记忆神经网络系列时序
#RBF径向基神经网络时序
