??syms?W?E?W1?R21?R2?f?s?R3?R4%蒸氣云爆炸模型VEC
X?=?[309.451300];%物質(zhì)的量數(shù)組（KG）
for?i=1:2%有幾組數(shù)據(jù)循環(huán)幾次
X?=?[Xsym（[‘x‘num2str（i）]）];
end
for?j=1:2
E（j）=1.8*0.04*X（j）*50000%50000為燃燒熱（kg/mol）
W1（j）=E（j）/4520
R1（j）=13.6*（（W1（j）/1000）^0.37）;%死亡半徑
RR1（j）=eval（R1（j））
R21（j）=1/（（E（j）/101325）^（1/3））
f1=inline（0.137*（（R21（j）*R2）^（-3））+0.119*（（R21