我正在实现隐马尔可夫模型,因此处理的概率非常小。我通过在对数空间中表示变量(所以x→log(x))来处理下溢,它的副作用是乘法现在被加法取代,加法是通过numpy.logaddexp或类似方法处理的。有没有一种简单的方法来处理对数空间中的矩阵乘法? 最佳答案 这是我能想到的最好的方法。fromscipy.specialimportlogsumexpdeflog_space_product(A,B):Astack=np.stack([A]*A.shape[0]).transpose(2,1,0)Bstack=np.stack([B]*