如何在Python中求矩阵对角线上方和下方元素的总和?

2023-12-09

我需要找到主对角线上方和下方元素的总和。我不知道如何调整算法以仅对这些数字求和。这是我到目前为止的代码,其中 A 是矩阵

A = []
N = int(raw_input("Input matrix size: "))
for i in range(0, N):
    row = []
    for j in range(0, N):
        row.append(int(raw_input("Input elements: ")))
    A.append(row)
sum = 0
for i in range(0, N):
    sum += A[i][i]
print sum       
sum2 = 0
for i in range(0, N):
    for j in range(i+1, N):
        sum2 += A[i][j]
print sum2

我想我应该使用更多的 for 语句。 谢谢


您可以使用np.triu, np.tril and np.trace计算这些总和(您的问题没有具体说明是否允许您利用杠杆numpy):

import numpy as np

np.random.seed(0)
A = np.random.randint(0,10,size=(5,5))

Gives:

[[5 0 3 3 7]
 [9 3 5 2 4]
 [7 6 8 8 1]
 [6 7 7 8 1]
 [5 9 8 9 4]]

Then:

upper_sum = np.triu(A).sum()-np.trace(A)
lower_sum = np.tril(A).sum()-np.trace(A)

Yields:

34
73
本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系:hwhale#tublm.com(使用前将#替换为@)

