单个方程的不动点迭代法
1. 不动点和不动点迭代法
设 f 是连续函数,为了解方程
把方程变换为等价的方程
其中φ 是连续函数。若x* 是f 的零点,即f(x*)=0 则有
x*称为函数φ 的一个不动点.构造的迭代法
这称为一种不动点迭代法, 称为迭代函数。由(1.3)式产生的序列 {xk}如果满足 ,
则 x*是φ的一个不动点。即方程(1.1)的一个根。
2.函数在区间上不动点的存在性和唯一性及推论
定理 2.1 设 φ∈C[a,b]且满足
则 φ 在[a,b] 上一定存在不动点。若φ 满足(2.1),又存在常数L∈(0,1) 使
则φ 在[a,b] 的不动点是唯一的。
条件 (2.2)称为Lipschitz条件,L称为Lipschitz常数。
定理 2.1的推论 设函数φ满足(2.1),又 φ'∈C[a,b],且存在常数L∈(0,1) ,使
则φ 在[a,b] 上存在唯一的不动点。
3.迭代法在区间[a,b]的收敛性
定理3.1 设φ∈C[a,b],满足条件(2.1)和(2.2),其中L∈(0,1)。则对任意的x0∈[a,b],由迭代法(1.3) 产生的序列{xk}收敛到φ在[a,b]上的不动点x* ,而且对整数P≥1,
定理3.1描述了对任意的x0∈[a,b],{xk}收敛到[a,b]唯一的不动点x*,这可以说是在区间[a,b]上的全局收敛性。
局部收敛性
定义 设函数φ在区间[a,b]上有不动点x* , 如果存在x*的一个邻域
对任意的 x0∈S,迭代法(1.3)产生的序列{xk}∈S,且{xk}收敛到x*,就称迭代法(1.3)局部收敛。
定理3.2 设 x*为函数φ 的不动点, 在φ'在x*的某个邻域上存在且连续,且
则迭代法(1.3)局部收敛。
4.牛顿迭代法
设x*是f的零点,并已有一个近似值xk≈x*,如果f''存在且连续,由Taylor展开式得
,
其中ξ在xk和x*之间。因为f(x*)=0,如果f'(xk)≠ 0,略去最后一项即得
.
这就是Newton迭代法,又称Newton-Raphson迭代法。
牛顿迭代法的局部收敛性
Newton迭代法对应的迭代函数是
定理4.1 设f(x*)=0 ,f'(x*)≠ 0,且f在包含x*的一个区间上有二阶连续导数,则牛顿迭代法(4.2)局部收敛到x* 。
练习题:设 a>0, 求平方根
可以转化为解方程
求迭代公式及迭代公式在(0,+∞)上的全局收敛性。并和二分法比较。
练习题解答请参考:点击打开链接
5.参考文献
最优化理论与算法(第二版).陈宝林编著.
单个方程的不动点迭代法相关推荐
- FPI(Fixed-point Iteration)不动点迭代法——迭代求方程的方法
一 不动点是什么? 不动点,其实定义比较简单,对于一些方程,例如f(x)=cosx,那么令cosx=x的点就是函数的不动点,说白了,就是y=x这条直线与函数曲线的交点.这个不动点有什么用呢?请继续往下 ...
- 计算方法实验:方程求根二分法、不动点迭代法、牛顿法
计算方法实验一,方程求根 分别用Matlab和C写了一下,初学Matlab,如有不足还请指正. 实验内容: 1)在区间[0,1]内用二分法求方程 e x + 10 ∗ x − 2 e^x+10*x-2 ...
- matlab求方程实根,简单迭代法求方程根的MATLAB程序
function [x_reality,n_reality] = Simple_stepit( f_name,x_start,tolerance,n_limit) %% % 简单迭代法(也叫不动点迭代 ...
- matlab牛顿法求区间根程序,MATLAB用二分法、不动点迭代法及Newton迭代(切线)法求非线性方程的根...
一.实验原理 二.实验步骤 三.实验过程 1.(程序) (1)二分法:求 在区间(1,2)之间的根,取 (a)bipart.m: function [x,m]=bipart(fun,a0,b0,to ...
- 不动点迭代法求函数根(非线性方程求解)
问题描述 在[a,b]区间内寻找方程x**5-2*x-1=0的根的初始近似值位置,确定不动点迭代的初始点(可能有多个),然后使用不动点迭代法求方程的根(可能有多个根).前后两次迭代的差的绝对值小于de ...
- 重根的二阶迭代法matlab,MATLAB用二分法、不动点迭代法及Newton迭代(切线)法求非线性方程的根...
作者:凯鲁嘎吉 - 博客园 http://www.cnblogs.com/kailugaji/ 一.实验原理 二.实验步骤 三.实验过程 1.(程序) (1)二分法:求 在区间(1,2)之间的根,取 ...
- 【MATLAB实验】数学实验实验求根问题三种方法、积分画图及古典密码设计(二分法、牛顿法、不动点迭代法)
目录 特征值与特征向量 求最简型 编程题 求积分问题(quad.trapz.int) 求根问题(二分法.不动点迭代法.牛顿法) 1.二分法: 2.不动点迭代法 3.牛顿切线法 画图题 画螺旋曲线图 画 ...
- [数值分析]不动点迭代法
1. 1. 1.不动点与不动点迭代法 对于方程 f(x)=0(1.1) (1.1) f ( x ) = 0 f(x) = 0\tag{1.1} 改写成 x=φ(x)(1.2) (1.2) x = φ ...
- 工程计算——不动点迭代法
不动点迭代法 理论基础: 基本思想 f(x)=0 等价变换 x=ф(x) f(x)的根 等价变换 ф(x)的不动点 从一个初值x0出发,计算x1=ф(x0),x2=ф(x1),--,x(k+1)=ф( ...
最新文章
- python中label组件参数_Python tkinter(六) 标签(Label)组件的属性说明及示例
- 信号完整性与电源完整性分析_「1月11日 免费视频直播」高速PCB信号完整性和电源完整性仿真实战分析、ADS创新解决方案分享-送给工程师的一份新年礼物...
- Python3 Scrapy爬取猫眼TOP100代码示例
- 从零开始搭建spring-cloud(2) ----ribbon
- 【Linux驱动】字符设备驱动
- 精简jdk包_在JDK 12精简数字格式中使用最小分数数字
- Get Requests with Json Data Get Requests with Url Parameters
- python剪切文件如何恢复_用python实现的可以拷贝或剪切一个文件列表中的所有文件...
- arm export 汇编_C/C++与汇编混合编程有什么好处?
- Mysql 一条SQL语句实现批量更新数据,update结合case、when和then的使用案例
- 黑猿大叔-译文 | TensorFlow实现Batch Normalization
- [转]22条经典的编程引言
- 编译原理-第一节(东南大学廖力老师)
- DIY智能小车篇(四):常见问题 BUG汇总
- 无线上网认证之Portal认证
- 可以实现树形结构的设计模式(组合模式)
- iOS 音乐 播放器 声音
- 你知道两台Linux之间如何传输文件吗?
- spring入门学习粗解(一)
- 屏蔽微信分享到朋友圈等按钮
热门文章
- h5打开麦克风权限录音_通过H5实现html页面的录音和播放 (Recorder用于html5录音)...
- 单臂路由配置实例学习记录
- Mokee android 10.0代码下载编译总结
- python单片机自动浇花_MicroPython实例之TurnipBit开发板自动浇水实验
- html实现价格随数量变化,一种商品价格智能分段递减定价方法与流程
- 喇叭自动化生产线点胶划线工艺
- 错题集一,基础错题集
- 四核处理器_攀升推出 MaxBook P1 笔记本:四核CPU+15.6英寸大屏,1899 元
- 【愚公系列】2023年06月 网络安全(交通银行杯)-APFS
- 2022升降机司机考题及模拟考试