基于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
投稿猜你喜欢
- '-----------------------------------------------------------
- 1.视图的概述 视图其实就是一条查询sql语句，用于显示一个或多个表或其他视图中的相关数据。视图将一个查询的结果作为一个表来使用，因此视图可
- 下拉框连动JavaScript代码，市区二级联动多级联动下拉选择框，动态获取下一级琥珀无限级联动菜单-JavaScript版 <htm
- 背景：pony是公司的首席体验官、首席产品经理。这次在产品峰会上pony将自己平时经验的积累与大家交流，体验较细。这次分享研发管理部，设计中
- 这个分页使用的是0游标，也就是Rs.Open Sql,Conn,0,1。但是感觉也快不了多少，10万条数据的分页时间300多豪秒之间。风格A
本文实例讲述了Python常用模块sys,os,time,random功能与用法。分享给大家供大家参考，具体如下：sys:介绍：主要包含涉及- 当鼠标滑过一个图片时，图片会变成另外一张图片的效果是怎么制作出来的呢？相对一些不熟悉HTML代码的朋友来说，可以使用网页编辑软件自带的动作来
- 假设mysql 安装在c：盘，mysql数据库的用户名是root，密码是123456，数据库名是database_name，在d：盘根目录下
- 如何在浏览器地址栏显示自己的Favicons小图标？显示效果如下列网站的图标：Google网易163上面两个网站都设置了自己的Favicon
- Linux环境MySQL服务器级优化讲解 摘要：本节简单介绍了如何在服务器级优化数据库的性能
- 阅读之前：在看文章具体内容之前，希望你可以 先打开IE8,打开http://www.taobao.com,然后在地址栏里输入：javascr
- 导言：本文一步一步手把手教你打造一个极酷的三层分离的标准滑动门导航菜单，从思路、原理、步骤，手段可谓“无所不用其极”，即便你是菜鸟，相信你看
- 今天做自定义的404页面，需要垂直水平居中，垂直居中不像水平居中那么简单，呵呵。这种效果会在企业站首页，或一些提示性页面经常用到，今天做了个
- PHP是一种面向对象的编程语言，它允许开发者使用面向对象的编程技术来构建复杂的应用程序。下面是一些关于PHP面向对象编程的讲解：类与对象类是
- 前言在应用程序开发的过程中，有的时候需要在代码中保存一些机密的信息，比如加密密钥，字符串，或者是用户名密码等。通常的做法是保存到一个配置文件
- 我们可使用Haskeys属性判别每个条目是否为一个集合，遍历完整的Request.Cookies集合，以来取得所有cookie的列表及其值：
- 在填写表单的时候为了让用户有更好的体验，有时需要根据用户的输入出现提示共用户选择，我这个就是根据输入的内容，从数据库取出相关内容以供选择，这
- 本文实例讲述了PHP实现ASCII码与字符串相互转换的方法。分享给大家供大家参考，具体如下：<?phpclass ascii { &n
- 1.现在我本机系统已内置python2.62.下载进行源码安装复制链接下载到/root/mypackage，解压接着mkdir /usr/l
- js汉字简繁转换源代码：<html> <head> <title>汉字简繁转换工具_asp之家</
