第4章–养老保险问题.pptVIP

因此仍然能保证在 领域内 ，使算法具有二阶收敛性。在实际应用中对 重根，迭代函数可改写成 ，但由于一般不能确定常数 ，则可考虑函数： ，如果 是正整数， 。 则 ，显然 是 的单重零点，故可得切线法（即Newton法）用于 ，得到二阶收敛的迭代函数： 或 因此作为迭代函数即可找到根 ，收敛阶为2阶。 例4.2.10 对于方程 是二重根用下面三种方法求根： Newton法： 即 有上面 所确定的修改方法化简为： 三种方法均取初值 。计算结果为： 表4.2.5 部分计算结果表 经3次迭代，方法2、3均达到精度。它们都是二阶收敛的方法。而方法1是一阶的，要进行30次迭代才能达到同样的精度。 方法1 方法2 方法3 1.5 1.5 1.5 1.458333333 1.416666667 1.411764706 1.4366o7143 1.414215686 1.414211438 1.425497619 1.414213562 1.414213562 4.3 养老保险模型的求解 对4.2中建立的养老保险模型，以25岁起保为例，假设男性平均寿命为75则 ，初始值为 ，我们可以得到： 在上面两式中，分别取 和 并利用 可以求出： 其中 利用MATLAB编程编写牛顿迭代法程序，令初值为 ，迭代最大次 数为10000，求出方程的根为： 同样方法可以求出：35岁和45岁起保所获得的月利率分别为 习 题 4 1、用迭代法求解如下方程在（1，2）内的实根 2、用Newton法求 的近似解 3、用弦截法求方程 在区间（1，2）内的实根。 4、利用二分法求方程 在 内的根 中南大学数学科学与计算技术学院 * 表4.2.4 部分计算结果表 0 0.0 1 02 014 1.1461931 15 1.1461932 (3)为了求[-1.9,-1]内方程的根。由迭代方程 ，显然 ，所以迭代过程 （初值 ）不能保证收敛于 。 (4)若将方程转化为等价方程 或 则 ，且 （ 时） 所以当选取 时迭代过程 收敛。如取 ，则迭代12次有 ，且 。 由此例可见，对于方程 ，迭代函数 取不同形式，相应的迭代法产生的 收敛情况也不一样，因此，我们应该选择迭代函数时构造的迭代过程 收敛，且收敛较快。 关于迭代公式的加工： 对于收敛的迭代过程，只要迭代足够多次，总可以使结果达到任意的精度。但有时迭代收敛缓慢，从而使计算量变得很大，因此迭代过程的加速是一个很重要的课题。 设 为根 的某个预测值，用迭代公式校正一次得： 由中值定理： ， 介于 之间，若 改变不大。近似地取某常数，则由 可以期望按上式右端求得的 是比更好的近似值。 若将每得到一次改进值算作一步，并用 和 分别表示第