如何在Python中求矩阵对角线上方和下方元素的总和? 的相关文章

  • 带括号的上下文管理器

    我试图了解新的新内容带括号的上下文管理器Python 3 10 中的功能 新功能中的顶部项目here https docs python org 3 10 whatsnew 3 10 html 我的测试示例是尝试编写 with open f
  • 帮助需要在可选条件下编写正则表达式[关闭]

    我有一个日志文件包含如下内容 log Using data from yyyy mm dd 2011 8 3 0 files queued for scanning Warning E test H ndler pdf File not F
  • Python 小数.InvalidOperation 错误

    当我运行这样的东西时 我总是收到此错误 from decimal import getcontext prec 30 b 2 3 Decimal b Error Traceback most recent call last File Te
  • 引发 RuntimeError(f"目录 '{directory}' 不存在") RuntimeError: 导入 fitz 时目录 'static/' 不存在

    当我运行 extract img py 文件时出现此错误 RuntimeError f 目录 directory 不存在 运行时错误 导入 fitz 时不存在目录 static 我不明白为什么这会给我发回此错误消息 我之前看到过关于这个话题
  • 在 Python 3 中动态导入模块的问题

    我遇到的情况是 在我的 Python 3 项目中 在运行时必须包含某些模块 我在用着importlib import module为了这 第二次更新 我确实找到了一种方法来做一些接近我想要的事情 一些额外的代码可能会使我的一些链接稍微偏离一
  • 使用字母而不是数字进行顺序计数[重复]

    这个问题在这里已经有答案了 我需要一种方法 将字符串 递增 到 z 然后将 aa 递增到 az 然后将 ba 递增到 bz 依此类推 就像 Excel 工作表中的列一样 我将向该方法提供前一个字符串 它应该增加到下一个字母 PSEUDO C
  • 如何从对应的数组值中获取数组键?

    您可以轻松地通过其键获取数组值 如下所示 value array key 但如果我有该值并且想要它的密钥怎么办 获得它的最佳方式是什么 你可以使用array search https www php net manual en functi
  • 矩阵向量变换

    我正在编写一个代码来制作软件蒙皮器 骨骼 皮肤动画 并且我正处于 优化 阶段 蒙皮器工作得很好 并且在 Core 上 1 09 毫秒内对 4900 个三角形网格与 22 个骨骼进行蒙皮Duo 2 Ghz 笔记本 我需要知道的是 1 有人可以
  • 如何使用python读取最后一行的特定位置

    我有一个太大的 txt 文件 并且有几行类似的行 如下所示 字1 字2 字3 字4 553 75 我对位置 4 值 感兴趣 即最后一行 553 75 我的文件文本 word1 word2 word3 word4 553 20 word1 w
  • 如何使用 jira-python 设置 fixVersions 字段

    我正在尝试使用 jira python 模块 http jira python readthedocs org en latest 更新现有的 JIRA 具体来说 我正在尝试设置问题的fixesVersion 列表 我已经尝试了一段时间但没
  • Floyd-Warshall 算法:获取最短路径

    假设一个图由一个表示n x n维数邻接矩阵 我知道如何获得所有对的最短路径矩阵 但我想知道有没有办法追踪所有最短路径 Blow是python代码实现 v len graph for k in range 0 v for i in range
  • 为 Python 2.4 改进“with”语句的直接替换

    您能否建议一种方法来编写可在 Python 2 4 中使用的 with 语句的直接替换代码 这将是一个 hack 但它可以让我更好地将我的项目移植到 Python 2 4 EDIT 删除了不相关的元类草图 只需使用 try finally
  • x11 - 导入错误:没有名为“kivy.core.window.window_x11”的模块

    目前我正在尝试构建一个我通过 buildozer 用 Python 和 Kivy 编写的应用程序 无论我在做什么 我都会遇到 window x11 的问题 即使我在代码中注释掉所有与 Windows 相关的内容或执行本文中描述的所有操作 这
  • 如何在 Python 中执行相当于预处理器指令的操作?

    有没有办法在 Python 中执行以下预处理器指令 if DEBUG lt do some code gt else lt do some other code gt endif There s debug 这是编译器预处理的特殊值 if
  • 请求response.iter_content()获取不完整的文件(1024MB而不是1.5GB)?

    您好 我一直在使用此代码片段从网站下载文件 到目前为止 小于 1GB 的文件都很好 但我注意到 1 5GB 文件不完整 s is requests session object r s get fileUrl headers headers
  • 从 Python 中编译的正则表达式中提取命名组正则表达式模式

    我有一个 Python 正则表达式 其中包含多个命名组 但是 如果先前的组已匹配 则可能会错过与一组匹配的模式 因为似乎不允许重叠 举个例子 import re myText sgasgAAAaoasgosaegnsBBBausgisego
  • 在 scrapy 中将基本 url 与结果 href 结合起来

    下面是我的蜘蛛代码 class Blurb2Spider BaseSpider name blurb2 allowed domains www domain com def start requests self yield self ma
  • Selenium Python 使用代理运行浏览器[重复]

    这个问题在这里已经有答案了 我正在尝试编写一个非常简单的脚本 该脚本从 txt 文件获取代理 不需要身份验证 并用它打开浏览器 然后沿着代理列表循环此操作一定次数 我确实知道如何打开 txt 文件并使用它 我的主要问题是让代理正常工作 我见
  • python中匹配3个或更多相同的字符

    我正在尝试使用正则表达式在字符串中查找三个或更多相同的字符 例如 你好 不匹配 噢 会的 我尝试过做类似的事情 re compile 1 3 a zA Z re compile w 1 5 但似乎都不起作用 w 1 2 是您正在寻找的正则表
  • Shap - 颜色条不显示在摘要图中

    显示summary plot时 不显示颜色条 shap summary plot shap values X train 我尝试过改变plot size 当绘图较高时 会出现颜色条 但它非常小 看起来不应该 shap summary plo

