已阅读5页,还剩19页未读, 继续免费阅读
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
线性方程与非线性方程的概述与运用,问题背景和研究目的,解方程(代数方程)是最常见的数学问题之一,也是众多应用领域中不可避免的问题之一。,求解一般非线性方程没有通用的解析方法,但如果 在任意给定的精度下,能够解出方程的近似解,则 可以认为问题已能够解决,至少可以满足实际需要。,本节主要介绍一些有效的求解方程的数值方法:二分法,迭代法 ( 牛顿法)。同时要求大家学会如何利用Matlab 来求方程的近似解。,2.6 非线性方程近似根,相关概念,如果 f(x) 是一次多项式,称上面的方程为线性方程; 否则称之为非线性方程。,线性方程 与 非线性方程,问题: 如何求连续的非线性方程 实根的近似值。,根的隔离,若函数 f(x) 在闭区间a,b上连续,且 f(a)f(b)0,则 f(x)在开区间 (a,b)内至少存在一个根。通过根的隔离,可假设此区间内存在唯一根 x*。,基本思想,二分法,将隔离区间进行对分,判断出解在某个子区间内,然后再对该子区间对分,依次类推,直到满足给定的精度为止。,算法,二分法,设方程在区间 a,b 内连续,且 f(a)f(b)0,给定精度要求 ,若有 |f(x)| ,则 x 就是f(x) 在区间 (a,b) 内的 近似根。,收敛性分析,二分法收敛性,设方程的根为 x* (an , bn ) ,又 ,所以,根据上面的算法,我们可以得到一个每次缩小一半的区间序列 an , bn ,在 (an, bn ) 中含有方程的根。,二分法总是收敛的,二分法的收敛速度较慢 通常用来给出根的一个 较为粗糙的近似。,简单迭代法, (x) 的不动点,f (x) = 0,x = (x),f (x) 的零点,(x) 称为迭代函数,若 收敛,即 ,假设 (x) 连续,则,收敛性分析,迭代法的收敛性,即,注:若得到的点列发散,则迭代法失效!,迭代法的收敛性判据,定理2.1:全局收敛,定理2.2:全局发散,定理2.3:局部收敛与发散,定理2.4:收敛速度,定义:,迭代法收敛性判断,如果存在 x* 的某个邻域 =(x*- , x* + ), 使得对 x0 开始的迭代 xk+1 = (xk) 都收敛, 则称该迭代法在 x* 附近局部收敛。,迭代法收敛性判断,L 越小,迭代收敛越快,收敛阶,为了进一步研究收敛速度问题,引入阶的概念: 记 ,如果 ( p=1时还要求01时称为超线性收敛。 p越大收敛越快。,牛顿迭代法,令:,设非线性方程 f (x)=0 , f (x) 在 xk 处作 Taylor 展开,牛顿迭代公式,k = 0, 1, 2, . .,牛顿迭代公式,牛顿法的优点,牛顿法是目前求解非线性方程 (组) 的主要方法,对于单重根迭代2阶收敛,收敛速度较快, 特别是当迭代点充分靠近精确解时。,在实际计算中,如果要求高精度,可以先用其它方法(如二分法)获得精确解的一个粗糙近似,然后再用牛顿法求解。,牛顿迭代法大范围收敛性,Matlab 解方程的函数,roots(p):多项式的所有零点,p 是多项式系数向量。,fzero(f,x0):求 f=0 在 x0 附近的根,f 可以使用 inline、字符串、或 ,但不能是方程或符号表达式!,solve(f,x):求方程关于指定自变量x的解, f 可以是用字符串表示的方程、符号表达式或符号方程; solve 也可解方程组(包含非线性); 得不到解析解时,给出数值解。,Ab:解线性方程组Ax=b。,其他 Matlab 相关函数,g=diff(f,x):求符号表达式 f 关于 x 的导数 g=diff(f):求符号表达式 f 关于默认变量的导数 g=diff(f,x,n):求 f 关于 x 的 n 阶导数,diff,f 是符号表达式,也可以是字符串,默认变量由 findsym(f,1) 确定, syms x f=sin(x)+3*x2; g=diff(f,x), g=diff(sin(x)+3*x2,x),作业,每题分别用两种一步迭代法(要求写出迭代格式): 1) Newton迭代法; 2)自己构造的非牛顿切线或割线法迭代格式(需讨论收敛性) 根据迭代格式用计算机(器)求下列非线性方程的根:,迭代法的加速,设迭代 xk+1 = (xk) ,第 k 步和第 k+1 步得到的 近似根分别为 xk 和 (xk) ,令,其中 wk 称为加权系数或权重。得新迭代 xk+1 = (xk),松弛迭代法,松弛法迭代公式:,松弛法具有较好的加速效果,甚至有些不收敛的迭代格式,通过加速后也能收敛。,缺点:每次迭代都需计算导数,Altken 迭代法,Al
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 2022-2023年人教版九年级物理上册期末考试卷及答案2
- 人教部编版五年级语文上册期末试卷【带答案】
- 2023年人教版九年级化学上册期末考试题(含答案)
- 新人教版七年级物理上册期末试卷【及答案】
- 苏教版七年级《生物》上册期末考试题及答案【新版】
- 2022-2023年人教版四年级数学上册期末考试题及答案【完美版】
- 人教部编版九年级化学(上册)期末试卷及答案(各版本)
- 2021-2022年部编版六年级语文上册期末考试卷(完整)
- 小学四年级数学上册期末考试题(必考题)
- 部编版四年级上册语文《期末》测试卷及答案一
- 《C++程序设计》分类练习题
- HJ 484-2009 水质 氰化物的测定 容量法和分光光度法(正式版)
- 2024届安徽省合肥市名校联考中考二模化学试题含解析
- 国开2024年春《形势与政策》大作业答案
- 新能源汽车火灾事故处置程序及方法
- 中国移动SEQ平台介绍
- 心电图导联连接快速确定位置及图示-常规胸导联电极位置
- 肺痿肺痈咳嗽上气病脉证并治第七.ppt
- 空调制冷第一讲制冷原理(压焓图).ppt
- ppt 各种人物图标ppt课件
- 无线电波的发射传播接收.ppt
评论
0/150
提交评论