MATLAB之用黄金分割法求函数极小点

ZeroZeroOneTwo

已于 2023-03-21 21:26:59 修改

阅读量1.5k

点赞数 4

文章标签： matlab

于 2023-03-21 20:26:51 首次发布

本文链接：https://blog.csdn.net/zerozeroonetwo/article/details/129696756

版权

函数： $f(X)=(x-3)\sqrt{x}$ ，初始点 $X^{(0)}=0$ ，精度 $\varepsilon =0.1$

主程序：

clc
clear
f=@(x) (x-3)*sqrt(x);
fplot(f)
a=-100;
b=100;
[fx,x]=huangjin(f,[a b],100)  % （a，b）为搜索区间
function [result,x]=huangjin(f,x0,n)
tol=0.1;%精度
x1=x0(1);
x2=x0(2);
i=1;
while i < n
    %取中间值
    a=x1+0.382*(x2-x1);
    b=x1+0.618*(x2-x1);
    fa=f(a);
    fb=f(b);
    % 判断fa  fb大小，缩小区间
    if fa < fb
        x2=b;
    else
        x1=a;
    end
    if  abs(x1-x2) < tol
        result=f((x1+x2)/2);
        x=(x1+x2)/2;
        break;
    end
    i=i+1;
end
end

计算结果：