注册 登录
编程论坛 Matlab

[求助]解常微分方程有赏!

hitlgq 发布于 2007-04-12 18:59, 1350 次点击
谁能帮我解一下这个常微分方程,感激不尽!1
只有本站会员才能查看附件,请 登录

12 回复
#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
#13
hitlgq2007-04-17 09:26
回复:(hitzhang)尽我所能吧[em25]hitzhangyanxiang...

我已经把附件发过去了!

1