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

《MATLAB 程序设计实践》课程考核 １、编程实现以下科学计算算法，并举一例应用之。（参考书籍《精通ＭＡＬＡＢ科学计算》，王正林等著，电子工业出版社，２００９年）“抛物线法非线性方程求解” 抛物线法 弦截法其实是用不断缩短的直线来近似函数f（x）的，而抛物线法则是采用三个点来近似函数f（x），并且用抛物线与横轴的交点来逼近函数f（x）的根。它的算法过程如下： 选定初始值x0，x1，x2，并计算f（x0），f（x1），f（x2）和以下差分： f[x2.x1]= f[x1,x0]= f[x2,x1,x0]= 一般取x0=a，x1=b，ax2b。注意不要使三点共线。 用牛顿插值法对三点（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)); tol=abs(root-t(3)); end end 例子：抛物线法求解非线性方程求方程lgx+=2在区间【1,4】上的一根。 流程图 实验2，编程解决以下科学计算问题。 用ode45，ode23和ode113解下列微分方程： y=x-y，y(0)=1,