如何用变步长的差分法解PDE？ - 未名空间MITBBS历史存档

国际科技财经博客移民网络热点娱乐民生时事公众号

Redian新闻

>未名空间

>Computation - 科学计算

如何用变步长的差分法解PDE？

如何用变步长的差分法解PDE？# Computation - 科学计算

b*h2004-03-03 08:03

1 楼

关键是如何确定步长。
没有学过计算数学，
还请大虾指点一下有什么网上资源或者书比较好。

h*o2004-03-03 08:03

2 楼

步长一般都是adaptive的吧。平滑的地方步长大，歪歪扭扭的地方步长小

【在 b******h 的大作中提到】

: 关键是如何确定步长。
: 没有学过计算数学，
: 还请大虾指点一下有什么网上资源或者书比较好。

b*h2004-03-03 08:03

3 楼

But how to do it?

【在 h***o 的大作中提到】

: 步长一般都是adaptive的吧。平滑的地方步长大，歪歪扭扭的地方步长小

a*s2004-03-03 08:03

4 楼

你怎么知道什么地方平滑，什么地方不平滑？如果在求解之前不知道的话，怎么知道呢？

【在 h***o 的大作中提到】

: 步长一般都是adaptive的吧。平滑的地方步长大，歪歪扭扭的地方步长小

h*o2004-03-03 08:03

5 楼

看看NR是怎么make adaptive step in Runge-Kutta method 的

【在 a******s 的大作中提到】

: 你怎么知道什么地方平滑，什么地方不平滑？如果在求解之前不知道的话，怎么知道呢？

h*o2004-03-03 08:03

6 楼

找个PDE solver..它一般都会自动take care这些咚咚

【在 b******h 的大作中提到】

: But how to do it?

b*h2004-03-03 08:03

7 楼

麻烦推荐一个，matlab据说有，但是很慢。

【在 h***o 的大作中提到】

: 找个PDE solver..它一般都会自动take care这些咚咚

h*o2004-03-03 08:03

8 楼

netlib.org里肯定有一大堆吧，自己慢慢刨

【在 b******h 的大作中提到】

: 麻烦推荐一个，matlab据说有，但是很慢。

g*e2004-03-03 08:03

9 楼

难道是解纯数学PDE?
如果是个数学模型,该知道解的大致情况,就是哪些地方平滑用长步长.

【在 a******s 的大作中提到】

: 你怎么知道什么地方平滑，什么地方不平滑？如果在求解之前不知道的话，怎么知道呢？

c*n2004-03-03 08:03

10 楼

In matlab code, for example RK45 is an adaptive solver.
The mechanisim is that every time, use RK4 get a solution, use RK5 to get
another, check if abs(RK4-RK5)> error1 that you set. if it is, cut interval
into half, if maintain the interval. Actually there are plenty of adaptive solvers, you
don't
need to program by yourself.

【在 b******h 的大作中提到】

: 关键是如何确定步长。
: 没有学过计算数学，
: 还请大虾指点一下有什么网上资源或者书比较好。