基于simulink的PEM燃料电池建模与控制仿真,对比PID,积分分离以及滑膜控制器,含仿真操作录像

clc; clear; close all; %% Common used physical and chemical parameters R=8.314; %单位J/(mol*K); F=96485; %法拉第常数，单位c/mol n=381; %电池片数 r=1.4; %空气比热比 Cp=1e3; %air specific heat p_atm=101325; p_aamb=101325; %环境干空气压力 psat_amb=0.4e5; %环境温度饱和蒸气压 yo2=0.21; %氧气体积分数 Mo2=0.032; %氧气摩尔质量（kg/mol） Mn2=0.028; %氮气摩尔质量（kg/mol） Mv=0.018; %水蒸气摩尔质量(kg/mol) Ma=0.0288; %干空气摩尔质量（kg/mol） Mh2=0.002; %氢气摩尔质量（kg/mol） T_atm=273+25; LHV_H2=119.64e6; %氢气低热值（J/kg） %% Compressor parameters Jcp=5e-5; %空压机转动惯量（kg*m2） Ra=286.9; %空气气体常数（J/(kg*K)） eta1=0.82; %电机效率 eta2=0.85; %空压机效率 kt=0.0153; %电机参数 Rcm=0.82; %电机阻抗 kv=0.0153; %电机参数 %% Supply Manifold Vol=0.02; %SM管道体积 k_nozzle=0.3629e-5; %% Return Manifold Crm=0.0124; %管道系数 Vrm=0.005; %RM管道体积 %% Membrane tm=1.25e-2; % thickness of membrane(cm) ep_m=0.2; % membrane porosity %EW=1.1; % equilient weight %rho_m=1981; % density of membrane lamda_ini=2; %initial value of water content rhom_dry=2000; Mm_dry=1.1; %% Stack gdl=300e-6; %GDL厚度 Afc=280e-4; %电堆表面积(m2) A_fc=280e-4; Cdl=2200; Ea=10000; Vca=0.01; %阴极道体积 Van=0.005; %阳极道体积 d_cl=30e-6; a_i0_c=2.5e3; Co2_ref=40.88; Do2_eff=3e-3; k_ca_out=0.2177e-5; %阴极输出空气系数 delta_p=0.001*101325; %堆内压降 imax=2.2; %最大电流密度 c3=2; c1=10; L_ch=0.1; d_gdl=300e-6; d_cl=35e-6; T_ref=298;