【PLA】基于Python实现的线性代数算法库之lu分解
2022/2/6 20:14:30
本文主要是介绍【PLA】基于Python实现的线性代数算法库之lu分解,对大家解决编程问题具有一定的参考价值,需要的程序猿们随着小编来一起学习吧!
【PLA】基于Python实现的线性代数算法库之lu分解
算法包下载链接:https://download.csdn.net/download/qq_42629529/79481514
from .Matrix import Matrix from .Vector import Vector from ._globals import is_zero#判断是否为0 #方阵的lu分解 def lu(matrix): assert matrix.row_num() == matrix.col_num(), "matrix must be a square matrix" n = matrix.row_num()#行数 A = [matrix.row_vector(i) for i in range(n)]#矩阵的行向量,上三角矩阵 L = [[1.0 if i == j else 0.0 for i in range(n)] for j in range(n)] #高斯消元的过程 for i in range(n): # 看A[i][i]位置是否可以是主元 if is_zero(A[i][i]): return None, None#不能分解两个都空 else: for j in range(i + 1, n): p = A[j][i] / A[i][i] A[j] = A[j] - p * A[i] L[j][i] = p return Matrix(L), Matrix([A[i].underlying_list() for i in range(n)])
这篇关于【PLA】基于Python实现的线性代数算法库之lu分解的文章就介绍到这儿,希望我们推荐的文章对大家有所帮助,也希望大家多多支持为之网!
- 2024-11-02Python编程基础
- 2024-11-01Python 基础教程
- 2024-11-01用Python探索可解与不可解方程的问题
- 2024-11-01Python编程入门指南
- 2024-11-01Python编程基础知识
- 2024-11-01Python编程基础
- 2024-10-31Python基础入门:理解变量与数据类型
- 2024-10-30Python股票自动化交易资料详解与实战指南
- 2024-10-30Python入行:新手必读的Python编程入门指南
- 2024-10-30Python入行:初学者必备的编程指南