第一次看见Python的运行感觉就让我想起了matlab,于是就上网嗖嗖他在矩阵方面的运算如何,如果不想安装Matlab那么大的软件,而你又只是想计算些矩阵,python绝对够用!尤其在Linux下太方便了

Python使用NumPy包完成了对N-维数组的快速便捷操作。使用这个包,需要导入numpy。SciPy包以NumPy包为基础,大大的扩展了numpy的能力。为了使用的方便,scipy包在最外层名字空间中包括了所有的numpy内容,因此只要导入了scipy,不必在单独导入numpy了!但是为了明确哪些是numpy中实现的,哪些是scipy中实现的,本文还是进行了区分。以下默认已经:import

numpy as np 以及 impor scipy as sp

下面简要介绍Python和MATLAB处理数学问题的几个不同点。1.MATLAB的基本是矩阵,而numpy的基本类型是多为数组,把matrix看做是array的子类。2.MATLAB的索引从1开始,而numpy从0开始。

1.建立矩阵

a1=np.array([1,2,3],dtype=int)  #建立一个一维数组,数据类型是int。也可以不指定数据类型,使用默认。几乎所有的数组建立函数都可以指定数据类型,即dtype的取值。

a2=np.array([[1,2,3],[2,3,4]])  #建立一个二维数组。此处和MATLAB的二维数组(矩阵)的建立有很大差别。

同样,numpy中也有很多内置的特殊矩阵:

b1=np.zeros((2,3))  #生成一个2行3列的全0矩阵。注意,参数是一个tuple:(2,3),所以有两个括号。完整的形式为:zeros(shape,dtype=)。相同的结构,有ones()建立全1矩阵。empty()建立一个空矩阵,使用内存中的随机值来填充这个矩阵。

b2=identity(n)  #建立n*n的单位阵,这只能是一个方阵。

b3=eye(N,M=None,k=0)  #建立一个对角线是1其余值为0的矩阵,用k指定对角线的位置。M默认None。

此外,numpy中还提供了几个like函数,即按照某一个已知的数组的规模(几行几列)建立同样规模的特殊数组。这样的函数有zeros_like()、empty_like()、ones_like(),它们的参数均为如此形式:zeros_like(a,dtype=),其中,a是一个已知的数组。

c1=np.arange(2,3,0.1)  #起点,终点,步长值。含起点值,不含终点值。

c2=np.linspace(1,4,10)  #起点,终点,区间内点数。起点终点均包括在内。同理,有logspace()函数

d1=np.linalg.companion(a)  #伴随矩阵

d2=np.linalg.triu()/tril()  #作用同MATLAB中的同名函数

e1=np.random.rand(3,2)  #产生一个3行2列的随机数组。同一空间下,有randn()/randint()等多个随机函数

fliplr()/flipud()/rot90()  #功能类似MATLAB同名函数。

xx=np.roll(x,2)  #roll()是循环移位函数。此调用表示向右循环移动2位。

2.数组的特征信息

先假设已经存在一个N维数组X了,那么可以得到X的一些属性,这些属性可以在输入X和一个.之后,按tab键查看提示。这里明显看到了Python面向对象的特征。

X.flags  #数组的存储情况信息。

X.shape  #结果是一个tuple,返回本数组的行数、列数、……

X.ndim  #数组的维数,结果是一个数

X.size  #数组中元素的数量

X.itemsize  #数组中的数据项的所占内存空间大小

X.dtype  #数据类型

X.T  #如果X是矩阵,发挥的是X的转置矩阵

X.trace()  #计算X的迹

np.linalg.det(a)  #返回的是矩阵a的行列式

np.linalg.norm(a,ord=None)  #计算矩阵a的范数

np.linalg.eig(a)  #矩阵a的特征值和特征向量

np.linalg.cond(a,p=None)  #矩阵a的条件数

np.linalg.inv(a)  #矩阵a的逆矩阵

3.矩阵分解

常见的矩阵分解函数,numpy.linalg均已经提供。比如cholesky()/qr()/svd()/lu()/schur()等。某些

4.矩阵运算

np.dot(a,b)用来计算数组的点积;vdot(a,b)专门计算矢量的点积,和dot()的区别在于对complex数据类型的处理不一样;innner(a,b)用来计算内积;outer(a,b)计算外积。

专门处理矩阵的数学函数在numpy的子包linalg中定义。比如np.linalg.logm(A)计算矩阵A的对数。可见,这个处理和MATLAB是类似的,使用一个m后缀表示是矩阵的运算。在这个空间内可以使用的有cosm()/sinm()/signm()/sqrtm()等。其中常规exp()对应有三种矩阵形式:expm()使用Pade近似算法、expm2()使用特征值分析算法、expm3()使用泰勒级数算法。在numpy中,也有一个计算矩阵的函数:funm(A,func)。

5.索引

numpy中的数组索引形式和Python是一致的。如:

x=np.arange(10)

print x[2]  #单个元素,从前往后正向索引。注意下标是从0开始的。

print x[-2]  #从后往前索引。最后一个元素的下标是-1

print x[2:5]  #多个元素,左闭右开,默认步长值是1

print x[:-7]  #多个元素,从后向前,制定了结束的位置,使用默认步长值

