🔢 mathematics

Analysis of Delay Differential Equations Using the Offset Linear Canonical Transform

本文提出了一种利用偏移线性规范变换（OLCT）将时滞微分方程转化为 Volterra 型积分方程的新型框架，并通过解析与数值结果证明，该方法相比传统的傅里叶变换方法具有更高的精度和有效性。

原作者： Gita Rani Mahato Manab Kundu

发布于 2026-09-30

📖 1 分钟阅读🧠 深度阅读

原作者： Gita Rani Mahato Manab Kundu

原始论文采用 CC BY 4.0 许可（http://creativecommons.org/licenses/by/4.0/）。 ✨ 这是对下方论文的AI生成解释。它不是由作者撰写或认可的。如需技术准确性，请参阅原始论文。阅读完整免责声明

自然界中的许多系统并不会对变化做出即时反应。植物不会在接收到水分的瞬间就开始生长；经济也不会在政策公布的瞬间发生转变。这些系统带有某种“记忆”，即当前的状态是由过去所发生的事情所塑造的。科学家们使用一种被称为延迟微分方程（delay differential equation）的特定数学方程来描述这种行为。这些方程对于模拟从传染病的传播到工程控制系统的稳定性等各种事物都至关重要，因为它们考虑了原因与结果之间的时间滞后。然而，求解这些方程是极其困难的。延迟引入了一种复杂性，使得标准的数学工具往往难以处理，经常导致难以预测或控制的振荡或不稳定性。

几十年来，研究人员一直依赖一种名为傅里叶变换（Fourier transform）的强大数学工具来分析此类系统。这个工具就像一个棱镜，将复杂的信号分解为更简单的、潜在的频率，从而更容易观察系统的随时间变化的行为。虽然有效，但傅里ย里变换是一个固定的仪器；它有一种特定的观察世界的方式，无法被轻易调整。近年来，数学家们开发了一种更灵活的工具版本，称为偏移线性正则变换（offset linear canonical transform）。可以将傅里叶变换想象成一个标准的相机镜头，而这个新工具则是一个带有可调节旋钮的镜头，允许用户在分析数据之前对其进行平移、缩放和倾斜处理。这种增加的灵活性提供了一种观察旧问题的新方式，可能揭示出标准方法可能会遗漏的细节。

在最近的一项研究中，研究人员 Gita Rani Mahato 和 Manab Kundu 探讨了这种更灵活的数学透镜是否能为求解延迟微分方程提供更好的框架。他们的目标是观察偏移线性正则变换是否不仅能求解这些方程，还能提供一种理解它们的统一方式，从而弥合不同数学方法之间的鸿沟。他们专注于一类其中延迟为常数的方程，这是现实世界建模中的常见场景。通过应用他们的新方法，他们将这种涉及变化率和过去值的困难微分方程，转化为另一种被称为积分方程（integral equation）的问题。这种新形式通常在数值处理上更容易处理，允许研究人员按时间间隔逐步构建解，就像走上台阶一样向前推进。

团队发现这种方法效果显著。他们证明了新方法可以产生与使用传统傅里叶变换获得的结果相匹配的精确解。事实上，他们表明传统的傅里叶方法实际上只是他们这个更广泛、更灵活框架中的一个特殊的、灵活性较低的特例。当他们利用不同的延迟值和其他参数测试其方法时，结果是一致的。他们观察到，随着系统中延迟的增加，解的行为变得更加平稳和稳定，振荡消退得更快。相反，当延迟较短时，系统表现出更剧烈且快速的变化。研究人员还注意到，方程中一个关键参数的正负决定了延迟效应是支持还是反对当前状态，从而影响系统趋于稳定的速度。

为了确保研究结果的可靠性，研究人员将他们的结果与其它科学家使用傅里叶变换得出的既定解进行了对比。对比显示出高度的一致性。两种方法的解的实部紧密重合，证实了新方法是可靠的。值得注意的唯一区别在于，他们的新方法产生了一个包含实部和虚部的复数结果，而传统方法仅产生实部。这种由新方法提供的额外信息是其更复杂数学结构的自然结果，在不违背已知结果的情况下，提供了对系统行为更丰富的描述。该研究还对误差行为进行了详细观察，显示出两种方法之间的差异很小，且主要发生在振荡的峰值附近，从而进一步验证了新技术的准确性。

最终，这项工作为解决这些复杂的时滞问题建立了极具前景的新基础。研究人员表明，通过使用偏移线性正则变换，他们可以将困难的延迟微分方程转换为一种在特殊情况下具有解析可解性、而在一般情况下具有数值鲁棒性的形式。虽然这项研究侧重于线性方程，但作者指出，他们方法的灵活性以及所得积分方程的特性，使其成为未来研究更困难的非线性问题（这些问题通常涉及目前极难求解的超越项）的有力候选方案。通过证明这种广义变换可以有效地处理延迟方程，这项研究为更精确地模拟那些“过去不断塑造现在”的系统打开了大门。

