[论文解读] Nonnegative Trigonometric Polynomials, Sturms Theorem, and Symbolic Computation

Man Kam Kwong|arXiv (Cornell University)|Feb 27, 2014

Mathematical functions and polynomials参考文献 12被引用 8

一句话总结

本文提出了一种基于Sturm定理的严格符号计算方法，用于确定给定三角多项式在指定区间（特别是[0, π]）上是否为非负。该方法将三角多项式转化为cos(x)的代数多项式，应用Sturm序列计算区间内实根数量，当无实根时确认非负性；该方法在数学上完全严格，并在MAPLE 13中实现以实现自动化与验证，从而能够发现并确认复杂三角和式及参数化族的非负性。

ABSTRACT

In this paper, we explain a procedure based on a classical result of Sturm that can be used to determine rigorously whether a given trigonometric polynomial is nonnegative in a certain interval or not. Many examples are given. This technique has been employed by the author in several recent works. The procedure often involves tedious computations that are time-consuming and error-prone. Fortunately, symbolic computation software is available to automate the procedure. In this paper, we give the details of its implementation in MAPLE 13. Some who are strongly attached to a more traditional theoretical research framework may find such details boring or even consider computer-assisted proofs suspicious. However, we emphasize again that the procedure is completely mathematically rigorous.

研究动机与目标

• 开发一种基于经典Sturm定理的数学上严格的方法，用于验证给定区间（特别是[0, π]）上三角多项式的非负性。
• 解决现有通用非负三角多项式判别准则的局限性，尤其是在低次多项式或单调性假设不成立的情况下。
• 利用MAPLE 13等符号计算软件自动化并简化三角多项式根计数与非负性验证中涉及的繁琐代数计算。
• 展示该方法在确认特定三角和式的非负性以及证明参数相关不等式方面的实用性。
• 为研究人员提供一个可验证、可复现的框架，通过计算机辅助但数学根基严谨的计算方式，验证三角不等式中的非负性声明。

提出的方法

• 利用三角恒等式（如cos(kx) = T_k(cos(x))，其中T_k为切比雪夫多项式），将三角多项式转化为关于y = cos(x)的代数多项式。
• 对所得代数多项式应用Sturm定理，精确计算在区间y ∈ [-1, 1]（对应x ∈ [0, π]）内的实根数量。
• 将[-1, 1]内无实根作为原三角多项式在[0, π]上非负的充分条件。
• 在MAPLE 13中使用内置函数（如'sturm'和'discrim'）实现Sturm序列计算，以自动化根计数与判别式分析。
• 对于参数化三角多项式族，分析边界情况（y = ±1处）及判别式为零的条件，以确定确保非负性的最大参数值区间。
• 验证多重根出现的关键情形（判别式为零），并通过Sturm序列确认边界参数值处的非负性，以确保正确性。

实验结果

研究问题

• RQ1 Sturm定理能否系统性地应用于以数学严格性与计算可行性确定三角多项式在[0, π]上的非负性？
• RQ2 如何有效利用MAPLE 13等符号计算软件自动化Sturm定理所需的根计数过程，以处理三角多项式？
• RQ3 参数a的最优区间为何值时，参数化三角多项式在[0, π]上保持非负？该区间能否通过算法确定？
• RQ4 不同边界行为（如端点处为负 vs. 内部存在根）如何影响参数化三角和式非负性区间的判定？
• RQ5 当已知通用准则（如Vietoris或Belov准则）不适用时（特别是低次或非单调系数序列情况下），该方法能否用于验证非负性？

主要发现

• Sturm过程成功验证了三角多项式$ C_1(x) = 5 + 4\text{cos}(x) + 3\text{cos}(2x) + 4\text{cos}(3x) $在$[0, \tau]$上的非负性，方法是将其转化为代数多项式，并确认在$[-1, 1]$内无实根。
• 对于$ C_2(x) = 7 + 6\text{cos}(x) + 5\text{cos}(2x) + 4\text{cos}(3x) + 3\text{cos}(4x) + 5\text{cos}(5x) $，尽管代数多项式$ P_2(y) $次数较高（为5次），该方法仍通过Sturm序列计算确认了其非负性。
• 参数$ a $确保$ S_8(x) = 2\text{sin}(x) + \text{sin}(2x) + a\text{sin}(3x) $在$[0, \tau]$上非负的最大区间为$[0, 1 + \frac{\tau}{2}]$，其中上界对应于关于$ y = \text{cos}(x) $的二次多项式存在重根的情况。
• 对于$ S_9(x) = 2\text{sin}(x) + \text{sin}(2x) + a\text{sin}(3x) + \text{sin}(4x) $，参数$ a $的最大非负区间为$[4/3, 1.881648914\text{...}]$，通过求解$ X(-1) \to 0 $并找到代数多项式判别式为零的位置确定。
• 在$ S_{10}(x) $中，确保非负性的$ a $的最大区间为$[0, 31.513\text{...}]$，上界对应于判别式的根，且Sturm过程在该边界处仍确认了非负性，尽管重根位于$[-1, 1]$之外。
• 在$ S_{11}(x) $于$[0, \tau/2]$上，最大区间为$[-4/3, 28.98537710\text{...}]$，右端点由判别式的根确定，且该方法通过Sturm序列在该边界处确认了非负性。