| 查看: 826 | 回复: 7 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
迭代法求助 已有2人参与
|
|||
| 请教大神一个迭代计算问题。由于对函数求偏导数不方便，所以对于列向量{x}采用固定点迭代法求解，{Xk+1}={f(Xk)}，收敛控制的条件是将列向量{Xk+1}和{Xk}作差，找到其中元素的最大值，判断此最大值是否满足误差要求，满足则收敛，不满足则把列向量{Xk+1}代入{Xk}重新计算。在计算过程中，发现无论初值取得差别有多大，比如10和10000000，计算速度都是一样的，迭代3次就算得同样的结果，而且结果也不太正确。想请教一下各位大神，这究竟是什么情况，方法有问题吗？ |
回复此楼» 猜你喜欢
305分求调剂（食品工程）
已经有7人回复
-
材料与化工考研调剂
已经有5人回复
-
299求调剂
已经有6人回复
-
材料292调剂
已经有5人回复
-
085404电子信息284分求调剂
已经有4人回复
-
收08调剂生
已经有3人回复
-
081700 调剂 267分
已经有9人回复
-
293求调剂
已经有5人回复
-
086003食品工程求调剂
已经有4人回复
-
0856材料专硕353求调剂
已经有5人回复
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
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
