d****o 发帖数: 32610 | 1 【 以下文字转载自 Military 讨论区 】
发信人: Xujiahui (南模是所好学校), 信区: Military
标 题: Re: 法国科学家
发信站: BBS 未名空间站 (Sun Apr 24 00:25:31 2016, 美东)
谢谢,我没干过数值计算的活,这次是赶鸭子上架。我大学是在法国读的,所以有些中
文术语不清楚,牛顿迭代是不是就是newton raphson?收敛应该是converge吧?
我的问题有点复杂,前面没有解释清楚。关键是逼近值一定要准,因为这只是一个单元
里的逼近值,然后输出到下一个单元作为输入。一共有几千个单元,所以如果不准的话
最后的误差会被放大。我最开始用的是brute forcing,一步一步向交点逼近,计算停止
在误差小于0.1%(举例),但是如果我的step取的不够小的话,可能最靠近交点的逼近
值都不能在0.1%以内,无法converge。但是step取的小的话,计算量又很大。
bisection我也用过,也不准尤其是在多个单元之后,误差会越来越大。
我还试过shooting method,其实差不多就是newton raphson了,两边取tangent line向
中间逼近,但是他的缺点是逼近点在不同的单元里有时在左侧,有时在右侧,不过倒是
有可能互相取消。但是横坐标的数值不是整数,给我后面的运算带来很多麻烦。 |
t****t 发帖数: 6806 | 2 这个水平果然写偏微方程是有难度的.
【在 d****o 的大作中提到】 : 【 以下文字转载自 Military 讨论区 】 : 发信人: Xujiahui (南模是所好学校), 信区: Military : 标 题: Re: 法国科学家 : 发信站: BBS 未名空间站 (Sun Apr 24 00:25:31 2016, 美东) : 谢谢,我没干过数值计算的活,这次是赶鸭子上架。我大学是在法国读的,所以有些中 : 文术语不清楚,牛顿迭代是不是就是newton raphson?收敛应该是converge吧? : 我的问题有点复杂,前面没有解释清楚。关键是逼近值一定要准,因为这只是一个单元 : 里的逼近值,然后输出到下一个单元作为输入。一共有几千个单元,所以如果不准的话 : 最后的误差会被放大。我最开始用的是brute forcing,一步一步向交点逼近,计算停止 : 在误差小于0.1%(举例),但是如果我的step取的不够小的话,可能最靠近交点的逼近
|
d****o 发帖数: 32610 | 3 按拉普拉斯算子的形状重新发明了个轮子
[在 thrust (祝阳阳早日康复) 的大作中提到:]
:这个水平果然写偏微方程是有难度的.
:
:........... |
d********f 发帖数: 43471 | 4 md,估计又是做纳米太阳能电池的
【在 d****o 的大作中提到】 : 按拉普拉斯算子的形状重新发明了个轮子 : [在 thrust (祝阳阳早日康复) 的大作中提到:] : :这个水平果然写偏微方程是有难度的. : : : :...........
|
a***e 发帖数: 27968 | 5 这个横坐标不是整数是咋回事?
★ 发自iPhone App: ChineseWeb 1.0.6
【在 d****o 的大作中提到】 : 按拉普拉斯算子的形状重新发明了个轮子 : [在 thrust (祝阳阳早日康复) 的大作中提到:] : :这个水平果然写偏微方程是有难度的. : : : :...........
|