how to express 1/(x(x+1)) in CVX, where x is a variable?
1/x(x+1)<=y is equivalent to -logx -log(x+1) <= logy
1 Like
I see CVX has prod_inv
so maybe prod_inv([x x+1])
should work. (Or whatever the correct Matlab syntax to make that list).
The formulation and syntax in @Michal_Adamaszek 's post is correct.
help prod_inv
prod_inv inverse of the product of a positive vector.
For a real vector, matrix, or X, prod_inv(X) returns 1.0 ./ PROD(X) if
the elements of X are all positive, and +Inf otherwise.For matrices, prod_inv(X) is a row vector containing the inverse product of each column of X. For N-D arrays, prod_inv(X) is an array of inverse products taken along the first non-singleton dimension of X. prod_inv(X,DIM) takes inverse products along the dimension DIM of X. prod_inv(X,DIM,P), where P is a positive real constant, computes prod_inv(X).^P. This is slightly more efficient than the equivalent POW_POS(prod_inv(X),P). Disciplined convex programming information: prod_inv(X) is convex and nonincreasing in X; therefore, when used in CVX specifications, its argument must be concave or affine.
The inequality to be modelled can be rewritten as
(y*(x)(x+1))^1/3 >= 1
The left hand side is the geometric mean of [y x x+1].
Most likely this is the idea used to implement prod_inv.
Yes, prod_inv
uses geo_mean
.