标题:[求助]解常微分方程有赏!
只看楼主
hitlgq
Rank: 1
等 级:新手上路
帖 子:43
专家分:0
注 册:2007-4-6
 问题点数:0 回复次数:12 
[求助]解常微分方程有赏!
谁能帮我解一下这个常微分方程,感激不尽!1
ch3FbEnl.rar (15.83 KB) [求助]解常微分方程有赏!


搜索更多相关主题的帖子: 常微分方程 感激不尽 
2007-04-12 18:59
hitlgq
Rank: 1
等 级:新手上路
帖 子:43
专家分:0
注 册:2007-4-6
得分:0 
怎么没有人恢复,是不是
这个题太难了 ?

2007-04-12 21:31
hitzhang
Rank: 20Rank: 20Rank: 20Rank: 20Rank: 20
等 级:版主
威 望:21
帖 子:369
专家分:52
注 册:2006-9-24
得分:0 
F =
 
1-exp(-Ka*exp(-Q/R/T)*(p-peq)*t)

2007-04-13 14:27
hitlgq
Rank: 1
等 级:新手上路
帖 子:43
专家分:0
注 册:2007-4-6
得分:0 
三楼的高手 能不能把程序过程写详细一点。因为我还好在在电脑上运行分析。谢谢!求你了!没有办法,主要是我的专业不擅长这方面的编程!

2007-04-13 15:03
hitzhang
Rank: 20Rank: 20Rank: 20Rank: 20Rank: 20
等 级:版主
威 望:21
帖 子:369
专家分:52
注 册:2006-9-24
得分:0 
>> 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)
 
>> 

2007-04-13 15:09
hitlgq
Rank: 1
等 级:新手上路
帖 子:43
专家分:0
注 册:2007-4-6
得分:0 

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


2007-04-13 15:51
hitzhang
Rank: 20Rank: 20Rank: 20Rank: 20Rank: 20
等 级:版主
威 望:21
帖 子:369
专家分:52
注 册:2006-9-24
得分:0 
η=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.

2007-04-13 16:12
冰天雪
Rank: 1
等 级:新手上路
威 望:1
帖 子:331
专家分:0
注 册:2007-2-27
得分:0 
高数微分想起就头大
2007-04-13 16:25
hitlgq
Rank: 1
等 级:新手上路
帖 子:43
专家分:0
注 册:2007-4-6
得分:0 
其实分别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,...) 针对于我这个题的参数怎样变动?  有劳您了,谢谢!

2007-04-13 16:37
hitzhang
Rank: 20Rank: 20Rank: 20Rank: 20Rank: 20
等 级:版主
威 望:21
帖 子:369
专家分:52
注 册:2006-9-24
得分:0 
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)

2007-04-13 17:04



参与讨论请移步原网站贴子:https://bbs.bccn.net/thread-131392-1-1.html




关于我们 | 广告合作 | 编程中国 | 清除Cookies | TOP | 手机版

编程中国 版权所有,并保留所有权利。
Powered by Discuz, Processed in 0.132530 second(s), 8 queries.
Copyright©2004-2024, BCCN.NET, All Rights Reserved