请高人指点这个微分方程组用matlab怎么解.dw/dt = -k1*w,ds/dt= k1*w-k2*s,w的积分范围是从w0 积分到w,w0是初始值,w0大于w需要求出w,s

来源:学生作业帮助网 编辑:作业帮 时间:2024/05/02 01:56:18
请高人指点这个微分方程组用matlab怎么解.dw/dt = -k1*w,ds/dt= k1*w-k2*s,w的积分范围是从w0 积分到w,w0是初始值,w0大于w需要求出w,s

请高人指点这个微分方程组用matlab怎么解.dw/dt = -k1*w,ds/dt= k1*w-k2*s,w的积分范围是从w0 积分到w,w0是初始值,w0大于w需要求出w,s
请高人指点这个微分方程组用matlab怎么解.
dw/dt = -k1*w,
ds/dt= k1*w-k2*s,
w的积分范围是从w0 积分到w,w0是初始值,w0大于w
需要求出w,s

请高人指点这个微分方程组用matlab怎么解.dw/dt = -k1*w,ds/dt= k1*w-k2*s,w的积分范围是从w0 积分到w,w0是初始值,w0大于w需要求出w,s
>> w=dsolve('Dw=-k1*w','w(0)=w0')
w =w0*exp(-k1*t)
>> dsdt= 'k1*w-k2*s'
dsdt =k1*w-k2*s
>> dsdt=subs(dsdt)
dsdt =k1*w0*exp(-k1*t)-k2*s
>> s=dsolve('Ds=k1*w0*exp(-k1*t)-k2*s')
s =(k1*w0/(-k1+k2)*exp(-t*(k1-k2))+C1)*exp(-k2*t)