python矩阵运算函数_Python之Numpy基础

以前总认为Numpy是渣渣，直到深入接触以后才知道功能这么强大。堪比Matlab啊。果然是人生苦短，我用Python。所以本文作为一个记录&笔记，文章内容大多数取自网络以&官网快速入门等(文末有参考链接，如有侵权请联系本人改正)，希望可以帮助大家快速入门Numpy。如果你有Matlab基础，那么你能很快看懂本文！！！

一个栗子

创建矩阵

对于Python中的numpy模块，一般用其提供的ndarray对象。 创建一个ndarray对象很简单，只要将一个list作为参数即可。 例如:

矩阵行数列数

矩阵按行列选取

矩阵的截取和list相同，可以通过[](方括号)来截取

矩阵按条件截取

按条件截取应用较多的是对矩阵中满足一定条件的元素变成特定的值。 例如将矩阵中大于6的元素变成0。

Stacking together different arrays

矩阵的合并可以通过numpy中的hstack方法和vstack方法实现:

矩阵的合并也可以通过concatenatef方法。

np.concatenate( (a1,a2), axis=0 ) 等价于 np.vstack( (a1,a2) )

np.concatenate( (a1,a2), axis=1 ) 等价于 np.hstack( (a1,a2) )

通过函数创建矩阵

arange

linspace/ logspace

ones、zeros、eye、empty

ones创建全1矩阵 ，zeros创建全0矩阵 ，eye创建单位矩阵 ，empty创建空矩阵(实际有值)

fromstring

fromstring()方法可以将字符串转化成ndarray对象，需要将字符串数字化时这个方法比较有用，可以获得字符串的ascii码序列。

random

fromfunction

fromfunction()方法可以根据矩阵的行号列号生成矩阵的元素。 例如创建一个矩阵，矩阵中的每个元素都为行号和列号的和。

矩阵的运算

常用矩阵运算符

Numpy中的ndarray对象重载了许多运算符，使用这些运算符可以完成矩阵间对应元素的运算。

运算符说明+矩阵对应元素相加-矩阵对应元素相减*矩阵对应元素相乘/矩阵对应元素相除，如果都是整数则取商%矩阵对应元素相除后取余数**矩阵每个元素都取n次方，如**2：每个元素都取平方

常用矩阵函数

同样地，numpy中也定义了许多函数，使用这些函数可以将函数作用于矩阵中的每个元素。 表格中默认导入了numpy模块，即 。a为ndarray对象。

常用矩阵函数说明np.sin(a)对矩阵a中每个元素取正弦,sin(x)np.cos(a)对矩阵a中每个元素取余弦,cos(x)np.tan(a)对矩阵a中每个元素取正切,tan(x)np.arcsin(a)对矩阵a中每个元素取反正弦,arcsin(x)np.arccos(a)对矩阵a中每个元素取反余弦,arccos(x)np.arctan(a)对矩阵a中每个元素取反正切,arctan(x)np.exp(a)对矩阵a中每个元素取指数函数,exnp.sqrt(a)对矩阵a中每个元素开根号

当矩阵中的元素不在函数定义域范围内，会产生RuntimeWarning，结果为nan(not a number)

矩阵乘法(点乘)

矩阵乘法必须满足矩阵乘法的条件，即第一个矩阵的列数等于第二个矩阵的行数。 矩阵乘法的函数为 dot 。

矩阵的转置 a.T

矩阵的转置还有更简单的方法，就是a.T。

矩阵的逆

设A是数域上的一个n阶方阵，若在相同数域上存在另一个n阶矩阵B，使得：AB=BA=E。 则我们称B是A的逆矩阵，而A则被称为可逆矩阵。

求矩阵的逆需要先导入，用的函数来求逆。矩阵求逆的条件是矩阵应该是方阵。

矩阵信息获取(如均值等)

最值

获得矩阵中元素最大最小值的函数分别是和，可以获得整个矩阵、行或列的最大最小值。

平均值

获得矩阵中元素的平均值可以通过函数。同样地，可以获得整个矩阵、行或列的平均值。

方差

方差的函数为,方差函数相当于函数,其中x为矩阵。

标准差

标准差的函数为。 相当于，或相当于。

中值

中值指的是将序列按大小顺序排列后，排在中间的那个值，如果有偶数个数，则是排在中间两个数的平均值。中值的函数是median()，调用方法为numpy.median(x,[axis])，axis可指定轴方向，默认axis=None，对所有数取中值。

求和

矩阵求和的函数是sum()，可以对行，列，或整个矩阵求和

累积和

某位置累积和指的是该位置之前(包括该位置)所有元素的和。例如序列[1,2,3,4,5]，其累计和为[1,3,6,10,15]，即第一个元素为1，第二个元素为1+2=3，……，第五个元素为1+2+3+4+5=15。矩阵求累积和的函数是cumsum()，可以对行，列，或整个矩阵求累积和。

极差

百分位数

序号参数及描述1. 输入数组2. 要计算的百分位数，在 0 ~ 100 之间3. 沿着它计算百分位数的轴

加权平均值

Shape Manipulation

Changing the shape of an array

数组的形状可以用以下方式改变。Note that the following three commands all return a modified array, but do not change the original array:

The function returns its argument with a modified shape, whereas the method modifies the array itself:

If a dimension is given as -1 in a reshaping operation, the other dimensions are automatically calculated:

Splitting one array into several smaller ones

Using , you can split an array along its horizontal axis, either by specifying the number of equally shaped arrays to return, or by specifying the columns after which the division should occur:

Copies and Views

When operating and manipulating arrays, their data is sometimes copied into a new array and sometimes not. This is often a source of confusion for beginners. There are three cases:

No Copy At All

a = b,改变b就相当于改变a,或者相反。

View or Shallow Copy

Different array objects can share the same data. The method creates a new array object that looks at the same data.

Slicing an array returns a view of it:

Deep Copy

The method makes a complete copy of the array and its data.

关于有三种情况，完全不复制、视图(view)或者叫浅复制()和深复制()。而就属于第二种，即视图，这本质上是一种切片操作()，所有的切片操作返回的都是视图。具体来说，会创建一个新的对象(所以说 id 和不一样)，但是的数据完全来自于，和保持完全一致，换句话说，b的数据完全由a保管，他们两个的数据变化是一致的，可以看下面的示例:

和 的差别就在于后者会创建新的对象，前者不会。两种方式都会导致和的数据相互影响。要想不让的改动影响到，可以使用深复制：

曼德勃罗