一、前言

1.在数值计算中，如三次样条插值或用差分法求解微分方程边值问题时会遇到三对角方程组。现在来讲一讲求解方法。

2.本文借鉴了https://www.bilibili.com/video/BV1aA41147Ab

二、方程组的特点

1.类型

2.要求

3.方法的优势

当矩阵比较大的时候用这种方法会快很多

三、算法步骤

四、例题+代码

1.用追赶法求解以下五阶方程组

2.代码

import numpy as npA = np.array([[4, 1, 0, 0, 0],[1, 4, 1, 0, 0],[0, 1, 4, 1, 0],[0, 0, 1, 4, 1],[0, 0, 0, 1, 4]])
f = np.array([2, 1, 1, 1, 2])
# 一、左乘A逆
print('传统方法：', np.linalg.inv(A) @ f)# 二、追赶法
a = np.array([0, 1, 1, 1, 1], dtype=float)
b = np.array([4, 4, 4, 4, 4], dtype=float)
c = np.array([1, 1, 1, 1, 0], dtype=float)
y = np.zeros(5)
B = np.zeros(5)
x = np.zeros(5)
y[0] = f[0] / b[0]
d = b[0]
for i in range(1, 5):B[i - 1] = c[i - 1] / dd = b[i] - a[i] * B[i - 1]y[i] = (f[i] - a[i] * y[i - 1]) / d
x[4] = y[4]for i in range(3, -1, -1):x[i] = y[i] - B[i] * x[i + 1]
print('追赶法：', x)

求解结果：

传统方法： [0.48076923 0.07692308 0.21153846 0.07692308 0.48076923]
追赶法： [0.48076923 0.07692308 0.21153846 0.07692308 0.48076923]

【数学建模之Python】12.追赶法求解三对角方程组相关推荐

追赶法求解三对角方程组
1. 来源和背景对于一个(主)三对角方程组,我们常用"追赶法"来进行求解. 而三对角方程组常常出现于微分方程的数值求解,例如热传导方程的边值问题 {y′′(x)=f(x,y,y′ ...
用追赶法求方程组c语言程序,数值计算——追赶法求解三对角方程组（附代码）...
目录追赶法基础理论追赶法c++程序代码程序运行结果源码文件下载地址追赶法基础理论在数值计算中,对三次样条曲线插值和用差分方法求解常微分方程边值问题时,通常会遇到Ax=d三对角形式的方程组: ...
Python解线性方程组的直接法(6)————求解三对角方程组的追赶法
求解三对角方程组的追赶法 import numpy as npdef zuiganfa(A, d):n = A.shape[0]l = np.mat(np.zeros(n, dtype=float
数学建模竞赛用python_2018全国中学生数学建模与Python编程夏令营
全国中学生数学建模与Python编程夏令营招生简章为什么要学习数学建模? 1)国务院关于新一代人工智能发展<规划>:在中小学阶段设置人工智能相关课程,实施全民智能教育项目,在中小学阶段设 ...
python数学建模大赛_2018全国中学生数学建模与Python编程冬令营
全国中学生数学建模与Python编程冬令营招生简章为什么要学习数学建模? 1)国务院关于新一代人工智能发展<规划>:在中小学阶段设置人工智能相关课程,实施全民智能教育项目,在中小学阶段设 ...
利用matlab求解线性规划,数学建模讲座之三利用matlab求解线性规划问题（linprog函数）...
数学建模讲座之三利用matlab求解线性规划问题(linprog函数) 利用利用 Matlab求解线性规划问题求解线性规划问题河北科技河北科技大学大学*第第 1页页利用 Matlab求解线性规划问题 ...
追赶法求解三对角线性方程组的MATLAB程序
function [M] = after_method(a,b,c,g,tolerance) %% 追赶法求解三对角矩阵 % a为三对角矩阵左下对角线上的值 % b为三对角矩阵中间对角线上的值 % c ...
备战数学建模（Python）
备战数学建模(Python) Python之建模规划篇 Python之建模数值逼近篇 Python之建模微分方程篇由于美国大学生数学建模大赛很快就要开赛了,所以我就打算在这几天内,好好的看看< ...
线性规划问题的数学建模matlab,数学建模讲座之三——利用Matlab求解线性规划问题（linprog函数）.ppt...
数学建模讲座之三--利用Matlab求解线性规划问题(linprog函数) 利用Matlab求解线性规划问题线性规划是一种优化方法,Matlab优化工具箱中有现成函数linprog对如下式描述的LP ...

【数学建模之Python】12.追赶法求解三对角方程组

一、前言

二、方程组的特点

1.类型

2.要求

3.方法的优势

三、算法步骤

四、例题+代码

1.用追赶法求解以下五阶方程组

2.代码

【数学建模之Python】12.追赶法求解三对角方程组相关推荐

最新文章

热门文章