| 查看: 763 | 回复: 7 | |||
[求助]
迭代法求助 已有2人参与
|
| 请教大神一个迭代计算问题。由于对函数求偏导数不方便，所以对于列向量{x}采用固定点迭代法求解，{Xk+1}={f(Xk)}，收敛控制的条件是将列向量{Xk+1}和{Xk}作差，找到其中元素的最大值，判断此最大值是否满足误差要求，满足则收敛，不满足则把列向量{Xk+1}代入{Xk}重新计算。在计算过程中，发现无论初值取得差别有多大，比如10和10000000，计算速度都是一样的，迭代3次就算得同样的结果，而且结果也不太正确。想请教一下各位大神，这究竟是什么情况，方法有问题吗？ |
回复此楼» 猜你喜欢
Materials Today Chemistry审稿周期
已经有4人回复
-
萌生出自己或许不适合搞科研的想法，现在跑or等等看？
已经有3人回复
-
参与限项
已经有3人回复
-
假如你的研究生提出不合理要求
已经有7人回复
-
实验室接单子
已经有4人回复
-
全日制（定向）博士
已经有4人回复
-
对氯苯硼酸纯化
已经有3人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有12人回复
-
所感
已经有4人回复
-
要不要辞职读博？
已经有7人回复
pippi6
铁杆木虫 (著名写手)
工程和科学数值计算咨询
- 数学EPI: 6
- 应助: 413 (硕士)
- 贵宾: 0.002
- 金币: 7116.5
- 散金: 15
- 红花: 63
- 帖子: 1639
- 在线: 798.9小时
- 虫号: 2469437
- 注册: 2013-05-14
- 专业: 计算数学与科学工程计算
2楼2015-05-22 10:51:02
3楼2015-05-22 15:04:10
4楼2015-05-22 16:40:09
feixiaolin
荣誉版主 (文坛精英)
-
专家经验: +518 - 应助: 942 (博后)
- 贵宾: 1.275
- 金币: 3430
- 散金: 58785
- 红花: 532
- 沙发: 11
- 帖子: 24215
- 在线: 2601.8小时
- 虫号: 2139575
- 注册: 2012-11-21
- 专业: 光学信息获取与处理
- 管辖: 数学
5楼2015-05-22 17:43:10
6楼2015-05-23 00:14:04
pippi6
铁杆木虫 (著名写手)
工程和科学数值计算咨询
- 数学EPI: 6
- 应助: 413 (硕士)
- 贵宾: 0.002
- 金币: 7116.5
- 散金: 15
- 红花: 63
- 帖子: 1639
- 在线: 798.9小时
- 虫号: 2469437
- 注册: 2013-05-14
- 专业: 计算数学与科学工程计算
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
kimileegdut: 金币+30, ★有帮助 2015-05-23 14:22:54
kimileegdut: 金币+30, ★有帮助 2015-05-23 14:22:54
|
嗯。如果是方程组，直接使用quasi-Newton确实有点问题。我做过2个变量的。假定方程是 f(x,y)=0 g(x,y)=0 你需要计算 f_x,f_y,g_x,g_y，用以下方式做 f_x = (f(x(k),y(k))-f(x(k-1),y(k)))/(x(k)-x(k-1)) f_y = (f(x(k-1),y(k))-f(x(k-1),y(k-1)))/(y(k)-y(k-1)) g_x = (g(x(k),y(k))-g(x(k-1),y(k)))/(x(k)-x(k-1)) g_y = (g(x(k-1),y(k))-g(x(k-1),y(k-1)))/(y(k)-y(k-1)) 比起通常的迭代，需要多算 f(x(k-1),y(k)) 和 g(x(k-1),y(k)) 两个值，成本稍微高些。更多的变量可以参照此法推广。 残差就是 f(x(k),y(k)) 和 g(x(k),y(k)) ，即方程满足的程度 |
7楼2015-05-23 11:41:06
8楼2015-05-23 14:23:28