?%test?計(jì)算PE（0）
Fb0=0.4412;
P0=0.897;
T0=650;?%初始溫度
beta=2000;?%升溫速率
%不穩(wěn)橋在全部鏈接中的概率
PB=Fb0*P0;
%計(jì)算反應(yīng)物中不穩(wěn)定橋的概率
w=[P0Fb0];?%w（1）---p（0）--完整鏈接在全煤中的分?jǐn)?shù)???w（2）--Fb（0）--不穩(wěn)定橋在完整鏈接中所占比例
Prb0=PRb（w）;
ww=[Prb0P0];
BR_0=BR（ww）;
%計(jì)算末端含有官能團(tuán)的概率
%x=[275101Fb0*P01.8560.9110.8387]?%x（1）--MWmon-單體單元分子量???x（2）--MWA--核的分子量??x（3）--pb--不穩(wěn)定橋占完整鏈接比例??x（4）--MWB/MWA???x（5）--p????x（6）--MWC/MWA%
%Pe0=pe（x）
Pe0=1;
%計(jì)算反應(yīng)物中的外圍結(jié)構(gòu)
y=[Pe0P0];?%y（1）--pe（0）??y（2）--p（0）
SR0=SR（y）;?
%計(jì)算后形質(zhì)體的摩爾數(shù)

m5_0=m5（P0）;
m4_0=m4（P0）;
m3_0=m3（P0）;
m2_0=m2（P0）;
m1_0=（1-P0）^2;

%計(jì)算后形質(zhì)體的不穩(wěn)定橋

z=[Fb0P0];
BM5_0=BM5（z）;
BM4_0=BM4（z）;
BM3_0=BM3（z）;
BM2_0=BM2（z）;
BM1_0=0;

%計(jì)算SM
SM0=SM（y）;

%計(jì)算中間體
n10_0=n10（P0）;
n9_0=n9（P0）;
n8_0=n8（P0）;
n7_0=n7（P0）;
n6_0=n6（P0）;

%計(jì)算中間體的不穩(wěn)定橋
BN10_0=BN10（z）;
BN9_0=BN9（z）;
BN8_0=BN8（z）;
BN7_0=BN7（z）;
BN6_0=BN6（z）;
%計(jì)算中間體的外圍結(jié)構(gòu)
SN_0=SN（y）;

%初始非冷凝氣體?表III??fuel?and?energy?55-1991663
YG_0=0;
%計(jì)算KB不穩(wěn)定橋斷裂速率常數(shù)的表達(dá)式及值?在T溫度的值
t=[0:0.014:0.275];
%T=t*beta+T0;
%kB=kB_splinet（T）;
%外圍結(jié)構(gòu)脫除反應(yīng)速率
%kG=1d15*exp（（-230000）*power（8.314*T-1））;?%kG--氣體逸出速率常數(shù)表達(dá)式