一元3次以上非线性方程的求解基于ExcelVBA求解法基于.pptVIP

基于ExceVBA的Newton迭代法基于ExceVBA的Newton迭代法*基于ExceVBA的Newton迭代法*一元3次以上非线性方程的求解

——基于ExcelVBA求解法*基于ExceVBA的Newton迭代法*牛顿法(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思考:曲线的斜率如何得到?*基于ExceVBA的Newton迭代法*牛顿法(Newton-Raphson)法一般步骤以一元3次非线性方程，如y=Ax3+Bx2+Cx+D为例：步骤1:选初值，如x0；步骤2:将初值分别代入原方程,求出原函数的值、一阶导数值,如:步骤3:由迭代公式计算出x1值；步骤4:计算两次计算结果的相对误差的绝对值(也有比较绝对误差的绝对值,但较少采用),若误差在允许范围内(通常称满足迭代精度,如手工计算取10?3,若采用计算机计算,一般取10?4),则可认为所得结果即为方程的解,否则将xi代替xi+1,重复步骤1~4,直到满足迭代精度为止。*基于ExceVBA的Newton迭代法*牛顿法(Newton-Raphson)法计算机程序框图*基于ExceVBA的Newton迭代法*应用牛顿迭代法计算原理，采用VBA自编迭代函数求解步骤：①打开Excel?“开发工具”?“VisualBasic”编辑器?插入?模块?过程,输入函数名“Newton3”,如右下图所示②在编辑窗口编写VBA代码,如下图所示*基于ExceVBA的Newton迭代法*应用牛顿迭代法计算原理，采用VBA自编迭代函数求解步骤：②在编辑窗口编写VBA代码,如下图所示PublicFunctionNewton3(V0,ESP,A,B,C,D)AsDouble‘V0-初值,ESP-迭代精度,A,B,C,D-分别代表方程的系数值Dimy0,Y1,X,x1AsDouble‘定义y0,Y1,X,X1四个实数型变量,分别代表方程的值,一阶导数值,方程的解x1=V0 ‘将方程的初值V0赋予变量x1Do ‘程序循环体开始语句X=x1 ‘将初值x1赋予中间变量xy0=A*X^3+B*X^2+C*X+D‘计算原函数的值Y1=3*A*X^2+2*B*X+C‘计算函数一阶导数值x1=X-y0/Y1‘由假设值通过牛顿迭代公式得到计算值LoopUntilAbs((x1-X)/X)ESP判断假设值与计算值之间的误差Newton3=x1‘当满足计算精度时,程序退出循环体,将计算结果赋予函数EndFunction*基于ExceVBA的Newton迭代法*应用牛顿迭代法计算原理，采用VBA自编迭代函数求解步骤：③在Excel表格中分别输入方程初值、迭代精度及方程中各系数之值;④在如下图单元格F4中输入：=newton3(E3,F3,G3,H3,I3,J3)，按“Enter”即得方程的解为：0.2749化1141聚物1321化1123单元一实训N文件,在VisualBasic编辑器界面,”工具宏逐语句(运行子程序”或”F5)F8”,单步运行,就可仔细观看计算按步的迭代过程.*单元一实训N文件,在VisualBasic编辑器界面,”工具宏逐语句(运行子程序”或”F5)F8”,单步运行,就可仔细观看计算按步的迭代过程.*对应的Excel文件单元一实训中有SubTEST3()子程度,可通过逐语句(按F8)说明牛顿迭代法.*基于ExceVBA的Newton迭代法基于ExceVBA的Newton迭代法化1141聚物1321化1123单元一实训N文件,在VisualBasic编辑器界面,”工具宏逐语句(运行子程序”或”F5)F8”,单步运行,就可仔细观看计算按步的迭代过程.*单元一实训N文件,在VisualBasic编辑器界面,”工具宏逐语句(运行