矩阵位乘是Python中非常常见的运算，它的实现也比较简单。在Python中，我们可以使用numpy库中的dot函数来进行矩阵位乘运算。

import numpy as np
# 定义两个矩阵
a = np.array([[1, 2], [3, 4]])
b = np.array([[5, 6], [7, 8]])
# 矩阵位乘
result = np.dot(a, b)
# 输出结果
print(result)

以上代码中，我们先导入了numpy库，然后定义了两个矩阵a和b。接着使用numpy中的dot函数进行矩阵位乘运算，并将结果存储在变量result中。

下面我们来解释一下代码中的numpy.dot函数，它的声明如下：

numpy.dot(a, b, out=None)

它接受两个矩阵作为参数a和b，然后将它们进行位乘运算。如果我们还想将结果存储在一个新的矩阵中，可以在调用dot函数时传入一个out参数，如下所示：

import numpy as np
# 定义两个矩阵
a = np.array([[1, 2], [3, 4]])
b = np.array([[5, 6], [7, 8]])
# 初始化一个新的矩阵
c = np.zeros((2, 2))
# 矩阵位乘并存储在新的矩阵中
np.dot(a, b, out=c)
# 输出结果
print(c)

以上代码中，我们在调用dot函数时传入了一个新的矩阵c作为out参数，dot函数将矩阵位乘的结果存储在了c中。最后我们输出了新矩阵c的值。