用进退法求函数的单峰区间

ZeroZeroOneTwo

已于 2023-03-21 20:30:33 修改

阅读量562

点赞数

文章标签：机器学习人工智能

于 2023-03-21 20:23:20 首次发布

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

版权

函数： $f(X)=x^{2}-6x+9$ ，初始点 $X^{(0)}=0$

定义函数：

function f = func(x)
    f = x^2-6*x+9;

进退法程序：

function [JT_a, JT_b, JT_c] = JinTuiFa(x0,a,d)  %x0为初值，a为初始步长，d为方向
    %计算x1
    x1 = x0 + a * d;
    %前进后退标志符，n=1意味着上一次是前进
    n=1;
    %判断前进还是后退运算
    if func(x1) > func(x0)
        JT_c = x1;
        [JT_a, JT_b]=Back(x0,x1,a,n,d);
    else
        JT_a = x0;
        [JT_b, JT_c]=Front(x0,x1,a,n,d);
    end

后退运算子函数：

function [B_a, B_b] = Back(x0,x1,a,n,d)
    if n == 1
        n = 0;
    else
        a = 2 * a;
    end
    x2 = x0 - a * d;
    if func(x2) < func(x0)
        x1 = x0;
        x0 = x2;
        [B_a, B_b] = Back(x0,x1,a,n,d);     %递归
    else
        B_a = x2;
        B_b = x1;
    end

前进运算子函数：

function [F_a, F_b] = Front(x0,x1,a,n,d)
    if n == 0
        n = 1;
    else
        a = 2 * a;
    end
    x2 = x1 + a * d;
    if func(x2) < func(x1)
        x0 = x1;
        x1 = x2;
        [F_a, F_b] = Front(x0,x1,a,n,d);
    else
        F_a = x0;
        F_b = x2;
    End

主程序：

%设置函数起始点
x0 = 2;
%调用进退法计算
[JT_a, JT_b, JT_c] = JinTuiFa(x0,0.1,1);      %x0为初值，a为初始步长，d=1为方向
fprintf('进退法求得单峰区间为：\n')
fprintf('左端点：%f\n',JT_a)
fprintf('右端点：%f\n',JT_c)

结果：