随机推荐

  • OpenCV视频采集和fps问题

    我正在 MacOSX 上使用 OpenCV 从网络摄像头捕获视频 它工作正常 但当我尝试在 QuickTime 上播放我捕获的视频时 它播放得太快了 即我从相机捕获 10 秒 但当我在 QuickTime 上播放时 视频为 2 秒 我尝试将
  • 如何使用Python创建Windows服务

    我编写了一个 python 脚本 它将作为 Windows 服务安装 下面是代码 import datetime import logging from logging handlers import RotatingFileHandler
  • ActionScript 中的普通对象和强类型对象哪个更快?

    So is public var user Object user firstName Bill user lastName Cosby 比我有像 User 这样的值对象慢吗 var user User new User user firs
  • 如何从服务器读取消息? (TCP)

    客户端写入服务器 服务器读取 服务器写入客户端 客户端不读取 Server using System using System Text using System Net using System Net Sockets using Sys
  • 为什么Java中的“this”关键字是final的? [复制]

    这个问题在这里已经有答案了 这似乎是一件几乎没有人意识到的事情 但Java中的 this引用 是最终的 在正常的编程日中 我认为我可以通过在自己的类中重新定义 this 引用来重新定义整个实例 public void method this
  • 1D numpy 数组,在新的 2D 数组中每个连续行向右移动

    我试图通过删除 for 循环并仅在处理大型数据集时使用 numpy 数组来优化一些代码 我想采用一维 numpy 数组 例如 a 1 2 3 4 5 并生成一个 2D numpy 数组 其中每列中的值沿着一个位置移动 例如在上面的情况下 我
  • 找不到 InsertOnSubmit() 方法

    我是实体框架的新手 我认为我在这里误解了一些东西 我试图在表中插入一行 在我找到的代码示例中 他们都调用方法 InsertOnSubmit 但问题是我在任何地方都找不到方法 InsertOnSubmit 或 SubmitChanges 错误
  • Asp.net Mvc Ajax Json(数组后)

    javascript jQuery var items new Array var obj Begin 444 End end items push obj items push obj var request DateStart Date
  • OpenGL 着色器版本编译错误

    我在 OpenGL 中编译简单的顶点着色器时遇到问题 收到以下错误消息 错误 106 GL2 不支持版本号 错误 279 布局限定符 位置 无效 我认为我一定使用了错误版本的 GL2 但我不知道如何找到我的版本号或去哪里升级 是的 我试图寻
  • 无法使用 VS11 创建 VC++ 项目

    我刚刚在 Windows 8 Consumer Preview 操作系统上安装了 VS11 Professional 版本 我尝试创建 win32 控制台应用程序 向导失败并出现以下错误 Microsoft Visual Studio C
  • cvFindContours 如何工作?

    我在 c 和 openCv 2 1 中使用 cvFindContours 函数 因为我必须编写我正在使用的函数的文档 有人可以告诉我这个函数是如何工作的 使用哪种查找轮廓的方法以及如果可能的话一些基础理论的功能 正如您在文档中所读到的 该函
  • 销毁地图实例的正确方法是什么?

    我最近开发了一个html5移动应用程序 该应用程序是一个单页面 其中导航哈希更改事件替换了整个 DOM 该应用程序的一部分是使用 API v3 的 Google 地图 在从 DOM 中删除地图 div 之前 我想删除所有事件处理程序 侦听器
  • java.lang.NoClassDefFoundError:无法初始化类 javax.crypto.JceSecurityManager

    java lang NoClassDefFoundError Could not initialize class javax crypto JceSecurityManager javax crypto Cipher getConfigu
  • 将正态分布拟合到加权列表

    我有一堆数据点 我想对数据进行正态分布 我看到 scipy 有stats norm fit方法 但这需要一个数据点列表 就像是 data 1 1 1 1 2 2 2 3 3 3 3 3 3 4 5 5 5 而我的数据包含在两个列表中 类似
  • 未解析的外部静态变量

    我有一个名为 Pub 的类 它具有以下标头 pragma once class Pub public static double X FACTOR static double Y FACTOR static const int INIT S
  • 可以采取哪些措施来加速同步 WCF 调用?

    我对 Silverlight 应用程序中的同步 WCF 调用的性能测量表明我可以7 次调用 秒在本地主机连接上 即非常慢 可以加快速度吗 或者这是正常现象吗 这是我的测试代码 const UInt32 nrCalls 100 ICalcul
  • Angular2 - 单击按钮时更新模型

    当使用 angular2 ngModel 进行双向数据绑定时
  • 主/从 Mysql 架构与服务器/读取数据库和仅写入的单独数据库[关闭]

    很难说出这里问的是什么 这个问题模棱两可 含糊不清 不完整 过于宽泛或言辞激烈 无法以目前的形式合理回答 如需帮助澄清此问题以便重新打开 访问帮助中心 主 从类型 mysql 架构与负载平衡 Web 服务器 每个服务器都使用单独的服务器保存
  • KEYCLOAK:请求中未提供客户端密钥

    keycloak js 版本 gt 10 0 2 角度 gt 7 2 16 keycloak json realm REALM auth server url
  • 如何在Python中求矩阵对角线上方和下方元素的总和?

    我需要找到主对角线上方和下方元素的总和 我不知道如何调整算法以仅对这些数字求和 这是我到目前为止的代码 其中 A 是矩阵 A N int raw input Input matrix size for i in range 0 N row