怎么没有人恢复,是不是
这个题太难了 ?
>> F=dsolve('DF=g*(1-F)','F(0)=0') F = 1-exp(-g*t) >> g='Ka*exp(-Q/R/T)*(p-peq)'; >> subs(F) ans = 1-exp(-Ka*exp(-Q/R/T)*(p-peq)*t) >>
呵呵,我比较苯啊,m文件怎样写?F=dsolve('DF=g*(1-F)','F(0)=0')里面的'DF=g*(1-F)'代表什么?还有哪个g呢?
我的m文件是这样写的,不知道对不对。
fuction f=xiqing(t,F)
η=n/2;
g= η(1-F){-Logexp(1-F)}^(η-1))/η
f=ka*exp(-Q/R/T)*g*(p-peq)
然后调用在窗口中键入:[t,F]=ode23(@xiqing,[0 500],0])
其中[0 500]是时间,0是初值 还有就是p,peq,T都是常数,但是要在同一个土中同时显示当为不同值时 的曲线!!谢谢,再麻烦您一下!!
η=n/2还是n 如果η=n=1 那么原方程可以通过分离变量求解 F = 1-exp(-Ka*exp(-Q/R/T)*(p-peq)*t) 如果η=n/2=1/2 你的g= η(1-F){-Logexp(1-F)}^(η-1))/η好像和原文不符 再及 执行>>hold on即可在同一个土中同时显示当为不同值时 的曲线 另外ode23及你的xiqing.m还需变动一下 [T,Y] = ODE23(ODEFUN,TSPAN,Y0,OPTIONS,P1,P2,...) passes the additional parameters P1,P2,... to the ODE function as ODEFUN(T,Y,P1,P2...), and to all functions specified in OPTIONS. Use OPTIONS = [] as a place holder if no options are set.
Ln对应log命令 ----------------------------------- function DF=qwe(t,F,n) Ka=0.0028; Q=47595; p=1; peq=2.5; T=250; R=8.314; yit=n/2; g= yit*(1-F)*(-log(1-F))^((yit-1)/yit) ; DF=Ka*exp(-Q/R/T)*g*(p-peq); ----------------------------------------- 对于n=2 >>ode23(@qwe,[0 500],0,[],2)