Optimization of a script

4 visualizaciones (últimos 30 días)
Alcide
Alcide el 24 de Ag. de 2023
Respondida: Bruno Luong el 24 de Ag. de 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
  2 comentarios
Dyuman Joshi
Dyuman Joshi el 24 de Ag. de 2023
Do you mean - you want to find value of b and c for which v will be maximum?
Torsten
Torsten el 24 de Ag. de 2023
b = c = Inf

Iniciar sesión para comentar.

Respuesta aceptada

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

Más respuestas (1)

Kevin Holly
Kevin Holly el 24 de Ag. de 2023
Editada: Kevin Holly el 24 de Ag. de 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

Categorías

Más información sobre Optimization en Help Center y File Exchange.

Etiquetas

Community Treasure Hunt

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

Start Hunting!

Translated by