[论文解读] Trigonometric Splines for Oscillator Simulation

Kai Bittner, Hans-Georg Brachtendorf|arXiv (Cornell University)|Apr 26, 2016

Advanced Measurement and Metrology Techniques被引用 5

一句话总结

本文提出使用三角函数样条来消除振荡器仿真中的数值阻尼，否则将导致能量人为损失并产生不可靠结果。通过用三角函数B样条替代多项式B样条，该方法能精确保持低频振荡，同时维持稳定性和自适应性，从而实现对像3MHz Colpitz-Quartz振荡器这样的自治振荡器的精确长期仿真。

ABSTRACT

We investigate the effects of numerical damping for oscillator simulation with spline methods. Numerical damping results in an artificial loss of energy and leads therefore to unreliable results in the simulation of autonomous systems, as e.g.\ oscillators. We show that the negative effects of numerical damping can be eliminated by the use of trigonometric splines. This will be in particular important for spline based adaptive methods.

研究动机与目标

• 解决基于样条的振荡器仿真中的数值阻尼问题，该问题会人为减少能量，损害自治系统中的精度。
• 研究经典多项式样条为何在基频处引起数值阻尼，特别是由于导数近似误差所致。
• 开发一种样条方法，能精确保持低频振荡，同时仍能抑制高频数值噪声。
• 证明三角函数样条可对低频实现精确的导数表示，从而消除人为能量损失。
• 通过确保在不同网格分辨率下保持一致的数值行为，实现可靠的自适应网格方法。

提出的方法

• 用基于正弦递归公式的三角函数B样条替代多项式B样条，以保持紧支撑性和光滑性。
• 在均匀网格上使用配置法，利用三角函数样条基函数在偏移的网格点处近似函数及其导数。
• 应用离散傅里叶变换分析导数近似方法的放大因子，定义为 ψ̃_m,h(x,ξ) = (∂/∂x φ̃_m,h(x,ξ)) / φ̃_m,h(x,ξ)。
• 利用三角函数样条在奇数阶 m=2μ+1 时能精确表示复指数 e^{2πikt}（当 |k|<μ 时）的性质，确保低频下导数计算的精确性。
• 通过选择配置参数 σ ∈ (−1/2, 0) 来确保稳定性和准确性，避免奇点并维持有利的阻尼行为。
• 通过模拟3MHz Colpitz-Quartz振荡器并将其结果与经典多项式样条进行比较，验证该方法。

实验结果

研究问题

• RQ1 三角函数样条能否在保持稳定性的前提下消除振荡器仿真中的数值阻尼？
• RQ2 为何经典多项式样条会导致振荡器中出现人为能量损失？这种现象如何依赖于网格分辨率和配置点位置？
• RQ3 三角函数样条能否精确表示基频及其导数，从而避免对物理信号的阻尼？
• RQ4 配置参数 σ 的选择如何影响三角函数样条配置法中的数值阻尼和稳定性？
• RQ5 三角函数样条能否扩展到非均匀自适应网格，而不会引入不同网格区域间不一致的阻尼效应？

主要发现

• 奇数阶 m=3,5,… 的三角函数样条能精确表示低频正弦波，确保基频无数值阻尼。
• 三角函数样条的导数放大因子 ψ̃_m,h(σ,ξ) 在 |k|<μ 时等于 2πik，证实了低频下导数计算的精确性。
• 当 σ = −1/4 时，经典多项式样条方法对基频表现出显著阻尼，而三角函数样条完全消除了该效应。
• 使用三角函数样条的仿真结果与预期振幅一致，仅存在近似误差，证实了能量守恒。
• 三角函数样条保持了稳定性，并避免了在多项式样条中因 σ>0 而出现的振铃伪影，后者缺乏高频阻尼。
• 该方法可扩展至自适应方法中使用的非均匀网格，防止在不同网格区域间出现不一致的阻尼效应。