24小时热门版块排行榜    

Znn3bq.jpeg
查看: 2611  |  回复: 12

hongxiaodouy

新虫 (初入文坛)

[求助] 用不动点迭代求非线性方程组,求得的不是想要的解,怎么办

我用不动点迭代法求解一个非线性方程组(在这个方程组里我会用到其他的子程序),求得的解不满足我的实际条件,里面有的是负数。例如我给的初始值x0=5,y0=0.1已经很接近我需要得到的解x*=5.1,y*=0.11可是迭代以后却远离这个解,得到x=2,y=0.03之类的,请问这个问题是我程序语言的问题还是计算方法的错误?
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

xiegangmai

版主 (职业作家)

我没头衔

优秀版主优秀版主优秀版主

【答案】应助回帖


感谢参与,应助指数 +1
臭水沟: 金币+1, 谢谢交流~~ 2012-05-31 19:15:58
你是自己写的迭代算法吗?建议仔细检查一下你的程序。不满足收敛条件,不会退出的。
另外你可以用matlab自带的函数试试,对比一下。
明德厚学、求是创新
3楼2012-05-30 22:57:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

guardian_qq

银虫 (小有名气)

【答案】应助回帖

★ ★
xiegangmai: 金币+2, 谢谢参与 2012-10-10 21:01:46
建议你用fsolve之类的求解,通常自带的都是优化过的,自己编的程序往往在考虑累积误差方面很粗糙,不容易发现问题。
如果要分析自己的程序也请你最好贴出代码来给大家看一下。
10楼2012-10-08 15:30:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

czxf

银虫 (正式写手)

【答案】应助回帖


感谢参与,应助指数 +1
xiegangmai: 金币+1, 谢谢参与! 2012-05-30 22:55:57
建议楼主看下不动点迭代收敛的条件,收敛区间应该是要求倒数绝对值小于1的吧,所以迭代初始点的选取还是有讲究的
2楼2012-05-30 12:57:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hongxiaodouy

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by czxf at 2012-05-30 12:57:04
建议楼主看下不动点迭代收敛的条件,收敛区间应该是要求倒数绝对值小于1的吧,所以迭代初始点的选取还是有讲究的

看了收敛条件 可是没有看懂呢   其实就像我上面的举例,如果我直接带入解x0=5.1,y0=0.11,那不就不用迭代,直接满足精度条件,得到解x*=5.1,y*=0.11,但是只要初始值偏离x*,y*,迭代结果就是x=2,y=0.03,不知道是怎么回事
4楼2012-05-31 13:27:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hongxiaodouy

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by czxf at 2012-05-30 12:57:04
建议楼主看下不动点迭代收敛的条件,收敛区间应该是要求倒数绝对值小于1的吧,所以迭代初始点的选取还是有讲究的

看了收敛条件 可是没有看懂呢   其实就像我上面的举例,如果我直接带入解x0=5.1,y0=0.11,那不就不用迭代,直接满足精度条件,得到解x*=5.1,y*=0.11,但是只要初始值偏离x*,y*,迭代结果就是x=2,y=0.03,不知道是怎么回事
5楼2012-05-31 13:27:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hongxiaodouy

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by czxf at 2012-05-30 12:57:04
建议楼主看下不动点迭代收敛的条件,收敛区间应该是要求倒数绝对值小于1的吧,所以迭代初始点的选取还是有讲究的

看了收敛条件 可是没有看懂呢   其实就像我上面的举例,如果我直接带入解x0=5.1,y0=0.11,那不就不用迭代,直接满足精度条件,得到解x*=5.1,y*=0.11,但是只要初始值偏离x*,y*,迭代结果就是x=2,y=0.03,不知道是怎么回事
6楼2012-05-31 13:28:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hongxiaodouy

新虫 (初入文坛)

