如何求分段函数断点最优化问题

以下还是不太成熟的练习。

首先通过gpt生成一个算例，如下：

% 定义分段函数
piecewise_function = @(x, b) (x < b) .* (x.^2) + (x >= b) .* (1./x);

% 定义目标函数
objective_function = @(b) -piecewise_function(0.5, b);  % 在这个示例中，我们希望最小化分段函数在x=0.5处的值，因此使用负号

% 定义约束条件
lb = 0;    % 断点的下界
ub = 1;    % 断点的上界
nonlcon = [];  % 非线性约束为空

% 求解最优断点
x0 = 0;  % 初始点
[b_opt, fval] = fmincon(objective_function, x0, [], [], [], [], lb, ub, nonlcon)

%输出结果为：

% b_opt = 0.4500
% fval = -2

% 接着，通过绘图检验b=0.45是否为最优

x = 0:0.001:1; b1 =0.45; b2 = 0.4; b3 = 0.6;

y1= (x < b1) .* (x.^2) + (x >= b1) .* (1./x);

y2= (x < b2) .* (x.^2) + (x >= b2) .* (1./x);

y3= (x < b3) .* (x.^2) + (x >= b3) .* (1./x);

plot(x,y1,'red',x,y2,'blue',x,y3,'green')

 

 

 关于分段函数的优化问题，我还有非常大的困惑。希望有兴趣的同学可以一起来探讨。