首页 理论教育 数值分析典型应用案例及理论分析:数值结果与收敛速度分析

数值分析典型应用案例及理论分析:数值结果与收敛速度分析

时间:2023-11-23 理论教育 版权反馈
【摘要】:图6.2y=f可看出该函数零点在(4,5)范围内,且在该区间单调。,且从图6.2中可知其一定收敛;弦截法是在牛顿迭代法的基础上,用两点所连成弦的斜率代替该点处导数:,…。MATLAB代码如下:程序输出结果见表6.1~表6.3。表6.2牛顿迭代法计算数据列表表6.3弦截法计算数据列表从本案例数据很容易看出,三种方法均收敛;仅从计算次数判断收敛速度来看,牛顿迭代法收敛速度最快,弦截法次之,二分法收敛速度最慢。

数值分析典型应用案例及理论分析:数值结果与收敛速度分析

设计一台三级圆锥圆柱齿轮减速器,已知总传动比i=50,低速级大齿轮上的工作扭矩T=45 N·m,高、中、低速级载荷系数K1=K2=K3=1.2,许用接触应力[σ]H1=[σ]H2=[σ]H3=675 MPa,圆锥齿轮齿宽系数φR=0.3,中、低速级齿轮齿宽系数分别为φd2=1,φd3=1.2,按最小长度条件分配传动比。

参考式(6.12),得

代入式(6.15),得

将i3代入式(6.15)第二个方程,得如下非线性方程:

为了更好地选用迭代方法,先画出了原方程f(x)=0中y=f(x)的图像,如图6.2所示。

图6.2 y=f(x)

可看出该函数零点在(4,5)范围内,且在该区间单调。因而本案例采用二分法牛顿迭代法和弦截法对该非线性方程进行求解。

二分法的初始区间即可设为[4,5],不断地取其中点,通过判断其和区间端点的函数值分别乘积得到的值的正负,来取下一个区间,给定精度范围,当满足后停止。(具体过程见MATLAB代码);牛顿迭代法的迭代公式为,k=0,1,…,且从图6.2中可知其一定收敛;弦截法是在牛顿迭代法的基础上,用两点所连成弦的斜率代替该点处导数,…。

MATLAB代码如下:(www.xing528.com)

程序输出结果见表6.1~表6.3。

表6.1 二分法计算数据列表

最终计算结果为4.200 2。

表6.2 牛顿迭代法计算数据列表

表6.3 弦截法计算数据列表

从本案例数据很容易看出,三种方法均收敛;仅从计算次数判断收敛速度来看,牛顿迭代法收敛速度最快,弦截法次之,二分法收敛速度最慢。下面从一般性角度来比较三种方法的收敛性和收敛的快速性。首先针对二分法来说,若原方程有解,初始范围选择正确且函数连续的话,则其一定收敛。收敛的速度和初始区间的选择有关,而且在程序运行中只涉及取中点、代入、比较的过程,即使迭代步数很多,运行速度依旧很快。对于牛顿迭代法来说,首先其收敛性和初值的选取有关,若初值x0与x*距离很远,则该方法有可能发散,其次其收敛性只能达到线性收敛,且在f′(xk)≠0时才能使用,当有重根情形时需要知道其根的重数才可使用该方法,且迭代公式改变为。针对其快速性,由于其要求出每个迭代点的导数,增加了计算繁复性,从该角度分析其快速性有提高空间;但其收敛速度较快,优于二分法。弦截法是牛顿迭代法的演变,其收敛性也和初值的选择有关;从快速性来看,该方法减少了求导数的部分,而是利用弦的斜率去代替,加快了计算速度,但是很可能出现越过精确值的情况,也就是说会出现从两侧逼近精确值的情形,因而从收敛的速度来看,仍然是牛顿迭代法更有优势。本案例中由于在取弦截法初值的第二个点时,运用了一次牛顿迭代法,多求了原函数的导函数,故速度较慢。

免责声明:以上内容源自网络,版权归原作者所有,如有侵犯您的原创版权请告知,我们将尽快删除相关内容。

我要反馈