bus_boy 2022-05-20 15:08 采纳率: 100%
浏览 419
已结题

如何用python构造未知数矩阵并计算

img

比如像图片这样的矩阵,构造好之后传输数据进行批量计算。

img


目前只会基础的sympy和numpy,如果有朋友知道学哪一个模块可以做出来也会采纳答案!十分感谢

  • 写回答

2条回答 默认 最新

  • 歇歇 2022-05-26 01:51
    关注

    python未知数的矩阵运算_python如何进行矩阵运算

    import numpy as np
    
    

    加法

    思路

    参考矩阵加法公式

    先判断两矩阵能否相加

    遍历两个矩阵

    各相同位置元素相加

    实现

    def add(a, b):
    
    if(a.shape != b.shape):
    
    print('两矩阵不为同型矩阵!')
    
    return
    
    c = np.zeros(a.shape)
    
    for i in range(a.shape[0]):
    
    for j in range(a.shape[1]):
    
    c[i][j] = a[i][j] + b[i][j]
    
    return c
    

    数乘

    思路

    参考矩阵数乘公式

    遍历矩阵

    相乘

    实现

    def numSub(a, b):
    
    c = np.zeros(a.shape)
    
    for i in range(a.shape[0]):
    
    for j in range(a.shape[1]):
    
    c[i][j] = a[i][j] * b
    
    return c
    

    Hadamard乘积

    思路

    参考Hadamard乘积公式

    与矩阵加法思路一致

    实现

    # 矩阵对应元素相乘 Hadamard乘积
    
    def had(a,b):
    
    if(a.shape != b.shape):
    
    print('两矩阵不为同型矩阵!')
    
    return
    
    c = np.zeros(a.shape)
    
    for i in range(c.shape[0]):
    
    for j in range(c.shape[1]):
    
    c[i][j] = a[i][j] * b[i][j]
    

    点积

    思路

    参考矩阵点积公式

    理解公式:乘积C的第m行第n列的元素等于矩阵A的第m行的元素与矩阵B的第n列对应元素乘积之和。

    判断两矩阵是否为同型矩阵

    构建m行n列的结果矩阵

    向结果矩阵中赋值:

    (1). 遍历结果矩阵

    (2). 求出矩阵A对应行与矩阵B对应列的元素乘积之和:

    根据公式,A矩阵的行与B矩阵的列均确定,因此根据累加条件遍历两个矩阵来确定另一个索引位置,并累加求和

    实现

    # 矩阵相乘
    
    def mul(a,b):
    
    if(a.shape[1] != b.shape[0]):
    
    print('两矩阵不能相乘')
    
    return
    
    m = a.shape[0]
    
    n = b.shape[1]
    
    s = a.shape[1]
    
    c = np.zeros((m,n))
    
    for i in range(m):
    
    for j in range(n):
    
    for k in range(s):
    
    c[i][j] += a[i][k]*b[k][j]
    
    pass
    
    pass
    
    return c
    
    本回答被题主选为最佳回答 , 对您是否有帮助呢?
    评论
查看更多回答(1条)

报告相同问题?

问题事件

  • 系统已结题 6月4日
  • 已采纳回答 5月27日
  • 创建了问题 5月20日

悬赏问题

  • ¥15 如何在scanpy上做差异基因和通路富集?
  • ¥20 关于#硬件工程#的问题,请各位专家解答!
  • ¥15 关于#matlab#的问题:期望的系统闭环传递函数为G(s)=wn^2/s^2+2¢wn+wn^2阻尼系数¢=0.707,使系统具有较小的超调量
  • ¥15 FLUENT如何实现在堆积颗粒的上表面加载高斯热源
  • ¥30 截图中的mathematics程序转换成matlab
  • ¥15 动力学代码报错,维度不匹配
  • ¥15 Power query添加列问题
  • ¥50 Kubernetes&Fission&Eleasticsearch
  • ¥15 報錯:Person is not mapped,如何解決?
  • ¥15 c++头文件不能识别CDialog