[论文解读] On positive functions with positive Fourier transforms

B. G. Giraud, R. Peschanski|ArXiv.org|Apr 5, 2005

Differential Equations and Boundary Problems参考文献 3被引用 8

一句话总结

本文提出一种构造性方法，通过使用埃尔米特-傅里叶本征态和斯图姆定理，利用多项式根约束来强制实现正性，从而生成其傅里叶变换也为正的实函数。关键贡献在于基于埃尔米特多项式和拉盖尔多项式的系统性代数框架，将正性条件简化为基于系数的多项式分析，实现了在1D和2D径向情形下此类函数的显式构造，并在物理学和概率论中有应用。

ABSTRACT

Using the basis of Hermite-Fourier functions (i.e. the quantum oscillator eigenstates) and the Sturm theorem, we derive the practical constraints for a function and its Fourier transform to be both positive. We propose a constructive method based on the algebra of Hermite polynomials. Applications are extended to the 2-dimensional case (i.e. Fourier-Bessel transforms and the algebra of Laguerre polynomials) and to adding constraints on derivatives, such as monotonicity or convexity.

研究动机与目标

• 为解决在物理学、概率论和晶体学中具有重要意义的公开问题——即显式构造其傅里叶变换也为正的函数——提出方法。
• 开发一种基于正交多项式的实用代数方法，以强制实现函数及其傅里叶变换的正性。
• 将该方法从1D傅里叶变换扩展至使用拉盖尔多项式的2D径向傅里叶-贝塞尔变换。
• 通过分析多项式导数，纳入单调性或凸性等附加约束。
• 提供一种系统性、数值可处理的框架，通过截断和斯图姆序列识别确保正性的参数区域。

提出的方法

• 该方法使用埃尔米特-傅里叶函数（量子谐振子本征态）作为基，因为其傅里叶变换是自身的相位平移版本。
• 将展开限制在偶数阶埃尔米特函数（本征值+1对应4的倍数，本征值-1对应4+2的倍数），以保持函数的实值性和变换的正性。
• 将函数及其变换的正性问题转化为检查相关多项式 P+(ρ) 和 P−(σ) 是否无正实根，其中 ρ = r² 且 σ = s²。
• 应用斯图姆定理以高效计数并定位这些多项式的实根，从而识别出无正实根存在的参数区域。
• 在2D径向问题中，通过在 r² 上使用拉盖尔多项式扩展形式体系，保持相同的代数结构。
• 截断至有限阶 N 允许实际计算，且通过分析多项式系数的导数，可识别凸性或单调性约束。

实验结果

研究问题

• RQ1 一个实函数及其傅里叶变换均正的充要条件是什么？此类函数如何被显式构造？
• RQ2 如何利用正交多项式展开和根求解定理来强制实现函数及其傅里叶变换的正性？
• RQ3 该方法能否通过在拉盖尔多项式中使用 r² 来扩展至2D径向傅里叶-贝塞尔变换？
• RQ4 在保持正性的同时，如何将导数约束（如单调性、凸性）纳入构造中？
• RQ5 截断至有限阶多项式如何影响该方法在逼近非高斯正函数时的灵活性与收敛性？

主要发现

• 该方法通过确保相关多项式 P+(ρ) 和 P−(σ) 无正实根，并利用斯图姆序列进行根计数，成功构造出其函数及其傅里叶变换均为正的函数。
• 高斯函数（n=0）是最简单的解，更高阶解可通过具有偶数索引的埃尔米特-傅里叶态的叠加构建，且在特定系数约束下保持正性。
• 在2D径向情形下，该方法使用 r² 上的拉盖尔多项式，ψ 和 φ 的正性区域被可视化，显示其与1D情形具有拓扑相似性，但具有不同的数值行为。
• 2D中 φ 正性的充分条件是三阶导数 d³ψ/dr³ 为负，但该条件可能仅在截断后 r ≥ rt 的有限区域成立。
• 该方法生成高度振荡的解（例如，最高达12阶多项式），其形态类似狄拉克梳，表明该基函数足够灵活，可逼近非高斯正函数。
• 数值示例表明，即使在有限截断下，该方法仍可产生具有长近似周期性和强正性的函数，表明其在截断下仍具有良好的收敛性。