今天小编就为大家分享一篇关于Python高斯消除矩阵,小编觉得内容挺不错的,现在分享给大家,具有很好的参考价值,需要的朋友一起跟随小编来看看吧

高斯消除矩阵

#! /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("输入数据")
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
  # 当前行交换最大的行
  for k in range(i, n + 1):
    data[maxRow][k], data[i][k] = data[i][k], data[maxRow][k]
  print("当前行交换最大的行 交换结果")
  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结果")
  pprint(data)
# 化简后结果
print("化简结果")
pprint(data)
# 结果计算
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)

总结

以上就是这篇文章的全部内容了,希望本文的内容对大家的学习或者工作具有一定的参考学习价值,谢谢大家对爱安网的支持。如果你想了解更多相关内容请查看下面相关链接

    无相关信息
最新资讯
Pinterest任命沃尔玛前CTO为工程主管 迎接IPO

Pinterest任命沃尔玛

图片社交网络Pinterest今日宣布,已任命沃尔玛前首席技
Facebook的电商雄心及商业机会

Facebook的电商雄心及

Facebook在电商领域已经进行了多次尝试,未来的商业机会
华为阎力大:去年企业BG收入近110亿美元 年增长40%

华为阎力大:去年企业BG

阎力大在演讲中透露,2018年华为企业BG销售收入接近110
美媒:Facebook明文存储数亿用户密码 被查看900万次

美媒:Facebook明文存储

在这些故障中,员工构建的应用程序记录了Facebook用户的
天文学家发现高速"飞奔"的脉冲星:速度近400万km/h

天文学家发现高速"飞

美国航天局20日发布公报说,天文学家发现一颗在银河系高
因不支持现金结账:Amazon Go在旧金山面临被禁风险

因不支持现金结账:Amaz

基于历史记录,对于这项新提案,亚马逊很可能会积极游说官
最新文章
漂亮的Django Markdown富文本app插件的实现

漂亮的Django Markdow

这篇文章主要介绍了漂亮的Django Markdown富文本app插
Python WEB应用部署的实现方法

Python WEB应用部署的

这篇文章主要介绍了Python WEB应用部署的实现方法,小编
Python中出现IndentationError:unindent does not match any outer indentation level错误的解决方法

Python中出现Indentat

今天小编就为大家分享一篇关于Python中出现Indentatio
浅谈python requests 的put, post 请求参数的问题

浅谈python requests

今天小编就为大家分享一篇浅谈python requests 的put,
Python多图片合并PDF的方法

Python多图片合并PDF

今天小编就为大家分享一篇关于Python多图片合并PDF的
Python关于excel和shp的使用在matplotlib

Python关于excel和shp

今天小编就为大家分享一篇关于Python关于excel和shp的