我的编程空间,编程开发者的网络收藏夹
学习永远不晚

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

短信预约 -IT技能 免费直播动态提醒
省份

北京

  • 北京
  • 上海
  • 天津
  • 重庆
  • 河北
  • 山东
  • 辽宁
  • 黑龙江
  • 吉林
  • 甘肃
  • 青海
  • 河南
  • 江苏
  • 湖北
  • 湖南
  • 江西
  • 浙江
  • 广东
  • 云南
  • 福建
  • 海南
  • 山西
  • 四川
  • 陕西
  • 贵州
  • 安徽
  • 广西
  • 内蒙
  • 西藏
  • 新疆
  • 宁夏
  • 兵团
手机号立即预约

请填写图片验证码后获取短信验证码

看不清楚,换张图片

免费获取短信验证码

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

文章目录

对角矩阵

scipy中的函数

scipy.linalg中,通过tri(N, M=None, k=0, dtype=None)可生成 N × M N\times M N×M对角矩阵,若M=None,则 M M M默认为 N N Nk表示矩阵中用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用于生成对角矩阵,diagdiagflat基础上,添加了提取对角元素的功能,例如

>>> 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用于调节01的分界线相对于对角线的位置,例如

>>> 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)而言,会生成如下形式矩阵

A 0 0 0 B 0 0 0 C \begin{matrix} A&0&0\\0&B&0\\0&0&C\\ \end{matrix} A000B000C

from scipy.linalg import *import numpy as npA = 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)

其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. ] \begin{bmatrix} 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. \end{bmatrix} 1.1.0.0.0.0.0.0.1.1.0.0.0.0.0.0.0.0.0.80.840.320.0.0.0.0.0.380.450.220.0.0.0.0.0.410.240.250.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.

其中

A= [ 1 1 1 1 ] B= [ 0.8 0.38 0.41 0.84 0.45 0.24 0.32 0.22 0.25 ] C= [ 1 0 0 0 2 0 0 0 3 ] A=\begin{bmatrix}1&1\\1&1\end{bmatrix}\quad B=\begin{bmatrix}0.8 &0.38&0.41\\0.84&0.45&0.24\\0.32&0.22&0.25\end{bmatrix}\quad C=\begin{bmatrix}1&0&0\\0&2&0\\0&0&3\end{bmatrix} A=[1111]B= 0.80.840.320.380.450.220.410.240.25 C= 100020003

来源地址:https://blog.csdn.net/m0_37816922/article/details/129936879

免责声明:

① 本站未注明“稿件来源”的信息均来自网络整理。其文字、图片和音视频稿件的所属权归原作者所有。本站收集整理出于非商业性的教育和科研之目的,并不意味着本站赞同其观点或证实其内容的真实性。仅作为临时的测试数据,供内部测试之用。本站并未授权任何人以任何方式主动获取本站任何信息。

② 本站未注明“稿件来源”的临时测试数据将在测试完成后最终做删除处理。有问题或投稿请发送至: 邮箱/279061341@qq.com QQ/279061341

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

下载Word文档到电脑,方便收藏和打印~

下载Word文档

猜你喜欢

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

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

怎么使用Python实现生成对角矩阵和对角块矩阵

这篇文章主要介绍了怎么使用Python实现生成对角矩阵和对角块矩阵的相关知识,内容详细易懂,操作简单快捷,具有一定借鉴价值,相信大家阅读完这篇怎么使用Python实现生成对角矩阵和对角块矩阵文章都会有所收获,下面我们一起来看看吧。对角矩阵s
2023-07-06

怎么生成对角矩阵numpy.diag

本文小编为大家详细介绍“怎么生成对角矩阵numpy.diag”,内容详细,步骤清晰,细节处理妥当,希望这篇“怎么生成对角矩阵numpy.diag”文章能帮助大家解决疑惑,下面跟着小编的思路慢慢深入,一起来学习新知识吧。生成对角矩阵 nump
2023-06-30

Python——生成矩阵

def makeMatrix(row_num, column_num, start=0.1, step=0): matrix = [] row = [] for i in range(row_num): fo
2023-01-31

Python 生成矩阵

限定步长,起始数字,然后生成x行,y列的矩阵>>> def range2rect(x,y,start=0,step=1):... N=[]... F=[]... for i in range(x):                  #等价
2023-01-31

Python numpy生成矩阵、串联矩阵代码分享

import numpy生成numpy矩阵的几个相关函数:numpy.array() numpy.zeros() numpy.ones() numpy.eye() 串联生成numpy矩阵的几个相关函数:numpy.array() numpy
2022-06-04

一个3x3矩阵对角线元素相加的题目

今天看完了数组,在做课后题的时候发现有这么一道题,让求出一个3x3矩阵对角线元素相加的问题,在网上大概的看了下感觉不是很全面,只是求了一个对角线元素的和,并没有求另外一个对角线的元素之和,当然我觉的我写的这个也不是很全面要是改成4x4矩阵或
2023-01-31

Python用scipy生成特殊矩阵

本文主要介绍了Pytho用scipy生成特殊矩阵详解,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友们可以参考一下
2023-05-16

python生成单位阵或对角阵的三种方式是什么

本篇内容介绍了“python生成单位阵或对角阵的三种方式是什么”的有关知识,在实际案例的操作过程中,不少人都会遇到这样的困境,接下来就让小编带领大家学习一下如何处理这些情况吧!希望大家仔细阅读,能够学有所成!python生成单位阵或对角阵前
2023-06-30

Python如何生成n行n列的矩阵

要生成n行n列的矩阵,可以使用嵌套的列表推导式。下面是一个示例代码:```pythonn = 3 # 矩阵的大小,这里为3x3# 生成n行n列的零矩阵matrix = [[0 for j in range(n)] for i in ran
2023-09-26

编程热搜

  • Python 学习之路 - Python
    一、安装Python34Windows在Python官网(https://www.python.org/downloads/)下载安装包并安装。Python的默认安装路径是:C:\Python34配置环境变量:【右键计算机】--》【属性】-
    Python 学习之路 - Python
  • chatgpt的中文全称是什么
    chatgpt的中文全称是生成型预训练变换模型。ChatGPT是什么ChatGPT是美国人工智能研究实验室OpenAI开发的一种全新聊天机器人模型,它能够通过学习和理解人类的语言来进行对话,还能根据聊天的上下文进行互动,并协助人类完成一系列
    chatgpt的中文全称是什么
  • C/C++中extern函数使用详解
  • C/C++可变参数的使用
    可变参数的使用方法远远不止以下几种,不过在C,C++中使用可变参数时要小心,在使用printf()等函数时传入的参数个数一定不能比前面的格式化字符串中的’%’符号个数少,否则会产生访问越界,运气不好的话还会导致程序崩溃
    C/C++可变参数的使用
  • css样式文件该放在哪里
  • php中数组下标必须是连续的吗
  • Python 3 教程
    Python 3 教程 Python 的 3.0 版本,常被称为 Python 3000,或简称 Py3k。相对于 Python 的早期版本,这是一个较大的升级。为了不带入过多的累赘,Python 3.0 在设计的时候没有考虑向下兼容。 Python
    Python 3 教程
  • Python pip包管理
    一、前言    在Python中, 安装第三方模块是通过 setuptools 这个工具完成的。 Python有两个封装了 setuptools的包管理工具: easy_install  和  pip , 目前官方推荐使用 pip。    
    Python pip包管理
  • ubuntu如何重新编译内核
  • 改善Java代码之慎用java动态编译

目录