print x[1:7:2]  #指定步长值

x.shape=(2,5)  #x的shape属性被重新赋值,要求就是元素个数不变。2*5=10

print x[1,3]  #二维数组索引单个元素,第2行第4列的那个元素

print x[0]  #第一行所有的元素

y=np.arange(35).reshape(5,7)  #reshape()函数用于改变数组的维度

print y[1:5:2,::2]  #选择二维数组中的某些符合条件的元素

python构建矩阵_python矩阵运算 | 学步园相关推荐

  1. python求解矩阵_python矩阵运算

    第一次看见Python的运行感觉就让我想起了matlab,于是就上网嗖嗖他在矩阵方面的运算如何,如果不想安装Matlab那么大的软件,而你又只是想计算些矩阵,python绝对够用!尤其在Linux下太 ...

  2. python分形算法_python 分形 | 学步园

    代码如下: #!/usr/bin/env python # -*-coding:UTF-8-*-# from common import * import Image, ImageFilter, ma ...

  3. python创建矩阵_Python创建对称矩阵的方法示例【基于numpy模块】

    本文实例讲述了Python创建对称矩阵的方法.分享给大家供大家参考,具体如下: 对称(实对称)矩阵也即: step 1:创建一个方阵 >>> import numpy as np & ...

  4. python开发服务程序_Python 编写Windows服务程序:将Python作为Windows服务启动 | 学步园...

    Python程序作为Windows服务启动,需要安装pywin32包.下载路径: #-*- coding:utf-8 -*- import win32serviceutil import win32s ...

  5. python numpy矩阵_python中numpy矩阵运算操作大全(非常全)!

    用vuejs实现一个todolist项目 用vue.js实现一个todolist项目:input输入框输入的值会呈现在下方,并且会保存在localStorage里面,而且下方的列表点击之后也会有变化: ...

  6. python 共轭转置_python矩阵运算,转置,逆运算,共轭矩阵实例

    我就废话不多说了,大家还是直接看代码吧! #先定义两个矩阵 X=np.array([[1,2104,5,1,45],[1,1416,3,2,40],[1,1534,3,2,30],[1,852,2,1 ...

  7. python构建矩阵 x y_生成Python函数一半没问题,当前最正统的代码生成是什么样的?...

    机器之心原创 参与:思源 大家都说深度神经网络能力很强,那么从函数注释生成函数代码,以及从函数代码总结函数注释这种最基础的代码任务到底能不能行?像 Python.Java 这样的通用高级语言,到底在代 ...

  8. python 轮廓矩阵_python – 在3D numpy矩阵中绘制给定值的曲面轮廓

    我有三个对应于xyz坐标空间的3D网格矩阵(X,Y,Z). 我还有一个3D Numpy矩阵A,其中A [i,j,k]包含一个与点(x,y,z)相关联的浮点数,其中x = X [i,j,k],y = Y ...

  9. python创建矩阵_Python创建二维数组的正确姿势

    List (列表)是 Python 中最基本的数据结构.在用法上,它有点类似数组,因为每个列表都有一个下标,下标从 0 开始.因此,我们可以使用 list[1] 来获取下标对应的值.如果我们深入下列表 ...

最新文章

  1. 记一次TIME_WAIT网络故障
  2. Python常见十六个错误集合,你知道那些?
  3. SpringCloud Zuul(九)之路由自动刷新原理
  4. 如何平衡存储系统的一致性和可用性?
  5. ubuntu mysql 二进制安装_ubuntu18.04手动安装二进制MySQL8.0
  6. 微信公布10月朋友圈十大谣言 包括牙膏能杀灭幽门螺杆菌等
  7. 第一个程序python-HelloWorld
  8. 关于vim的误删除经历
  9. tshark/wireshark抓包小结
  10. eclipse对mysql进行增删改查_用eclipse编写的 能对MySQL数据库中表实现增删改查功能...
  11. 小米 samba linux,U-BOOT移植准备之(一)Samba服务器
  12. python nlpir_NLPIR(ICTCLAS 2013)分词工具Python封装
  13. 【Excel】数据的输入
  14. 童年计算机课企鹅游戏,这些童年游戏,你一定玩过
  15. EF Core学习之路02
  16. 毛孔很大很难看该怎么处理
  17. open3d画3D矩形框
  18. 蓝牙耳机哪款性价比高?2023年高性价比蓝牙耳机盘点
  19. 寄存器总结之通用寄存器
  20. 咸鱼ZTMS实例—加速传感器制作平衡仪

热门文章

  1. 大数据基础(5) - Raft协议
  2. python - dtype
  3. 中国消费电子高频主轴市场趋势报告、技术动态创新及市场预测
  4. upload类型题目总结
  5. 怎样直接下载网页中的视频
  6. java无法解析符号什么意思,Intellij:无法解析符号'springframework'
  7. 智能问答客服系统的设计和开发(ssm(spring mvc+mybatis)+netty4)
  8. 解决vlookup函数查找参数不在第一列_if({1,0},,)用法
  9. matlab状态空间模型,Matlab做状态空间辨识2
  10. 关于TB67S109AFNAG的应用分享及应用