追赶法\Thomas算法在常微分方程两点边值问题的差分解法，紧差分解法中的应用（附python代码实现及作图）

最新推荐文章于 2025-03-23 22:18:42 发布

感受野_

最新推荐文章于 2025-03-23 22:18:42 发布

阅读量2k

点赞数

分类专栏：题目记录文章标签：算法 python 矩阵线性代数

本文链接：https://blog.youkuaiyun.com/carpediem_sisy/article/details/123490399

版权

本文介绍了追赶法和Thomas算法在求解常微分方程两点边值问题的差分解法中的应用，详细分析了相关数学模型，并提供了Python代码实现及图形化展示。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

偏微分数值解法课程作业系列

Chapter 1 常微分方程两点边值问题的差分解法

题目引入：

也可为
注：差分法紧差分法最终实现形式都为追赶法求解 代码通用

分析：

6题分析

第7题紧差分法分析

下附上第七题紧差分法代码实现：

import matplotlib.pyplot as plt
import numpy as np

def thomas(N,a,b,c,d):
    l = np.zeros((N+1,1))
    u = np.zeros((N+1,1))
    x = np.zeros((N+1,1))
    y = np.zeros((N+1,1))
    u[0] = b[0]
    for i in range(N):
        l[i+1] = a[i+1]/u[i]
        u[i+1] = b[i+1