24小时热门版块排行榜    

查看: 2590  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 274求调剂 +7 S.H1 2026-03-18 7/350 2026-03-20 11:11 by Delta2012
[考研] 317求调剂 +4 申子申申 2026-03-19 9/450 2026-03-20 11:08 by 申子申申
[考研] 271材料工程求调剂 +7 .6lL 2026-03-18 7/350 2026-03-20 09:10 by xingguangj
[考博] 申博26年 +3 八6八68 2026-03-19 3/150 2026-03-19 19:43 by nxgogo
[考研] 321求调剂 +8 何润采123 2026-03-18 10/500 2026-03-19 16:46 by 何润采123
[考研] 085601材料工程专硕求调剂 +10 慕寒mio 2026-03-16 10/500 2026-03-19 15:26 by 丁丁*
[考研] 266求调剂 +5 阳阳哇塞 2026-03-14 10/500 2026-03-19 15:08 by 阳阳哇塞
[考研] 化学求调剂 +3 临泽境llllll 2026-03-17 4/200 2026-03-19 13:59 by houyaoxu
[考研] 085600材料与化工求调剂 +6 绪幸与子 2026-03-17 6/300 2026-03-19 13:27 by houyaoxu
[考研] 一志愿中海洋材料工程专硕330分求调剂 +7 小材化本科 2026-03-18 7/350 2026-03-19 10:46 by Linda Hu
[考研] 材料工程专硕调剂 +5 204818@lcx 2026-03-17 6/300 2026-03-18 22:55 by 204818@lcx
[考研] 085601专硕,总分342求调剂,地区不限 +5 share_joy 2026-03-16 5/250 2026-03-18 14:48 by haxia
[考研] 0854可跨调剂,一作一项核心论文五项专利,省、国级证书40+数一英一287 +8 小李0854 2026-03-16 8/400 2026-03-18 14:35 by 搏击518
[考研] 280求调剂 +6 咕噜晓晓 2026-03-18 7/350 2026-03-18 11:25 by 无际的草原
[考研] 085601求调剂 +4 Du.11 2026-03-16 4/200 2026-03-17 17:08 by ruiyingmiao
[考研] 290求调剂 +3 p asserby. 2026-03-15 4/200 2026-03-17 16:35 by wangkm
[考研] 东南大学364求调剂 +5 JasonYuiui 2026-03-15 5/250 2026-03-16 21:28 by 木瓜膏
[考研] 一志愿211 0703方向310分求调剂 +3 努力奋斗112 2026-03-15 3/150 2026-03-16 16:44 by houyaoxu
[考研] 321求调剂 +5 大米饭! 2026-03-15 5/250 2026-03-16 16:33 by houyaoxu
[考研] 本科南京大学一志愿川大药学327 +3 麦田耕者 2026-03-14 3/150 2026-03-14 20:04 by 外星文明
信息提示
请填处理意见