擅长:python、mysql、java
<p>好的,我找到解决办法了。你知道吗</p>
<p>正确的方法是完全避免计算逆。你知道吗</p>
<p>如果我们有一个线性问题<code>Ax = b</code>,解决方法是:<code>x = A^(-1)*b</code></p>
<p>例如,Python已经有了计算这个解决方案的函数scipy.sparse.linalg公司斯索尔夫先生。你知道吗</p>
<p>如果我需要计算<code>(AHA.T)^(-1) * b</code>,我只需要调用<code>spsolve(A * H * A.T, b)</code>。你知道吗</p>
<p>我对t的计算如下:<code>t_k = (c - A.T * y_k)/(A.T * spsolve(A * H * A.T, b))</code></p>
<p>这给了我一个稳定的解决方案,我的整个算法在20次迭代中收敛。你知道吗</p>