UC知道谁能用matlab帮我解个微分方程组啊？详细点啊，我是新手啊！急啊！！
来源：百度知道 编辑：UC知道 时间：2024/05/18 12:51:58
d(T)/d(Z)=(RA*(-DHR)-4*U*(T-TC)/dt)/(G*CPT)
RA=C1*C2*C3*KE*((PE-PS*PH/KP)+PE*(KT+KB))/(1+LEME*PE+LEMS*PS)^2*ROU
其中：Z范围[0,40] XA初值=0,T初值=793.15K
PE=34.1952*(1-XA)/(356.452+37.477*XA)*P
PS=(31.801*XA+0.0208)/(356.452+37.477*XA)*P
PH=38.982*XA/(356.452+37.477*XA)*P
KE=5.7e8*exp(-34420/(R*T))
KT=1.3e3*exp(-18880/(R*T))
KB=9.1e7*exp(-38190/(R*T))
KP=exp(16.12-15350/T)
LEME=3.6e-3*exp(-40780/(R*T))
LEMS=1.5e-3*exp(-16430/(R*T))
C1=((3624.641-61.549*XA)/2265.401*(1-XA)/2.45)^0.25
C2=0.247
C3=0.95
R=1.987
G=9430.217/(250*3.14159/4*0.123^2)
WAO=0.3884
ROU=1100
MA=106
CPT=2.22
U=20.02
TC=923.15
DHR=30850
P=0.155/0.101325
dt=0.123
关键是要解100组啊!!!
终于完成了。
function hahha
%global XA T G WAO MA DHR U TC dt G CPT RA
syms XA T
C2=0.247
C3=0.95
R=1.987
G=9430.217/(250*3.14159/4*0.123^2)
WAO=0.3884
ROU=1100
MA=106
CPT=2.22
U=20.02
TC=923.15
DHR=30850
P=0.155/0.101325
dt=0.123
PE=34.1952*(1-XA)/(356.452+37.477*XA)*P
PS=(31.801*XA+0.0208)/(356.452+37.477*XA)*P
PH=38.982*XA/(356.452+37.477*XA)*P
KE=5.7e8*exp(-34420/(R*T))
KT=1.3e3*exp(-18880/(R*T))
KB=9.1e7*exp(-38190/(R*T))
KP=exp(16.12-15350/T)
LEME=3.6e-3*exp(-40780/(R*T))
LEMS=1.5e-3*exp(-16430/(R*T))
C1=((3624.641-61.549*XA)/2265.401*(1-XA)/2.45)^0.25
RA=C1*C2*C3*KE*((PE-PS*PH/KP)+PE*(KT+KB))/(1+LEME*PE+LEMS*PS)^2*ROU
dXAdZ=RA/(G*WAO/MA)
dTdZ=(RA*(-DHR)-4*U*(T-TC)/dt)/(G*CPT)
[Z,XA_T] = ODE45(@myfun,[0:40],[0,793.15])
plotyy(Z,XA_T(:,1),Z,XA_T(:,2))
function dy=myfun(t,x)
%global XA