CSU_科学计算与数学建模_学习笔记（6）

科学计算与数学建模，郑洲顺老师
平台——学堂在线 https://www.xuetangx.com/course/CSU07011000630/5882763?channel=learn_title

非线性方程的数值解法

养老保险问题分析

-时间序列实际月收益率，保险公司每月应至少有多少获益率假设投保过程{F1, F2,…F_k.F_k+1,…, F_N, F_N+1} *
+自己推导笔记

+问题，非线性方程很难解
科学工程中，很多问题归结为解方程的问题
而我们知道n>=5次的一般代数方程式是不能用代数公式求解的，1，2，3，4都有相应的求根公式（16世纪）

+采用，逐步搜索的思路
足够小的步长h，得到足够精度的近似点
像电视节目里猜产品价格，低价a和高价b之间

二分法

二分法—— 一种收敛的方法
缺点，不能求复根及偶数重根，只能求实根；收敛速度慢

迭代思路

f(x) = x -sinx - 0.5
迭代函数
构造数列，希望数列的极限存在—— 基本迭代思想

第一步，首先将方程转化为等价方程： x=g(x)
f(x) = x -sinx - 0.5
+课堂截图，设x=g(x) 逐次代入，目标构造近似解的序列

g(x)迭代过程，称为迭代函数
+重点，迭代函数多个，不同的函数产生不同的序列{xk}
并且不同迭代函数序列的收敛情况不相同

迭代函数收敛的条件

（1） 设g(x)于[a,b]一阶导数存在
（2）当x属于[a,b]时，有g(x)属于[a,b] #保证迭代过程可以进行
（3）导函数在[a,b]的绝对值有上确界，且上确界小于1
+课堂截图，迭代函数收敛的条件

+了解，把可导条件放宽——
（1）连续函数（不一定可导）
（2）|任意两个自变量对应的函数值的差| <=常数 * |两个自变量的差|

迭代解法的收敛性及改进

已知误差|x* - xk|有如上公式
+问题，什么样的算法好？
L需要足够的小，如何改进？
+了解，Atken加速算法
两项加速的新序列思想：

用两次迭代的结果做一个新的线性组合
更靠近近似解

Atken加速算法（三项加速）

用慢的收敛序列的三项，来加速得到一个新的收敛序列的方法

Newton迭代法

基本思想，将非线性方程f(x) = 0逐步线性化求解
几何意义：用切线的交点作为曲线（与x轴）交点近似的方法
用f（x）本身来构造迭代函数的思想
在什么条件下，牛顿迭代收敛？
+思考，如果x0取得很远，牛顿法还有效吗？
—— 牛顿法收敛依赖于初值的选取，要求初值在根的附近

+实例，用Newton法求解

初值，6次，精度足够高

牛顿下山法

如果x0偏离x*较远，牛顿法可能发散

保证函数下降条件

+概念，迭代收敛的阶

序列P阶收敛：
保证迭代函数在精确解附近，有直到p阶的连续导数，1阶导数到（p-1）阶导数均为零，而p阶导不为零

证明，牛顿迭代收敛阶为2阶（此处略）

弦截法与抛物法

构造插值多项式
知道f(x)的表达式，知道r-1个点，构造r次插值多项式
用r次插值多项式的根近似原函数的根

弦截法—— 用一次插值函数
抛物法—— 用二次插值函数
+了解，单点弦截法，始终固定x0和其他点做x轴交点
+分析，单点、双点、抛物线方法的收敛和收敛阶

+重点，弦截法与Newton切线法有本质区别
-牛顿切线法，计算时用到前一步的节点和导数信息
-弦截法，计算式用到前面两步的节点信息，不须计算导数

收敛

双点弦截的方法

+认识，Atken加速算法的几何意义—— 扭转了发散的问题

+了解，抛物线法按阶p=1.840收敛于根x*，要求领域内有三阶连续导数
以上方法要求一阶导数不为零，因此不适合重根

+总结，
抛物线法（1.840阶）比弦截法（1.618阶）的收敛性更接近于Newton法（2阶），二分法是线性收敛（1阶）的
以上，牛顿法的收敛是四者里最快的
这四种方法，都是用线性的思想，求解非线性方程组的根