# Optimization of a script

5 views (last 30 days)
Alcide on 24 Aug 2023
Answered: Bruno Luong on 24 Aug 2023
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 on 24 Aug 2023
Do you mean - you want to find value of b and c for which v will be maximum?
Torsten on 24 Aug 2023
b = c = Inf

Bruno Luong on 24 Aug 2023
The v you compute is simply
132*b*(c-1)
You want to optimiza it? Take b=Inf and c=Inf

Kevin Holly on 24 Aug 2023
Edited: Kevin Holly on 24 Aug 2023
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

### Categories

Find more on Solver Outputs and Iterative Display in Help Center and File Exchange

### Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!