一元3次以上非线性方程的数学求解法逐步扫描法迈步法二分法牛.pptVIP

一元非线性方程的数学求解法一元非线性方程的数学求解法*一元非线性方程的数学求解法*一元3次以上非线性方程的数学求解法逐步扫描法(迈步法)二分法牛顿法弦截法*一元非线性方程的数学求解法*一、概述思考:以前所学知识的哪方面遇到过求解一元三次以上非线性问题?ax3+bx2+cx+d=0,如何求解?因式分解!?(x±A)(x±B)(x±C)=0?x=±A、±B、±C☆3x3-4x2-43x+84=0?(x+4)(x-3)(x-2.3333)=0?x=-4,3,2.333☆真实气体p-V-T之间的关系,如:vanderWaalsEq☆已知供热量、气体的初温T1、恒压摩尔热容cp=f(T)=a+bT+cT2，求出口温度T2*一元非线性方程的数学求解法*二、逐步扫描法(迈步法)☆函数f(x)在区间(a,b)连续,且f(a)与f(b)异号,则区间(a,b)内必定至少有一个实根;☆函数f(x)在区间(a,b)连续并单调(上升或下降),则在区间(a,b)必定只有一个实根;☆选一步长h,计算函数f(a),f(a+h),f(a)与f(a+h)乘积,若乘积大于零,说明该区间内无实根?计算f(a+2h),f(a+h)与f(a+2h)乘积,直至两函数值的乘积小于或等于零,即直至相邻两个函数值异号,此时可把此区间的始点作为方程式根的近似值。这个方法叫迈步法或逐步扫描法。返回*一元非线性方程的数学求解法*三、二分法设方程f(x)=0通过逐步扫描法已知有根区间为(x1,x2)①取x1与x2的中点x0，即:x0=(x1+x2)/2②检查f(x0)与f(x1)的符号是否相同☆如为同号,根必在x0与x2之间,如图(A)。令:x1=x0☆如为异号,则根必在x1与x0之间,如图(B)。令:x2=x0③再取新区间(x1,x2)的中点,重复以上步骤,直至x1与x2之间的距离小于某指定值?为止。

只要区间(x1,x2)内有根,此法定能求出结果,但收敛速度较慢;

此法只能得到一个实根,需增加一个终值B,当有根区间的终点x2小于终值B,再施行迈步法求下一个实根的近似值,进而求出精确值,如此反复,直至达到或超过终值B为止。返回*一元非线性方程的数学求解法*四、牛顿法(Newton-Raphson)法假设方程f(x)=0有一个实根x*①取一初值x0,过x0作垂线交曲线f(x)于点P0,过P0点作曲线切线交x轴于x1点,显然x1点较x0点更接近于根x*;②若?(x1－x0)/x0?＜?,则方程根x*=x1,否则按上述同样方法过x1作x轴垂线交于曲线f(x)于点P1,过P1点作曲线f(x)的切线并于x轴相交于x2……，直到?(xk+1－xk)/xk?＜?(迭代精度,一般手工计算取10-3,计算机计算取10-4)为止,方程的根为:x*=xk+1思考:曲线的斜率如何得到?返回*一元非线性方程的数学求解法*四、牛顿法(Newton-Raphson)法以一元3次非线性方程，如y=Ax3+Bx2+Cx+D为例：步骤1:选初值,如x0,如何选择需根据具体情况决定;步骤2:将初值分别代入原方程,求出原函数的值、一阶导数值,如:步骤3:由迭代公式计算出x1值；步骤4:计算两次计算结果的相对误差的绝对值(也有计算绝对误差的绝对值的情况,但很少),若误差在允许范围内(通常称满足迭代精度,如手工计算取10?3,计算机计算取10?4),则可认为得结果即为方程的解,计算结束;否则将xi代替xi+1,重复步骤1~4,直到满足迭代精度为止。*一元非线性方程的数学求解法*☆采用Newton迭代法求解一元3次非线性方程举例例1-8常压气相反应体系某一组分平衡分压p(单位atm)符合以下方程:4p3－1.640p2+1.640p－0.410=0,试求其分压的数值?思考:p的大致范围?解:步骤1：令： y=4p3－1.640p2+1.640p－0.410 (1)则: y?=12p2－3.280p+1.640 (2)步骤2:确定p的初值,取为p0=0.3步骤3:将p0=0.3代入式(1),(2)计算得:y0=0.04240、y0?=1.73600步骤4:步骤5:计算p0、p1之间相对误差?1步骤6：判断?1,满足计算结果,否则重复步骤3~6,直到满足要求。*一元非线性方程的数学求解法*☆采用Newton迭代法求解一元3次非线性方程举例例1-8的求解步骤用下表表示例1-8Newton迭代法计算数据汇总表返回迭代次序p初值y计算值y′计算值p计算值迭代误差计算判断10.30