linear inequality constrains based on absolute values
2 vues (au cours des 30 derniers jours)
Afficher commentaires plus anciens
Please help me to define the inequality constrains for quadprog in the below scenario
x+y <= 0.1*abs(x)
x+y >= -0.1*abs(x)
0 commentaires
Réponse acceptée
Matt J
le 6 Oct 2022
Modifié(e) : Matt J
le 6 Oct 2022
The constraints correspond to a non-convex region in (as Walter's second plot shows). You would have to break it into two regions and optimize over each one separately:
Region 1:
x<=0
x+y <= 0.1*(-x)
x+y >= -0.1*(-x)
Region 2:
x>=0
x+y <= 0.1*(x)
x+y >= -0.1*(x)
3 commentaires
Walter Roberson
le 7 Oct 2022
There is no point which is not in one of the regions or the other, so solving separately and looking for the best between the two is going to get you the same result as if you had no constraint.
It would make more sense if the conditions were "and" and you processed the intersection of the constraints in two pieces, one for negative x and the other for non-negative x, and took the best between the two of those.
Plus de réponses (1)
Walter Roberson
le 6 Oct 2022
x = linspace(-0.005, 0.005, 100);
y = linspace(-0.005, 0.005, 101).';
M1 = x + y <= 0.1*abs(x);
M2 = x + y >= -0.1*abs(x);
[r1, c1] = find(M1);
[r2, c2] = find(M2);
plot(x(c1), y(r1), 'k.', x(c2), y(r2), 'ro');
As you can see from the plot, there is nowhere which is not part of one of the regions or the other, so nothing is constrained out.
Matters would be different if the constraints were "and".
x = linspace(-0.0003, 0.0003, 1000);
y = linspace(-0.0003, 0.0003, 1001).';
M1 = x + y <= 0.1*abs(x);
M2 = x + y >= -0.1*abs(x);
[r1, c1] = find(M1 & M2);
plot(x(c1), y(r1), 'k.');
I think the small gap is a matter of resolution.
0 commentaires
Voir également
Catégories
En savoir plus sur Solver Outputs and Iterative Display dans Help Center et File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!