技术摘要：利用偏移线性正则变换分析时滞微分方程

问题陈述
时滞微分方程（DDEs）用于模拟当前状态依赖于过去状态的系统，这种现象在生物学、经济学和工程学中广泛存在。与常微分方程不同，DDEs 会表现出复杂的动力学特性，包括由时间延迟引起的稳定性问题和振荡行为。虽然拉普拉斯变换和傅里叶变换（FT）已被应用于求解线性 DDEs，但它们往往面临局限性。具体而言，由于计算闭式形式的反变换具有挑战性，通过拉普拉斯变换获取复杂时滞系统的显式解往往十分困难；此外，时滞项经常导致超越特征方程的出现。尽管 Ohira 等人 [13,14] 已成功利用傅里叶分析求解了特定的 DDEs，但对于这类问题，具有更多自由度的更广义积分变换的潜力在很大程度上仍未得到充分探索。

方法论
作者提出了一个基于偏移线性正则变换（OLCT）的框架，OLCT 是一种具有六个参数的积分变换，它推广了傅里叶变换（FT）、分数傅里叶变换（FrFT）和线性正则变换（LCT）。该方法论通过两种主要途径进行：

直接 OLCT 方法：
- 作者将 OLCT 应用于一类形如 $\frac{d}{dt}X(t) + \alpha t X(t) = \beta X(t - \tau)$ 的线性时滞微分方程。
- 通过利用 OLCT 的运算性质（特别是平移和微分性质），将原始 DDE 转换为变换域中的泛函微分方程。
- 随后，利用积分因子将该方程转换为 Volterra 型时滞积分方程。
- 由于存在平移项 $OM(u - a\tau)$ ，通常无法获得闭式解析解。因此，作者采用 Brunner 的步进法（method of steps） 在连续的区间内迭代构建数值解。
- 最终在时域内的解通过逆 OLCT 恢复。
基于傅里叶变换的 OLCT 方法：
- 为了验证直接法并建立统一的视角，作者通过利用 OLCT 与 FT 之间的已知关系推导出了解。
- 他们在傅里叶域中求解 DDE（遵循 Ohira 的方法论），并将结果映射到 OLCT 域。
- 对于参数 $a=0$ 的特殊情况，他们在 OLCT 域内推导出了一个显式解析解，该解随后表示为涉及历史函数的积分。

核心贡献

框架开发： 本文建立了第一个基于 OLCT 的求解时滞微分方程的框架，将其转化为 Volterra 型时滞积分方程。
统一表述： 文中证明了经典的基于傅里叶的表述是所提 OLCT 框架的一个特例，从而提供了一个统一的变换域视角。
解析与数值解： 作者针对特殊情况（ $a=0$ ）推导了显式解析解，并针对一般情况提供了稳健的数值方案（步进法）。
参数分析： 研究调查了关键参数——特别是时滞 $\tau$ 、系数 $\beta$ 和阻尼参数 $\alpha$ ——对解的行为（包括振荡、衰减和稳定性）的影响。

结果

一致性： 数值模拟和图形对比显示，在各种时延值（ $\tau = 4, 10, 16$ ）下，直接 OLCT 解与基于傅里叶的 OLCT 解在实部和虚部上几乎完全一致。仅在振荡峰值附近观察到微小差异，且这些差异随着解的衰减而消失。
参数效应：
- 时延 ( $\tau$ )： 增加时延参数会降低振荡幅度，并产生更稳定、更平滑的解。
- 系数 ( $\beta$ )： $\beta$ 的符号显著影响小时延下的解（负 $\beta$ 与解相反，导致更快衰减；正 $\beta$ 则支持解的发展）。对于大时延，解的幅度对 $\beta$ 符号的敏感度降低。
- 阻尼 ( $\alpha$ )： 较大的 $\alpha$ 值会导致更快的衰减和更平滑的解，而较小的 $\alpha$ 值则保留时延诱导的振荡。
与 Ohira 研究的比较： 与 Ohira 等人 [14] 的经典傅里叶变换解相比，OLCT 解在实部上高度吻合。然而，由于其复值核的特性，OLCT 解还提供了虚部成分，从而提供了纯实数傅里叶解所不具备的额外相位和偏移信息。

意义与主张
本文声称，所提出的 OLCT 框架通过利用 OLCT 增强的灵活性（时间平移、频率平移、缩放和啁啾调制），为分析时滞微分方程提供了极具前景的基础。作者断言，由于误差分析和与现有基于傅里叶结果的一致性验证，该方法对于线性 DDEs 是准确且可靠的。

至关重要的是，作者指出这项工作是实现更广泛目标的一个初步步骤：即求解包含超越项的非线性时滞微分方程。他们认为，OLCT 的卷积和相关性质，结合所得 Volterra 积分方程的结构，具备处理此类非线性的必要属性。然而，本文明确限定了目前的研究范围仅限于线性方程，并未展示非线性系统的结果，并指出将该方法扩展到此类问题将是未来研究的重点。

类似论文