[论文解读] Extraction of harmonics from trigonometric polynomials by amplitude and phase operators

V. I. Danchenko, Darya Vasilchenkova|arXiv (Cornell University)|Jun 28, 2016

Electromagnetic Scattering and Analysis参考文献 6被引用 5

一句话总结

本文提出了一种新颖的解析方法，通过振幅和相位算子（AFOs）精确提取实三角多项式中的特定谐波分量。该方法通过原多项式的有限线性组合（经相位移位和缩放）构造目标谐波，在条件 $n = s\mu - 1$ 下获得所需振幅和相位的精确闭式公式，其中 $s \in \mathbb{N}$，从而克服了传统谱方法和Prony方法的局限性。

ABSTRACT

Extraction of harmonics of a given order from real trigonometric polynomials (signals) is one of the main problems in harmonic analysis. It has many applications in physics, radio and electrical engineering, in particular, in filtration of harmonic signals of different nature. There exist many methods (mainly, approximative) for solution of this problem. The most common ones are spectral methods based on Fourier transform and other resonance principles. In this paper we propose a new method for extracting harmonics by amplitude and phase transformation of trigonometric polynomials. These transformations use the two simplest operations -- multiplication by a real constant and phase shift -- to obtain polynomials similar to the initial ones. A harmonic is extracted by an amplitude and phase operator that simply overlays (sums up) a finite number of such similar polynomials. The overlay method enables us to obtain precise analytical formulas for calculating harmonics of a given order.

研究动机与目标

• 解决从实三角多项式 $T_n(x)$ 中提取特定阶数 $\mu$ 谐波的根本问题。
• 开发一种可生成谐波提取精确解析公式的算法，避免谱方法和共振方法中常见的近似。
• 解决构造一个阶数 $\leq n$ 的实值振幅和相位算子（AFO）以分离单个谐波分量的挑战。
• 解决标准Prony方法在条件 $n = s\mu - 1$ 下因非正则解而失效的问题。
• 设计正则化技术以处理非正则解，从而恢复有效且实值的AFO参数。

提出的方法

• 该方法使用定义为 $H_m(T_n, \{X_j\}, \{\lambda_j\}; x) = \sum_{j=1}^n X_j \cdot T_n(x - \lambda_j)$ 的振幅和相位算子（AFOs），其中 $X_j$ 为振幅，$\lambda_j$ 为相位移。
• 问题被重新表述为涉及复变量 $Y_j = X_j z_j$ 和 $z_j = e^{-i\lambda_j}$ 的离散矩问题，导出与谐波提取条件等价的方程组。
• 通过添加辅助方程并使用极限过程对非正则解进行正则化，特别处理 $n = s\mu - 1$ 的情况。
• 正则化过程涉及构建一系列 $({\mathcal{B}}^{\prime\prime})$-正则系统，其生成多项式为 $\hat{G}_k$，确保收敛至有效解。
• 通过在正则化系统上应用经典Prony方法求解，且当 $k \to \infty$ 时，根 $z_j^{(k)}$ 收敛至期望值。
• 最终解确保AFO为实值且阶数恰好为 $n$，满足条件 $\tau_\mu(x) = H_n(\cdots)$。

实验结果

研究问题

• RQ1 能否构造一个实值振幅和相位算子，以从次数为 $n$ 的三角多项式中精确提取阶数为 $\mu$ 的谐波？
• RQ2 在 $n$ 和 $\mu$ 满足何种条件下，谐波提取问题存在精确解析解？
• RQ3 为何标准Prony方法在 $n = s\mu - 1$ 条件下失效，以及如何克服这一失效？
• RQ4 如何将正则化技术应用于非正则解，以恢复有效且实值的AFO参数？
• RQ5 正则化系统中根和系数的极限行为如何，以及其如何确保收敛至正确谐波？

主要发现

• 当 $n = s\mu - 1$ 且 $s \in \mathbb{N}$ 时，本文推导出提取 $\mu$-阶谐波所需振幅 $X_j$ 和相位移 $\lambda_j$ 的显式精确解析公式。
• 解通过应用于非正则系统的正则化过程构造，且正则化系统的极限给出有效且实值的AFO。
• 正则化系统根 $z_j^{(k)}$ 收敛至正确值：当 $k \to \infty$ 时，$z_j^{(k)} \to 0$ 对于 $j = 1,\dots,\mu-1$，且 $z_j^{(k)} \to z_j$ 对于 $j = \mu,\dots,n$。
• 该方法确保AFO为实值且阶数恰好为 $n$，满足条件 $\tau_\mu(x) = H_n(\cdots)$，最终解由正则化系统的极限导出。
• 系数 $Z_j^{(k)}$ 收敛至有限极限，证实了在给定条件下原系统（1.3）存在解。
• 该方法成功克服了经典Prony方法在 $n = s\mu - 1$ 条件下因非正则解而失效的问题。