()非线性方程的弦割法,抛物线法..docxVIP

精品文档 西京学院数学软件实验任务书 课程名称 数学软件实验 班级 数 0901 学号 0912020119 姓名 王震 实验课题 非线性方程的弦割法，抛物线法 实验目的 熟悉非线性方程的弦割法，抛物线法 运用 Matlab/C/C++/Java/Maple/Mathematica 等其中 实验要求 一种语言完成 非线性方程的弦割法 实验内容 非线性方程的抛物线法 成绩 教师 精品文档 精品文档 【实验课题】 非线性方程的弦割法，抛物线法 【实验目的】 熟悉非线性方程的弦割法，抛物线法 【实验内容】 1、弦割法 弦割法是函数逼近法的一种，基本思想是用用区间 [ xk 1， xk ] 上的割线近似代替目标函数的导函数的曲线。并用割 线与横轴交点的横坐标作为方程式的根的近似。 对于牛顿迭代式 xk 1 xk f xk k 0,1,2,... f xk 由于上式存在微分，为了避免计算导数值，我们用差商来代替导数。即： xk 1 f ( xk ) ( xk xk 1 ) xk f ( xk ) f ( xk 1 ) 算法过程： 1、定义函数 f (x) ，输入控制精度 和迭代初始值 x0 , x1 ，并求 出 f ( x0 ), f ( x1 ) 2、进行循环迭代，若相邻的两个数值之差的绝对值小于所给 定的 , 则跳出循环，输出 x k 2、抛物线法 精品文档 精品文档 抛物线法是在函数的最优点附近取三个构造点， 然后用这 三个点构造一条抛物线， 把这条抛物线的极值点作为函数的极 值点的近似。每次构造一条抛物线后，抛物线的极值点就可作 为一个新的构造点， 新的构造点与原来的三个构造点经过某种 算法，得到下一步抛物线逼近的三个构造点，这就是抛物线法 的算法过程。 算法过程： 1、选定初始值 x0 , x1, x2 ，并计算 f ( x0 ), f (x1), f ( x2 ) 和差分 f ( x2 f ( x2 ) f ( x1 ) , x1 ) x1 x2 f ( x1 f ( x1 ) f ( x0 ) , x0 ) x0 x1 f (x2 , x1 ) f ( x1 , x0 ) f ( x2 , x1, x0 ) x0 x2 注意， x0, x1, x2 三点不能共线。 2、用牛顿插值法对三点 ( x0 , f ( x0 )) 、 ( x1, f ( x1 )) 、 ( x2 , f (x2 )) 进 行插值得到一条抛物线，存在两个根： x3 BB2 4 AC x2 2C 其中， A f ( x2 ) ， C f ( x2 , x1, x0 ) B f ( x2 , x1 ) f (x2 , x1 , x0 )( x2 x1 ) 两个根中取其中靠近 x2 的那个根。 3、重复以上步骤，并控制精度。 精品文档 精品文档 【程序】 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% function [x1,n]=gexian() x0=input( 请输入初始值 x0=); x1=input( 请输入初始值 x1=); err=input( 请输入误差极限 err=); N=input( 请输入最大迭代次数 N=); x2=x1-fun(x1)*(x1-x0)/(fun(x1)-fun(x0)); n=0; while abs(x1-x0)err x0=x1; x1=x2; x2=x1-fun(x1)*(x1-x0)/(fun(x1)-fun(x0)); n=n+1; if nN disp( 迭代次数过大！ ); return end end end %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% function root=paowuxian() x0=input( 请输入初值 x0=); x1=input( 请输入初值 x1=); x2=input( 请输入初值 x2=); N=input( 请输入最大迭代次数 N=); f1=fun(x0); f2=fun(x1); if f1==0 root=x0; end if f2==0 root=x1; end if f1*f20 disp( 两点之积大于 0！ ); 精品文档 精品文档 return else err=input( 请输入误差极限 err=); f1=fun(x0); f2=fun(x1); f3=fun(x2); d1=(f2-f1)/(x1-x0); d2=(f3-f2)/(x2-x1); d3=(d2-d1)/(x2-x0); B=d2+d3*(x2-x1); root=x2-2*f3/(B+sign(B)*sqrt(B^2-4*f3*d3)); t(1