Python不使用库进行矩阵运算详解

 更新时间:2023年08月23日 09:37:58   作者:小小草帽  
这篇文章主要介绍了Python不使用库进行矩阵运算详解,矩阵乘法中,需要判断两个矩阵是否可以进行相乘,即前一个矩阵的列是否等于后一个矩阵的行,,需要的朋友可以参考下

Python矩阵运算 

矩阵乘法

矩阵乘法中,需要判断两个矩阵是否可以进行相乘,即前一个矩阵的列是否等于后一个矩阵的行。具体代码如下:

class Solution:
    def multiplyofmatrix(self, A, f):
        """
        Compute matrix A * matrix f.
        For example:
        Input:
        A: [[1, 2], [3, 4]]    Row of A: 2 Col of A: 2
        f: [[1, 2], [3, 4]]    Row of f: 2 Col of f: 2
        Output:
        ans: [[7, 10], [15, 22]]
        s.t. A.col == f.col
        """
        Arow, Acol = len(A), len(A[0])
        frow, fcol = len(f), len(f[0])
        if Acol != frow:
            return 'False!'
        ans = [[0] * fcol for _ in range(Arow)]
        for i in range(fcol):
            for j in range(Arow):
                for k in range(frow):
                    ans[j][i] += A[j][k]*f[k][i]
        return ans
if __name__=='__main__':
    a = Solution().multiplyofmatrix([[1, 2], [3, 4]], [[1, 2], [3, 4]])
    print(a)

转置矩阵

    def transposeofmatrix(A):
        """
        Matrix transpose.
        For example:
        Input:
        A = [[1, 2],
             [3, 4]]
        Output:
        A_transpose = [[1, 3],
                       [2, 4]]
        """
        return [list(raw) for raw in zip(*A)]

矩阵余矩阵

    def submatrix(A, i, j):
        # 矩阵A第i行第j列元素的余矩阵
        m, n = len(A), len(A[0])
        C = [[A[x][y] for y in range(n) if y != j] for x in range(m) if x != i]  # 列表推导式
        return C

矩阵行列式

需要用到矩阵余矩阵函数,代码如下:

    def detofmatrix(A):
        m = len(A)  # 矩阵的行数
        n = len(A[0])  # 矩阵的列数
        if (m == 1 and n == 1):
            return A[0][0]
        else:
            value = 0
            for j in range(n):
                value += ((-1) ** (j + 2)) * A[0][j] * detofmatrix(submatrix(A, 0, j))
            return value

矩阵的逆

需要用到矩阵行列式函数和余矩阵函数,代码如下:

    def inverseofmatrix(A):
        m = len(A)  # 矩阵的行数
        n = len(A[0])  # 矩阵的列数
        C = [[0] * n for _ in range(m)]
        d = self.detofmatrix(A)
        for i in range(m):
            for j in range(n):
                C[i][j] = ((-1) ** (i + j + 2)) * detofmatrix(submatrix(A, j, i))
                C[i][j] = C[i][j] / d
        return C

求方程组根

求∣ A x − b ∣ = 0即为求 x = A-1 ∗b(A可逆的情况下),调用求逆函数和矩阵相乘函数即可。

求一元n次方程正实数根

例如:  x+ ax+bx = c,其中a, b, c, 均为正整数,求正实数根。

到此这篇关于Python不使用库进行矩阵运算详解的文章就介绍到这了,更多相关Python矩阵运算内容请搜索脚本之家以前的文章或继续浏览下面的相关文章希望大家以后多多支持脚本之家!

相关文章

  • NetWorkX使用方法及nx.draw()相关参数解读

    NetWorkX使用方法及nx.draw()相关参数解读

    这篇文章主要介绍了NetWorkX使用方法及nx.draw()相关参数解读,具有很好的参考价值,希望对大家有所帮助。如有错误或未考虑完全的地方,望不吝赐教
    2022-12-12
  • python中快速进行多个字符替换的方法小结

    python中快速进行多个字符替换的方法小结

    最近在用python给自己的seo工作提高效率和节省时间,发现python真的很不错,可以完成很多事情。多个字符替换是大家可能都会遇到的一个问题,昨天在工作中就碰到了这么一个问题,所以想着记录一下解决方案及其过程,方便以后参考。下面来一起看看吧。
    2016-12-12
  • pytorch--之halfTensor的使用详解

    pytorch--之halfTensor的使用详解

    这篇文章主要介绍了pytorch--之halfTensor的使用详解,具有很好的参考价值,希望对大家有所帮助。如有错误或未考虑完全的地方,望不吝赐教
    2021-05-05
  • Django中URL视图函数的一些高级概念介绍

    Django中URL视图函数的一些高级概念介绍

    这篇文章主要介绍了Django中URL视图函数的一些高级概念,Django是Python重多人气框架中最为著名的一个,需要的朋友可以参考下
    2015-07-07
  • 如何用python给数据加上高斯噪声

    如何用python给数据加上高斯噪声

    这篇文章主要介绍了如何用python给数据加上高斯噪声问题,具有很好的参考价值,希望对大家有所帮助。如有错误或未考虑完全的地方,望不吝赐教
    2023-05-05
  • python从入门到精通(DAY 2)

    python从入门到精通(DAY 2)

    本文是python从入门到精通系列文章的第二篇,主要给大家简单讲解下Python的字典复制和内置函数,非常简单实用,需要的小伙伴可以参考下
    2015-12-12
  • python隐藏终端执行cmd命令的方法

    python隐藏终端执行cmd命令的方法

    今天小编就为大家分享一篇python隐藏终端执行cmd命令的方法,具有很好的参考价值,希望对大家有所帮助。一起跟随小编过来看看吧
    2019-06-06
  • python自动化测试用例全对偶组合与全覆盖组合比较

    python自动化测试用例全对偶组合与全覆盖组合比较

    这篇文章主要为大家介绍了python自动化测试用例全对偶组合与全覆盖组合比较,有需要的朋友可以借鉴参考下,希望能够有所帮助,祝大家多多进步,早日升职加薪
    2022-06-06
  • Pandas 如何筛选包含特定字符的列

    Pandas 如何筛选包含特定字符的列

    这篇文章主要介绍了Pandas 如何筛选包含特定字符的列,具有很好的参考价值,希望对大家有所帮助。一起跟随小编过来看看吧
    2021-03-03
  • python实现学生信息管理系统(精简版)

    python实现学生信息管理系统(精简版)

    这篇文章主要为大家详细介绍了python实现学生信息管理系统的精简版,文中示例代码介绍的非常详细,具有一定的参考价值,感兴趣的小伙伴们可以参考一下
    2020-11-11

最新评论