您好,登錄后才能下訂單哦!
我就廢話不多說了,直接上代碼大家一起看吧!
#Gauss迭代法 輸入系數(shù)矩陣mx、值矩陣mr、迭代次數(shù)n(以list模擬矩陣 行優(yōu)先) def Gauss(mx,mr,n=100): if len(mx) == len(mr): #若mx和mr長度相等則開始迭代 否則方程無解 x = [] #迭代初值 初始化為單行全0矩陣 for i in range(len(mr)): x.append([0]) count = 0 #迭代次數(shù)計數(shù) while count < n: for i in range(len(x)): nxi = mr[i][0] for j in range(len(mx[i])): if j!=i: nxi = nxi+(-mx[i][j])*x[j][0] nxi = nxi/mx[i][i] x[i][0] = nxi count = count + 1 return x else: return False #調(diào)用 Gauss(mx,mr,n=100) 示例 mx = [[8,-3,2],[4,11,-1],[6,3,12]] mr = [[20],[33],[36]] print(Gauss(mx,mr,20))
以上這篇python實現(xiàn)高斯(Gauss)迭代法的例子就是小編分享給大家的全部內(nèi)容了,希望能給大家一個參考,也希望大家多多支持億速云。
免責(zé)聲明:本站發(fā)布的內(nèi)容(圖片、視頻和文字)以原創(chuàng)、轉(zhuǎn)載和分享為主,文章觀點不代表本網(wǎng)站立場,如果涉及侵權(quán)請聯(lián)系站長郵箱:is@yisu.com進行舉報,并提供相關(guān)證據(jù),一經(jīng)查實,將立刻刪除涉嫌侵權(quán)內(nèi)容。