Optimization of a script

Question

Alcide 2023-8-24

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2012392-optimization-of-a-script

回答： Bruno Luong 2023-8-24

采纳的回答： Bruno Luong

在 MATLAB Online 中打开

Hi,

I'm working on a script like this:

a=[ 1 1 0 1 1 1 1 0 1]
b=2
c=3.2
v=0
for i=1: 9
    f=b*2^(i-1)
        if a(i)==0
           g=f*c-f
           v=v+g
        end
end

I would like to maximize the value of v using b and c as variables.

How can i do that?

TY

2 个评论
显示无隐藏无

Dyuman Joshi 2023-8-24

Do you mean - you want to find value of b and c for which v will be maximum?

Torsten 2023-8-24

b = c = Inf

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Bruno Luong 2023-8-24

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2012392-optimization-of-a-script#answer_1293217

在 MATLAB Online 中打开

The v you compute is simply

132*b*(c-1)

You want to optimiza it? Take b=Inf and c=Inf

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

Answer 2

Kevin Holly 2023-8-24

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2012392-optimization-of-a-script#answer_1293212

编辑：Kevin Holly 2023-8-24

在 MATLAB Online 中打开

You can use fmincon.

a = [1 1 0 1 1 1 1 0 1];
b0 = 2; % Initial guess for b
c0 = 3.2; % Initial guess for c
% Define the objective function to maximize
objective = @(x) -calculateV(a, x(1), x(2));
% Set the lower and upper bounds for the variables
lb = [0, 0]; % Lower bounds for b and c
ub = [Inf, Inf]; % Upper bounds for b and c
% Perform constrained optimization to maximize v
[x, fval] = fmincon(objective, [b0, c0], [], [], [], [], lb, ub);
Problem appears unbounded.

fmincon stopped because the objective function value is less than
the value of the objective function limit and constraints 
are satisfied to within the value of the constraint tolerance.
% Retrieve the optimized values
b_optimized = x(1);
c_optimized = x(2);
v_maximized = -fval;
% Display the optimized values and the corresponding v
disp("Optimized Values:")
Optimized Values:
disp("b = " + b_optimized)
b = 10167158107334.54
disp("c = " + c_optimized)
c = 10639554538105.3
disp("Maximized v = " + v_maximized)
Maximized v = 1.4278972379828e+28
% Function to calculate v for given a, b, and c
function v = calculateV(a, b, c)
    v = 0;
    for i = 1:numel(a)
        f = b * 2^(i-1);
        if a(i) == 0
            g = f * c - f;
            v = v + g;
        end
    end
end