l*1
2 楼
在matlab中求解一个常微分方程组(ODEs),使用ode45时,求解速度非常的慢。
MATLAB帮助中提示如果ode45的求解速度非常慢,说明该微分方程组是stiff的,可以尝
试ode15s。
我的问题是:是什么造成了方程组出现stiff的情况? 如何避免方程组出现stiff的情
况?
请高手指点,不胜感激!
MATLAB帮助中提示如果ode45的求解速度非常慢,说明该微分方程组是stiff的,可以尝
试ode15s。
我的问题是:是什么造成了方程组出现stiff的情况? 如何避免方程组出现stiff的情
况?
请高手指点,不胜感激!
f*l
3 楼
那里可以看见用机会成本理由来分析案例的。
h*c
5 楼
singular, the rank of Jacobian is not n;
use coninuation.
AUTO
http://indy.cs.concordia.ca/auto/
it is written by my boss,
you can contact me via PM.
use coninuation.
AUTO
http://indy.cs.concordia.ca/auto/
it is written by my boss,
you can contact me via PM.
h*w
6 楼
key word?
l*1
7 楼
多谢!
AUTO可能很好,只是它是用Fortran写的,而我的那个微分方程组(一个动态系统)是
用MATLAB写的。这么看来,我暂时还无法借用AUTO了。
AUTO可能很好,只是它是用Fortran写的,而我的那个微分方程组(一个动态系统)是
用MATLAB写的。这么看来,我暂时还无法借用AUTO了。
b*e
9 楼
清华徐士良的算法书里有一节是关于刚性方程的,解法叫吉尔(GILL)方法。
h*c
11 楼
correction:
stiff problem:
x'=\lambda x, when Re(\lambda)-> -\infty
x could be extended to n-vector, \lambda thus could be a matrix D, thus one
of the eigenvalues of D has its real part approaches negative infinity.
This is a totally different thing from singular Jaconbian of ODEs.
Stiff problem can be attacked by stability analysis. Normally a multi-step
implicit mesh scheme should help.
Just for discussion.
stiff problem:
x'=\lambda x, when Re(\lambda)-> -\infty
x could be extended to n-vector, \lambda thus could be a matrix D, thus one
of the eigenvalues of D has its real part approaches negative infinity.
This is a totally different thing from singular Jaconbian of ODEs.
Stiff problem can be attacked by stability analysis. Normally a multi-step
implicit mesh scheme should help.
Just for discussion.
h*w
12 楼
讲什么,大家没人买我帖的东西吗?奇怪,头一遭,看来稀有的不受欢迎.
j*p
13 楼
dude, good to see that you figured it out at last :-)
one
【在 h**********c 的大作中提到】![](/moin_static193/solenoid/img/up.png)
: correction:
: stiff problem:
: x'=\lambda x, when Re(\lambda)-> -\infty
: x could be extended to n-vector, \lambda thus could be a matrix D, thus one
: of the eigenvalues of D has its real part approaches negative infinity.
: This is a totally different thing from singular Jaconbian of ODEs.
: Stiff problem can be attacked by stability analysis. Normally a multi-step
: implicit mesh scheme should help.
: Just for discussion.
one
【在 h**********c 的大作中提到】
![](/moin_static193/solenoid/img/up.png)
: correction:
: stiff problem:
: x'=\lambda x, when Re(\lambda)-> -\infty
: x could be extended to n-vector, \lambda thus could be a matrix D, thus one
: of the eigenvalues of D has its real part approaches negative infinity.
: This is a totally different thing from singular Jaconbian of ODEs.
: Stiff problem can be attacked by stability analysis. Normally a multi-step
: implicit mesh scheme should help.
: Just for discussion.
l*1
17 楼
What is stiffness?
A stiff ODE is an ordinary differential equation that has a transient region
whose behavior is on a different scale from that outside this transient
region. A physical example of a stiff system involves chemical reaction
rates, where typically the convergence to a final solution can be quite
rapid.
A stiff ODE is an ordinary differential equation that has a transient region
whose behavior is on a different scale from that outside this transient
region. A physical example of a stiff system involves chemical reaction
rates, where typically the convergence to a final solution can be quite
rapid.
h*c
19 楼
I believe this is rather an explanation ; not a concrete definition or the
qualification.
I personally believe in the "transient region", during integration, at some
step the system A X= f , A is very ill-conditioned.
I did take a course covering this problem two years ago, by my supervisor.
And I happened to audit that course again this term.
The definition is as per my previous post.
region
【在 l********1 的大作中提到】![](/moin_static193/solenoid/img/up.png)
: What is stiffness?
: A stiff ODE is an ordinary differential equation that has a transient region
: whose behavior is on a different scale from that outside this transient
: region. A physical example of a stiff system involves chemical reaction
: rates, where typically the convergence to a final solution can be quite
: rapid.
qualification.
I personally believe in the "transient region", during integration, at some
step the system A X= f , A is very ill-conditioned.
I did take a course covering this problem two years ago, by my supervisor.
And I happened to audit that course again this term.
The definition is as per my previous post.
region
【在 l********1 的大作中提到】
![](/moin_static193/solenoid/img/up.png)
: What is stiffness?
: A stiff ODE is an ordinary differential equation that has a transient region
: whose behavior is on a different scale from that outside this transient
: region. A physical example of a stiff system involves chemical reaction
: rates, where typically the convergence to a final solution can be quite
: rapid.
j*1
20 楼
因为水平有限,不敢买稀有的啊,材料价的放心点,亏不到哪儿去.
f*l
21 楼
为什么要关注这个定义呢? 关于stiff的定义可以有很多种. 但关键还是在于
数值稳定性的考虑. 你所说的都是把方程线性化之后的分析. 这些分析只是为了
给人一个直觉的基础罢了. 数学要是只关心定义, 就失去其应用的一面了.
对你的问题, 你直接用ode15s就是了. 一般用matlab来解的都不会太复杂.
ode15s是建立在BDF公式上的变阶变步长算法, 由C. W. Gear早先提出的.
如果问题比较大, 建议使用CVODE, LSODE, DASSL, DASPK等成熟的软件. 严肃的
计算一般是不用matlab搞的.
消除stiffness取决于问题. 如果你知道是什么部分导致一些变量很快地进入
动态平衡, 把这部分用代数的方法取代, 就是了. 这个具体到问题可能很复杂.
the
some
【在 h**********c 的大作中提到】![](/moin_static193/solenoid/img/up.png)
: I believe this is rather an explanation ; not a concrete definition or the
: qualification.
: I personally believe in the "transient region", during integration, at some
: step the system A X= f , A is very ill-conditioned.
: I did take a course covering this problem two years ago, by my supervisor.
: And I happened to audit that course again this term.
: The definition is as per my previous post.
:
: region
数值稳定性的考虑. 你所说的都是把方程线性化之后的分析. 这些分析只是为了
给人一个直觉的基础罢了. 数学要是只关心定义, 就失去其应用的一面了.
对你的问题, 你直接用ode15s就是了. 一般用matlab来解的都不会太复杂.
ode15s是建立在BDF公式上的变阶变步长算法, 由C. W. Gear早先提出的.
如果问题比较大, 建议使用CVODE, LSODE, DASSL, DASPK等成熟的软件. 严肃的
计算一般是不用matlab搞的.
消除stiffness取决于问题. 如果你知道是什么部分导致一些变量很快地进入
动态平衡, 把这部分用代数的方法取代, 就是了. 这个具体到问题可能很复杂.
the
some
【在 h**********c 的大作中提到】
![](/moin_static193/solenoid/img/up.png)
: I believe this is rather an explanation ; not a concrete definition or the
: qualification.
: I personally believe in the "transient region", during integration, at some
: step the system A X= f , A is very ill-conditioned.
: I did take a course covering this problem two years ago, by my supervisor.
: And I happened to audit that course again this term.
: The definition is as per my previous post.
:
: region
l*1
23 楼
to fool,
ode15s我尝试过,速度一样很慢。
现在我用c编写部分程序,求解ode时采用了dassl算法,编译后在matlab下进行调用,
同其他程序部分协调进行,效果不错,速度很快:)
ode15s我尝试过,速度一样很慢。
现在我用c编写部分程序,求解ode时采用了dassl算法,编译后在matlab下进行调用,
同其他程序部分协调进行,效果不错,速度很快:)
相关阅读
转贴:激辩广东模式:广东模式是否该被终结?[合集] 通货膨胀请教计量高手:人民币的valuation问题请教几个有关回复审稿报告的问题~THX[合集] 光华金融硕士炙手可热,08毕业生就业再创佳绩灾区重建的城市规划,有搞房地产金融以及区域经济学的参与吗?上财与图大IDEI达成合作协议[合集] 发paper的时机和tenure的关系?[合集] 近十年来Harvard培养了哪些中国学术强人?搞经济的朋友,这篇文章说的有道理吗?An old Chinese myth[合集] George Mason University怎么样呢[合集] 可能跟LZ的问题不是很相关,只是澄清一下 (转载)Information structure and firm[合集] 请教FINANCE FIELD EXAM[合集] 请pro们分析一下俺的原始问题郎咸平:08年的8大危机[合集] 一篇 JPE 可以值 10万刀?请教panel data: T >> N刚有一篇JFE被接受,自己cong一下