| 查看: 967 | 回复: 7 | |||
[求助]
迭代法求助 已有2人参与
|
| 请教大神一个迭代计算问题。由于对函数求偏导数不方便，所以对于列向量{x}采用固定点迭代法求解，{Xk+1}={f(Xk)}，收敛控制的条件是将列向量{Xk+1}和{Xk}作差，找到其中元素的最大值，判断此最大值是否满足误差要求，满足则收敛，不满足则把列向量{Xk+1}代入{Xk}重新计算。在计算过程中，发现无论初值取得差别有多大，比如10和10000000，计算速度都是一样的，迭代3次就算得同样的结果，而且结果也不太正确。想请教一下各位大神，这究竟是什么情况，方法有问题吗？ |
回复此楼» 猜你喜欢
能否退出参与的面上项目解除限项
已经有15人回复
-
2026国自然函评费到账
已经有18人回复
-
建议基金发布提前给出明确的时间点
已经有15人回复
-
范进中举一文的中心思想
已经有7人回复
-
朋友圈看到的
已经有9人回复
-
让我中一个面上吧！
已经有16人回复
-
放榜前的不淡定
已经有17人回复
-
科研孤儿太难了
已经有19人回复
-
明天应该可查了！？
已经有4人回复
-
什么时候开奖？
已经有11人回复
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
- 金币: 3881
- 散金: 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
20