24小时热门版块排行榜    

查看: 889  |  回复: 6

关月阑珊

木虫 (小有名气)

[求助] 数值分析 迭代算法 急! 已有2人参与

见图 22.jpg
求m;谢谢!
楼主疑问:让迭代函数的导数值为0就是加速迭代了吧?

数值分析 迭代算法  急!
22.jpg

[ Last edited by feixiaolin on 2014-3-21 at 22:29 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

关月阑珊

木虫 (小有名气)

2楼2014-03-21 21:56:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

feixiaolin

荣誉版主 (文坛精英)

优秀版主

取f(x)=1+x-sin(x)
证明对后者|f(x_k+1)|<= 前者的 |f(x_k+1)|,即可; 从而获得m的范围。
3楼2014-03-21 22:39:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mathstudy

金虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
关月阑珊: 金币+35, ★★★★★最佳答案 2014-03-23 21:16:36
关月阑珊: 金币+5, ★★★★★最佳答案 2014-03-23 21:17:35
引用回帖:
3楼: Originally posted by feixiaolin at 2014-03-21 22:39:08
取f(x)=1+x-sin(x)
证明对后者|f(x_k+1)|<= 前者的 |f(x_k+1)|,即可; 从而获得m的范围。

1.一般情况下 可以按照 2楼说的做,但是针对这个题目应该不太适合。
2.但是 楼主的问题是不是有误,因为在第二个格式x_(k+1)=(m x_k+1-sin x_k)/(1+m)对不动点不成立,所以我怀疑楼主漏掉或者输入错了  (应该是sin x_k -1吧)

3.可以试试这样,设不动点是x*, 两中方法均能收敛,则
lim(k->无穷) (x_k_1-x*)/(x_k_2-x*)=0 一定说明x_k_1代表的迭代法收敛速度比x_k_2对应的收敛方法快(这只是一个充分条件)。
4. 如果将楼主给出的 加速迭代中 (应该是sin x_k-1)   则3中分析导出
lim (m+cos x_k)/cos x_k/(m+1)=0 进一步导出 m=-cos x*
4楼2014-03-22 00:43:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wurongjun

专家顾问 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
关月阑珊: 金币+10, ★★★很有帮助 2014-03-23 21:16:15
这是一个典型的加权提速法!
范围:-1<m<0
首先容易证明:解的大概范围为[-pi,-pi/2];
然后在这一区间里面讨论即可!
显然原迭代函数为:h(x)=sin(x)-1
修改后的代函数为:g(x)=(mx-sin(x))/(1+m)
要使得迭代加速,其实就是选择合适的m
使得|g'(x)|<|h'(x)|=-cos(x)
再由这个条件就可推出:-1<m<0
不过讨论复杂一点而已,就留给楼主吧!
善恶到头终有报,人间正道是沧桑.
5楼2014-03-22 08:23:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wurongjun

专家顾问 (职业作家)

【答案】应助回帖

3楼的2说得很对啊!要保证是不动点才行!
刚才忽略了这一点!
修改后再用刚才的办法,那么结果就不是那样了!
善恶到头终有报,人间正道是沧桑.
6楼2014-03-22 08:34:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

关月阑珊

木虫 (小有名气)

送红花一朵
非常感谢大家的热心帮助!
7楼2014-03-23 21:18:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 关月阑珊 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 s3fFTmArrBt6 2026-09-14 4/200 2026-09-14 17:25 by Aj1rhIDL5ixY
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 LwdutQ8HoqWP 2026-09-13 5/250 2026-09-14 16:25 by Aj1rhIDL5ixY
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 LwdutQ8HoqWP 2026-09-13 4/200 2026-09-14 16:13 by Aj1rhIDL5ixY
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 LwdutQ8HoqWP 2026-09-13 4/200 2026-09-14 16:12 by Aj1rhIDL5ixY
[找工作] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 3n8v2C8RimXI 2026-09-13 4/200 2026-09-14 16:00 by Aj1rhIDL5ixY
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 3n8v2C8RimXI 2026-09-13 4/200 2026-09-14 15:48 by Aj1rhIDL5ixY
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 3n8v2C8RimXI 2026-09-13 4/200 2026-09-14 15:45 by Aj1rhIDL5ixY
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 3n8v2C8RimXI 2026-09-13 4/200 2026-09-14 15:37 by Aj1rhIDL5ixY
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 QUjhNVAcOSff 2026-09-13 6/300 2026-09-14 15:21 by Aj1rhIDL5ixY
[考博] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 QUjhNVAcOSff 2026-09-13 4/200 2026-09-14 15:13 by Aj1rhIDL5ixY
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 QUjhNVAcOSff 2026-09-13 3/150 2026-09-14 15:09 by Aj1rhIDL5ixY
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 QUjhNVAcOSff 2026-09-13 4/200 2026-09-14 15:01 by Aj1rhIDL5ixY
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 QUjhNVAcOSff 2026-09-13 3/150 2026-09-14 14:57 by Aj1rhIDL5ixY
[基金申请] 求教各位大神:2026教育部人文社科青年基金项目何时公示呀? +5 云宴山人 2026-09-10 7/350 2026-09-14 12:11 by ecnu2013
[找工作] 浙江师范大学是怎么坑我的 +6 recruit123 2026-09-09 6/300 2026-09-14 12:00 by 萨仁其其格
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 QUjhNVAcOSff 2026-09-13 3/150 2026-09-14 10:11 by 23jxep3nCNZb
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 QUjhNVAcOSff 2026-09-13 4/200 2026-09-14 10:01 by 23jxep3nCNZb
[考博] 西北工业大学材料学院电化学传感与催化课题组招收2027年推免研究生、博士生、博士后 +3 马方园 2026-09-10 4/200 2026-09-14 09:29 by Chem张zz
[基金申请] 面上有专家说实验设备不是我们单位的 +6 zj_muchong 2026-09-08 8/400 2026-09-11 06:56 by suifengdao
[有机交流] 科研求助 10+4 我吃一只鱼 2026-09-07 4/200 2026-09-10 08:32 by 0940628
信息提示
请填处理意见