- 1、本文档共48页,可阅读全部内容。
- 2、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
7非线性方程求根7非线性方求根7非线性方程求根7非线性方程求根
Numerical Analysis 第七章 非线性方程(组)的数值解法 本章内容 本讲内容 非线性方程数值解法 非线性方程数值解法 二分法(对分法) 二分法 误差分析 不动点迭代 不动点迭代 连续性分析 解的存在唯一性 收敛性分析 收敛性分析 举例 局部收敛 收敛速度 收敛速度 Aitken 加速 Aitken 加速 Steffenson 加速 Steffenson 加速 举例 本讲内容 Newton 法 Newton 法 Newton 法 收敛性 举例 举例 牛顿法 简化的Newton法 Newton下山法 重根情形 重根情形 举例 弦截法与抛物线法 弦截法 收敛性 弦截法几何含义 抛物线法(略) 抛物线法 抛物线法 抛物线法 收敛性 例:用 Newton 法求 f(x) = xex – 1=0 的解 % Newton 迭代 clear; clc f = inline(x*exp(x)-1,x); df = inline(exp(x)*(1+x),x); n = 10; tol = 1e-6; x0 = 0.5; % 迭代初始值 for k = 1 : n x = x0 - f(x0)/df(x0); fprintf(k=%2d, x=%.7f\n,k,x); if abs(x-x0)tol, break, end x0 = x; end 例(P224):用 Newton 法求 f(x) = x2 – C=0 的正根 解: 对任意 x00, 总有 |q|1, 即牛顿法收敛 牛顿的优点 牛顿法是目前求解非线性方程 (组) 的主要方法 至少二阶局部收敛,收敛速度较快,特别是当迭代点充分靠近精确解时。 牛顿的缺点 对重根收敛速度较慢(线性收敛) 对初值的选取很敏感,要求初值相当接近真解 先用其它算法获取一个近似解,然后使用牛顿法 需要求导数! 线性收敛 简化的 Newton 法 基本思想:用 f’(x0) 替代所有的 f’(xk) 下山因子的取法: 从 ?=1 开始,逐次减半,直到满足下降条件 基本思想:要求每一步迭代满足下降条件 具体做法:加下山因子 ? Newton下山法 保证全局收敛 且 解法一:直接使用 Newton 法 线性收敛 解法二:改进的 Newton 法 二阶收敛 缺点:需要知道 m 的值 重根情形 令 x* 是 ? (x)=0 的单重根 解法三:用 Newton 法解 ? (x) = 0 迭代格式: 例:求 x4 - 4x2 + 4=0 的二重根 (1) 普通 Newton 法 (2) 改进的 Newton 法 (3) 用 Newton 法解 ? (x) = 0 % 重根 clear; clc f = inline(x^4 -4*x^2 + 4,x); g1 = inline(x-(x^2-2)/(4*x),x); g2 = inline(x-(x^2-2)/(2*x),x); g3 = inline(x-x*(x^2-2)/(x^2+2),x); fprintf(True solution: x=%.7f\n,sqrt(2)); n = 10; x0 = 1.5; % 迭代初始值 x1 = x0; x2 = x0; x3 = x0; for k = 1 : n x1 = g1(x1); x2 = g2(x2); x3 = g3(x3); fprintf(k=%2d, x1=%.7f, x2=%.7f, x3=%.7f\n,k,x1,x2,x3); end 弦截法与抛物线法 目的:避免计算 Newton 法中的导数,且具有较高的收敛性(超线性收敛) 弦截法(割线法):用差商代替微商 抛物线法:用二次多项式近似 f(x) 弦截法迭代格式: k = 1, 2, 3, . . . 注:弦截法需要提供两个迭代初始值 定理:设 x* 是 f(x) 的零点, f(x) 在 x* 的某邻域 U(x,?) 内有二阶连续导数,且 f’(x)?0,若初值 x0,x1 ?U(x,?),则当 U(x,?) 充分小时,弦截法具有 p 阶收敛性,其中 x y x* xk-1 xk xk+1 基本思想: 用二次曲线与 x 轴的交点作为 x* 的近似值 抛物线法 y xk-2 xk-1 xk xk+1 * * 数值分析 非线性方程求解 二分法 不动点迭代法及其加速 牛顿法、弦截法、抛物线法 求根问题的敏感性与多项式的零点 非线性方程组的数值求解 迭代格式 加速算法 收敛性 非线性方程求解介绍 二分法及其收敛性 不动点迭代及其加速 考虑方程 若 f(x) 是一次多项式,则称为线性方程;
您可能关注的文档
- 6.1.4使用plsql语对表中数据进行更新6.1.4使用plsql语句对表中数据进行更新6.1.4使用plsql语句对表中数据进行更新6.1.4使用plsql语句对表中数据进行更新.ppt
- 06-09市政计价试题0609市政计价试题06-09市政计价试题06-09市政计价试题.doc
- 6.《台阶》ppt6.《台》ppt阶》ppt.ppt
- 6.1线性空间的概念.6.线性空间的概念.6.1线性空间的概念.6.1线性空间的概念..ppt
- 6.1力课件第一节6.1力件第一节课件第一节.ppt
- 6. 食物在口腔里的变化 件16. 食物在口腔里的变化 课件16. 食物在口腔里的变化 课件16. 食物在口腔里的变化 课件1.ppt
- 6.2(新课标)同角三角函基本关系式及诱导公式(典型例题+习题+答案)老师10页6.2(新课标)同角三角函数基本关系式及诱导公式(典型例题+习题+答案)老师10页6.2(新课标)同角三角函数基本关系式及诱导公式(典型例题+习题+答案)老师10页6.2(新课标)同角三角函数基本关系式及诱导公式(典型例题+习题+答案)老师10页.doc
- 6. 基本工作底稿简介6.基本工作底稿简介6. 基本工作底稿简介6. 基本工作底稿简介.doc
- 6.1让我们选择坚强6.1我们选择坚强6.1让我们选择坚强6.1让我们选择坚强.ppt
- 5外国媒体研究(二)5外国体研究(二)5外国媒体研究(二)5外国媒体研究(二).ppt
- 第5课 控制系统的三个环节(教学设计)2023-2024学年五年级下册信息科技浙教版.docx
- 2024-2025学年高中历史必修三北师大版教学设计合集.docx
- 2024-2025学年初中音乐九年级下册苏少版(2024)教学设计合集.docx
- 2024-2025学年初中信息技术(信息科技)七年级上册川教版(2019)教学设计合集.docx
- 2024-2025学年高中地理必修二鲁教版教学设计合集.docx
- 【新课标】Unit 2 What time is it整体单元教学设计.docx
- 2024-2025学年初中信息技术(信息科技)八年级上册滇人版(2016)教学设计合集.docx
- 2024-2025学年小学信息技术(信息科技)四年级上册西师大版教学设计合集.docx
- 2024-2025学年高中地理选择性必修2中图中华地图版教学设计合集.docx
- 2024-2025学年初中信息技术(信息科技)八年级上册湘电子版(2019)教学设计合集.docx
文档评论(0)