RegularChains[ParametricSystemTools]
RealRootClassification
compute a real root classification of a parametric semi-algebraic system
Calling Sequence
Parameters
Description
Options
Examples
References
RealRootClassification(F, N, P, H, d, n, R)
RealRootClassification(F, N, P, H, d, r, R)
RealRootClassification(F, N, P, H, d, n, R, 'output'='formula')
RealRootClassification(F, N, P, H, d, n, R, 'output'='samples')
R
-
polynomial ring
F
list of polynomials of R
N
P
H
d
positive integer
n
non-negative integer
r
non-negative integer range or an integer range
For a semi-algebraic system with parametric variables (parameters), RealRootClassification computes conditions on the parameters for the system to have a given number of real solutions.
More precisely, consider a semi-algebraic system S with polynomial equations, non-negative polynomial inequalities, (strictly) positive polynomial inequalities, and polynomial inequations given respectively by F, N, P, and H, and with the last d variables of R as parameters.
This command returns a pair C,BP, such that, provided BP does not vanish, the semi-algebraic system S has exactly n distinct real solutions (or, the number of distinct real solutions of S falls into the range r) if and only if the parameters are taking values in C.
In the returned pair, C,BP, C is of type semi_algebraic_set (a list of regular_semi_algebraic_sets) representing the conditions, and BP is of type border_polynomial excluding some exceptional parameter values. Note that all the polynomials in C and BP involve only parameters.
Note that, in special circumstances, the returned pair C,BP is not well defined. See BorderPolynomial or the related explanation below.
r is an integer range of the form a..k where a is a non-negative integer and k is either a non-negative integer satisfying a≤k or a name without value standing for infinity.
There are two special cases where the computed condition C,BP is not well defined or does not return the expected answer. One is when the system S is positive dimensional for almost all parameter values. In this case some variables will be treated as parameters, and the command returns conditions for the system to have real solutions w.r.t. the new parameters set. The second case is when the input system is inconsistent.
There are 2 representations of regular_semi_algebraic_sets in the current implementation according to the different nature of the input systems. The related representations are quantifier-free formula, quantifier-free formula with n-th root indices, and isolating box. For the latter one, see RealRootIsolate.
One can get a more detailed description on the result of the command RealRootClassification by setting the infolevel variable to be greater than 0 (see the examples below).
Because of the inner representation of data in Maple and the strategy of formula simplification in RealRootClassification, the output conditions C for mathematically equivalent input (e.g. all are same but different variable order in ring R) may be different in their form but are equivalent to each other, of course.
The algorithm is described in the paper by Yang, L., Hou, X., Xia, B. "A complete algorithm for automated discovering of a class of inequality-type theorems." Sci. China Ser. F, Vol. 44 (2001): 33--49.
The algorithm is complete in theory, however, the command RealRootClassification may quit and output a necessary condition for some examples due to too heavy computation.
The option 'output'='samples' modifies the algorithm as follows: instead of returning a pair, C,BP, where C is of type semi_algebraic_set (a list of regular_semi_algebraic_sets) the command returns sample points in the parameter space. All these sample points belong to the complement of the hypersurface defined by BP in the parameter space. Moreover there is at least sample point per connected component satisfying the following property: the input parametric semi-algebraic system S has n distinct real solutions (or, the number of distinct real solutions of S falls into the range r) above this component. In addition, there is no sample points belonging to the connected components not satisfying the previous property.
The option 'output'='samples' also modifies the output format as follows. It is a record with five fields. DecomposedSystems is a list of basic semi-algebraic systems (each of them consisting of a list of equational constraints and a list positive inequality constraints) such that the union of their solution sets equals the solution set of the input semi-algebraic system S outside of the hypersurface defined by its border polynomial. BorderPolynomial is the border polynomial of S as computed by this command. Parameters lists the parameters in the projection order used by the underlying OPENCAD sampling algorithm. ProjectionPolynomials is a list of lists of polynomials, where each inner list consists of the polynomials used in each OPENCAD projection. SamplePoints is the list of the sample points described above.
The option 'output'='formula' is the default option of this command.
The computations performed for the option 'output'='samples' are a first step in the computations performed for option 'output'='formula'. This latter option is often more computationally intensive than the former one while the information provided by this former may often be sufficient in practice.
with⁡RegularChains:
with⁡ParametricSystemTools:
with⁡SemiAlgebraicSetTools:
R≔PolynomialRing⁡x,a,b,c
R≔polynomial_ring
F≔a⁢x2+b⁢x+c;N≔;P≔;H≔
F≔a⁢x2+b⁢x+c
N≔
P≔
H≔
We compute the condition on a,b,c for the system to have no real solutions.
rrc≔RealRootClassification⁡F,N,P,H,3,0,R
rrc≔regular_semi_algebraic_set,border_polynomial
We first inspect the border polynomial.
Info⁡rrc2,R
a,−4⁢a⁢c+b2
We then inspect the regular semi-algebraic set.
rc≔RepresentingChain⁡rrc11,R
rc≔regular_chain
Info⁡rc,R
rsas≔rrc11
rsas≔regular_semi_algebraic_set
pbx≔RepresentingBox⁡rsas,R
pbx≔parametric_box
qff≔RepresentingQuantifierFreeFormula⁡pbx
qff≔quantifier_free_formula
Info⁡qff,R
4⁢a⁢c−b2,1
We now pretty print the quantifier-free formula.
DisplayQuantifierFreeFormula⁡qff:
4⁢a⁢c−b2>0
We compute the condition on a,b,c for the system to have real solutions.
rrc≔RealRootClassification⁡F,N,P,H,3,1..k,R
4⁢a⁢c−b2<0
Below we make use of the option 'output'='samples' for the cases where the input system has exactly one positive solution when a is positive.
rrc≔RealRootClassification⁡F,N,a,x,H,3,1,R,output=samples
rrc≔Parameters:⁢c,b,aBorderPolynomial:⁢a,c,−4⁢a⁢c+b2ProjectionPolynomials:⁢c,b,a,−4⁢a⁢c+b2SamplePoints:⁢c=−12,b=−12,a=12,c=−12,b=12,a=12DecomposedSystems:⁢x2⁢a+x⁢b+c,a,x
Using the info printing an alternative and more verbose output is illustrated hereafter.
infolevelRegularChains≔1:
We consider another parametric semi-algebraic system.
R≔PolynomialRing⁡z,y,x,b,a
F≔x2+y2−x⁢y−1,y2+z2−y⁢z−a2,z2+x2−z⁢x−b2
F≔x2−x⁢y+y2−1,−a2+y2−y⁢z+z2,−b2+x2−z⁢x+z2
N≔a−1,b−a;P≔x,y,z,a,b,a+1−b;H≔
N≔a−1,b−a
P≔x,y,z,a,b,a+1−b
We compute the conditions on b,a for the system to have exactly two distinct real solutions.
rrc≔RealRootClassification⁡F,N,P,H,2,2,R
RealRootClassification: FINAL RESULT: RealRootClassification: The system has given number of real solution(s) IF AND ONLY IF RealRootClassification: [R[1] < 0, 0 < R[2], 0 < R[3]] RealRootClassification: where RealRootClassification: R[1] = a^2-b^2+b-1 RealRootClassification: R[2] = a^2-b^2+a+1 RealRootClassification: R[3] = 81*a^12-216*a^10*b^2+144*a^8*b^4-18*a^6*b^6+144*a^4*b^8-216*a^2*b^10+81*b^12-216*a^10+414*a^8*b^2-204*a^6*b^4-204*a^4*b^6+414*a^2*b^8-216*b^10+144*a^8-204*a^6*b^2+241*a^4*b^4-204*a^2*b^6+144*b^8-18*a^6-204*a^4*b^2-204*a^2*b^4-18*b^6+144*a^4+414*a^2*b^2+144*b^4-216*a^2-216*b^2+81 RealRootClassification: PROVIDED THAT RealRootClassification: a-1 <> 0 RealRootClassification: a-b <> 0 RealRootClassification: b-1 <> 0 RealRootClassification: b+1 <> 0 RealRootClassification: a+1 <> 0 RealRootClassification: a+b <> 0 RealRootClassification: a^2-b^2-b-1 <> 0 RealRootClassification: a^2-b^2+b-1 <> 0 RealRootClassification: a^2-b^2-a+1 <> 0 RealRootClassification: a^2-b^2+a+1 <> 0 RealRootClassification: a^2-a*b+b^2-1 <> 0 RealRootClassification: a^2+a*b+b^2-1 <> 0 RealRootClassification: 81*a^12-216*a^10*b^2+144*a^8*b^4-18*a^6*b^6+144*a^4*b^8-216*a^2*b^10+81*b^12-216*a^10+414*a^8*b^2-204*a^6*b^4-204*a^4*b^6+414*a^2*b^8-216*b^10+144*a^8-204*a^6*b^2+241*a^4*b^4-204*a^2*b^6+144*b^8-18*a^6-204*a^4*b^2-204*a^2*b^4-18*b^6+144*a^4+414*a^2*b^2+144*b^4-216*a^2-216*b^2+81 <> 0 RealRootClassification: 1.861*seconds
a−1,a+1,a+b,a−b,b−1,b+1,1−29⁢b6+169⁢b4+169⁢b8+169⁢a4−6827⁢a6⁢b2+24181⁢a4⁢b4−6827⁢a2⁢b6−6827⁢a4⁢b2−6827⁢a2⁢b4+469⁢a2⁢b2−83⁢a2+a12+b12−83⁢a10−83⁢b10−83⁢b2−29⁢a6−6827⁢a4⁢b6+469⁢a2⁢b8−83⁢a10⁢b2+169⁢a8⁢b4−29⁢a6⁢b6+169⁢a4⁢b8−83⁢a2⁢b10+469⁢a8⁢b2−6827⁢a6⁢b4+169⁢a8,a2−b2−b−1,a2−b2+b−1,a2−b2−a+1,a2−b2+a+1,a2−a⁢b+b2−1,a2+a⁢b+b2−1
pbx≔RepresentingBox⁡rrc11,R
a2−b2+b<1anda2−b2+a+1>0and81⁢a12−216⁢a10⁢b2+144⁢a8⁢b4−18⁢a6⁢b6+144⁢a4⁢b8−216⁢a2⁢b10+81⁢b12−216⁢a10+414⁢a8⁢b2−204⁢a6⁢b4−204⁢a4⁢b6+414⁢a2⁢b8−216⁢b10+144⁢a8−204⁢a6⁢b2+241⁢a4⁢b4−204⁢a2⁢b6+144⁢b8−18⁢a6−204⁢a4⁢b2−204⁢a2⁢b4−18⁢b6+144⁢a4+414⁢a2⁢b2+144⁢b4−216⁢a2−216⁢b2+81>0
We discuss the situation when the parameter satisfies a^2+a+1=0.
RealRootClassification⁡a2+a+1,op⁡F,N,P,H,2,2,R
RealRootClassification: The system does NOT have the given number of real solution(s)! RealRootClassification: .72e-1*seconds
,border_polynomial
In the printed information below, (1)S[1] means that if b is specified and satisfies the output inequality, a should be the least real root of S[1].
In the following example, we compute the real root classification of x4+x2⁢a+x⁢b+c=0 when the parameters satisfy b=0 and −2⁢a3+8⁢a⁢c−9⁢b2=0.
F≔b,−2⁢a3+8⁢a⁢c−9⁢b2,x4+a⁢x2+b⁢x+c;N≔;P≔;H≔
F≔b,−2⁢a3+8⁢a⁢c−9⁢b2,x4+x2⁢a+x⁢b+c
RealRootClassification⁡F,N,P,H,3,1..4,R
RealRootClassification: The system has given number of real solution(s) IF AND ONLY IF RealRootClassification: [R[1] < 0, `(1)S`[1], `(1)S`[2]] RealRootClassification: where RealRootClassification: R[1] = a RealRootClassification: and RealRootClassification: S[1] = b RealRootClassification: S[2] = a^2-4*c RealRootClassification: PROVIDED THAT RealRootClassification: a <> 0 RealRootClassification: .12e-1*seconds
regular_semi_algebraic_set,border_polynomial
Yang, L.; Hou, X.; and Xia, B. "A complete algorithm for automated discovering of a class of inequality-type theorems." Science in China Ser F. Vol. 44, (2001): 33-49.
See Also
BorderPolynomial
ComplexRootClassification
DiscriminantSequence
RealRootCounting
RealRootIsolate
RegularChains
Download Help Document