[論文レビュー] 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.

研究の動機と目的

• 与えられた区間、特に[0, π]上で三角多項式の非負性を確認する数学的に厳密な手続きを開発すること。
• 特に低次のケースや単調性の仮定が成り立たない場合に限界を示す既存の一般的基準の限界を克服すること。
• MAPLE 13などの記号計算ソフトウェアを用いて、根の数え上げと非負性の検証に伴う煩雑な代数的計算を自動化・簡素化すること。
• 特定の三角和やパラメータ依存不等式の証明における本手法の有効性を示すこと。
• 研究者が非負性の主張をコンピュータ支援だが数学的に根拠のある計算に基づいて再現可能に検証できるフレームワークを提供すること。

提案手法

• 三角多項式を、T_k(cos(x)) として表されるチェビシェフ多項式 T_k を用いた三角恒等式 cos(kx) = T_k(cos(x)) を利用して、y = cos(x) に関する代数的多項式に変換する。
• 得られた代数的多項式に対してSturmの定理を適用し、区間 y ∈ [-1, 1] 内に存在する実数解の正確な個数を計算する。これは x ∈ [0, π] に対応する。
• [-1, 1] 内に実数解が存在しないことを条件として、元の三角多項式が[0, π]上で非負であることを確認する。
• MAPLE 13の組み込み関数 'sturm' および 'discrim' を用いてSturm列の計算を実装し、根の数え上げと判別式解析を自動化する。
• パラメータ依存の三角多項式族に対しては、境界ケース（y = ±1 での端点）と判別式が0になる条件を分析し、非負性を保証するパラメータ値の最大区間を求める。
• 重解が存在する場合（判別式が0）の臨界ケースを検証し、Sturm列を用いて境界パラメータ値における非負性を確認することで正しさを保証する。

実験結果

リサーチクエスチョン

• RQ1 Sturmの定理を[0, π]上での三角多項式の非負性を数学的に厳密かつ計算的に実行可能に判定するために体系的に適用可能か？
• RQ2 MAPLE 13 などの記号計算ソフトウェアを用いて、Sturmの定理に必要な根の数え上げプロセスを三角多項式に対して効果的に自動化できるか？
• RQ3 パラメータ a に対して、三角多項式 S_8(x) = 2sin(x) + sin(2x) + a sin(3x) が[0, τ]上で非負であるための最大区間は何か？ そして、この区間をアルゴリズム的にどのように特定できるか？
• RQ4 異なる境界挙動（例：端点で負になる、内部に解が存在する）は、パラメータ依存の三角和の非負性区間の決定にどのように影響するか？
• 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列の計算により非負性が確認された。
• S_8(x) = 2\text{sin}(x) + \text{sin}(2x) + a\text{sin}(3x) が[0, \tau]上で非負であるためのパラメータ a の最大区間は $[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 となる条件と、代数的多項式の判別式が0になる条件を解くことで得られた。
• S_{10}(x) では、非負性を保証する a の最大区間は $[0, 31.513\text{...}]$ であり、上界は判別式の根に対応する。Sturm手順により、重解が [-1, 1] 外に存在するが、境界においても非負性が確認された。
• S_{11}(x) は $[0, \tau/2]$ 上で定義され、a の最大非負性区間は $[-4/3, 28.98537710\text{...}]$ である。右端は判別式の根によって決定され、Sturm列を用いて境界での非負性が確認された。