Python 實(shí)現(xiàn)順序高斯消元法示例
更新時(shí)間：2019年12月09日 09:28:19 作者：倚樓灬風(fēng)細(xì)
今天小編就為大家分享一篇Python 實(shí)現(xiàn)順序高斯消元法示例，具有很好的參考價(jià)值，希望對(duì)大家有所幫助。一起跟隨小編過(guò)來(lái)看看吧
我就廢話不多說(shuō)，直接上代碼吧！
# coding: utf8
import numpy as np
# 設(shè)置矩陣
def getInput():
matrix_a = np.mat([[2, 3, 11, 5],
[1, 1, 5, 2],
[2, 1, 3, 2],
[1, 1, 3, 4]],dtype=float)
matrix_b = np.mat([2,1,-3,-3])
#答案：-2 0 1 1
return matrix_a, matrix_b
def SequentialGauss(mat_a):
for i in range(0, (mat_a.shape[0])-1):
if mat_a[i, i] == 0:
print("終斷運(yùn)算：")
print(mat_a)
break
else:
for j in range(i+1, mat_a.shape[0]):
mat_a[j:j+1 , :] = mat_a[j:j+1,:] - \
(mat_a[j,i]/mat_a[i,i])*mat_a[i, :]
return mat_a
def revert(new_mat):
#創(chuàng)建矩陣存放答案 初始化為0
x = np.mat(np.zeros(new_mat.shape[0], dtype=float))
number = x.shape[1]-1
# print(number)
b = number+1
x[0,number] = new_mat[number,b]/new_mat[number, number]
for i in range(number-1,-1,-1):
try:
x[0,i] = (new_mat[i,b]-np.sum(np.multiply(new_mat[i,i+1:b],x[0,i+1:b])))/(new_mat[i,i])
except:print("錯(cuò)誤")
print(x)
if __name__ == "__main__":
mat_a, mat_b = getInput()
# 合并兩個(gè)矩陣
print("原矩陣")
print(np.hstack((mat_a, mat_b.T)))
new_mat = SequentialGauss(np.hstack((mat_a, mat_b.T)))
print("三角矩陣")
print(new_mat)
print("方程的解")
revert(new_mat)
運(yùn)行結(jié)果如下
以上這篇Python 實(shí)現(xiàn)順序高斯消元法示例就是小編分享給大家的全部?jī)?nèi)容了，希望能給大家一個(gè)參考，也希望大家多多支持腳本之家。
相關(guān)文章
python OpenCV的imread不能讀取中文路徑問(wèn)題及解決
這篇文章主要介紹了python OpenCV的imread不能讀取中文路徑問(wèn)題及解決方案，具有很好的參考價(jià)值，希望對(duì)大家有所幫助。如有錯(cuò)誤或未考慮完全的地方，望不吝賜教2022-07-07
Python實(shí)現(xiàn)學(xué)校管理系統(tǒng)
這篇文章主要為大家詳細(xì)介紹了Python實(shí)現(xiàn)學(xué)校管理系統(tǒng)，具有一定的參考價(jià)值，感興趣的小伙伴們可以參考一下2018-01-01
python對(duì)文件目錄的操作方法實(shí)例總結(jié)
這篇文章主要介紹了python對(duì)文件目錄的操作方法,結(jié)合實(shí)例形式總結(jié)分析了Python針對(duì)文件目錄相關(guān)的遍歷、刪除、移動(dòng)、查找等操作技巧,需要的朋友可以參考下2019-06-06
今天小編就為大家分享一篇PyQt5圖形界面播放音樂(lè)的實(shí)例，具有很好的參考價(jià)值，希望對(duì)大家有所幫助。一起跟隨小編過(guò)來(lái)看看吧2019-06-06
這篇文章主要介紹了Python運(yùn)行DLL文件的方法，本文給大家介紹的非常詳細(xì)，具有一定的參考借鑒價(jià)值,需要的朋友可以參考下2020-01-01
今天小編就為大家分享一篇python 視頻逐幀保存為圖片的完整實(shí)例，具有很好的參考價(jià)值，希望對(duì)大家有所幫助。一起跟隨小編過(guò)來(lái)看看吧2019-12-12
這篇文章主要介紹了python中*與**用法的相關(guān)資料,Python中的*和**用于數(shù)值計(jì)算、序列解包、函數(shù)定義和調(diào)用,*用于乘法和解包序列,**用于乘方、解包字典和擴(kuò)展字典,需要的朋友可以參考下2025-01-01
Python seek()和tell()函數(shù)的具體使用
本文主要介紹了Python seek()和tell()函數(shù)的具體使用，文中通過(guò)示例代碼介紹的非常詳細(xì)，對(duì)大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價(jià)值，需要的朋友們下面隨著小編來(lái)一起學(xué)習(xí)學(xué)習(xí)吧2023-02-02
這篇文章主要介紹了Python反轉(zhuǎn)序列的方法,結(jié)合實(shí)例形式分析了列表、元組及字符串反轉(zhuǎn)的具體實(shí)現(xiàn)技巧,需要的朋友可以參考下2018-03-03
Python語(yǔ)言中的重要函數(shù)對(duì)象用法小結(jié)
Python作為一種強(qiáng)大的編程語(yǔ)言,提供了多種高級(jí)函數(shù)對(duì)象,如lambda匿名函數(shù)、map()、reduce()函數(shù),以及迭代器和生成器的使用,本文給大家介紹Python語(yǔ)言中的重要函數(shù)對(duì)象用法,感興趣的朋友跟隨小編一起看看吧2024-09-09