热电联产CGAM问题建模

%decision variable
T4=1492.63;
T3=914.28;
efficiency_ac=0.8468;
efficiency_gt=0.8786;
r_cp=8.5234;


% initialization
p0=1.013; 
p1=1.013;
p5=1.099;
p6=1.066;
p7=1.013;
p8=20;
p8p=20;
p9=20;


T0=298.15;
T1=298.15;
T8=298.15;
T8p=470.52;
T9=485.52;

Wnet=30000;
delta_Pap=0.05;
delta_Pcb=0.05;
efficiency_cb=0.98;

Cpa=1.004;
r_air=1.4;
LHV=50000;
e_fuel=51850;
Cpg=1.17;
r_gas=1.33;
Msteam=14;
h98=2690;
h98p=1956;
e98=910;


% Thermodynamic model 
T2=T1*（1+ （r_cp^（（r_air-1）/r_air）-1）/efficiency_ac）
p2=r_cp*p1
p3=（1-delta_Pap）*p2
p4=（1-delta_Pcb）*p3
T5=T4*（1-efficiency_gt*（1-（p4/p5）^（（1-r_gas）/r_gas）））
T6=T5-（Cpa*（T3-T2）*（Cpg*（T4-T0）-efficiency_cb*LHV））/（Cpg*（Cpa*（T3-T0）-efficiency_cb*LHV））

% （Cpg*（T4-T0）-efficiency_cb*LHV）/（Cpa*（T3-T0）-efficiency_cb*LHV）
Mgas=Wnet/（Cpg*（T4-T5）-Cpa*（T2-