2021
04-07
04-07
python共轭梯度法特征值迭代次数讨论
共轭梯度法,特征值聚堆情况下迭代次数讨论输入各种特征值聚堆与分散时的矩阵,并应用共轭梯度法,观察迭代次数与聚堆情况的关系。因为对角矩阵的对角线元素为其特征值,则用对角矩阵讨论较为方便代码importnumpyasnpdefcg(x0,A,b):r0=np.dot(A,x0)-bp0=-r0rk=r0pk=p0xk=x0t=0#记录迭代次数whilenp.linalg.norm(rk)>=1e-6:rr=np.dot(rk.T,rk)ak=rr/np.dot(np.dot(pk.T,A),pk)xk=...
继续阅读 >