[论文解读] Identification of Fractional-Order Dynamical Systems

Ľ. Dorčák, V. Lesko+1 more|ArXiv.org|Apr 15, 2002

Advanced Control Systems Design参考文献 1被引用 29

一句话总结

本文提出了一种新颖的方法，通过结合分数阶导数的数值逼近与传递函数微分技术，识别分数阶动力系统。该方法能够利用测量的输入-输出数据，同时估计系统参数与导数阶次，实现了在仿真系统和真实系统（包括实验室加热炉）中均具有高精度的参数识别。

ABSTRACT

This contribution deals with identification of fractional-order dynamical systems. We consider systems whose mathematical description is a three-member differential equation in which the orders of derivatives can be real numbers. We give a discretization method and a numerical solution of differential equations of this type. An experimental method of identification is given which is based on evaluation of transfer characteristics. This is a combination of the method of derivatives of transfer characteristics and of the method of passive search. The verification was performed on systems with known parameters and also on a laboratory object.

研究动机与目标

• 开发一种可靠的方法，用于识别导数阶次为实数（而非整数）的分数阶动力系统。
• 解决传统整数阶建模忽略真实系统分数阶动力特性的局限性。
• 实现从实验输入-输出数据中准确识别系统参数与导数阶次（α, β）。
• 在具有已知参数的系统及一个真实实验室对象上验证该方法，证明其在拟合精度上优于整数阶近似。

提出的方法

• 采用三元分数阶微分方程模型：$ a_2 y^{(\alpha)}(t) + a_1 y^{(\beta)}(t) + a_0 y(t) = u(t) $，其中 α 和 β 为实数。
• 应用格伦沃尔德-莱特尼科夫公式对分数阶导数进行数值逼近：$ y^{(\alpha)}(t) \approx h^{-\alpha} \sum_{j=0}^{N(t)} b_j y(t-jh) $，其中二项式系数 $ b_j = (-1)^j \binom{\alpha}{j} $。
• 推导出一个显式递推公式（公式5），用于利用离散时间值与时间步长 $ h $ 数值求解系统。
• 采用最小二乘误差准则（公式6），最小化模型预测与实测数据在时间上的差异。
• 通过求解由梯度条件 $ \partial E/\partial \bar{a} = 0 $ 导出的三元线性方程组（公式8），利用求和近似离散积分。
• 采用基于二分法的迭代过程（公式10–11），在指定区间内通过最小化误差泛函，对估计的阶次 α 和 β 进行优化。

实验结果

研究问题

• RQ1 当系统参数与导数阶次均未知时，能否准确识别分数阶系统？
• RQ2 与经典整数阶建模相比，该方法在真实动力系统中的拟合精度与准确性如何？
• RQ3 测量精度与时步对分数阶系统识别性能有何影响？
• RQ4 该方法能否从实验数据中可靠估计非整数阶导数阶次（α, β），即使真实阶次初始未知？

主要发现

• 对于一个整数阶系统（α=2，β=1），该方法恢复的参数为 $ a_2=1.0001 $，$ a_1=2.99987 $，$ a_0=1.99998 $，估计阶次为 α=1.99993，β=0.99996，与真实值高度一致。
• 对于一个分数阶系统（α=2.2，β=0.9），该方法估计的参数为 $ a_2=0.80005 $，$ a_1=0.49996 $，$ a_0=0.99998 $，估计阶次为 α=2.19996，β=0.89989，与真实参数非常接近。
• 当将同一分数阶系统近似为二阶整数阶模型时，该方法得到 $ a_2=0.76639 $，$ a_1=0.23184 $，$ a_0=1 $，其响应明显不同，表明整数阶近似存在模型缺陷。
• 对于一个真实加热炉，该方法识别出分数阶模型参数为 $ a_2=-14994.3 $，$ a_1=6009.52 $，$ a_0=1.69 $，α=1.31，β=0.97，准则值为 $ 2.7 \times 10^{-4} $，显著优于整数阶近似（$ 1.02 \times 10^{-3} $）。
• 该方法表现出鲁棒性，成功识别出二元模型 $ a_1=788.35 $，$ a_0=1.39 $，β=0.73，准则值为 $ 6.3 \times 10^{-4} $，仍优于整数阶拟合结果。