用抛物线与横轴的交点来逼近函数f(x)的根.docxVIP

用抛物线与横轴的交点来逼近函数 f （x）的根 《 MATLAB程序设计实践》课程考核 , 、编程实现以下科学计算算法～并举一例应用之。 , 参考书籍《精通 ,,,,, 科 学计算》～王正林等著～电子工业出版社～ ,,,, 年,?抛物线法非线性方程求解 ? 抛物线法 弦截法其实是用不断缩短的直线来近似函数 f,x, 的～而抛物线法则是采用三个 点来近似函数 f,x, ～并且用抛物线与横轴的交点来逼近函数 f,x, 的根。它的算法 过程如下 : ,1, 选定初始值 x0～x1～ x2～并计算 f,x0, ～f,x1, ～f,x2, 和以下差 分: f[x2.x1]= f[x1,x0]= f[x2,x1,x0]= 一般取 x0=a～x1=b～ax2b。注意不要使三点共线。 ,2, 用牛顿插值法对三点 ,x0 ～f,x0,, ～,x1 ～f,x1,, 。 ,x2 ～f ,x2,, 进行插值得一条抛物线～它有两个根 : x3=x2+ 其中 A=f(x2),C=f[x2,x1,x0] ～ B=f[x2,x1]+f[x2,x1,x0](x2-x1). 两个根中只取靠近 x2 的那个根～即号取与 , 同号～ 即 x3=x2-- ,？ , 用 x1～x2～ x3 代替 ,, ～ ,, ～,, ～重复以上步骤～并有以下递推公式 : xn+1=xn- 其中 =f(xn),Cn=f[xn,xn-1,xn-2] ～ =f[xn,xn-1]+f[xn,xn-1,xn-2](xn-xn-1) 。 ,,, 进行精度控制。 在 MATLAB中编程实现的抛物线法的函数为 :Parabola 。 功能: 用抛物线法求函 数在某个区间上的一个零点。 调用格式 :root=Parabola,f ～a～ b～x～ eps,。 其 中～f 为函数名 : a为区间左端点 : b为区间右端点 : x 为初始迭代点 : eps 为根的精度 : root 为求出的函数零点。 抛物线法的 MATLAB程序代码如下 : function root=Parabola,f ～a～ b～x～ eps, %抛物线法求函数 f 在区间【 a～ b】上的一个零点 %函数名 :f %区间左端点 :a %区间右端点 :b %初始迭代点 :x %根的精度 :eps %求出的函数零点 :root if(nargin==4) eps=1.0e-4; end f1=subs(sym(f),findsym(sym(f)),a); f2=subs(sym(f),findsym(sym(f)),b); if(f1==0) root=a; end if(f2==0) root=b; end if(f1*f20) disp( ‘两端点函数值乘积大于 0! ’); return; else tol=1; fa=subs(sym(f),findsym(sym(f)),a); fb=subs(sym(f),findsym(sym(f)),a); fx=subs(sym(f),findsym(sym(f)),x); d1=(fb-fa)/(b-a); d2=(fx-fb)/(x-b); d3=(f2-f1)/(x-a); B=d2+d3*(x-b); root=x-2*fx/(B+sign(B)*sqrt(B^2-4*fx*d3)); t=zeros(3); t(1)=a; t(2)=b; t(3)=x; while(toleps) t(1)=t(2); % 保存 3 个点 t(2)=t(3); t(3)=root; f1=subs(sym(f),findsym(sym(f)),t(1)); % 计算 3 个点的函数值 f2=subs(sym(f),findsym(sym(f)),t(2)); f3=subs(sym(f),findsym(sym(f)),t(3)); d1=(f2-f1)/(t(2)-t(1)); % 计算 3 个 差分 d2=(f3-f2)/(t(3)-t(2)); d3=(d2-d1)/(t(3)-t(1)); B=d2+d3*(t(3)-t(2)); % 计算算法中的 B root=t(3)- 2*f3/(B+sign(B)*sqrt(B^2-4*f3*d3));