引用回帖:
3楼: Originally posted by xiegangmai at 2012-05-30 22:57:37
你是自己写的迭代算法吗?建议仔细检查一下你的程序。不满足收敛条件,不会退出的。
另外你可以用matlab自带的函数试试,对比一下。

是我自己写的,自带的函数需要定义自变量什么的,比较麻烦,我就没有用
7楼2012-05-31 13:31:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hongxiaodouy

新虫 (初入文坛)

引用回帖:
3楼: Originally posted by xiegangmai at 2012-05-30 22:57:37
你是自己写的迭代算法吗?建议仔细检查一下你的程序。不满足收敛条件,不会退出的。
另外你可以用matlab自带的函数试试,对比一下。

电脑有点慢,就点了好几次,结果就回复了好几次~删不去~嘻~
8楼2012-05-31 13:33:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mataohk

银虫 (初入文坛)

【答案】应助回帖


xiegangmai: 金币+1, 谢谢参与 2012-10-05 20:45:53
估计是你自己程序的问题。其实一般问题都可以用Matlab自带的函数求解,比如:fsolve,fzero
9楼2012-09-29 15:25:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 hongxiaodouy 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料复试求调剂 +16 xhhdjdjsjks 2026-04-09 16/800 2026-04-09 20:56 by bljnqdcc
[考研] 296求调剂 +11 汪!?! 2026-04-08 11/550 2026-04-09 19:30 by 小小树2024
[考研] 一志愿华工085600 331分 +6 天下ww 2026-04-09 6/300 2026-04-09 18:59 by l_paradox
[考研] 085600材料与化工专硕329 求调剂 +24 额cc 2026-04-06 25/1250 2026-04-09 16:01 by wp06
[考研] 299求调剂 +6 ZVVZ13 2026-04-08 6/300 2026-04-09 15:48 by wp06
[考研] 材料专硕322 +14 哈哈哈吼吼吼哈 2026-04-05 14/700 2026-04-09 13:25 by 5268321
[考研] 288求调剂 +15 没有答案_ 2026-04-05 15/750 2026-04-09 10:22 by 5268321
[考研] 085404,334分,求调剂 +5 sunjie8888 2026-04-08 8/400 2026-04-09 07:26 by sunjie8888
[考研] 320分人工智能调剂 +9 振—TZ 2026-04-03 10/500 2026-04-08 19:56 by 振—TZ
[考研] 274求调剂求调剂 +10 Jachenbingoo 2026-04-06 13/650 2026-04-08 14:25 by zhq0425
[考研] 一志愿生物与医药,296分,求调剂 +14 66鹿 2026-04-03 16/800 2026-04-08 10:38 by tjzhao
[考研] 求考研材料调剂 +3 材化李可 2026-04-07 3/150 2026-04-08 00:21 by JourneyLucky
[考研] 计算机408|在校多次国家级竞赛获奖|申请调剂 +4 东山大白鹅 2026-04-05 4/200 2026-04-08 00:18 by chongya
[考研] 331求调剂 +5 张元一 2026-04-07 6/300 2026-04-07 22:13 by hemengdong
[考研] 307求调剂 +3 Youth@@ 2026-04-07 3/150 2026-04-07 09:25 by 小黑不怕难
[考研] 一志愿武汉理工大学080200机械工程308分,求调剂 +4 终不似从前 2026-04-05 4/200 2026-04-06 11:46 by 考研学校招点人
[考研] 求调剂 +7 张.1 2026-04-05 7/350 2026-04-05 20:40 by 啵啵啵0119
[考研] 311分 22408 求调剂 +3 bing_bot 2026-04-03 3/150 2026-04-05 00:43 by chongya
[考研] 化工求调剂 +11 荔香芝士椰奶 2026-04-03 11/550 2026-04-03 22:06 by 啵啵啵0119
[考研] 285求调剂 +5 AZMK 2026-04-03 8/400 2026-04-03 18:17 by AZMK
信息提示
请填处理意见