1 基本用法

scipy.linalg.solve_triangular(a, b, trans=0, lower=False, unit_diagonal=False, overwrite_b=False, debug=None, check_finite=True)

解决ax=b方程中的x,(假定a是一个上/下三角矩阵)

2 参数介绍

a

(M, M) array_like

三角矩阵a

b (M,) or (M, N) array_like
lower

bool, optional

如果是True的话,那么a是一个下三角矩阵;否则是一个上三角矩阵(默认也是一个上三角矩阵)

trans

{0, 1, 2, ‘N’, ‘T’, ‘C’}, optional

T是转置,H是共轭转置

unit_diagonal

bool, optional

如果为True,则a的对角元素被假定为1并且不被引用。

overwrite_b

bool, optional

如果为True,结果写入原来的b

check_finite

bool, optional

是否检查输入矩阵只包含有限的数。禁用可能会提高性能,但如果输入确实包含无限或nan,则可能导致问题(崩溃、未终止)。

3 使用举例

from scipy.linalg import solve_triangular
a = np.array([[3, 0, 0, 0], [2, 1, 0, 0], [1, 0, 1, 0], [1, 1, 1, 1]])
b = np.array([4, 2, 4, 2])
x=solve_triangular(a,b,lower=True)
print(x)
#[ 1.33333333 -0.66666667  2.66666667 -1.33333333]
print(a @ x)
#[4. 2. 4. 2.]

如果我们设置unit_diagonal=True的话,即使原来对角线有值,也会变成1来进行计算

from scipy.linalg import solve_triangular
a = np.array([[3, 0, 0, 0], [2, 1, 0, 0], [1, 0, 1, 0], [1, 1, 1, 1]])
b = np.array([4, 2, 4, 2])
x=solve_triangular(a,b,lower=True,unit_diagonal=True)
print(x)
#[ 4. -6.  0.  4.]
print(a @ x)
#[12.  2.  4.  2.]
a1 = np.array([[1, 0, 0, 0], [2, 1, 0, 0], [1, 0, 1, 0], [1, 1, 1, 1]])
print(a1  @ x)
#[4. 2. 4. 2.]

scipy 笔记:solve_triangular相关推荐

  1. scipy笔记:wishart (威沙特分布)

    威沙特分布理论部分:概率统计笔记:威沙特分布(Wishart Distribution)_UQI-LIUWJ的博客-CSDN博客 1 使用方法 scipy.stats.wishart(df=None, ...

  2. scipy笔记:scipy.sparse

    1 稀疏矩阵介绍 在networkx包中,很多运算返回的是sparse matrix(如nx.laplacian_matrix),这是稀疏矩阵格式.隶属于scipy.sparse import net ...

  3. scipy笔记—scipy.misc.imresize用法(方便训练图像数据)

    scipy.misc.imresize 不同于普通的reshape, imresize不是单纯的改变图像矩阵的维度,而是能将图片重采样为指定像素,这样给深度学习中训练图像数据带来方便. import ...

  4. numpy笔记整理 multivariate_normal(多元正态分布采样)

    1 基本用法 np.random.multivariate_normal(mean, cov, size=None, check_valid=None, tol=None) 根据均值和协方差矩阵的情况 ...

  5. DeepSORT 多目标跟踪算法笔记

    SORT 是一种实用的多目标跟踪算法,然而由于现实中目标运动多变且遮挡频繁,该算法的身份转换(Identity Switches)次数较高.DeepSORT 整合外观信息使得身份转换的数量减少了45% ...

  6. pytorch笔记:torch.sparse类

    PyTorch 提供了 torch.Tensor 来表示一个包含单一数据类型元素的多维数组. 默认情况下,数组元素连续存储在内存中,从而可以有效地实现各种数组处理算法,这些算法依赖于对数组元素的快速访 ...

  7. 从零开始学习Deep SORT+YOLO V3进行多目标跟踪(附代码)

    点击上方"小白学视觉",选择加"星标"或"置顶" 重磅干货,第一时间送达 1.跟踪基础知识简介 首先要说明一点,现在多目标跟踪算法的效果,与 ...

  8. 【项目实践】从零开始学习Deep SORT+YOLO V3进行多目标跟踪(附注释项目代码)...

    点击上方"小白学视觉",选择加"星标"或"置顶" 重磅干货,第一时间送达本文转自|集智书童 1.跟踪基础知识简介 首先要说明一点,现在多目标 ...

  9. 【汇总记忆:命令、语法、函数、等】

    目录 Linux常用操作 Python语法 numpy&scipy语法 Python常用 函数 和 操作 Linux常用操作 查看硬件信息: 磁盘空间:df -h 显卡型号和个数(没装驱动时) ...

最新文章

  1. java 怎么让一个jlabel和一个jtextfield在一行_今天小哥大家分享一下Java编程语言的第一个程序应该怎么样去写...
  2. python3的pyqt5 qtablewidget按数据大小排列_InnoDB为什么要选择B+树来存储数据?
  3. 用公式求e^x的近似值
  4. “基因编辑婴儿”惹争议,那你知道机器学习在脱靶效应中的作用吗?
  5. vscode 模拟本地服务器打开文件,设置默认浏览器为chrome
  6. 剑指offer JZ02/05 替换空格 Python
  7. pythonjson格式化输出_pythonjson格式化输出_Python json格式化打印实现过程解析
  8. 环境土壤物理模型HYDRUS1D/2D/3D实践技术
  9. 心理学与生活 - 文化与社会
  10. 使用vue-router+vuex进行导航守卫(转)
  11. 基于Mendix的云上全流程透明性备品备件协同管理
  12. 高版本node,无法编译程序,删除所有文件
  13. 怎样解决ABBYY FineReader 12识别字体字符错误
  14. starting Intent from ProcessRecord with revoked permission android.permission.CALL_PHONE 的错误原因
  15. 步进电机驱动程序,5线
  16. 毕设 电脑维修_入库出库结算论文
  17. Monkey常见面试问题
  18. Silk这个网站是用来在线绘制对称图形的
  19. 学习通考试刷题(ocs网课助手)
  20. OpenStack ironic裸金属部署(裸金属作为独立服务)

热门文章

  1. Delete Volume 操作 - 每天5分钟玩转 OpenStack(57)
  2. 进程外Session和进程内Session存储
  3. 自定义的 ListBoxItem 自适应ListBox的宽度
  4. 关于“习惯”的精彩分析
  5. GridView的多主键(Key)取值问题
  6. 如何解决Beyond Compare内容相同仍然标示红色
  7. 写在开年:移植wolfssl4.3.0到w60x_sdk_3.04时的一点问题
  8. lwip协议栈中超时定时器实现原理
  9. 通信原理2习题课汇总(随机信号、信道、最佳接收)
  10. 吴恩达机器学习Week4神经网络表述