Optimization of a script

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

[Dyuman Joshi]

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

[Torsten]

b = c = Inf

Answer 1

[Bruno Luong]

The v you compute is simply

132*b*(c-1)

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

Answer 2

[Kevin Holly]

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