做作业要用迭代计算湿球温度Tw，刚开始的时候把值输错了但是可以出结果，第二次把输错的值改了之后算不出结果，得到的结果一直在两个值之间振荡。我猜可能是因为不收敛，但是前后两次的计算式都完全相同，仅仅改了几个值就会不收敛吗？以下是编的程序和结果及报错的截图，希望能得到各位大神的指导和帮助，小菜鸟感激不尽 T=296.5; L=2445.3*10^3; cp=1005; eps=0.622; p=91; A=2.53*10^8; B=5.42*10^3; w=0.01123; Tw(1)=273; Tw(2)=T-(L/cp)*(eps/p*A*exp(-B/Tw(1))-w); n=1; while abs(Tw(n+1)-Tw(n))>=1e-10     Tw(n+2)=T-(L/cp)*(eps/p*A*exp(-B/Tw(n+1))-w);     n=n+1; end

matlab, 算法

微信图片_20181113111819.png (471.04 KB, 下载次数: 1)

报错截图

微信图片_20181113111758.png (559.54 KB, 下载次数: 0)

结果截图

是的,就跟牛顿法有时候无法收敛一样,这个是迭代算法本身的问题
要么改初值要么换其他迭代法

举个例子,牛顿法求f = -x^4 + 3*x^2 + 2 ==0的根时,取初值x(1) = 1的迭代会陷入1,-1,1,-1,...的震荡死循环
要么换初值比如x(1) = 1.2
要么换其他迭代法

地狱异类 发表于 2018-11-20 11:40
是的,就跟牛顿法有时候无法收敛一样,这个是迭代算法本身的问题
要么改初值要么换其他迭代法

谢谢您的回复，感谢！
还有点问题，希望您不要介意。我是个新手，对迭代算法也是要做作业才开始接触，请问有没有推荐的讲迭代算法比较基础的教程之类的？还有就是每种算法都会出现这种不收敛的问题吗？您调整初值的时候有没有什么经验可以分享一下？
再次感谢！

书染桔灯 发表于 2018-11-20 11:41
谢谢您的回复，感谢！
还有点问题，希望您不要介意。我是个新手，对迭代算法也是要做作业才开始接触，请 ...

迭代法一般都是包含在讲数值方法的课里面,搜numerical methods有一大堆文档视频,可以挑着看看这方面的公开课

初值设定这个是个玄学的活儿,总之没有一定的方法来保证一个必定收敛的初值

收敛性跟算法有关
二分法收敛速度慢,但保证收敛;
盈不足术法也保证收敛,不过有的书里,不把这两种方法归在迭代法类目之下;
牛顿法对可解的方程收敛速度快,但不保证收敛;
割线法对可解的方程收敛速度也较快,,但同样不保证收敛;
Ridders'法保证收敛;
Muller法收敛性忘了,就记得这玩意儿好像容易算出来复数;
Brent法算是把几种方法混合,最差情形就是退化成二分法,所以也保证收敛.

总而言之,收敛速度快的难保证收敛,收敛速度慢的计算次数或耗时会显著高于快的.

地狱异类 发表于 2018-11-20 11:42
迭代法一般都是包含在讲数值方法的课里面,搜numerical methods有一大堆文档视频,可以挑着看看这方面的公 ...

谢谢您！