刚性常微分方程的有限元方法

《AppliedMath》：The Finite Element Method for Stiff Ordinary Differential Equations Yanhui Ding, Qiong Tang and Sijia Tang

　　

来自世界各地大学的教授、学生和研究人员使用根据许可分发的软件进行数值模拟，这需要一台具有相当硬件能力的计算机。这意味着在需要动态建模的工程应用（尤其是机器人和非线性控制）中，仿真成本高昂。本文比较并分析了一种基于低成本、免费开源技术的节俭仿真方案的性能，具体使用低功耗单板微型计算机（Raspberry Pi）配合GNU-Octave软件。基准测试是对选择性顺应性装配机器人手臂（Selective Conformal Assembly Robot Arm, SCARA）关节空间轨迹跟踪控制的数值模拟。为执行此任务，研究人员使用称为常微分方程求解器（ODE solver）的数值方法求解矩阵形式的耦合非线性微分方程组，该解包含控制律和来自欧拉-拉格朗日形式（Euler-Lagrange formalism）的动态系统模型。研究人员分析时间复杂度和精度，以将节俭仿真工具与由个人计算机和运行相同仿真代码的MATLAB^TM组成的常规仿真设置进行比较。分析显示数值解的偏差最小且时间复杂度合理。此外，该方法的节俭分数和仿真工具的低采购成本使得预算有限的高校能够创建用于教育和研究的仿真实验室。

**论文解读文章**

**研究背景与问题**

在工程教育中，仿真软件常用于辅助学习和研究，尤其是在机器人学、非线性控制等需要动态建模的领域。然而，传统仿真方案依赖高昂的许可软件（如MATLAB^TM）和高性能计算机，导致许多学生和研究者难以承担，尤其在疫情后全球经济危机背景下，数字工具缺乏的问题更加突出。为此，学术界引入节俭创新（Frugal Innovation, FI）概念，强调“以更少做更多”，通过低成本技术实现可接受的仿真质量。先前的工作从开源软件（Open-Source Software, OSS，如GNU-Octave、Python）和开源硬件（Open-Source Hardware, OSH，如Raspberry Pi）两个角度探索解决方案，但缺乏将两者结合并系统比较其与常规方案的性能。本研究旨在提出一种基于Raspberry Pi（RPi）与GNU-Octave的节俭仿真方案，并以选择性顺应性装配机器人手臂（SCARA）的轨迹跟踪控制为基准，分析其与传统MATLAB^TM加笔记本方案的性能差异，验证其在教育和研究中的可行性。论文发表在《AppliedMath》。

**关键技术与方法**

研究人员采用以下关键技术：（1）开源硬件与软件组合：使用搭载基于Debian操作系统的低功耗单板微型计算机Raspberry Pi 3 B+（型号RPi 3B+），配合开源数值计算软件GNU-Octave 4.0.3。该软件具有社区维护的包（Octave Packages），支持控制、通信等领域的计算，但需自主实现非线性控制律的数值逼近算法。（2）基准测试案例：以4自由度（degrees of freedom, dof）的SCARA机器人为对象，通过欧拉-拉格朗日形式建立其动力学模型，得到惯性矩阵（M矩阵）、科里奥利力和离心力矩阵（C矩阵）及重力向量（G向量），并使用基于李雅普诺夫方法的计算力矩控制律实现关节空间轨迹跟踪。（3）数值求解算法：采用ode45求解器，基于龙格-库塔和Dormand-Prince方法，将二阶微分方程组改写为一阶耦合常微分方程组（ODE）并迭代求解。代码文件中定义了动力学参数、控制器增益（kp=50，kd=6）、期望轨迹（余弦形式）及初始条件，并设置相对容差（RelTol=1e-4）和绝对容差（AbsTol=1e-3）。（4）性能评估指标：计算时间复杂度（以tic/toc函数测量10次连续执行的平均墙钟时间）和均值标准误差（Standard Error of the Mean, α），并与常规方案（MATLAB^TM加笔记本，配置见文末说明）进行对比。此外，提出节俭分数（Frugality Score）公式，综合精度和时间复杂度，动态权重w可调（0~1），用于量化节俭方案的优劣。样本队列来源：仿真参数中机器人参数来自文献中的SCARA模型，无真实样本。

**研究结果**

**5.1. 时间复杂度与均值标准误差α**
通过对比实验，研究人员测得节俭方案（GNU-Octave + RPi 3B+）的平均时间复杂度为7.32秒，而常规方案（MATLAB^TM + 笔记本）为0.86秒，GNU-Octave + 笔记本为1.081秒。节俭方案的速度约为常规方案的8.5倍慢。但通过计算各变量（关节位置和速度）的均值标准误差α，发现节俭方案与参考方案（MATLAB^TM或GNU-Octave笔记本）之间差异极小（例如关节位置变量的α值在10^-4至10^-5量级），说明数值精度高。

**5.2. 节俭分数**
研究人员利用公式Frugality_lc = α_all - w + 1 / (t_lc - t_hc)计算节俭分数，其中α_all为所有变量α的均值，t_lc为节俭方案时间复杂度，t_hc为常规方案时间复杂度。通过改变权重w（0到1），绘制节俭线。结果显示：当w=0时，节俭分数仅由精度决定；当w=1时，分数权衡时间差。假设案例表明，若节俭方案时间小于常规方案，节俭线斜率为正（如案例5），而实际节俭方案（案例1和2）斜率为负，但精度损失可忽略。节俭分数提供了一个量化比较的基线。

**讨论与结论**

在讨论部分，研究人员指出速度误差的α值高于位置误差的原因在于ode45求解器将二阶系统降为一阶时，数组元素顺序（速度在前，加速度在后）导致内存管理问题，尤其在RPi上有限的内存资源下效率降低。尽管节俭方案时间成本是常规方案的7倍，但其硬件成本（RPi 3B+约35美元）和软件成本（GNU-Octave免费）远低于常规方案（笔记本约669美元，MATLAB^TM许可约2250美元），且代码无需额外学习，因此在预算受限的大学或家庭场景中具有显著社会价值。

**研究结论翻译**
在这项研究中，研究人员将基于低成本开源技术的节俭仿真工具与其他相关技术进行了比较，结果验证了其在涉及机器人系统非线性控制的教育或研究应用中的可用性。此外，仿真性能分析表明各仿真方案之间差异极小。这证明了所提出方法的节俭性，即在限制性环境中利用技术“以少博多”来获得灵活性。因此，所提出的仿真工具具有社会影响，因为它可以在预算有限的高校中实施，并且每个学生或研究人员都能在家中建造自己的仿真实验室，用于在线教育、补充研究或科学研究。本研究受限于所用计算机设备的技术规格。Raspberry Pi微型计算机与其他开发板一样在持续演进；更强的硬件能力将提升性能，减少硬件在环应用的运算时间。因此，后续工作将继续通过实现一个在Raspberry Pi上运行的非线性控制律，来设计用于SCARA样机轨迹跟踪的低成本实验室设备。此外，预计将设计一个连接互联网的图形用户界面（Graphical User Interface, GUI），用于实时远程控制应用。