非线性方程的弦割法,抛物线法.docVIP

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