[求助]解常微分方程有赏！ - 编程论坛

#2

hitlgq2007-04-12 21:31

怎么没有人恢复，是不是
这个题太难了？

#3

hitzhang2007-04-13 14:27

F =
 
1-exp(-Ka*exp(-Q/R/T)*(p-peq)*t)

#4

hitlgq2007-04-13 15:03

三楼的高手能不能把程序过程写详细一点。因为我还好在在电脑上运行分析。谢谢！求你了！没有办法，主要是我的专业不擅长这方面的编程！

#5

hitzhang2007-04-13 15:09

>> 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)
 
>>

#6

hitlgq2007-04-13 15:51

呵呵，我比较苯啊，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都是常数，但是要在同一个土中同时显示当为不同值时的曲线！！谢谢，再麻烦您一下！！

#7

hitzhang2007-04-13 16:12

η=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.

#8

冰天雪2007-04-13 16:25

高数微分想起就头大

#9

hitlgq2007-04-13 16:37

其实分别n=1和 n=3，以上方程是当n=1的时候。g= η(1-F){-Logexp(1-F)}^（η-1))/η 就是原文中的f(η,F)，可以把这个式子提出来吗？还有就是Ln（1-F）可不可以写成Logexp(1-F)？另外：[T,Y] = ODE23(ODEFUN,TSPAN,Y0,OPTIONS,P1,P2,...) 针对于我这个题的参数怎样变动？有劳您了，谢谢！

#10

hitzhang2007-04-13 17:04

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)

#11

hitlgq2007-04-13 19:03

回复：（hitzhang）Ln对应log命令---------------...

谢谢您的讲解，不过和理论有很大的差别（我记得Ln的底数好象是e，不过我带进去又不能运行了），可能是我给你的方程推导不够明确，这样吧，把您的邮箱给我，我把整个过程都给您，不过这个文章您不能向别人透露，因为这关系到专利问题。呵呵，不管怎样，我还是相信您的！

#12

hitzhang2007-04-15 11:38

尽我所能吧

hitzhangyanxiang@sohu.com

#13

hitlgq2007-04-17 09:26

回复：（hitzhang）尽我所能吧[em25]hitzhangyanxiang...

我已经把附件发过去了！