Python高斯消除矩陣
更新時(shí)間：2019年01月02日 15:19:12 作者：staHuri
今天小編就為大家分享一篇關(guān)于Python高斯消除矩陣，小編覺得內(nèi)容挺不錯(cuò)的，現(xiàn)在分享給大家，具有很好的參考價(jià)值，需要的朋友一起跟隨小編來看看吧
高斯消除矩陣
#! /usr/bin/env python
# -*- coding: utf-8 -*-
#
def pprint(A):
for i in A:
print(i)
print("")
data = [
[1, 2, 1, 2],
[3, 8, 1, 12],
[0, 4, 1, 2]
]
n = len(data)
print("輸入數(shù)據(jù)")
pprint(data)
for i in range(n):
print("第{}次操作".format(i))
maxE = abs(data[i][i])
# 最大值在第幾行
maxRow = i
for k in range(i + 1, n):
if abs(data[k][i] > maxE):
maxE = abs(data[k][i])
maxRow = k
# 當(dāng)前行交換最大的行
for k in range(i, n + 1):
data[maxRow][k], data[i][k] = data[i][k], data[maxRow][k]
print("當(dāng)前行交換最大的行 交換結(jié)果")
pprint(data)
# 改0
for k in range(i + 1, n):
c = -data[k][i] / data[i][i]
print(-data[k][i] ,"/", data[i][i])
for j in range(i, n + 1):
if i == j:
data[k][j] = 0
else:
data[k][j] += c * data[i][j]
print("改0結(jié)果")
pprint(data)
# 化簡后結(jié)果
print("化簡結(jié)果")
pprint(data)
# 結(jié)果計(jì)算
result = [0 for i in range(n)]
for i in range(n - 1, -1, -1):
result[i] = data[i][n] / data[i][i]
for k in range(i - 1, -1, -1):
data[k][n] -= data[k][i] * result[i]
print(result)
總結(jié)
以上就是這篇文章的全部內(nèi)容了，希望本文的內(nèi)容對(duì)大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價(jià)值，謝謝大家對(duì)腳本之家的支持。如果你想了解更多相關(guān)內(nèi)容請(qǐng)查看下面相關(guān)鏈接
相關(guān)文章
Scrapy框架爬取Boss直聘網(wǎng)Python職位信息的源碼
今天小編就為大家分享一篇關(guān)于Scrapy框架爬取Boss直聘網(wǎng)Python職位信息的源碼，小編覺得內(nèi)容挺不錯(cuò)的，現(xiàn)在分享給大家，具有很好的參考價(jià)值，需要的朋友一起跟隨小編來看看吧2019-02-02
這篇文章主要介紹了怎樣確定K-means算法中的k值問題,具有很好的參考價(jià)值,希望對(duì)大家有所幫助,如有錯(cuò)誤或未考慮完全的地方,望不吝賜教2024-06-06
今天小編就為大家分享一篇用Python一鍵搭建Http服務(wù)器的方法，具有很好的參考價(jià)值，希望對(duì)大家有所幫助。一起跟隨小編過來看看吧2018-06-06
使用django的ORM框架按月統(tǒng)計(jì)近一年內(nèi)的數(shù)據(jù)方法
今天小編就為大家分享一篇使用django的ORM框架按月統(tǒng)計(jì)近一年內(nèi)的數(shù)據(jù)方法，具有很好的參考價(jià)值，希望對(duì)大家有所幫助。一起跟隨小編過來看看吧2019-07-07
python socket模塊創(chuàng)建和使用套接字示例詳解
這篇文章主要為大家介紹了python socket模塊來創(chuàng)建和使用套接字示例詳解，有需要的朋友可以借鑒參考下，希望能夠有所幫助，祝大家多多進(jìn)步，早日升職加薪2023-06-06
Python采集某評(píng)論區(qū)內(nèi)容的實(shí)現(xiàn)示例
本文主要介紹了Python采集某評(píng)論區(qū)內(nèi)容的實(shí)現(xiàn)示例，文中通過示例代碼介紹的非常詳細(xì)，對(duì)大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價(jià)值，需要的朋友們下面隨著小編來一起學(xué)習(xí)學(xué)習(xí)吧2023-04-04
今天小編就為大家分享一篇基于Python中的yield表達(dá)式詳解，具有很好的參考價(jià)值，希望對(duì)大家有所幫助。一起跟隨小編過來看看吧2019-11-11
pyshp創(chuàng)建shp點(diǎn)文件的方法
今天小編就為大家分享一篇pyshp創(chuàng)建shp點(diǎn)文件的方法，具有很好的參考價(jià)值，希望對(duì)大家有所幫助。一起跟隨小編過來看看吧2018-12-12
python調(diào)用API實(shí)現(xiàn)智能回復(fù)機(jī)器人
這篇文章主要為大家詳細(xì)介紹了python調(diào)用API實(shí)現(xiàn)智能回復(fù)機(jī)器人，具有一定的參考價(jià)值，感興趣的小伙伴們可以參考一下2018-04-04
這篇文章主要介紹了python爬蟲之bs4數(shù)據(jù)解析,文中有非常詳細(xì)的代碼示例,對(duì)正在學(xué)習(xí)python爬蟲的小伙伴們有非常好的幫助,需要的朋友可以參考下2021-04-04