python实现高斯(Gauss)迭代法的例子

yipeiwu_com5年前Python基础

我就废话不多说了,直接上代码大家一起看吧!

#Gauss迭代法 输入系数矩阵mx、值矩阵mr、迭代次数n(以list模拟矩阵 行优先)
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 #迭代次数计数
  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
 
#调用 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实现高斯(Gauss)迭代法的例子就是小编分享给大家的全部内容了,希望能给大家一个参考,也希望大家多多支持【听图阁-专注于Python设计】。

相关文章

python3连接kafka模块pykafka生产者简单封装代码

1.1安装模块 pip install pykafka 1.2基本使用 # -* coding:utf8 *- from pykafka import KafkaClient...

python实现windows倒计时锁屏功能

python实现windows倒计时锁屏功能 # 倒计时锁屏 import time from ctypes import * def closewindows(closetime)...

python 自动去除空行的实例

code 原文档 1.txt : Hello Nanjing 100 实现代码: file_ = "1.txt" r_file = open(file_, "r"...

分析Python中解析构建数据知识

Python 可以通过各种库去解析我们常见的数据。其中 csv 文件以纯文本形式存储表格数据,以某字符作为分隔值,通常为逗号;xml 可拓展标记语言,很像超文本标记语言 Html ,但主...

python多任务之协程的使用详解

1|0使用yield完成多任务 import time def test1(): while True: print("--1--") time.sleep(0.5)...