[论文解读] Optimal stabilization of a cycle in nonlinear discrete systems

Dmitriy Dmitrishin, Anna Khamitova+2 more|arXiv (Cornell University)|Jul 28, 2013

Chaos control and synchronization参考文献 5被引用 4

一句话总结

本文提出了一种最优延迟反馈控制（DFC），用于在最小预历史深度下稳定非线性离散系统中的两周期轨道。通过将控制表示为过去状态差分的线性组合，并利用费耶尔多项式，作者推导出唯一系数，以最小化所需记忆深度，确保在乘数范围为 (−μ*, −1) 内的负乘数周期轨道的渐近稳定性。

ABSTRACT

A method of stabilizing 2-cycles in discrete dynamic systems by Delayed Feedback Control is developed by using classic Harmonic Analysis.

研究动机与目标

• 解决在小而有界的控制下，稳定非线性离散系统中不稳定两周期轨道的挑战。
• 在确保周期轨道渐近稳定的同时，最小化延迟反馈控制（DFC）所需的预历史深度。
• 将先前关于平衡点稳定性的结果扩展至周期轨道，特别是那些具有负乘数（表明混沌行为）的周期轨道。
• 推导出显式的最优控制系数，以最小化控制复杂度（即最小化 N*），针对给定的周期乘数上界 μ*。
• 建立控制设计与极值三角多项式及费耶尔核之间理论框架的联系。

提出的方法

• 将系统建模为非线性离散映射 $x_{n+1} = f_h(x_n)$，其具有一个非稳定的两周期轨道 $(\eta_1, \eta_2)$ 和乘数 $\mu = f_h'(\eta_1)f_h'(\eta_2) \in (-\mu^*, -1)$。
• 设计一种延迟反馈控制形式 $u = -\sum_{j=1}^{N-1} \varepsilon_j (f_h(x_{n-2j+2}) - f_h(x_{n-2j}))$，简化为线性组合 $x_{n+1} = \sum_{j=1}^N a_j f_h(x_{n-2j+2})$，其中系数为 $a_j$。
• 在线性化周期轨道周围后，推导出关于 $\lambda$ 的特征方程，得到 $\lambda - \mu p(\lambda)^2 = 0$，其中 $p(\lambda) = \sum_{j=1}^N a_j \lambda^{-(j-1)}$。
• 通过令 $\nu = \lambda^{1/2}$，将特征方程转换为三角形式，得到一个关于 $\nu$ 的多项式，其系数与 $\mu$ 相关。
• 利用费耶尔多项式的极值性质，最小化控制谱的最大偏差，从而得到最优系数 $a_j^0 = \frac{2(N-j)+1}{N^2}$。
• 通过 $\varepsilon_j = \sum_{k=j+1}^N a_k^0$ 从 $a_j^0$ 推导出最优 $\varepsilon_j$，确保最小预历史深度 $N^* = 2(N_0 - 1)$。

实验结果

研究问题

• RQ1 在非线性离散系统中，对于乘数位于 $(−\mu^*, −1)$ 内的负乘数两周期轨道，实现最优稳定所需的最小预历史深度是多少？
• RQ2 如何选择控制系数，以在确保周期轨道渐近稳定的同时最小化内存需求？
• RQ3 当周期乘数远离零且为负时，延迟反馈控制的最优结构是什么？
• RQ4 能否利用极值三角多项式推导出最优控制？其与费耶尔核有何关联？
• RQ5 在何种条件下控制无法稳定周期轨道？为何在逻辑斯蒂映射的 $h=4$ 时控制失效？

主要发现

• 最优控制系数由 $a_j^0 = \frac{2(N-j)+1}{N^2}$ 唯一确定，该表达式源自平均值为 $1/N$ 的极值费耶尔多项式。
• 最小预历史深度为 $N^* = 2(N_0 - 1)$，其中 $N_0$ 是满足 $N_0 > \sqrt{\mu^*}$ 的最小整数，确保对所有 $\mu \in (-\mu^*, -1)$ 实现稳定性。
• 对于 $\mu^* = 4$ 的逻辑斯蒂映射，$N_0 = 3$，因此 $N^* = 4$，最优系数为 $\varepsilon_1 = \frac{4}{9}, \varepsilon_2 = \frac{1}{9}$。
• 在 $h=4$ 时控制失效，是因为乘数 $\mu = -4$ 恰好位于稳定区域的边界，且控制谱触及单位圆。
• 当 $\mu^* = 1$ 时，$N_0 = 2$，因此 $N^* = 2$，最优控制为 $u = -\frac{1}{4}(f_h(x_n) - f_h(x_{n-2}))$，可在 $h < 4$ 时稳定周期轨道，但在 $h=4$ 时不成立。
• 数值模拟表明，控制在 $h=3.95$ 时可稳定周期轨道，解收敛至周期轨道的邻域，验证了理论结果。