回答此问题可获得 20 贡献值,回答如果被采纳可获得 50 分。
<p>问题:</p>
<p>我有一个3*3的转移矩阵和另一个3*1矩阵。把这两个相乘,我得到另一个3*1矩阵。我必须将新的3*1矩阵与原始的3*3矩阵相乘,这将得到另一个3*1矩阵。这必须持续,直到连续步骤中获得的3*1矩阵足够接近,或者如果完成10000次此类乘法(以先到者为准)。要查看它们是否足够接近,我有一个函数:</p>
<p>np.allclose(上一个,新的)#确定此过程是否必须继续或停止</p>
<p>previous_x表示之前的3*1矩阵,new_x表示当前的新3*1矩阵</p>
<p>最大步数=10000</p>
<p>为了进入else循环,下面代码中的“res”肯定会在10k步数结束之前的某个时间点变为“True”。但由于某些原因,它不会发生。你能帮我解决这个问题吗</p>
<p>[守则]</p>
<pre><code>def random_walk(P,x_0,max_steps = 10000):
n_steps = 0
def matrix_mult(P,x_0,n_steps):
x = np.dot(P,x_0)
res = np.allclose(x_0,x)
n_steps = n_steps+1
return x,res,n_steps
x,res,n_steps = matrix_mult(P,x_0,n_steps)
print(res)
for i in range(max_steps-1):
if(res==False):
print('h')
x_0 = x
x=np.dot(P,x_0)
else:
print('g')
x_0 = x
x,res,n_steps = matrix_mult(P,x_0,n_steps)
n_steps =i
return x, n_steps
random_walk(np.array([[0,1,0.5],[1,0,0],[0,0,0.5]]),np.array([1,0,0]),max_steps = 10000)
</code></pre>