基于python解线性矩阵方程(numpy中的matrix类)
作者：yaominghui 发布时间：2023-11-03 06:54:12
标签：python,解,线性,矩阵,方程,numpy,matrix
这学期有一门运筹学，讲的两大块儿：线性优化和非线性优化问题。在非线性优化问题这里涉及到拉格朗日乘子法，经常要算一些非常变态的线性方程，于是我就想用python求解线性方程。查阅资料的过程中找到了一个极其简单的解决方式，也学到了不少东西。先把代码给出。
import numpy as np
# A = np.mat('1 2 3;2 -1 1;3 0 -1')
A = np.array([[1, 2, 3], [2, -1, 1], [3, 0, -1]])
b = np.array([9, 8, 3])
x = np.linalg.solve(A, b)
print(x)
是不是很简洁？因为调用了强大的包numpy~ 我们想解决的问题是求解矩阵方程Ax=bAx=b。在这里调用numpy中的线性代数包np.linalg，使用其中的function->solve(A, b)。几行代码就解决了问题。在这里solve函数有两个输入，第一个输入是矩阵，可以采用numpy里的矩阵数据类型或者最常用的数组数据类型。第二个输入是右端项b，一个一维numpy数组即可。函数返回方程的解，shape和b是相同的。如果矩阵A是奇异的或者不是方阵，函数就会报错。
0
投稿猜你喜欢
- 一、基本形式sorted(iterable[, cmp[, key[, reverse]]])iterable.sort(cmp[, key
- Python实战系列用于记录实战项目中的思路，代码实现，出现的问题与解决方案以及可行的改进方向本文为第2篇–200行Python代码实现20
Python数据库编程之pymysql学习之前务必安装MySQL并已启动相关服务。一、pymsql的安装在python3的环境中直接使用以下- 1、仅 Safari 和 Opera 识别的 Hack@media all and (min-width:&n
- 看代码：Vue提供了强大的前端开发架构，很多时候我们需要判断数据对象是否为空，使用typeof判断是个不错选择，具体代码见图。补充知识：vu
- Python HTTP客户端自定义Cookie实现实例几乎所有脚本语言都提供了方便的 HTTP 客户端处理的功能，Python 也不例外，使
- 索引初识最普通的情况，是为出现在where子句的字段建一个索引。为方便讲述，我们先建立一个如下的表。CREATE TABLE mytable
- ASP中转换unicode编码为GB编码方法：<％ function urldecode(encodestr)
- 这个问题是如何在一些场景下使用断言表达式，通常会有人误用它，所以我决定写一篇文章来说明何时使用断言，什么时候不用。为那些还不清楚它的人，Py
- 如下所示：stopword = ''str = ''for line in iter(raw_input,
- 上篇文章给大家介绍了 在 webpack 中使用 ECharts的实例详解 ，可以点击查看。1. 使用NPM安装（全局引入）执行下面的命令：
- 前言接口在软件工程扮演重要角色，随着应用程序的功能不断扩展，代码库的更新和改变也难以管理。在许多情况下，会发现有一些看起来非常相似，但却不相
- xml问价中查询数据中包含list，需要使用collection<resultMap id="XX" type=&
- 通常报表都存储在一个地方，但是在单独的文件中。随着文件数量的增加在结构和搜索上也困难多多，针对这个问题，我遇到了一个非常有趣的报表对象属性&
- “站内信”不同于电子邮件，电子邮件通过专门的邮件服务器发送、保存。而“站内信”是系统内的消息，说白了，“站内信”的实现，就是通过数据库插入记
- Django 自称是“最适合开发有限期的完美WEB框架”。本文参考《Django web开发指南》，快速搭建一个blog 出来，在中间涉及诸
- 通常情况下我们在更新数据时需要先从数据库里将原数据取出后放在内存里，然后编辑某些字段或属性，最后提交更新数据库。使用F方法则可以帮助我们避免
- 这个东西很重要，可以经常用在项目当中，所以我们单独拿出来进行讲解。在使用它之前我们需要导包： go get golang.org/
- 本文主要介绍了Python利用numpy实现三层神经网络的示例代码，分享给大家，具体如下：其实神经网络很好实现，稍微有点基础的基本都可以实现
- nginx简单配置php服务（多个）摘要：大部分网站开发语言都要运行在服务器，比如主流的nginx、apache等等，部署服务器环境对于大部
