MATLAB二次插值应用程序(共1页).doc

精选优质文档-倾情为你奉上function ecczsyms x;%定义f=25.*(5-4*x+x2)+38.5*(1+x2); %目标函数a1=0.5;a2=1.5;a3=1; %插值区间 k=0;while (abs(a3-a1)1.0e-7) %精度 f1=subs(f,x,a1);f2=subs(f,x,a2);f3=subs(f,x,a3);C1=(f3-f1)/(a3-a1);C2=(f2-f1)/(a2-a1)-C1)/(a2-a3);ap=0.5*(a1+a3-C1/C2);fp=subs(f,x,ap); if apa2; if f2=fp a1=a2; f1=f2; a2=ap; f2=fp; else a3=ap;