首页 > 其他分享 >【数值分析】非线性方程求根

【数值分析】非线性方程求根

时间:2023-11-01 14:23:25浏览次数:36  
标签:phi 非线性 方程 sqrt leq 数值 不动点 求根

第7章-非线性方程求根

不动点:对于\(f(x)\),若存在\(a\)使得\(f(a)=a\),则称 \(x=a\)为\(f(x)\)的不动点。

参考链接§1.2.6 不动点

7.1.2 简单迭代法(Jacobi迭代)

\[f(x)=0 \iff x = \phi(x) \]

利用\(x_{k+1} = \phi(x_k)\)迭代求解不动点,即得方程的根。

【例】求\(f(x)=x^3-x-1=0\),在 \(x=1.5\) 附近的根。

\[x^3-x-1=0 \Rightarrow x=\sqrt[3]{x+1} \Rightarrow \left\{ \begin{array}{l} x_0 = 初始值 \\ x_{k+1} = \sqrt[3]{x_k+1}, k =0,1,2,\dots \end{array} \right. \]

计算:
$ x_0 = 1.5$
$ x_1 = \sqrt[3]{x_0+1} = \sqrt[3]{1.5+1} = 1.35721 $
$ x_2 = \sqrt[3]{x_1+1} = 1.33086 $
\(\dots\)
$ x_6 = 1.32473$
$ x_7 = 1.32472$
$ x_8 = 1.32472$

\[x^* \approx 1.32472 即为方程的根 \]

\(\left\{ x \right\}\)

迭代收敛判定:

对于方程 $x=\phi(x), \phi(x) \in C[a,b] $ ,若

(1)当$x \in [a,b] $ 时,\(\phi(x) \in [a,b]\)
(2)\(\exists 0 \leq L < 1\),使得$|\phi'(x)| \leq L < 1 $ 对 $ \forall x \in [a,b]$ 成立。

则任取\(x_0 \in [a,b]\) ,由 $x_{k+1}=\phi(x_k) $ 得到的序列\(\left\{ x \right\}_{k=0}^{\infty}\) 收敛于 \(\phi(x)\) 在 \([a,b]\) 上的唯一不动点。

并且有误差估计:

\[\begin{align} & ① |x^*-x_k| \leq \frac{1}{1-L}|x_{k+1}-x_k|\nonumber\\ & ②|x^*-x_k| \leq \frac{L^k}{1-L}|x_{1}-x_0| , \ (k=1,2,\dots) \nonumber \end{align} \]

\[\]

且存在极限 \(\lim_{k\rightarrow \infty}{ \frac{x^* - x_{k+1} }{x^*-x_k}} = \phi'(x^*)\)

上面的判定是不动点存在的充分条件

标签:phi,非线性,方程,sqrt,leq,数值,不动点,求根
From: https://www.cnblogs.com/code-pigeon/p/17802658.html

相关文章

  • 【数值分析笔记】
    第3章-函数逼近3.1内积空间3.1.1内积设$f(x),g(x)\inC[a,b],\rho(x)$是\([a,b]\)上的权函数,积分\[(f,g)=\int_a^b\rho(x)f(x)g(x)dx\]称为函数\(f(x)\)与\(g(x)\)在\([a,b]\)上的内积。\(C[a,b]\)表示在区间\([a,b]\)内连续的全体函数组成的集合。满......
  • MySQL主键使用数值型和字符型的区别
    在了解数值型主键和字符型主键之前,需要先谈谈什么是主键、主键的特性、以及字段类型的存储空间,掌握了这些才能更好的理解数值型主键和字符型主键的区别。什么是主键在数据库中,主键是用于唯一标识数据库表中的每一行数据。主键能够确保数据的唯一性,并在表中建立索引,以提高查询效率。......
  • Xilinx VIvado学习-01 数值处理之乘法(有符号)
    Verilog数值处理,在处理减法的时候,需要注意溢出问题。实例:a*b=c 1modulesi_product(2inputsigned[9:0]a,3inputsigned[7:0]b,4outputsigned[17:0]product5);6assignproduct=a*b;7endmodule仿真代码:1modulesi_product_tb;2regsys_......
  • Xilinx VIvado学习-01 数值处理之乘法(无符号)
    Verilog数值处理,在处理减法的时候,需要注意溢出问题。实例:a*b=c 1`timescale1ns/1ps2//////////////////////////////////////////////////////////////////////////////////3//Company:4//Engineer:5//6//CreateDate:2023/10/2323:33:077//......
  • 斐波那契数列&数值传递
    #include<iostream>usingnamespacestd;intp1(inta){ if(a<=2){ return1; }else{ returnp1(a-1)+p1(a-2); }}intmain(){ intn; cin>>n; cout<<p1(n); return0;} #include<iostream>usingnamespacestd;int......
  • 金蝶云星空业务单据审核校验根据条件分录数值必录
     一、打开业务单据的审核操作 二、新增规则 设置 三、测试未录入质保期限,审核时 录录入质保期限后,审核 END......
  • Python中数值类型、布尔类型、字符串类型的特征及方法的使用
    数值类型包括int(整形)、float(浮点型)、complex(复数型)1.整型(int)整数就是没有小数部分的数字,Python中的整数包括正整数、0和负整数。为任意大小的整形数,可以表示成十进制八进制、十六进制和二进制。例:>>>0b1101+100 113>>>0x10+218>>>10+212>>>1_2_3_41234 2.......
  • 关于高级定时器 重复计数值寄存器的使用介绍
    来源:https://www.cnblogs.com/liaigu/p/17782198.html在使用高级定时器进行初始化的时候,相较于通用定时器,在初始化的时候会有一个重复计数的配置,如下图:该位主要是对重复计数值寄存器进行配置,如下图:关于该配置的使用说明,具体如下:以定时器中断为例:1、一般默认情况下,将重复计......
  • C++前缀和算法应用:矩形区域不超过 K 的最大数值和
    题目给你一个mxn的矩阵matrix和一个整数k,找出并返回矩阵内部矩形区域的不超过k的最大数值和。题目数据保证总会存在一个数值和不超过k的矩形区域。示例1:输入:matrix=[[1,0,1],[0,-2,3]],k=2输出:2解释:蓝色边框圈出来的矩形区域[[0,1],[-2,3]]的数值和是......
  • Xilinx VIvado学习-01 数值处理之减法器
    Verilog数值处理,在处理减法的时候,需要注意溢出问题。实例:a-b=c moduleun_sub(inputunsigned[7:0]a,inputunsigned[7:0]b,output[7:0]sub,outputcarry);assign{carry,sub}=a-b;endmoduleViewCode仿真代码:`timescale1ns/1ps////////////......