动力学分析基础--雅克比矩阵.pdfVIP

动力学分析基础——雅可比矩阵

代码编写,资料整理——ZH1110

动力学仿真计算归结为对典型的常微分方程组的初值问题。在解上述的初值问题时，除了应用常微分方程初值问题的数值积分

外，还将用到求解线性代数方程组的数值方法,所以首先我们必须先研究这两个常用的计算机算法,已便于后面的计算.

高斯消去法求解线性代数方程组(直接法,即消去法)，已在线性代数课程中有详细的讨论，在此给出些说明以及具体的算法描述。

大致可以分为以下两步。

1．将系数矩阵经过一系列的初等行变换（归一化）

在变换过程中，采用原地工作，即经变换后的元素仍放在原来的位置上。

2.消去。它的作用是将主对角线以下的均消成0，而其它元素与向量中的元素也应作相应的变换

最后，进行回代依次解出

如:我们要解如下方程组：

初等行变换: