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的爬虫框架scrapy用21行代码写一个爬虫
Apr 24 Python
Python实现全排列的打印
Aug 18 Python
Python关于excel和shp的使用在matplotlib
Jan 03 Python
Django 实现图片上传和显示过程详解
Jul 18 Python
Python使用Slider组件实现调整曲线参数功能示例
Sep 06 Python
sklearn-SVC实现与类参数详解
Dec 10 Python
python计算二维矩形IOU实例
Jan 18 Python
python将dict中的unicode打印成中文实例
May 11 Python
keras 读取多标签图像数据方式
Jun 12 Python
python利用os模块编写文件复制功能——copy()函数用法
Jul 13 Python
5道关于python基础 while循环练习题
Nov 27 Python
基于Python实现将列表数据生成折线图
Mar 23 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
jQuery EasyUI API 中文文档 - ComboGrid 组合表格
2011/10/13 Javascript
用客户端js实现带省略号的分页
2013/04/27 Javascript
jQuery不使用插件及swf实现无刷新文件上传
2014/12/08 Javascript
JavaScript控制两个列表框listbox左右交换数据的方法
2015/03/18 Javascript
分享我对JS插件开发的一些感想和心得
2016/02/04 Javascript
JS组件系列之Bootstrap table表格组件神器【二、父子表和行列调序】
2016/05/10 Javascript
Angular.js之作用域scope'@','=','&'实例详解
2017/02/28 Javascript
vue-resourse将json数据输出实例
2017/03/08 Javascript
JS栈stack类的实现与使用方法示例
2019/01/31 Javascript
js逆向解密之网络爬虫
2019/05/30 Javascript
了解JavaScript函数中的默认参数
2019/05/30 Javascript
简单了解JavaScript中常见的反模式
2019/06/21 Javascript
Quasar Input:type="number" 去掉上下小箭头 实现加减按钮样式功能
2020/04/09 Javascript
记一次用ts+vuecli4重构项目的实现
2020/05/21 Javascript
vue+elementui实现点击table中的单元格触发事件--弹框
2020/07/18 Javascript
[09:40]DAC2018 4.5 SOLO赛 MidOne vs Miracle
2018/04/06 DOTA
python模块简介之有序字典(OrderedDict)
2016/12/01 Python
Python遍历numpy数组的实例
2018/04/04 Python
Tensorflow实现卷积神经网络的详细代码
2018/05/24 Python
Python Pywavelet 小波阈值实例
2019/01/09 Python
Python 转换RGB颜色值的示例代码
2019/10/13 Python
Python range与enumerate函数区别解析
2020/02/28 Python
关于Python字符串显示u...的解决方式
2020/03/06 Python
html5 学习简单的拾色器
2010/09/03 HTML / CSS
HTML5 canvas基本绘图之图形组合
2016/06/27 HTML / CSS
HTML5中使用json对象的实例代码
2018/09/10 HTML / CSS
英国领先的瓷砖专家:Walls and Floors
2018/04/27 全球购物
捷克鲜花配送:Florea.cz
2018/10/29 全球购物
描述RIP和OSPF区别以及特点
2015/01/17 面试题
环保项目建议书
2014/08/26 职场文书
2014年检验科工作总结
2014/11/22 职场文书
工作失误检讨书范文
2015/01/26 职场文书
社会主义核心价值观主题教育活动总结
2015/05/07 职场文书
车间安全生产管理制度
2015/08/06 职场文书
2016年社区中秋节活动总结
2016/04/05 职场文书
电脑开机弹出documents文件夹怎么回事?弹出documents文件夹解决方法
2022/04/08 数码科技