如何使用 Numpy\Python 求解三对角矩阵

Question

我正在做一个项目，我想在 20 个步骤中求解一个包含 9 个方程和 9 个未知数的三对角矩阵系统。如何使用 numpy\python 处理代码。

这是一个例子 Tridiagonal Matrix

我正在求解 x11、x12、x13、x21、x22、x23、x31、x32、x33 我根本没有尝试过任何代码。我刚开始学习python。如果我能得到上面矩阵的代码，我就可以为我正在做的项目调出代码。

Answer 1

第一步：前向扫描消除次对角线；

第 2 步：执行反向替换以消除上对角线。

这是托马斯算法。