Python实现生成对角矩阵和对角块矩阵

更新时间：2023年04月13日 08:31:36 作者：微小冷

这篇文章主要为大家详细介绍了如何利用Python实现生成对角矩阵和对角块矩阵，文中的示例代码讲解详细，感兴趣的小伙伴可以跟随小编一起了解一下

对角矩阵

scipy中的函数

在scipy.linalg中，通过tri(N, M=None, k=0, dtype=None)可生成N×M对角矩阵，若M=None，则M默认为N。k表示矩阵中用1填充的次对角线个数。

print(tri(3,5,2,dtype=int))
'''
[[1 1 1 0 0]
 [1 1 1 1 0]
 [1 1 1 1 1]]
'''

在numpy中也提供了多种对角矩阵生成函数，包括diag, diagflat, tri, tril, triu等，

numpy.diagflat

diagflat用于生成对角矩阵，diag在diagflat基础上，添加了提取对角元素的功能，例如

>>> np.diagflat([1,2,3])
array([[1, 0, 0],
       [0, 2, 0],
       [0, 0, 3]])
>>> np.diag([1,2,3])
array([[1, 0, 0],
       [0, 2, 0],
       [0, 0, 3]])
>>> np.diag(np.ones([3,3])) #提取对角元素
array([1., 1., 1.])

numpy.tri

tri(M,N,k)用于生成M行N列的三角阵，其元素为0或者1，k用于调节0和1的分界线相对于对角线的位置，例如

>>> np.tri(3,5,1)
array([[1., 1., 0., 0., 0.],
       [1., 1., 1., 0., 0.],
       [1., 1., 1., 1., 0.]])
>>> np.tri(3,5,2)
array([[1., 1., 1., 0., 0.],
       [1., 1., 1., 1., 0.],
       [1., 1., 1., 1., 1.]])
>>> np.tri(3,5,3)
array([[1., 1., 1., 1., 0.],
       [1., 1., 1., 1., 1.],
       [1., 1., 1., 1., 1.]])

tril, triu可用于提取出矩阵的左下和右上的三角阵，其输入参数除了待提取矩阵之外，另一个参数与tri中的k相同。

x = np.arange(12).reshape(4,3)
>>> np.tril(x,-1)
array([[ 0,  0,  0],
       [ 3,  0,  0],
       [ 6,  7,  0],
       [ 9, 10, 11]])
>>> np.triu(x,-1)
array([[ 0,  1,  2],
       [ 3,  4,  5],
       [ 0,  7,  8],
       [ 0,  0, 11]])

对角块矩阵

对于scipy.linalg.block_diag(A,B,C)而言，会生成如下形式矩阵

from scipy.linalg import *
import numpy as np
A = np.ones([2,2])
B = np.round(np.random.rand(3,3),2)
C = np.diag([1,2,3])
bd = block_diag(A,B,C)
print(bd)
'''
[[1.   1.   0.   0.   0.   0.   0.   0.  ]
 [1.   1.   0.   0.   0.   0.   0.   0.  ]
 [0.   0.   0.8  0.38 0.41 0.   0.   0.  ]
 [0.   0.   0.84 0.45 0.24 0.   0.   0.  ]
 [0.   0.   0.32 0.22 0.25 0.   0.   0.  ]
 [0.   0.   0.   0.   0.   1.   0.   0.  ]
 [0.   0.   0.   0.   0.   0.   2.   0.  ]
 [0.   0.   0.   0.   0.   0.   0.   3.  ]]
'''

其中

到此这篇关于Python实现生成对角矩阵和对角块矩阵的文章就介绍到这了,更多相关Python矩阵内容请搜索脚本之家以前的文章或继续浏览下面的相关文章希望大家以后多多支持脚本之家！

您可能感兴趣的文章:

python如何实现DES加密
这篇文章主要介绍了python如何实现DES加密，帮助大家更好的理解和学习密码学，感兴趣的朋友可以了解下
2020-09-09
python3 破解 geetest（极验）的滑块验证码功能
这篇文章主要介绍了python3 破解 geetest（极验）的滑块验证码功能,本文通过实例代码给大家介绍的非常详细，对大家的学习或工作具有一定的参考借鉴价值，需要的朋友可以参考下
2018-02-02
Python如何通过内存管理提升程序执行效率
Python提供了自动内存管理的功能，但是如果不小心使用，可能会导致内存泄漏和性能问题，所以巧妙使用内存管理是提高Python执行效率的关键，下面就来和大家仔细讲讲Python的内存管理技巧吧
2023-06-06
用Python实现石头剪刀布游戏
大家好，本篇文章主要讲的是用Python实现石头剪刀布游戏，感兴趣的同学赶快来看一看吧，对你有帮助的话记得收藏一下
2022-01-01
OpenCV角点检测的实现示例
角点通常被定义为两条边的交点，本文主要介绍了OpenCV角点检测的实现示例，文中通过示例代码介绍的非常详细，具有一定的参考价值，感兴趣的小伙伴们可以参考一下
2022-03-03
windows下python使用ffmpeg实现rtsp推流
这篇文章主要为大家详细介绍了在windows环境下python如何使用ffmpeg实现rtsp推流,文中的示例代码讲解详细,有需要的小伙伴可以了解一下
2023-09-09
python常用数据结构字典梳理
这篇文章主要介绍了python常用数据结构字典梳理，文章围绕主题展开详细的内容介绍，具有一定的参考价值，需要的小伙伴可以参考一下
2022-08-08
Python Pygame实现俄罗斯方块
这篇文章主要为大家详细介绍了Python Pygame实现俄罗斯方块，文中示例代码介绍的非常详细，具有一定的参考价值，感兴趣的小伙伴们可以参考一下
2021-02-02
解析numpy中的iscomplex方法及实际应用
NumPy 的 iscomplex 方法为检查数组中的元素是否为复数提供了一种高效且易于使用的接口,本文介绍了 iscomplex 方法的基本概念、使用方法以及它在解决实际问题中的应用,需要的朋友可以参考下
2024-06-06
Python中使用print函数进行不换行打印问题
这篇文章主要介绍了Python中使用print函数进行不换行打印问题,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教
2024-02-02