实验四 非线性方程的求根

(2) 如果将原方程化为等价方程 x = 仍取初值
x0 = 0
3
3
x+1 2
x1 =
x0 + 1 = 2
3
3
1 ≈ 0 .7937 2
1 .7937 ≈ 0 .9644 2
x2 =
依此类推,得 依此类推, x3 = 0.9940 同样的方程 x4 = 0.9990 不同的迭代格式 x5 = 0.9998 x6 = 1.0000 有不同的结果 x7 = 1.0000 已经收敛, 已经收敛,故原方程的解为 x = 1.0000
用一般迭代法求方程x lnx＝ 在区间（ 例：用一般迭代法求方程x-lnx＝2在区间（2，∞） 内的根，要求|x 内的根，要求|xk-xk-1|/|xk|<=10-8 f(x)=x-lnx解：令f(x)=x-lnx-2 f(2)<0,f(4)>0,故方程在 2,4） 故方程在（ f(2)<0,f(4)>0,故方程在（2,4）内至少有一个根
y p1 p0
y=x y=g(x)
y p0
y=x
p1 x x0 y y=g(x) x1 x* y=x y y=g(x) p0 x0 x* x1
y=g(x)
x
y=x
p0 p1 x x1 x0 x* x0 x* x1 p1 x
f(x)=0化为等价方程x=g(x)的方式是不唯 f(x)=0化为等价方程x=g(x)的方式是不唯 化为等价方程x=g(x) 一的,有的收敛,有的发散. 一的,有的收敛,有的发散. example： For example：2x3-x-1=0
1 3 ( x + 1 ) = g ( x ) 进行迭代，则在 进行迭代，则在(1, 3

Matlab软件实习论文非线性方程求根系别信息与计算科学专业信息与计算科学学号姓名指导教师2008年8 月10 日非线性方程求根摘要随着科学技术，生产力经济的发展，在科学与工程计算中存在着大量方程求根问题，例如贷款购房问题，工厂的最佳订货问题等都需要求解一类非线性方程的根，而本文就针对这些求根问题提出了解决方案，本文利用牛顿迭代法来结决方程的求根问题．首先根据实际问题列出数学模型，确定变量，给出各个条件及相关函数；然后对建立的模型进行具体分析和研究，选择合适的求解方法；编写函数的程序，用计算机求出方程的解，通过所求解分析具体情况．关键词：非线性方程，牛顿迭代法，Matlab目录摘要 (I)1 绪论 (1)1.1非线性方程求根的背景 (1)1.2非线性方程求根的目的： (1)1.3非线性方程求根的内容： (1)2 牛顿迭代法的实现及应用 (3)2.1N EWTON迭代法具体例子的实现 (3)2.2应用牛顿法解决购房贷款利率问题 (4)2.3应用牛顿迭代法计算最佳订货量 (6)结论 (8)参考文献 (9)1 绪 论1.1 非线性方程求根的背景随着社会的进步，科学技术的快速发展，各种工程等也如雨后春笋一般破土而出，对我们的日常生活产生了巨大的影响如天气预报、石油的勘探、地质灾害的预报等．牛顿迭代法是牛顿在17世纪提出的一种求解方程()0f x =．多数方程不存在求根公式，从而求精确根非常困难，甚至不可能，从而寻找方程的近似根就显得特别重要．而在各种科学和工程计算中往往要用到非线性方程组的求解，而牛顿法又是最基础的迭代法，在各种计算力学、控制工程等领域中发挥了不可代替的作用．而在数值计算中，非线性方程组的求解同样具有重要意义．随着计算机技术的成熟和高速发展，对于非线性方程求根问题出现了大量的数学软件（如MATLAB,Matheamatica,Maple,SAS,SPSSD 等），计算机已经成为工程师应用数学解决工程问题的主要运算工具．同时，工程专业的学生对数学教育的需求重点正在从手工演绎和运算能力的培养转变到结合计算机软件进行建模、求解和论证能力的培养[]1．本文采用Matlab 作为软件平台，介绍了非线性方程求根的内容．1.2 非线性方程求根的目的为了推动科学的进步，能够很简便的完成各种工程计算，非线性方程组的求解方法以其独有的方法解决了各种计算，为今天以及将来的应用打下了坚实的基础．非线性方程组的求解正是为了这个目的才广泛被人们应用，此文也将给出非线性方程组求解的实际应用．1.3 非线性方程求根的内容解非线性方程()0f x =的主要算法是迭代法，如fsolve 、二分法、牛顿迭代法等．迭代法是从已知的解的初始近似值0 x (简称初值)开始，利用某种迭代格式( )x g x =求得一近似值序列121,,,,,k k x x x x +逐步逼近于所求的解α(称为不动点)．这一方法是否成功取决于三个因素，首先 ( )x g x =应与()0f x =同解，其次初值0 x 的选取是否合适，一般要与真解靠近，最后也是最关键的是迭代序列是否收敛，为了保证收敛性，在真解附近应有'| ()|1g x <否则迭代序列可能发散．最基本的迭代法是Newton 迭代法，其迭代格式为1'()()k k k k f x x x f x +=-. 从几何上说1k x +为用 ()f x 在 k x 出切线代替 ()f x 求得的解，所以也称为切线法，当初值0 x 与真值α足够靠近，Newton 迭代法收敛．对于单根，Newton 法收敛速度很快；对于重根，收敛较慢．牛顿迭代法的大概算法为：给定初始值0x ，ε为根的容许误差，η为()f x 的容许误差，N 为迭代次数的容许值.① 如果'0()0f x =或迭代次数大于N ，则算法失败，结束；否则执行② ② 计算010'0()()f x x x f x =-③ 若10x x ε-<或1()f x η<，则输出1x ，程序结束；否则执行④ ④ 令01x x =，转向①下面给出了Newton 迭代法的计算程序． function x=newton(fname,dfname,x0,e) %用途：Newton 迭代法解非线性方程f(x)=0%格式：x=nanewton( fname,dfname,x0,e) x 返回数值解， %fname 和dfname 分别表示f(x)及其导函数 %f'(x),x0为迭代初值，e 精度要求（默认为1e-4) if nargin<4,e=1e-4; %精度默认为1e-4 endx=x0;x0=x+2*e; %使while 成立，进入while 后x0得到赋值 while abs(x0-x)>e x0=x;x=x0-feval(fname,x0)/feval(dfname,x0); end2 牛顿迭代法的实现及应用2.1 Newton 迭代法具体例子的实现用Newton 迭代法解方程 32() -3-30 f x x x x =+= 在1.5附近的根．解：当2x >时，()0, ()0f x f x >>，即()f x 恒正，所以根在[0,2].我们先用图解法找初值，在用Newton 法程序newton.m 求解．fun= inline('x^3+x^2-3*x-3'); fplot(fun,[0,2]); grid on;图 2.1 ()f x 的函数图像由图可知方程有唯一正根在[1.6,1.8]之间，我们取初值1.5代入Newton 程序中． dfun=inline('3*x^2+2*x-3'); format long;newton(fun,dfun,1.5,1e-4); format short; ans =1.73205080756888而用Matlab 本身的函数fzero 求出来的结果为： format long;fzero(inline('x^3+x^2-3*x-3'),1.5); format short ans =1.73205080756888下面用牛顿迭代法解决一些实际问题 2.2 应用牛顿法解决购房贷款利率问题住房是居民消费的一个主要部分，大部分人选择银行按揭贷款，然后在若干年内逐月分期还款．如果你借了10万，还款额一定超过10万．设贷款总额为0x ,贷款期限为N 个月，采取逐月等额方式偿还本息．若k x 为第k 个月的欠款数，a 为月还款，r 为月利率．我们得到些列迭代关系式1 (1) - ,k k x r x a +=+那么1(1)k k x r x a -=+-22(1)(1)k r x r a a -=+-+-= 210(1)[1(1)(1)]k r x a r r -=+-+++++0(1)[(1)1]/k k r x a r r =+-+-, 由此可以得到月还款计算公式(1)(1)1N Nr x a r +=+-下面是《新民晚报》2000年3月30日第七版上的一则房产广告：不难算出，你向银行总共借了25.2万，30年内共要还51.96万，约为当初借款的两倍．这个案例中的贷款年利率是多少呢？我们根据a =0.1436, 0x =25.2, 360N =,由以上a 的求解公式得到：360360++=.25.2(1)-0.1436[(1)-1]0r r r我们令360360=++,( ) 25.2(1)-0.1436[(1)-1]f r r r r则次问题就转化成非线性方程求根的问题，令( )0,f r=求出r．我们先用Newton函数求解．在Matlab中输入如下程序：常识上，r应比当时活期存款月利率略高一些．我们用当时的活期存款月利率0.0198/2作为迭代初值，为了剔除0r=这个没有意义的根，我们对( )f r稍作变形：clear;fun=inline('25.2*(1+r)^360/0.1436-((1+r)^360-1)/r','r')fun =Inline function:fun(r) = 25.2*(1+r)^360/0.1436-((1+r)^360-1)/rdfun=inline('25.2*360*(1+r)^359/0.1436-(360*(1+r)^359*r-((1+r)^360-1))/(r^2)');r=newton(fun,dfun,0.0198/2,1e-4);R=12*r然后求得结果：R =0.0553于是得出年利率为 5.53%．下面我们用Matlab中的fzero函数检验一下：clear;fun=inline('25.2*(1+r)^360-((1+r)^360-1)/r*0.1436','r')fun =Inline function:fun(r) = 25.2*(1+r)^360-((1+r)^360-1)/r*0.1436r=fzero(fun,0.0198/2);R=12*rR =0.0553结果相同，可见牛顿迭代法的正确性． 2.3 应用牛顿迭代法计算最佳订货量汽车工厂为了保证生产的正常运作，配件供应一定要由保障．这些配件并不是在市场上随时可以买到的，所以往往要预先从配件供应商那里定货．由于配件供应商并不是生产单一产品，为你的定货必须要在流水线上作出调整，所以每次定货需要收取一定量的生产准备费．配件供应商的生产能力很大，开工后很快可以生产许多配件，但是你的汽车工厂并不是立即需要这么多，往往要在仓库里储存一段时间，为此你要付出储存费．如果订货量很小，必然需要频繁定货，造成生产准备费的增加；反之，若订货量很大，定货周期必然延长，生产准备费下降，但这样会造成储存费的增加．如何确定合适的订货量？实践中，这是一个相当复杂的问题，因为市场波动的影响是多方面的．我们先作一些必要的假设将问题简化．1) 汽车工厂对配件的日需量是恒定的，每日为r 件； 2) 所订配件按时一次性交货，生产准备费每次1 k 元； 3) 储存费按当日实际储存量计算，储存费每日每件2k 元； 4) 你的工厂不许缺货．设一次定货x 件，由于工厂不允许缺货，而为了节省存储费，交货日期应定为恰好用完时，所以定货日期/.T x r = (1)由于日需求量是恒定的，可以计算出第t 天的存储量为( )-, 0.q t x rt t T =<< (2)由于第t 天的储存费为2k q( t ),一个周期的总储存费为2201()()TTt t k q t k q t d -≈∑⎰. (3)根据（1），（2），（3）得到一个周期总费用212() 2x C x k k r=+,优化目标是使单位产品费用12()()2k k xC x f x x x r==+, 达到最小．由'()0f x =即122-02k k x r+=， 可直接解得x =这就是著名的经济批量定货公式． 当我们给出具体值时，非线性方程就可以求解了，由于具体的值不定，在此就不给出具体程序了．结论通过以上的论述我们可以知道计算机在现代生活中的应用已经如此普及，尤其是在数学计算当中，Matlab软件更是发挥了不可替代的作用．Matlab以其强大的功能，方便了当今数值计算，数学教程，及工程计算等众多领域．本文在以Matlab软件为平台的基础上，给出了非线性方程的一般解法，非线性方程的求解有二分法，牛顿迭代法，简单牛顿法，牛顿下山法，弦截法，抛物线法等．二分法的优点是算法简单，且总是收敛的，但由于二分法的收敛速度太慢，故一般不单独将其用于求根，只用其为根求得一个较好的近似[]2值．其他的求根方法各有优缺点，这里就不一一赘述．本文主要介绍了牛顿迭代法及其在现实生活中的应用．牛顿迭代法为平方收敛，故其收敛速度较快，但对初值的选取需要谨慎，如果初值选取错误，则可能导致方程迭代发散，最终不能求解出正确解．在计算一些对精度要求特别苛刻时，最好给出较高的精度输入及输出，防止因为精度问题导致误差过大，最终影响结果．牛顿迭代法可以应用于分形理论．分形理论是近二、三十年才发展起来的一门新的学科,其主要描述自然界和非线性系统中不光滑和不规则的几何形体.在地质、材料科学、物理学、计算机科学、艺术设计等方面有着十分广阔的应用前景. 利用牛顿迭代的数学原理和方法,实现牛顿迭代法的分形图形生成算法,将分形理论应用于计算机图形设计中,形．利用VC++6.0开发工具,实现了生成绚丽多彩的分形图[]3非线性方程的求根问题在计算机发展的基础上，被广泛应用于各种工程计算，大大方便了工程师们的计算过程，在现代的工业发展中发挥了重要作用．牛顿迭代法简便易学，为今后的学子们提供了更多的学习内容．参考文献[1] 胡良剑,孙晓君.Matlab数学实验[M]. 高等教育出版社,2006.[2] 李庆扬，王能超，易大义.数值分析（第4版）[M].清华大学出版社.施普林格出版社.2001.[3] 吴运兵，李勇.《西安科技大学学报》2005年03期.。

计算方法4方程求根的迭代法四方程求根的迭代法是一种用于解决非线性方程的数值方法。

在计算方法中，非线性方程指的是形如f(x)=0的方程，其中f(x)包含x的非线性项。

在实际中，非线性方程的求解是非常常见的问题，因此有很多不同的迭代法可以用于解决这些问题。

以牛顿迭代法为例，它是一种基于线性近似的迭代方法。

该方法的基本思想是将非线性方程转化为线性方程，通过不断迭代来逼近方程的根。

具体而言，牛顿迭代法的步骤如下：1.选择初始估计值x0作为方程的根，并计算f(x0)的值。

2.计算f(x)的导数f'(x)，并计算方程的线性近似式x-x0=-f(x0)/f'(x0)。

3.计算下一个近似值x1，即x1=x0-f(x0)/f'(x0)。

4.判断，x1-x0，是否小于给定的收敛条件，如果是则停止迭代，否则转到步骤55.将x1作为新的近似值x0，转到步骤2牛顿迭代法具有快速收敛的特点，尤其适用于具有单根的方程。

然而，该方法也存在一些限制，如在计算f'(x)时需要知道方程的导数，当方程的导数不易计算时，该方法可能不适用。

除了牛顿迭代法，还有其他一些常用的四方程迭代方法，如割线法、弦截法等。

每种方法都有其特点和适用范围，选择合适的方法对于求根问题的解决至关重要。

总结起来，四方程求根的迭代法是一种用于解决非线性方程的数值方法。

牛顿迭代法是其中一种常用的方法，通过不断迭代来逼近方程的根。

根据方程的特点和计算条件，选择合适的迭代方法是解决求根问题的关键。

希望以上的介绍可以帮助您更好地理解和应用这一方法。

Print["第", k, "次迭代x=", xx];
While[Abs[x0 - xx] > e,
If[k > max,
x0 = xx; Print["达到最大迭代次数,可能为死循环!"],
x0 = xx;
xx = g[x0];
k = k + 1;
Print["第", k, "次迭代x= ", xx];
x0为迭代初值,e为允许误差,k为计数器*)
2、（1）*输入数据,max为最大迭代次数,g[x]为迭代格式,
x0为迭代初值,e为允许误差,k为计数器*)
Clear[x];
max = 10;
g[x_] := (x^2 + x + 1)^(1/3);
x0 = 1.5;
e = 10^(-3);
k = 1;
(*利用迭代公式形成迭代序列*)
xx = g[x0];
第1次二分区间得:[-2.5,-2],区间中点为:-2.25, f[x]=0.0427101
第2次二分区间得:[-2.25,-2],区间中点为:-2.125, f[x]=-0.0295192
第3次二分区间得:[-2.25,-2.125],区间中点为:-2.1875, f[x]=0.00500699
第4次二分区间得:[-2.1875,-2.125],区间中点为:-2.15625, f[x]=-0.0126648
第5次二分区间得:[-2.1875,-2.15625],区间中点为:-2.17188, f[x]=-0.00392973
第6次二分区间得:[-2.1875,-2.17188],区间中点为:-2.17969, f[x]=0.000513602