您好,登錄后才能下訂單哦!
計算:Ax-b
A: 2*2
x: 2*1
b: 2*1
so, Ax-b: 2*1
if __name__ == "__main__": A = np.array([[4.0, 1.0], [1.0, 3.0]]) b = np.array([[1.0], [2.0]]) x_0 = np.array([[2.0], [1.0]]) r_k = A * x_0 - b print(r_k)
錯誤!?。?/p>
修改:
if __name__ == "__main__": A= mat([[4.0, 1.0], [1.0, 3.0]]) b = mat([[1.0], [2.0]]) x_0 = mat([[2.0], [1.0]]) r_k = A * x_0 - b print(r_k)
if __name__ == "__main__": A= mat([[4.0, 1.0], [1.0, 3.0]]) b = mat([[1.0], [2.0]]) x_k = mat([[2.0], [1.0]]) p_k = -x_k # 2行1列 r_k = A * x_k - b # 2行1列 alpha_k = (np.transpose(r_k) * r_k) / (np.transpose(p_k) * A * p_k) # 1行1列 ----------------------------------------------- x_k = x_k + alpha_k * p_k #2行1列 ?。。?!這里報錯 ----------------------------------------------- print(x_k)
修改:
if __name__ == "__main__": A= mat([[4.0, 1.0], [1.0, 3.0]]) b = mat([[1.0], [2.0]]) x_k = mat([[2.0], [1.0]]) p_k = -x_k # 2*1 r_k = A * x_k - b # 2*1 alpha_k = (np.transpose(r_k) * r_k) / (np.transpose(p_k) * A * p_k) # 1*1 ----------------------------------------------- x_k = x_k + p_k *alpha_k ----------------------------------------------- print(x_k)
以上這篇解決Python計算矩陣乘向量,矩陣乘實數(shù)的一些小錯誤就是小編分享給大家的全部內(nèi)容了,希望能給大家一個參考,也希望大家多多支持億速云。
免責(zé)聲明:本站發(fā)布的內(nèi)容(圖片、視頻和文字)以原創(chuàng)、轉(zhuǎn)載和分享為主,文章觀點不代表本網(wǎng)站立場,如果涉及侵權(quán)請聯(lián)系站長郵箱:is@yisu.com進(jìn)行舉報,并提供相關(guān)證據(jù),一經(jīng)查實,將立刻刪除涉嫌侵權(quán)內(nèi)容。