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正则表达式抓取成语网站
Nov 20 Python
Python实现的数据结构与算法之链表详解
Apr 22 Python
Scrapy的简单使用教程
Oct 24 Python
K-近邻算法的python实现代码分享
Dec 09 Python
Python验证文件是否可读写代码分享
Dec 11 Python
python实现图书管理系统
Mar 12 Python
Python进度条的制作代码实例
Aug 31 Python
python 爬虫百度地图的信息界面的实现方法
Oct 27 Python
python用requests实现http请求代码实例
Oct 31 Python
Python爬取豆瓣视频信息代码实例
Nov 16 Python
使用Python求解带约束的最优化问题详解
Feb 11 Python
配置python的编程环境之Anaconda + VSCode的教程
Mar 29 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
用ADODB来让PHP操作ACCESS数据库的方法
2006/12/31 PHP
PHP Pear 安装及使用
2009/03/19 PHP
php数组中删除元素的实现代码
2012/06/22 PHP
记录mysql性能查询过程的使用方法
2013/05/02 PHP
PHP内核探索:变量存储与类型使用说明
2014/01/30 PHP
php二分查找二种实现示例
2014/03/12 PHP
PHP实现简单ajax Loading加载功能示例
2016/12/28 PHP
PHP文件管理之实现网盘及压缩包的功能操作
2017/09/20 PHP
Yii框架连表查询操作示例
2019/09/06 PHP
javascript中不提供sleep功能如何实现这个功能
2014/05/27 Javascript
直接拿来用的15个jQuery代码片段
2015/09/23 Javascript
基于jQuery实现Ajax验证用户名是否存在实例
2016/03/30 Javascript
深入理解JavaScript中的浮点数
2016/05/18 Javascript
JS 数字转换为大写金额的简单实例
2016/08/04 Javascript
jQuery如何解决IE输入框不能输入的问题
2016/10/08 Javascript
jquery.Jcrop结合JAVA后台实现图片裁剪上传实例
2016/11/05 Javascript
Angular2 自定义validators的实现方法
2017/07/05 Javascript
ES6学习笔记之map、set与数组、对象的对比
2018/03/01 Javascript
axios拦截设置和错误处理方法
2018/03/05 Javascript
JavaScript+Canvas实现彩色图片转换成黑白图片的方法分析
2018/07/31 Javascript
webpack4手动搭建Vue开发环境实现todoList项目的方法
2019/05/16 Javascript
在vue中使用image-webpack-loader实例
2020/11/12 Javascript
[03:01]DOTA2英雄基础教程 露娜
2014/01/07 DOTA
[05:08]2014DOTA2国际邀请赛 Hao专访复仇的胜利很爽
2014/07/15 DOTA
Java中重定向输出流实现用文件记录程序日志
2015/06/12 Python
python实现批量修改文件名代码
2017/09/10 Python
用python的requests第三方模块抓取王者荣耀所有英雄的皮肤实例
2017/12/14 Python
python实现图书管理系统
2018/03/12 Python
Selenium鼠标与键盘事件常用操作方法示例
2018/08/13 Python
Python3.6中Twisted模块安装的问题与解决
2019/04/15 Python
纯CSS3实现圆角效果(含IE兼容解决方法)
2014/05/07 HTML / CSS
L’AGENCE官网:加州女装品牌
2018/06/03 全球购物
经济学人订阅:The Economist
2018/07/19 全球购物
中专毕业自我鉴定
2013/10/16 职场文书
position:sticky 粘性定位的几种巧妙应用详解
2021/04/24 HTML / CSS
15个值得收藏的JavaScript函数
2021/09/15 Javascript