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

相关文章

python的依赖管理的实现

主流开发语言的包管理工具一般都是支持依赖管理的,比如PHP的composer、Java的mvn。 对于python来说又该如何管理依赖呢? pip基本用法 python还不错,它提供了...

Python:Numpy 求平均向量的实例

如下所示: >>> import numpy as np >>> a = np.array([[1, 2, 3], [3, 1, 2]]) >...

浅谈django rest jwt vue 跨域问题

django rest framework 使用 router 注册url时,访问接口 包302错误 可能是因为请求url 写法有问题, 如请求 /api/login/ 报302 ,需要...

一些Python中的二维数组的操作方法

需要在程序中使用二维数组,网上找到一种这样的用法:   #创建一个宽度为3,高度为4的数组 #[[0,0,0], # [0,0,0], # [0,0,0], # [0,...

python操作kafka实践的示例代码

python操作kafka实践的示例代码

1、先看最简单的场景,生产者生产消息,消费者接收消息,下面是生产者的简单代码。 #!/usr/bin/env python # -*- coding: utf-8 -*- impor...