Reducing system of 10 inequalities
up vote
1
down vote
favorite
I have the following system of inequalities of real variables
sys = 0.383706 x1 + 0.720204 x2 + 1.4568 x3 - 0.000244875 y > 0 && -0.0941312 x1 - 0.176681 x2 - 0.357383 x3 + 0.0000592689 y > 0 && 1.41819 x1 + 3.079 x2 + 2.53801 x3 - 0.00182772 y > 0 && -0.00258971 x1 - 0.00562247 x2 - 0.00463458 x3 + 3.93843*10^-6 y > 0 && 0.129511 x1 + 0.214698 x2 + 0.286739 x3 - 0.0000795243 y > 0 && -0.660462 x1 - 1.09489 x2 - 1.46227 x3 + 0.000405426 y > 0 && 0.351966 x1 + 3.44019 x2 + 7.59863 x3 + 0.00257072 y > 0 && -0.00265509 x1 - 0.0259514 x2 - 0.0573209 x3 - 0.0000199557 y > 0 && 1.34471 x1 + 2.58639 x2 + 3.39561 x3 + 0.0000742173 y > 0 && -0.0314846 x1 - 0.0605571 x2 - 0.079504 x3 - 1.69328*10^-6 y > 0
I want to Reduce
the system in order to see if it admits solutions or not. However, the command Reduce
takes too much time and gets stuck. How can you check if this system admits solutions?
I have tried with FindInstance[sys,{x1,x2,x3,y},Reals]
which returns {}
. However, I am not sure this is a fully proof that the system does not admit solution. Indeed, If I run the command on a subsystem of inequalities I experience the following behaviour
FindInstance[sys[[1 ;; 4]], {x1, x2, x3, y}, Reals]
(* {} *)
but If I ask for more instances, Mathematica found the points
FindInstance[sys[[1 ;; 4]], {x1, x2, x3, y}, Reals, 2]
(*{{x1 -> -33., x2 -> 13.5619, x3 -> 1.98727, y -> 0.255564}, {x1 -> -51., x2 -> 20.9595, x3 -> 3.07117, y -> 0.601467}}*)
So, the subsystem [[1;;4]] admits at least a solution. I have run FindInstance[sys, {x1, x2, x3, y}, Reals, 2]
for all the night and got the output {}
. Still, is this a proof that the system does not admits solution?
equation-solving inequalities findinstance
add a comment |
up vote
1
down vote
favorite
I have the following system of inequalities of real variables
sys = 0.383706 x1 + 0.720204 x2 + 1.4568 x3 - 0.000244875 y > 0 && -0.0941312 x1 - 0.176681 x2 - 0.357383 x3 + 0.0000592689 y > 0 && 1.41819 x1 + 3.079 x2 + 2.53801 x3 - 0.00182772 y > 0 && -0.00258971 x1 - 0.00562247 x2 - 0.00463458 x3 + 3.93843*10^-6 y > 0 && 0.129511 x1 + 0.214698 x2 + 0.286739 x3 - 0.0000795243 y > 0 && -0.660462 x1 - 1.09489 x2 - 1.46227 x3 + 0.000405426 y > 0 && 0.351966 x1 + 3.44019 x2 + 7.59863 x3 + 0.00257072 y > 0 && -0.00265509 x1 - 0.0259514 x2 - 0.0573209 x3 - 0.0000199557 y > 0 && 1.34471 x1 + 2.58639 x2 + 3.39561 x3 + 0.0000742173 y > 0 && -0.0314846 x1 - 0.0605571 x2 - 0.079504 x3 - 1.69328*10^-6 y > 0
I want to Reduce
the system in order to see if it admits solutions or not. However, the command Reduce
takes too much time and gets stuck. How can you check if this system admits solutions?
I have tried with FindInstance[sys,{x1,x2,x3,y},Reals]
which returns {}
. However, I am not sure this is a fully proof that the system does not admit solution. Indeed, If I run the command on a subsystem of inequalities I experience the following behaviour
FindInstance[sys[[1 ;; 4]], {x1, x2, x3, y}, Reals]
(* {} *)
but If I ask for more instances, Mathematica found the points
FindInstance[sys[[1 ;; 4]], {x1, x2, x3, y}, Reals, 2]
(*{{x1 -> -33., x2 -> 13.5619, x3 -> 1.98727, y -> 0.255564}, {x1 -> -51., x2 -> 20.9595, x3 -> 3.07117, y -> 0.601467}}*)
So, the subsystem [[1;;4]] admits at least a solution. I have run FindInstance[sys, {x1, x2, x3, y}, Reals, 2]
for all the night and got the output {}
. Still, is this a proof that the system does not admits solution?
equation-solving inequalities findinstance
1
Reduce[Rationalize[sys, 2^-12], {x1, x2, x3, y}, Reals]
givesFalse
.
– kglr
yesterday
1
andRegionIntersection @@ (ImplicitRegion[#, {x1, x2, x3, y}] & /@ (List @@ sys))
givesEmptyRegion[4]
.
– kglr
yesterday
add a comment |
up vote
1
down vote
favorite
up vote
1
down vote
favorite
I have the following system of inequalities of real variables
sys = 0.383706 x1 + 0.720204 x2 + 1.4568 x3 - 0.000244875 y > 0 && -0.0941312 x1 - 0.176681 x2 - 0.357383 x3 + 0.0000592689 y > 0 && 1.41819 x1 + 3.079 x2 + 2.53801 x3 - 0.00182772 y > 0 && -0.00258971 x1 - 0.00562247 x2 - 0.00463458 x3 + 3.93843*10^-6 y > 0 && 0.129511 x1 + 0.214698 x2 + 0.286739 x3 - 0.0000795243 y > 0 && -0.660462 x1 - 1.09489 x2 - 1.46227 x3 + 0.000405426 y > 0 && 0.351966 x1 + 3.44019 x2 + 7.59863 x3 + 0.00257072 y > 0 && -0.00265509 x1 - 0.0259514 x2 - 0.0573209 x3 - 0.0000199557 y > 0 && 1.34471 x1 + 2.58639 x2 + 3.39561 x3 + 0.0000742173 y > 0 && -0.0314846 x1 - 0.0605571 x2 - 0.079504 x3 - 1.69328*10^-6 y > 0
I want to Reduce
the system in order to see if it admits solutions or not. However, the command Reduce
takes too much time and gets stuck. How can you check if this system admits solutions?
I have tried with FindInstance[sys,{x1,x2,x3,y},Reals]
which returns {}
. However, I am not sure this is a fully proof that the system does not admit solution. Indeed, If I run the command on a subsystem of inequalities I experience the following behaviour
FindInstance[sys[[1 ;; 4]], {x1, x2, x3, y}, Reals]
(* {} *)
but If I ask for more instances, Mathematica found the points
FindInstance[sys[[1 ;; 4]], {x1, x2, x3, y}, Reals, 2]
(*{{x1 -> -33., x2 -> 13.5619, x3 -> 1.98727, y -> 0.255564}, {x1 -> -51., x2 -> 20.9595, x3 -> 3.07117, y -> 0.601467}}*)
So, the subsystem [[1;;4]] admits at least a solution. I have run FindInstance[sys, {x1, x2, x3, y}, Reals, 2]
for all the night and got the output {}
. Still, is this a proof that the system does not admits solution?
equation-solving inequalities findinstance
I have the following system of inequalities of real variables
sys = 0.383706 x1 + 0.720204 x2 + 1.4568 x3 - 0.000244875 y > 0 && -0.0941312 x1 - 0.176681 x2 - 0.357383 x3 + 0.0000592689 y > 0 && 1.41819 x1 + 3.079 x2 + 2.53801 x3 - 0.00182772 y > 0 && -0.00258971 x1 - 0.00562247 x2 - 0.00463458 x3 + 3.93843*10^-6 y > 0 && 0.129511 x1 + 0.214698 x2 + 0.286739 x3 - 0.0000795243 y > 0 && -0.660462 x1 - 1.09489 x2 - 1.46227 x3 + 0.000405426 y > 0 && 0.351966 x1 + 3.44019 x2 + 7.59863 x3 + 0.00257072 y > 0 && -0.00265509 x1 - 0.0259514 x2 - 0.0573209 x3 - 0.0000199557 y > 0 && 1.34471 x1 + 2.58639 x2 + 3.39561 x3 + 0.0000742173 y > 0 && -0.0314846 x1 - 0.0605571 x2 - 0.079504 x3 - 1.69328*10^-6 y > 0
I want to Reduce
the system in order to see if it admits solutions or not. However, the command Reduce
takes too much time and gets stuck. How can you check if this system admits solutions?
I have tried with FindInstance[sys,{x1,x2,x3,y},Reals]
which returns {}
. However, I am not sure this is a fully proof that the system does not admit solution. Indeed, If I run the command on a subsystem of inequalities I experience the following behaviour
FindInstance[sys[[1 ;; 4]], {x1, x2, x3, y}, Reals]
(* {} *)
but If I ask for more instances, Mathematica found the points
FindInstance[sys[[1 ;; 4]], {x1, x2, x3, y}, Reals, 2]
(*{{x1 -> -33., x2 -> 13.5619, x3 -> 1.98727, y -> 0.255564}, {x1 -> -51., x2 -> 20.9595, x3 -> 3.07117, y -> 0.601467}}*)
So, the subsystem [[1;;4]] admits at least a solution. I have run FindInstance[sys, {x1, x2, x3, y}, Reals, 2]
for all the night and got the output {}
. Still, is this a proof that the system does not admits solution?
equation-solving inequalities findinstance
equation-solving inequalities findinstance
asked yesterday
apt45
506211
506211
1
Reduce[Rationalize[sys, 2^-12], {x1, x2, x3, y}, Reals]
givesFalse
.
– kglr
yesterday
1
andRegionIntersection @@ (ImplicitRegion[#, {x1, x2, x3, y}] & /@ (List @@ sys))
givesEmptyRegion[4]
.
– kglr
yesterday
add a comment |
1
Reduce[Rationalize[sys, 2^-12], {x1, x2, x3, y}, Reals]
givesFalse
.
– kglr
yesterday
1
andRegionIntersection @@ (ImplicitRegion[#, {x1, x2, x3, y}] & /@ (List @@ sys))
givesEmptyRegion[4]
.
– kglr
yesterday
1
1
Reduce[Rationalize[sys, 2^-12], {x1, x2, x3, y}, Reals]
gives False
.– kglr
yesterday
Reduce[Rationalize[sys, 2^-12], {x1, x2, x3, y}, Reals]
gives False
.– kglr
yesterday
1
1
and
RegionIntersection @@ (ImplicitRegion[#, {x1, x2, x3, y}] & /@ (List @@ sys))
gives EmptyRegion[4]
.– kglr
yesterday
and
RegionIntersection @@ (ImplicitRegion[#, {x1, x2, x3, y}] & /@ (List @@ sys))
gives EmptyRegion[4]
.– kglr
yesterday
add a comment |
1 Answer
1
active
oldest
votes
up vote
3
down vote
accepted
Reduce[Rationalize[sys, 2^-12], {x1, x2, x3, y}, Reals] // RepeatedTiming
{0.24878, False}
RegionIntersection @@ (ImplicitRegion[#, {x1, x2, x3, y}] & /@ (List @@ sys)) //
RepeatedTiming
{0.0152, EmptyRegion[4]}
add a comment |
1 Answer
1
active
oldest
votes
1 Answer
1
active
oldest
votes
active
oldest
votes
active
oldest
votes
up vote
3
down vote
accepted
Reduce[Rationalize[sys, 2^-12], {x1, x2, x3, y}, Reals] // RepeatedTiming
{0.24878, False}
RegionIntersection @@ (ImplicitRegion[#, {x1, x2, x3, y}] & /@ (List @@ sys)) //
RepeatedTiming
{0.0152, EmptyRegion[4]}
add a comment |
up vote
3
down vote
accepted
Reduce[Rationalize[sys, 2^-12], {x1, x2, x3, y}, Reals] // RepeatedTiming
{0.24878, False}
RegionIntersection @@ (ImplicitRegion[#, {x1, x2, x3, y}] & /@ (List @@ sys)) //
RepeatedTiming
{0.0152, EmptyRegion[4]}
add a comment |
up vote
3
down vote
accepted
up vote
3
down vote
accepted
Reduce[Rationalize[sys, 2^-12], {x1, x2, x3, y}, Reals] // RepeatedTiming
{0.24878, False}
RegionIntersection @@ (ImplicitRegion[#, {x1, x2, x3, y}] & /@ (List @@ sys)) //
RepeatedTiming
{0.0152, EmptyRegion[4]}
Reduce[Rationalize[sys, 2^-12], {x1, x2, x3, y}, Reals] // RepeatedTiming
{0.24878, False}
RegionIntersection @@ (ImplicitRegion[#, {x1, x2, x3, y}] & /@ (List @@ sys)) //
RepeatedTiming
{0.0152, EmptyRegion[4]}
answered yesterday
kglr
171k8194399
171k8194399
add a comment |
add a comment |
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
StackExchange.ready(
function () {
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmathematica.stackexchange.com%2fquestions%2f186279%2freducing-system-of-10-inequalities%23new-answer', 'question_page');
}
);
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
1
Reduce[Rationalize[sys, 2^-12], {x1, x2, x3, y}, Reals]
givesFalse
.– kglr
yesterday
1
and
RegionIntersection @@ (ImplicitRegion[#, {x1, x2, x3, y}] & /@ (List @@ sys))
givesEmptyRegion[4]
.– kglr
yesterday