快捷導(dǎo)航

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)

轉(zhuǎn)置矩陣

    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]  # 列表推導(dǎo)式
        return C

矩陣行列式

需要用到矩陣余矩陣函數(shù)，代碼如下：

    def detofmatrix(A):
        m = len(A)  # 矩陣的行數(shù)
        n = len(A[0])  # 矩陣的列數(shù)
        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

矩陣的逆

需要用到矩陣行列式函數(shù)和余矩陣函數(shù)，代碼如下：

    def inverseofmatrix(A):
        m = len(A)  # 矩陣的行數(shù)
        n = len(A[0])  # 矩陣的列數(shù)
        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可逆的情況下），調(diào)用求逆函數(shù)和矩陣相乘函數(shù)即可。

求一元n次方程正實數(shù)根

例如： x³+ ax²+bx = c，其中a, b, c, 均為正整數(shù)，求正實數(shù)根。

到此這篇關(guān)于Python不使用庫進行矩陣運算詳解的文章就介紹到這了,更多相關(guān)Python矩陣運算內(nèi)容請搜索腳本之家以前的文章或繼續(xù)瀏覽下面的相關(guān)文章希望大家以后多多支持腳本之家！

您可能感興趣的文章: