有机合成工职业技能培训有机合成工职业技能培训单元1任务4求解一元非线性方程.pptVIP

单元1 数据处理 * 单元1 数据处理 * 任务4 求解一元非线性方程 求解一元非线性方程 求根方法 手工求解 Excel求解 自编软件求解 逐步扫描法(迈步法) 二分法 牛顿法 弦截法 * 单元1 数据处理 * 任务4 求解一元非线性方程 思考:以前所学知识的哪方面遇到过求解一元三次以上非线性问题? ☆3x3-4x2-43x+84 = 0 ☆真实气体p-V-T之间的关系 如:van der Waals Eq (p+a/V2)(V-b)=RT?V3-(b+RT/p)V2+aV/p-ab/p=0 ☆已知供热量、气体的恒压摩尔热容cp=f(T)=a+bT+cT2，求出口温度T2 * 单元1 数据处理 * 任务4 求解一元非线性方程 一、一元非线性方程求根方法 (一)逐步扫描法(迈步法) ☆函数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) 乘积,直至两函数值的乘积小于或等于零,即直至相邻两个函数值异号,此时可把此区间的始点作为方程式根的近似值。 这个方法叫迈步法或逐步扫描法。 返回 * 单元1 数据处理 * 任务4 求解一元非线性方程 一、一元非线性方程求根方法 (二)二分法 设方程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为止。 返回 * 单元1 数据处理 * 任务4求解一元非线性方程 一、一元非线性方程求根方法 (三)牛顿法 (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?＜?为止,方程的根为:x*=xk+1 思考: 曲线的斜率如何得到? 返回 * 单元1 数据处理 * 任务4 求解一元非线性方程 一、一元非线性方程求根方法 (四)弦截法 思考: f ?(x)很难求得? 1)若已知非线性函数f(x)的根区间(x0, x1),过x0、x1作垂线交函数f(x)于P0、P1点,连接P0、P1交x轴于x2点 ; 2) ①若f(x2)=0,则方程的解为x*= x2; ② 若f(x2)?f(x1)＞0,如图(A),则用x2代替x1; ③ 若f(x2)?f(x1)＜0,如图(B),则用x2代替x0; 3)用新得到的两个点用以上方法继续迭代,直到相邻两次值满足?(xk+1－xk)/xk)?＜?为止,弦截法的迭代格式为 返回 * 单元1 数据处理 * 任务4 求解一元非线性方程 二、手工求解一元非线性方程 例1-8 常压气相反应体系某一组分平衡分压p(单位atm)符合以下方程:4p3－1.640p2+1.640p－0.410=0,试求其分压的数值? 思考: p的大致范围? 解:以Newton迭代法为例 步骤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)计算得：y=0.04240、y?=1.73600 步骤4: 步骤5:计算p0、p1之间相对误差?1 步骤6：判断?1,满足计算结果,否则重复步骤3~6,直到满足要求。 * 单元1 数据处理 * 任务4 求解一元非线性方程 二、手工求解一元非线性方程 例1-8 常压气相反应体系某一组分平衡分压p(单位atm)符合以下方程:4p3－1.640p2+1.640p－0.410=0,