[b1, b2]            
[a1, a2, a3]  *  [b3, b4]   =  [a1*b1+a2*b2+a3*a3]
[a4, a5, a6]     [b5, b6]       [a4*b2+a5*b4+a6*b6]

[a11, a12, a13]   [b11, b12, b13]     [a11*b11, a12*b12, a13*b13]
[a21, a22, a23]  *  [b21, b22, b23] =  [a21*b21, a22*b22, a23*b23]
[a31, a32, a33]   [b31, b32, b33]     [a31*b31, a32*b32, a33*b33]

# 使用多维数组定义矩阵
matrix_a = [[1,2,3],
            [4,1,2],
            [2,3,1]]

matrix_b = [[1,2,3,2],
            [2,3,6,3],
            [3,1,4,2]]

# 定义显示矩阵的函数
def display(matrix):
   for row in matrix:
      print(row)
   print()

# 显示两个输入矩阵
print('第一个矩阵定义如下：') 
display(matrix_a)
print('第二个矩阵定义如下：')
display(matrix_b)

# 初始化所有为0的矩阵
result = [[0, 0, 0, 0],[0, 0, 0, 0],[0, 0, 0, 0]]

# 两个矩阵相乘
for i in range(len(matrix_a)):

   # 遍历矩阵A的行
   for j in range(len(matrix_b[0])):

      # 遍历矩阵B的列
      for k in range(len(matrix_b)):        
         result[i][j] = matrix_a[i][k] * matrix_b[k][j]

print('两个矩阵的乘法结果为：')
display(result)

第一个矩阵定义如下：
[1, 2, 3]
[4, 1, 2]
[2, 3, 1]

第二个矩阵定义如下：
[1, 2, 3, 2]
[2, 3, 6, 3]
[3, 1, 4, 2]

两个矩阵的乘法结果为：
[9, 3, 12, 6]
[6, 2, 8, 4]
[3, 1, 4, 2]

import numpy as np

# 使用numpy array定义矩阵
matrix_a = np.array([[1,2,5], [1,0,6], [9,8,0]])
matrix_b = np.array([[0,3,5], [4,6,9], [1,8,0]])

# 显示两个输入矩阵
print('第一个矩阵定义如下：') 
print(matrix_a)

print('第二个矩阵定义如下：')
print(matrix_b)

# 两个矩阵相乘
result = matrix_a @ matrix_b

print('两个矩阵的乘法结果为：')
print(result) 

第一个矩阵定义如下：
[[1 2 5]
[1 0 6]
[9 8 0]]
第二个矩阵定义如下：
[[0 3 5]
[4 6 9]
[1 8 0]]
两个矩阵的乘积为：
[[13 55 23]
[6 51 5]
[32 75 117]]

import numpy as np

# 使用numpy array定义矩阵
matrix_a = np.array([[1,2,5], [1,0,6], [9,8,0]])
matrix_b = np.array([[0,3,5], [4,6,9], [1,8,0]])

# 显示两个输入矩阵
print('第一个矩阵定义如下：') 
print(matrix_a)

print('第二个矩阵定义如下：')
print(matrix_b)

# 两个矩阵的乘积
result = matrix_a * matrix_b

print('两个矩阵逐元素相乘的结果为：')
print(result)

第一个矩阵定义如下：
[[1 2 5]
 [1 0 6]
 [9 8 0]]
第二个矩阵定义如下：
[[0 3 5]
 [4 6 9]
 [1 8 0]]
两个矩阵逐元素相乘的结果为：
[[ 0  6 25]
 [ 4  0 54]
 [ 9 64  0]]