黄金分割法及MATLAB函数GSM（免费下载）

https://blog.csdn.net/luzhanbo207/article/details/121559905

syms X Y a b c e m x y;
%X为自变量，Y为函数，a为单峰区间左端点，b为单峰区间右端点,e为满足要求的区间长度,x为函数极值点，y为函数极值
c=0.618;%c为黄金分割系数
e=0.001;%e为满足要求的区间长度
Y=@(X) (X-1) *(X-3) *(X-9) *(X-12);
a=10.3;b=11.5;%函数单峰区间对应的左右端点
m=b-a;
x1=b-(c*(b-a));
f1=Y(x1);
x2=a+(c*(b-a));
f2=Y(x2);
while(m >= e)
    if(f1 >= f2)%区间压缩为[x1,b],将x1变为a，x2变为x1
        a=x1;
        x1=x2;
        f1=f2;
        x2=a+(c*(b-a));%重新计算x2
        f2=Y(x2);%重新计算f2
        m=b-a;%重新计算区间长度
    else%区间压缩为[a,x2],将x2变为b，x1变为x2
        b=x2;
        x2=x1;
        f2=f1;
        x1=b-(c*(b-a));%重新计算x1
        f1=Y(x1);%重新计算f1
        m=b-a;%重新计算区间长度
    end
end
x=(a+b)/2;
y=Y(x);
x
x =
   10.7528

 y
y =
 -165.2937

结果如下图所示：

matlab函数百度云链接：
https://pan.baidu.com/s/1Ib4dOv7wVhJTAnQaJJb8jA
提取码：0jc5