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

yipeiwu_com6年前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设计】。

相关文章

Centos5.x下升级python到python2.7版本教程

首先到官网下载python2.7.3版本,编译安装 复制代码 代码如下: $wget http://www.python.org/ftp/python/2.7.3/Python-2.7....

Python第三方库h5py_读取mat文件并显示值的方法

mat数据格式是Matlab默认保存的数据格式。在Python中,我们可以使用h5py库来读取mat文件。 >>> import h5py >>>...

Python学习pygal绘制线图代码分享

Python学习pygal绘制线图代码分享

pygal的安装大家可以参阅:pip和pygal的安装实例教程 线图: import pygal line_chart = pygal.Line() line_chart.titl...

python实现倒计时的示例

复制代码 代码如下:import timecount = 0 a = input('time:') b = a * 60 while (count <...

Django 路由层URLconf的实现

分组 分组的目的:让服务端获得url中的具体数据,通过分组,把需要的数据按函数传参的方式传递给服务器后台 1-无名分组 若要从URL 中捕获一个值,只需要在它周围放置一对圆括号 #...