用二次插值法求解函数极小值

最新推荐文章于 2024-05-10 17:23:50 发布

ZeroZeroOneTwo

最新推荐文章于 2024-05-10 17:23:50 发布

阅读量837

点赞数 1

文章标签： matlab

本文链接：https://blog.csdn.net/ZeroZeroOneTwo/article/details/129697040

版权

函数： $f(X)=x+20/x$ 的极小点。要求 $X\geq 0$ ，初始点 $X^{^{(0)}}=0$ ，精度 $\varepsilon =0.1$

主程序：

clc
clear
f=@(x) x+ (20./x);
ezplot(f,[0 100])
[x,fx]=Min_erci(f,[1 21],100)

子程序：

function [x,result]=Min_erci(f,x0,k) %x0为初始区间端点，k为迭代次数
 
x1=x0(1);
x3=x0(2);
x2=(x1+x3)/2;
 
n=1;
while n < k
    % 确定抛物线的系数
    f1=f(x1);
    f2=f(x2);
    f3=f(x3);
    A=[x1^2    x1   1;
       x2^2    x2   1;
       x3^2    x3   1;];
   b=[f1;f2;f3];
   XS=A\b;  %求出抛物线系数a b c 存放在xs中
   xp=-XS(2)/(2*XS(1)); %二次多项式的极值点在x=-b/2a
   fp=f(xp); %求出该点函数值
   if abs(xp-x2) < 1e-8  %该点满足极值点条件
       x=xp;        %输出极值点
       result=f(x); %输出函数值
       return;
   end
   if fp > f2       %判断新的迭代区间
       if xp < x2
           x1=xp;
       else
           x3=xp;
       end
   else
       if xp < x2
           x3=x2;
           x2=xp;
       else
           x1=x2;
           x2=xp;
       end
   end
    n=n+1;   %迭代次数+1
end
 
if n == k        %如果超出迭代次数
    x=[];        %输出空    
    result=[];
    disp('超过迭代次数');
end
end

求解结果：