Mixed trigonometric-polynomials (MTPs) are functions of the form f(x,sinx,cosx) where f is a trivariate polynomial with rational coefficients, and the argument x ranges over the reals. In this paper, an algorithm “isolating” all the real roots of an MTP is provided and implemented. It automatically divides the real roots into two parts: one consists of finitely many roots in an interval [μ−,μ+] while the other consists of countably many roots in R﹨[μ−,μ+]. For the roots in [μ−,μ+], the algorithm returns isolating intervals and corresponding multiplicities while for those greater than μ+, it returns finitely many mutually disjoint small intervals Ii⊂[−π,π], integers ci>0 and multisets of root multiplicity {mj,i}j=1ci such that any root >μ+ is in the set (∪i∪k∈N(Ii+2kπ)) and any interval Ii+2kπ⊂(μ+,∞) contains exactly ci distinct roots with multiplicities m1,i,...,mci,i, respectively. The efficiency of the algorithm is shown by experiments. The method used to isolate the roots in [μ−,μ+] is applicable to any other bounded interval as well. The algorithm takes advantages of the weak Fourier sequence technique and deals with the intervals period-by-period without scaling the coordinate so to keep the length of the sequence short. The new approaches can easily be modified to decide whether there is any root, or whether there are infinitely many roots in unbounded intervals of the form (−∞,a) or (a,∞) with a∈Q.
Read full abstract