Python高斯消除矩阵


Posted in Python onJanuary 02, 2019

高斯消除矩阵

#! /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)

总结

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

Python 相关文章推荐
Python基于分水岭算法解决走迷宫游戏示例
Sep 26 Python
Python实现扩展内置类型的方法分析
Oct 16 Python
Python编程pygal绘图实例之XY线
Dec 09 Python
Windows下的Jupyter Notebook 安装与自定义启动(图文详解)
Feb 21 Python
pandas.DataFrame.to_json按行转json的方法
Jun 05 Python
使用Python微信库itchat获得好友和群组已撤回的消息
Jun 24 Python
详解Python的hasattr() getattr() setattr() 函数使用方法
Jul 09 Python
APIStar:一个专为Python3设计的API框架
Sep 26 Python
python并发和异步编程实例
Nov 15 Python
对Python 除法负数取商的取整方式详解
Dec 12 Python
应用OpenCV和Python进行SIFT算法的实现详解
Aug 21 Python
python Canny边缘检测算法的实现
Apr 24 Python
python遍历小写英文字母的方法
Jan 02 #Python
Python最小二乘法矩阵
Jan 02 #Python
Centos部署django服务nginx+uwsgi的方法
Jan 02 #Python
Python 一句话生成字母表的方法
Jan 02 #Python
使用python将请求的requests headers参数格式化方法
Jan 02 #Python
浅谈python requests 的put, post 请求参数的问题
Jan 02 #Python
Python中出现IndentationError:unindent does not match any outer indentation level错误的解决方法
Apr 18 #Python
You might like
php输出xml格式字符串(用的这个)
2012/07/12 PHP
PHP实现对文本数据库的常用操作方法实例演示
2014/07/04 PHP
PHP开发中解决并发问题的几种实现方法分析
2017/11/13 PHP
PHP安装扩展mcrypt以及相关依赖项深入讲解
2021/03/04 PHP
form中限制文本字节数js代码
2007/06/10 Javascript
JavaScript对Json的增删改属性详解
2016/06/02 Javascript
jQuery获取浏览器类型和版本号的方法
2016/07/05 Javascript
JS转换HTML转义符的方法
2016/08/24 Javascript
基于JavaScript实现前端文件的断点续传
2016/10/17 Javascript
JS根据生日月份和日期计算星座的简单实现方法
2016/11/24 Javascript
Jquery EasyUI $.Parser
2017/06/02 jQuery
JS Input里添加小图标的两种方法
2017/11/11 Javascript
VuePress 静态网站生成方法步骤
2019/02/14 Javascript
详解Vuex下Store的模块化拆分实践
2019/07/31 Javascript
关于vue里页面的缓存详解
2019/11/04 Javascript
vue transition 在子组件中失效的解决
2019/11/12 Javascript
400多行Python代码实现了一个FTP服务器
2012/05/10 Python
Python常见加密模块用法分析【MD5,sha,crypt模块】
2017/05/24 Python
基于Python闭包及其作用域详解
2017/08/28 Python
Python+Pandas 获取数据库并加入DataFrame的实例
2018/07/25 Python
python实现微信小程序自动回复
2018/09/10 Python
pyqt5 实现在别的窗口弹出进度条
2019/06/18 Python
Python PyCharm如何进行断点调试
2019/07/05 Python
Python input函数使用实例解析
2019/11/22 Python
python 实现检验33品种数据是否是正态分布
2019/12/09 Python
python3.8.1+selenium实现登录滑块验证功能
2020/05/22 Python
matplotlib.pyplot.plot()参数使用详解
2020/07/28 Python
h5封装下拉刷新
2020/08/25 HTML / CSS
Lands’ End官网:经典的美国生活方式品牌
2016/08/14 全球购物
英国异国风情旅游网站:Travel Talk Tours(团体旅游、探险旅游、帆船假期)
2018/07/26 全球购物
庆元旦文艺演出主持词
2014/03/27 职场文书
群众路线教育实践活动整改落实情况汇报
2014/10/28 职场文书
教师求职信怎么写
2015/03/20 职场文书
小学秋季运动会通讯稿
2015/11/25 职场文书
《颐和园》教学反思
2016/02/19 职场文书
MySQL数据库安装方法与图形化管理工具介绍
2022/05/30 MySQL