matlab - Matlab 乘法矩阵 M * P(n-1)

Question

假设我有一个 3x3 矩阵和一个 3x1 向量，我想将我的矩阵A与向量P多次相乘。

A=[0 0.3 0.5; 0.8 0 0.5; 0.2 0.7 0] * P=[1; 1; 1] 

我想将这个答案乘以PN多次。

现在我可以手动完成，但我被迫使用 Matlab，这让我很头疼。

干杯

Answer 1

这只是一个猜测，但我认为你想要的是

PN = mpower(A,N)*P0

这里N是世代数，P0是初始种群向量。

正如@LuisMendo 和@DennisJaheruddin 所建议的，这相当于

PN = A^N*P0

稍微详细说明一下：在每一代中，您使用 计算旧人口的新人口Pnew = A*Pold，即：

 P1 = A*P0                % generation 1
 P2 = A*P1 = A*A*P0       % generation 2
 P3 = A*P2 = A*A*A*P0     % generation 3

等等，这样

 PN = (A*A*A*A...*A)*P0     % generation N
    = A^N*P0