ds-k1t673m

欧曲滢3183实验测得反应:2A+B→2C+D的速率方程为:r=k[A][B],反应历程为:A...
伍咸彪15140007238 ______ 自己设个方程 代入初值开始解,每次把步长减半,看误差的减小是不是4阶的 例如 y'=f(t,y) y=e^(t)*sin t y'=e^tsint+e^tcost =y+e^tcost 代入初值y(0)=0 把f(t,y)=y+e^tcost代入你的程序 然后算到T=1比如 dt=1/4,1/8,1/16...然后算误差,即数值解-精确解e^tsint的离散2范数 然后看误差在每次缩小一半步长后是不是缩小成之前的1/16.

欧曲滢3183请高人指点这个微分方程组用matlab怎么解.谢谢 -
伍咸彪15140007238 ______ >> 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)