不动点迭代法
- 格式:doc
- 大小:142.00 KB
- 文档页数:13
非线性方程求根——不动点迭代法一、迭代法的基本思想迭代法是一种逐次逼近的方法,用某个固定公式反复校正根的近似值,使之逐步精确化,最后得到满足精度要求的结果。
例:求方程x 3-x -1=0 在x =1.5 附近的一个根。
解:将所给方程改写成31x x =+假设初值x 0=1.5是其根,代入得33101 1.51 1.35721x x =+=+=x 1≠x 0,再将x 1代入得33211 1.357211 1.33086x x =+=+=x 2≠x 1,再将x 2代入得33321 1.330861 1.32588x x =+=+=如此继续下去,结果如下:k x kk x k 01234 1.51.357211.330861.325881.324945678 1.324761.324731.324721.32472仅取六位数字,x 7与x 8相同,即认为x 8是方程的根。
x *≈x 8=1.32472这种逐步校正的过程称为迭代过程。
这里用的公式称为迭代公式,即311k k x x +=+k =0,1,2,……若x *满足f (x*)=0,称x *为ϕ(x )的一个不动点。
将连续函数方程f (x )=0改写为等价形式:x=ϕ(x ),其中ϕ(x )也是连续函数。
1()k k x x ϕ+=(k =0,1,……)不动点迭代法就是指以迭代格式二、不动点迭代法进行迭代求解的方法。
其中ϕ(x )称为迭代函数。
三、不动点迭代法的实现——MATLAB程序function[root,n]=stablepoint_solver(phai,x0,tol) if(nargin==2)tol=1.0e-5;enderr=1;root=x0;n=0;while(err>tol)n=n+1; %迭代次数r1=root;root=feval(phai,r1); %计算函数值err=abs(root-r1);end程序应用示例:function testmain% x^3-x-1=0% =>x^3=1+x% =>x=(1+x)^(1/3)ph=inline(‘(1+x)^(1/3)’,’x’);[root,n]=stablepoint_solver(ph,1)运行结果:root=1.3247n=8若对任意x 0∈[a , b ],由不动点迭代格式lim *k k x x →∞=则称迭代过程收敛,且x *=ϕ(x *)即f (x*)=0,x *为不动点。
函数迭代与不动点迭代法函数迭代和不动点迭代法是数值分析中常用的数值迭代方法,用于求解方程或优化问题。
它们在不同的应用领域都有广泛的应用,并且具有简单易懂、易于实现等优点。
本文将介绍函数迭代的基本原理和步骤,并详细介绍不动点迭代法的定义、性质以及求解过程。
函数迭代函数迭代是一种基本的数值迭代方法,用于求解非线性方程或优化问题。
它的基本思想是通过多次迭代,使得每次迭代得到的结果趋近于方程的根或优化问题的极值点。
函数迭代的基本步骤如下:1.选择一个初始值x0作为迭代的起点。
2.根据迭代公式x n+1=f(x n),计算出下一个迭代点x n+1。
3.判断是否达到迭代的停止条件。
如果满足停止条件,则输出近似解x n+1;否则,返回第2步。
函数迭代的收敛性与迭代函数f(x)的选择密切相关。
如果函数迭代收敛,即x n收敛于方程的根或优化问题的极值点,那么我们可以通过多次迭代得到近似解。
反之,如果函数迭代发散或者收敛速度非常慢,那么我们需要考虑其他的数值方法。
不动点迭代法不动点迭代法是函数迭代的一种特殊形式,它通过将方程转化为f(x)=x的形式,求解方程的根或优化问题的极值点。
不动点迭代法的基本思想是选择一个适当的迭代函数g(x),通过迭代公式x n+1=g(x n),不断迭代,直到找到满足f(x)=x的不动点。
不动点迭代法的步骤如下:1.将方程f(x)=x转化为g(x)=x的形式,即f(x)=x等价于g(x)−x=0。
2.选择一个初始值x0作为迭代的起点。
3.根据迭代公式x n+1=g(x n),计算出下一个迭代点x n+1。
4.判断是否达到迭代的停止条件。
如果满足停止条件,则输出近似解x n+1;否则,返回第3步。
不动点迭代法的关键是选择合适的迭代函数g(x)。
迭代函数g(x)应该满足以下条件:1.在方程f(x)=x的根或优化问题的极值点附近,迭代函数g(x)的导数g′(x)存在且连续。
2.在方程f(x)=x的根或优化问题的极值点附近,满足|g′(x)|<1。
浅谈用不动点迭代法求解非线性方程解的教学方法
不动点迭代法是求解非线性方程的一种有效的解法,它的核心思想是:通过迭代,使得方程的某个特定解变得更加接近,最终达到收敛的程度,从而得到最终的解。
一、教学目标
1、让学生了解不动点迭代法的基本原理;
2、让学生能够熟练使用不动点迭代法求解非线性方程;
3、让学生能够分析不动点迭代法的收敛性,并能够给出合理的迭代步骤。
二、教学步骤
1、讲解不动点迭代法的基本原理:
首先,教师应给学生讲解不动点迭代法的基本原理,包括它的定义、其基本思想、应用场景等。
2、提供实例:
然后,教师可以提供一些简单的实例,让学生熟悉不动点迭代法的使用,并能够熟练求解。
3、分析收敛性:
最后,教师可以让学生分析不动点迭代法的收敛性,并给出合理的迭代步骤,以便更好地求解非线性方程。
三、教学评价
教学评价主要是通过实例检验学生对不动点迭代法的掌握情况,以及分析收敛性的能力,以便更好地指导学生掌握不动点迭代法。
第6章方程与方程组的迭代解法§6.2不动点迭代法及其收敛定理将非线性方程/(X)= 0化为一个同解方程X = 0(兀) (2)并且假设旗兀)为连续函数任取一个初值兀。
,代入(2)的右端,得= 0(兀0 )继续兀2=0(兀1)X* = <p(x k)(k =0,12…•>⑶称0(兀)为迭代函数,称叫为第&步迭代值 如果存在一点兀*,使得迭代序列{乜}满足(4)例1 •用迭代法求解方程2x 3-x-l = O 解:(1)将原方程化为等价方程limx^ =兀 *ksx = 2x3— 1如果取初值兀。
=0,由迭代法(3),得(2)如果将原方程化为等价方程V 妙仍取初值兀0=。
書 «0.79371'7^37依此类推,得 x2 = 0.9644x3 = 0.9940x4 = 0.9990 x5 = 0.9998 x6= 1.0000 x7 = 1.0000迭代函数的构造有关 已经收敛,故原方程的解为x —1.0000什么形式的迭代法 能够收敛呢?迭代过程的收敛性定理1设迭代函数卩⑴在SQ上连续,且满足(1)当兀〃耐,a <卩(兀)V";(2)存在一正数L,满足Ovtvl,且有I 03 E L------- (5)则1。
・方程r=0(x)在也上呐有唯一解兀*2。
.对于任意初值To e[a,b],迭代法%1 =卩(比)均收敛于**3°. x k—x^ <X, — X *L1-L< ---------一(局部收敛性)--(6)77)设f (x) = x - ^>(x),则/(兀)在⑷勿上连续可导由条件(1) /(a) = a- 0(Q)< 0/(Z?) = b —(p(b) > 0由根的存在定理,方程/(兀)=0在[°力]上至少有一个根证:由| 0‘(x) \< L<1f\x)^l-cp\x) >0贝在[%]上单调递增,/(%) = 0在[%]上仅有一个根所以1°.方程% =(p{x)在[%]内有唯一解X*2°.对于迭代法x k+1 =(p(x k\由微分中值定理| x k+1— X*— 0(耳)—0(h) =(p\^)(x k—x*)兀上+1 —X尸0(耳)—0(忑一1) =0(己)(耳—耳一1)由于|0(兀)\<L耳+1—血I — 4耳—耳一1V丨妙X k+l ~ X k 3 LXk—耳-1-x*< Lx k—x^-LXk+i _ 兀 * 一(忑+1 —<L|x,+1-x*| + L|(x,+1忑+1—兀法X k+1 ~ X k3匸步1-对由于乙 < 1, lim (忑- %*) = 0ks因此对任意初值兀0/迭代法忑+1 = 0(耳)均收敛于兀*Z?< X, — x n 1-L 1°证毕.< ------ 1-LL 21-LXk - X k-1X k-1 ~ Xk-2Xk —X 決L< ------1-LXk - ^-1定理1指出,只要构造的迭代函数满足I 0(%) 乙v 1 迭代法忑+产0(忑)就收敛对于预先给定的误差卩匪即要求I忑-,只要L1-L因此,当丨忑—耳- 迭代就可以终止,忑可以作为方程的近似解x k ~ X k-1,' —L〜< ----------- S〜定义1:如果存在/的某个邻域R.x-x <,5使迭代过程x k+l =(P(")对于任意初值X。
李雅普诺夫第一方法和第二方法李雅普诺夫第一方法,也称为不动点迭代法或迭代法。
这种方法基于一个重要的定理,即如果函数g(x)在给定区间[a,b]上连续,并且对于这个区间上任意的x,都有g(x)在[a,b]内的闭区间上有g(x)∈[a,b],那么在这个区间上存在唯一的不动点c,满足c=g(c)。
所谓不动点,即函数g(x)的值等于其自变量x的值。
李雅普诺夫第一方法的核心思想是通过迭代计算不动点c的近似值,即x_n+1=g(x_n),不断逼近真实的解。
迭代的过程中,从一个初始值x_0开始,通过将x_0代入g(x)得到新的近似值x_1,再将x_1代入g(x)得到新的近似值x_2,以此类推,直到达到预定的精度要求或者迭代次数时停止迭代。
具体地,李雅普诺夫第一方法的步骤如下:1.选择一个合适的初始值x_0。
2.根据迭代公式x_n+1=g(x_n)计算x_1,x_2,...,x_n,直到满足停止条件。
3.如果满足停止条件,则该迭代过程收敛于不动点c,并且c是非线性方程的实根的一个近似解。
李雅普诺夫第二方法,也称为牛顿迭代法,是一种更加高效的求解非线性方程的数值计算方法。
牛顿迭代法的基本思想是利用函数的局部线性近似来逼近非线性方程的根。
根据泰勒级数展开,可以将非线性方程f(x)=0在一些近似解x_n的邻域中展开成一个一次项和高阶项的级数。
利用一次项的值来逼近非线性方程的解,可以得到迭代公式x_n+1=x_n-f(x_n)/f'(x_n),其中f'(x_n)代表函数f(x)在x_n处的导数值。
具体地,李雅普诺夫第二方法的步骤如下:1.选择一个合适的初始值x_0。
2.根据迭代公式x_n+1=x_n-f(x_n)/f'(x_n)计算x_1,x_2,...,x_n,直到满足停止条件。
3.如果满足停止条件,则该迭代过程收敛于非线性方程的实根,并且x_n是方程的一个近似解。
与李雅普诺夫第一方法相比,李雅普诺夫第二方法的优点在于收敛速度更快。
泛函分析中的不动点迭代方法泛函分析是数学中的一个重要分支,研究的是函数空间中的元素及其性质。
在泛函分析领域中,不动点迭代方法是一种重要的解题工具。
本文将介绍不动点的基本定义和性质,并详细阐述如何应用不动点迭代方法解决泛函分析中的问题。
一、不动点的定义在泛函分析中,给定一个函数或者算子,我们称其不动点是指在函数或者算子作用下仍保持不变的点。
具体而言,对于一个函数f,若存在一个点x使得f(x) = x成立,则称x为f的不动点。
对于一个算子T,若存在一个向量v使得T(v) = v成立,则称v为T的不动点。
二、不动点的性质不动点在泛函分析中具有一些重要的性质,这些性质是不动点迭代方法能够顺利应用的基础。
以下是其中几个主要性质的介绍:1. 唯一性:不动点不一定是唯一的。
一个函数或算子可能有多个不动点。
2. 不动点集:对于一个函数或算子的所有不动点构成了一个集合,称为不动点集。
3. 不动点不变性:若一个函数或算子的不动点集是它自身的子集,则该函数或算子称为具有不动点不变性。
三、不动点迭代方法的应用不动点迭代方法是一种常用的解决泛函分析问题的方法,其基本思想是通过迭代过程逐步逼近不动点。
下面将介绍不动点迭代方法的一般步骤:1. 确定函数或算子:首先需要确定一个函数或算子,假设为f或T,且希望求解的不动点为x。
2. 构造递推公式:根据函数或算子的不动点定义,可以构造一个递推公式来逐步逼近不动点。
3. 进行迭代计算:根据递推公式进行迭代计算,直到满足停止准则或达到预设的精度要求。
4. 检验收敛性:对于得到的逼近结果,需要进行收敛性分析,确保逼近的不动点与真实不动点的误差在可接受范围内。
四、案例应用以一个简单的线性方程组为例来说明不动点迭代方法的应用。
设有一线性方程组Ax=b,其中A为已知系数矩阵,x为未知变量向量,b为已知常数向量。
我们可以将该方程组转化为不动点问题进行求解。
首先,根据线性方程组的定义,有Ax=b,可以将其转化为x=Bx+c的形式,其中B为系数矩阵,c为常数向量。
不动点迭代法c语言不动点迭代法是一种求解非线性方程的数值方法。
在这个过程中,通过迭代公式不断逼近方程的不动点(即方程的根)。
以下是一个使用C语言实现不动点迭代法的示例代码:```c#include <stdio.h>#include <math.h>// 定义迭代函数double iterate(double x){return (x - (x * x * x - 1) / (3 * x * x - 1));}// 计算不动点迭代法的迭代次数int iterations(double x0, double tol, double max_iter){int n = 0;double x = x0;double x_prev;do{x_prev = x;x = iterate(x);n++;} while (fabs(x - x_prev) > tol && n < max_iter);return n;}int main(){double x0, tol, max_iter;printf("请输入初始值x0、收敛容差tol和最大迭代次数max_iter:");scanf("%lf %lf %lf", &x0, &tol, &max_iter);int n = iterations(x0, tol, max_iter);printf("迭代次数:%d\n", n);printf("不动点近似值:%lf\n", x0);return 0;}```这个示例代码中,我们首先定义了迭代函数`iterate()`,它用于计算方程的迭代值。
然后,我们实现了`iterations()`函数,该函数计算不动点迭代法的迭代次数。
在`main()`函数中,我们获取用户输入的初始值、收敛容差和最大迭代次数,并调用`iterations()`函数求解不动点。
不动点迭代算法
以下是一份简单的不动点迭代算法的示例:
```
1. 初始化一个初始猜测值 x0。
2. 设置一个迭代次数的上限 max_iterations。
3. 设置一个收敛准则 epsilon,用于确定算法何时停止迭代。
4. 初始化一个计数器 iterations,用于记录迭代次数。
5. 重复以下步骤直到满足停止条件:
a. 计算下一个迭代值 x_{i+1},根据具体迭代公式。
b. 增加计数器 iterations 的值。
c. 如果达到最大迭代次数 max_iterations 或满足收敛准则(|x_{i+1} - x_i| < epsilon),则停止迭代。
6. 输出最终的迭代结果 x。
注意:具体的迭代公式需要根据问题的具体情况来确定,此处省略。
```
以上是一个通用的不动点迭代算法框架,你可以根据具体的问题,确定使用何种迭代公式。
在实际问题中,需要根据问题的特点选择合适的初始猜测值、收敛准则和迭代次数上限来保证算法的有效性和效率。
matlab不动点迭代法
Matlab不动点迭代法是一种求解非线性方程的方法,其原理是将原方程变形为x=f(x),然后通过对f(x)进行迭代,不断逼近方程的解。
具体实现时,我们可以先给出一个初值x0,然后通过x1=f(x0),x2=f(x1),x3=f(x2)……不断迭代,直到满足某个精度要求为止。
在Matlab中,实现不动点迭代法的代码非常简单,只需要编写一个函数来表示要求解的方程f(x),然后循环迭代即可。
这种方法的优点是收敛速度较快,但需要注意的是,对于不同的初始值,可能会得到不同的解,因此需要谨慎选择初始值。
- 1 -。
一、概述不动点迭代法是数值分析中常用的一种数值求解方法,广泛应用于求解方程及优化问题。
而判断不动点迭代法的收敛速度,对于有效地应用该方法具有重要意义。
本文将针对不动点迭代法的收敛速度判断准则展开讨论,以期为相关领域的研究和应用提供一定的参考。
二、不动点迭代法概述不动点迭代法(Fixed-Point Iteration Method)是一种通过不断迭代逼近解的常见求解方法。
其基本思想是利用迭代公式不断更新初始值,直至满足一定的条件为止。
其迭代公式通常具有如下形式:xn+1 = g(xn)其中,xn为第n次迭代的近似解,xn+1为第n+1次迭代的近似解,g(x)为迭代函数。
不动点迭代法的核心在于选择合适的迭代函数g(x),并通过迭代逼近不动点,即满足x = g(x)的点,从而得到近似解。
三、不动点迭代法的收敛速度不动点迭代法的收敛速度是指在迭代过程中,解逼近真实解的速度。
通常情况下,我们希望迭代能够快速收敛,即在迭代次数较少的情况下就能得到满足精度要求的近似解。
判断不动点迭代法的收敛速度是一个至关重要的问题。
四、判断不动点迭代法收敛速度的准则判断不动点迭代法收敛速度的准则有多种,下面将介绍几种常用且较为实用的方法:1. 利普希茨常数条件利普希茨常数条件是判断不动点迭代法收敛速度的重要准则之一。
对于迭代函数g(x),如果存在一个常数L,满足对于任意x1和x2有: |g(x1) - g(x2)| <= L|x1 - x2|则称L为迭代函数g(x)的利普希茨常数。
此时,如果L < 1,则不动点迭代法收敛速度较快,反之则收敛速度较慢。
2. 收敛域分析收敛域分析是判断不动点迭代法收敛速度的另一种常用准则。
通过对迭代函数的性质进行分析,可以确定不动点迭代法的收敛速度。
对于某些特定的函数形式,可以利用收敛域的性质来判断不动点迭代法的收敛速度。
3. 收敛速度估计收敛速度估计是通过对迭代过程中的误差进行分析,从而估计不动点迭代法的收敛速度。
二分法、不动点迭代法和牛顿迭代法是数值计算中常用的三种求根方法。
它们在不同的数学领域及实际问题中有着广泛的应用。
本文将对这三种方法进行介绍和比较。
一、二分法1. 原理二分法是一种基于区间不断缩小的求根方法。
其原理是通过在函数值的两个不同点处得到异号的情况下缩小区间来逼近实根。
具体过程为:首先确定一个区间[a,b],使得f(a)和f(b)异号,然后将区间一分为二,取中点c=(a+b)/2,若f(c)为零或在一定误差范围内,则c即为所求的根;否则,根据f(a)和f(c)的符号确定新的区间[a,c]或[c,b],重复上述步骤,直到满足要求。
2. 特点二分法的优点是简单易实现,对于连续且单调函数一定能收敛。
但其缺点是收敛速度较慢,尤其在根附近时迭代次数较多。
二、不动点迭代法1. 原理不动点迭代法是求解方程f(x)=0的一种迭代方法,通过将方程变换为x=g(x),其中g(x)为连续函数,然后通过不断地迭代计算得到方程的根。
具体过程为:给定一个初始近似值x0,通过不断迭代计算xn+1=g(xn)来逼近实根。
2. 特点不动点迭代法的优点是迭代过程简单,不需要对函数进行求导。
但其缺点是收敛性有一定要求,不是所有的g(x)函数都能得到收敛结果。
三、牛顿迭代法1. 原理牛顿迭代法是一种通过不断线性化函数来逼近方程根的方法。
其原理是通过对函数f(x)进行泰勒展开,并取展开式的线性部分来进行迭代计算。
具体过程为:给定一个初始近似值x0,通过不断迭代计算xn+1=xn-f(xn)/f'(xn)来逼近实根。
2. 特点牛顿迭代法的优点是收敛速度较快,在根附近有二次收敛性。
但其缺点是需要对函数进行求导,且初始值的选取对迭代结果有一定影响。
二分法、不动点迭代法和牛顿迭代法都是求解方程根的有效方法,各有其优缺点和适用范围。
在实际应用中,根据问题的特性和计算要求来选择适当的方法,以达到准确和高效的求解目的。
4. 比较与应用通过对二分法、不动点迭代法和牛顿迭代法的介绍,我们可以对它们进行比较与应用。
不动点迭代法举例
嘿,朋友们!今天咱就来讲讲不动点迭代法。
比如说,就像你要找到一个隐藏在迷宫里的宝藏!
你瞧啊,假设有个方程 f(x)=x^3+2x-1,咱就想办法用不动点迭代法
来找它的解。
这就好比是在茫茫大海中寻找那根能让你靠岸的锚索!咱先随便选一个初始值,比如说 x0=0。
然后呢,通过这个迭代公式 xn+1=f(xn)
不断地计算下去。
就好像你在一步一步地摸索着走向宝藏的道路!
“哎呀,这能找得到吗?”你可能会这样问。
哈哈,别着急呀!就这么一步步来,你会发现真的能逐渐接近那个神秘的解。
想象一下,就如同你在黑暗中渐渐看到了一丝亮光!
再比如说,在解决一个复杂的工程问题时。
哇,那问题就像是一座难以翻越的高山,而不动点迭代法就是你的登山杖呀!大家一起讨论着,“咱用这个方法试试吧!”“行啊,那就干呗!”然后就齐心协力地开始计算。
当每一次迭代后看到结果一点点变化,那兴奋感,就像发现了宝藏的一角!
“这也太神奇了吧!”有人会忍不住惊叹。
可不是嘛,不动点迭代法就是这么厉害。
它就像一把神奇的钥匙,能打开那些看似难以解开的谜题之门。
我觉得啊,不动点迭代法就是我们在数学世界里探索的有力工具。
它虽然有点复杂,但只要我们用心去理解、去运用,就能发现它的巨大魅力和强大力量!所以呀,别害怕去尝试,勇敢地利用不动点迭代法去攻克那些难题吧,朋友们!。
中国矿业大学(北京)理学院
数值分析实验报告
实验名称 不动点迭代法求方程的近似根 实验时间 201年3月20日
组长签名
龙纯鹏 班级 信息与计算
科学(1)班
学号
11107200110 成绩
组员签名
11107200101
11107200102
11107200103 11107200119
11107200120
一、实验目的,内容 二、相关背景知识介绍 三、代码
四、数值结果 五、计算结果的分析 六、计算中出现的问题,解决方法及体会 一、实验目的,内容
5、先确定方程x e x 51=
的一个收敛的有根区间[a,b], 然后用不动点迭代法求x
e x 5
1=在此有根区间的近似根,初值0x 自己确定,要求根精确到510-,并求迭代次数. ,掌握二分法,会编写相关代码。
二、相关背景知识介绍
(1)算法原理或计算公式
二分法的基本思路是通过计算隔根区间的中点,逐步将隔根区间缩小,从而可得方程
的近似根数列{}n x 。
ε<-≤-+1*2
k k a
b x x
(2)程序设计思路 二分法原理:
不妨设方程0)(=x f 在隔根区间],[b a 上.0)(,0)(><b f a f . (1)将区间],[b a 二分得中点
2
b a +,计算)(x f 在中点的函数值)2(b a f +,若0)2(=+b
a f ,则2
b
a x +=
*就是方程的根; 否则,若0)2(<+b a f ,由于)(x f 在]2,[b a a +内不变号,所以隔根区间变为],2
[b b
a +.同理,若0)2(>+
b a f , 隔根区间变为]2
,[b a a +.将新区间记为[]11,b a .
(2)将[]11,b
a重复上述步骤(1).得到一系列隔根区间:
[][][]⋯⋯⊃⊃⋯⋯⊃⊃k k b b a b ,a ,,a 11
并有),(,0)(f )(f *k k k k b a x b a ∈<*,且后一区间的长度都是前一区间长度的一半,所以[]k k b ,a 的长度为)(,2
∞→-=
-k a
b a b k k k ,区间[]k k b ,a 的长度趋向于零,即这些区间最终收缩于一点*x ,
显然*x 就是方程0)(=x f 的根。
三、代码(Matlab ) clear x0=1
e=10^(-5) k=1
x1=sqrt(0.2*e^x0)
while (abs(x0-x1)>e) k=k+1 x0=x1
x1=sqrt(0.2*e^x0) end x0
x0 = 0.1753 x1 = 0.1630 k = 100
x0 = 0.1630 x1 = 0.1750 k = 101
x0 = 0.1750 x1 = 0.1633 k = 102
x0 = 0.1633 x1 = 0.1747 k = 103 x0 =
0.1747 x1 = 0.1636 k = 104
x0 = 0.1636 x1 = 0.1744 k = 105 x0 = 0.1744 x1 = 0.1639 k = 106
x0 = 0.1639 x1 = 0.1741
k = 107
x0 = 0.1741 x1 = 0.1642 k = 108
x0 = 0.1642 x1 = 0.1738 k = 109
x0 = 0.1738 x1 = 0.1644 k = 110
x0 = 0.1644 x1 = 0.1735 k = 111
x0 = 0.1735 x1 = 0.1647 k = 112
x0 = 0.1647 x1 = 0.1733 k = 113
x1 =
0.1653
k =
118
x0 =
0.1653
x1 =
0.1727
k =
119
x0 =
0.1727 x1 =
0.1655 k =
120
x0 =
0.1655 x1 =
0.1725 k =
121
x0 =
0.1725 x1 =
0.1657 k =
122
x0 =
x1 =
0.1666 k =
134
x0 =
0.1666 x1 =
0.1714 k =
135
x0 =
0.1714 x1 =
0.1668 k =
136
x0 =
0.1668 x1 =
0.1712 k =
137
x0 =
0.1712 x1 =
0.1669 k =
138
x0 =
0.1669 x1 =
0.1711 k =
139
x0 =
0.1711
x1 =
0.1670 k =
140
x0 =
0.1670 x1 =
0.1710 k =
141
x0 =
0.1710 x1 =
0.1671 k =
142
x0 =
0.1671
x1 =
0.1709 k =
143
x0 =
0.1709 x1 =
0.1672 k =
144
x0 =
0.1672 x1 =
0.1708 k =
145
x0 =
0.1708
x1 =
0.1673 k =
146
x0 =
0.1673 x1 =
0.1707 k =
147
x0 =
0.1707 x1 =
0.1674 k =
148
x0 =
0.1674 x1 =
x1 =
0.1676
k =
154
x0 =
0.1676
x1 =
k =
158
x0 =
0.1678
x1 =
0.1702
k =
159
x0 =
0.1702
x1 = 0.1679
k =160
x0 = 0.1691 x1 = 0.1689 k = 250
x0 = 0.1689 x1 = 0.1691
k = 251
x0 = 0.1691
x1 = 0.1689
k = 252
x0 = 0.1689
x1 = 0.1691
k = 253
x0 = 0.1691
x1 = 0.1689
k = 254
x0 = 0.1689
x1 = 0.1691
k = 255
x0 = 0.1691
x1 = 0.1689
k = 256
x0 = 0.1689
x1 = 0.1691
k = 257
x0 = 0.1691
x1 = 0.1689
k = 258
x0 = 0.1689
x1 = 0.1691
k = 259
x0 = 0.1691
x1 = 0.1689
k = 260
x0 = 0.1689
x1 = 0.1691
k = 261
x0 = 0.1691
x1 = 0.1690
k = 262
x0 = 0.1690
x1 = 0.1691
四、数值结果
k 1 2 3 4 5 6 7 8 9 x 0.50000.75000.62500.68750.65630.64060.64840.64450.6426
k 10 11 12 13 14 15 16 17
x 0.6416 0.6411 0.6414 0.6412 0.6412 0.6412 0.6412 0.6412
f =
-5.4870e-006
k =
17
五、计算结果的分析
这个方程运用二分法求得四位有效数字的解为x=0.6414。
图一
六、计算中出现的问题,解决方法及体会
1.刚开始设计程序,出现死循环,无法得出结果,要经过不断调试。
2.不熟悉Matlab的用法,导致解决问题效率下降,要学好相关基础知识。
3.由图一可以看出二分在区间很大时,收敛接近根的速度很快,而当区间较小时,靠近要求精度根速度变得相当缓慢,计算量也很大,对满足高精度的要求比较困难。
我们想象此法可以和别的后半段接近根较快的方法结合,先求出适用于别的方法的初值,应当会优化算法。
教
师
评
语指导教师:
2013 年 3 月 12 日。