python-numpy一些方法总结

2023-11-03

1、multiply
例子:

x1=[1,2,3];x2=[4,5,6]
print multiply(x1,x2)

输出:

[ 4 10 18]

multiply函数得到的结果是对应位置上面的元素进行相乘。

2、std 标准方差 ,var 方差
例子:

b=[1,3,5,6]
print var(b)
print power(std(b),2)
ll=[[1,2,3,4,5,6],[3,4,5,6,7,8]]
print var(ll[0])
print var(ll,0)#第二个参数为0,表示按列求方差
print var(ll,1)#第二个参数为1,表示按行求方差

输出:

3.6875
3.6875
2.91666666667
[ 1.  1.  1.  1.  1.  1.]
[ 2.91666667  2.91666667]

var输出向量的方差,std输出向量的标准方差。

3、mean
例子:

b=[1,3,5,6]
print mean(b)
print mean(ll) #全部元素求均值
print mean(ll,0)#按列求均值
print mean(ll,1)#按行求均值

结果:

3.75
4.5
[ 2.  3.  4.  5.  6.  7.]
[ 3.5  5.5]

mean函数得到向量的均值。

4、sum
求和
例子

x=[[0, 1, 2],[2, 1, 0]]
b=[1,3,5,6]
print sum(b)
print sum(x)

输出:

15
6

sum求向量的和。也可以求矩阵所有元素的和

5、cov()

例子:

b=[1,3,5,6]
print cov(b)
print sum((multiply(b,b))-mean(b)*mean(b))/3

输出:

4.91666666667
4.91666666667

cov求的是样本协方差
公式如下:
这里写图片描述
公式一是样本均值
公式二是样本方差
公式三是样本协方差。
上面的函数cov就是使用样本协方差求得的。

例子:

x=[[0, 1, 2],[2, 1, 0]]
print cov(x)
print sum((multiply(x[0],x[1]))-mean(x[0])*mean(x[1]))/2

输出结果:

[[ 1. -1.]
 [-1.  1.]]
-1.0

cov的参数是矩阵,输出结果也是矩阵!输出的矩阵为协方差矩阵!
计算过程如下:
这里写图片描述
截图来源维基百科:协方差矩阵

6、corrcoef
该函数得到相关系数矩阵。
例子:

    vc=[1,2,39,0,8]
    vb=[1,2,38,0,8]
    print mean(multiply((vc-mean(vc)),(vb-mean(vb))))/(std(vb)*std(vc))
    #corrcoef得到相关系数矩阵(向量的相似程度)
    print corrcoef(vc,vb)

输出结果:

0.999986231331
[[ 1.          0.99998623]
 [ 0.99998623  1.        ]]

相关系数公式:
这里写图片描述

对应着公式理解上面的代码,应该是很容易的。

7、vdot 向量的点积
例子:

    #vdot 返回两向量的点积
    l1=[1,2,3]
    l2=[4,5,6]
    print  vdot(l1,l2)

结果:

32

点积没什么说的,对应位置相乘求和。

8、mat
例子:

    l1=[1,2,3]
    l2=[4,5,6]
    ll=[l1,l2]
    print  vdot(l1,l2)
    print mat(l1)*mat(l2).T
    print mat(ll)

结果:

32
[[32]]
[[1 2 3]
 [4 5 6]]

mat函数把列表转换成矩阵形式。这在矩阵运算中不可缺少。

9、shape
例子:

    #矩阵有一个shape属性,是一个(行,列)形式的元组
    a = array([[1,2,3],[4,5,6]])
    print a.shape

结果:

(2, 3)

shape返回矩阵的行列数

10、ones
例子:

    #返回按要求的矩阵
    ones = ones((2,1))
    print ones

结果:

[[ 1.]
 [ 1.]]

ones返回指定行列数的全一矩阵

11、xrange
例子:

    for i in xrange(3):
        print i
    test=[1,2,3,4]
    print test[:]
    print test[2:3]
    for i in xrange(2,5):
        print i

结果:

0
1
2
[1, 2, 3, 4]
[3]
2
3
4

xrange用于循环中,参数为一个整数的话,可循环遍历小于该参数的值。两个参数,则循环遍历两个整数之间的值。
test[:]则表示获取test列表中的所有元素。
test[2:3]则表示获取从第2个位置到第三个位置间的元素。

12、strptime
例子:

import time
from datetime import datetime,date
dd = datetime.strptime('2014-04-03T10:53:49.875Z', "%Y-%m-%dT%H:%M:%S.%fZ")
    print time.mktime(dd.timetuple())#1396493629.0

结果:

1396493629.0

strptime把时间按照指定的格式转换。处理时间的时候会用到这个函数。

13、tuple和数组转换成字符串
例子:

    tuple =(1,2,3)
    print tuple[len(tuple)-1]
    print tuple[-1]
    print 9.99.__repr__()
    print str(9.99)

结果:

3
3
9.99
9.99

上面的示例中tuple是一个元组,访问元素的时候,可以通过[index]这种方式访问。
访问最后一个元素的时候,可以通过[-1]访问.
那么访问倒数第二个元素:

print tuple[-2]
输出结果:2

另外,数字转换成字符串有两种方式:

_ _ repr _ _()
str()

14、transpose和.T
例子:

    aa =[[1],[2],[3]]
    aa= mat(aa)#将列表转换成矩阵,并存放在aa中
    print aa
    print aa.transpose()#将矩阵进行转置
    print aa #transpose()进行矩阵的转置,aa并没有改变
    print '*'*20
    print aa.T
    print aa #T转置矩阵,也没有发生改变

输出结果:

[[1]
 [2]
 [3]]
[[1 2 3]]
[[1]
 [2]
 [3]]
********************
[[1 2 3]]
[[1]
 [2]
 [3]]

有两种方式实现矩阵的转置。

15、zeros() ones()
例子:

    print zeros((2,1))
    print ones((2,3))#参数指明了矩阵的行列数

结果:

[[ 0.]
 [ 0.]]
[[ 1.  1.  1.]
 [ 1.  1.  1.]]

zeros返回指定行列全零矩阵
ones返回指定行列的全一矩阵

16、列表 数组 linspace
例子:

    #列表和数组的区别:
    #列表: [1, 2, 3, 4]
    #数组: [1 2 3 4]
    print '*'*20
    ll =[1,2,3,4]
    print '列表:',ll
    arr = array(ll)
    print '数组:',arr

    print '*'*20
    print linspace(0,3,6) #返回的是数组

输出结果:

********************
列表: [1, 2, 3, 4]
数组: [1 2 3 4]
********************
[ 0.   0.6  1.2  1.8  2.4  3. ]

数组中间元素没有分隔符。列表逗号分割。
linspace返回的是指定的开始结束位置的指定个数的数。

linspace(0,3,6)返回0到3之间6个数字,且间隔均匀

17、argsort 排序索引
例子:

    print '**************数组排序问题****************'
    #数组的构建问题,初始化使用array()
    ary=array(zeros(4))
    ary[0]=0.1
    ary[1]= 0.6
    ary[2]= 0.5
    ary[3]= 0.7
    #有-号,降序排列
    #无-号,升序排列
    sortindex = argsort(ary)
    for id in sortindex:
        print '索引:',id
    for i in ary:
        print i

结果:

**************数组排序问题****************
索引: 0
索引: 2
索引: 1
索引: 3
0.1
0.6
0.5
0.7

argsort函数返回数组按照从小到大排序的位置的索引。

sortindex = argsort(-ary)
for id in sortindex:
print ‘索引:’,id
for i in sortindex:
print ary[i]

输出结果:
索引: 1
索引: 2
索引: 0
0.7
0.6
0.5
0.1

18、 [:,:]矩阵元素切片

    #矩阵元素的获取
    ll = [[1,2,3],[4,5,6],[7,8,9]]
    #获取第二行第0个元素
    print mat(ll)[2,0]
    #第一个冒号代表获取行的起止行号
    #第二个冒号代表获取列的起止行号
    print mat(ll)[:,:]

结果:

7
[[1 2 3]
 [4 5 6]
 [7 8 9]]

19、diag
例子:

    #构建对角矩阵 
    #diag()参数为列表即可
    dd = [1,2,3]
    dilogg = diag(dd)
    print 'diag=',dilogg

结果:

diag= [[1 0 0]
     [0 2 0]
     [0 0 3]]

diag构建对角矩阵

20、linalg.inv和.I 求逆矩阵
例子:

    dd = [1,2,3]
    dilogg = diag(dd)
    print 'diag=',dilogg
    print 'dd:',linalg.inv(dilogg)
    print 'I=',mat(dilogg).I
    ll = [[1,2,3],[4,5,6],[7,8,9]]
    #求逆矩阵
    lv = linalg.inv(mat(ll))
    print 'inv:',lv
    print 'I:',mat(ll).I

结果:

diag= [[1 0 0]
 [0 2 0]
 [0 0 3]]
dd: [[ 1.          0.          0.        ]
 [ 0.          0.5         0.        ]
 [ 0.          0.          0.33333333]]
I= [[ 1.          0.          0.        ]
 [ 0.          0.5         0.        ]
 [ 0.          0.          0.33333333]]
inv: [[ -4.50359963e+15   9.00719925e+15  -4.50359963e+15]
 [  9.00719925e+15  -1.80143985e+16   9.00719925e+15]
 [ -4.50359963e+15   9.00719925e+15  -4.50359963e+15]]
I: [[ -4.50359963e+15   9.00719925e+15  -4.50359963e+15]
 [  9.00719925e+15  -1.80143985e+16   9.00719925e+15]
 [ -4.50359963e+15   9.00719925e+15  -4.50359963e+15]]

上面求逆矩阵有两种方式。
linalg.inv()和 矩阵.I
两种方式均可。

21、dot矩阵点积
例子:

    ll = [[1,2,3],[4,5,6],[7,8,9]]
    ld = dot(ll,ll)
    print 'dot:',ld
    print mat(ll)*mat(ll)

结果:

dot: [[ 30  36  42]
 [ 66  81  96]
 [102 126 150]]
[[ 30  36  42]
 [ 66  81  96]
 [102 126 150]]

22、eye 单元矩阵
例子:

print 'eye:',eye(2)#单元矩阵

结果:

eye: [[ 1.  0.]
 [ 0.  1.]]

23 eig 矩阵的特征值和特征向量
例子:

    A=mat([[1,0,0,0,2],[0,0,3,0,0],[0,0,0,0,0],[0,4,0,0,0]])
    U =A*A.T
    lamda,hU=linalg.eig(U)
    print 'hU:',hU
    print lamda

结果:

hU: [[ 1.  0.  0.  0.]
 [ 0.  1.  0.  0.]
 [ 0.  0.  1.  0.]
 [ 0.  0.  0.  1.]]
 [  5.   9.   0.  16.]

24 sorted 排序
例子:

    ll=[8,0,3,6,1,0,5,3,8,9]
    print sorted(ll,reverse=True) #降序
    print sorted(ll,reverse=False) #升序

结果:

[9, 8, 8, 6, 5, 3, 3, 1, 0, 0]
[0, 0, 1, 3, 3, 5, 6, 8, 8, 9]

25、linalg.svd 奇异值分解
例子:

    A=mat([[1,0,0,0,2],[0,0,3,0,0],[0,0,0,0,0],[0,4,0,0,0]])
    U,S,VT =linalg.svd(A)
    print 'U:',U
    print 'V:',VT
    print 's:',S
    print '===',U*(mat(diag(S))*eye(4,5))*VT

结果:

U: [[ 0.  0.  1.  0.]
 [ 0.  1.  0.  0.]
 [ 0.  0.  0. -1.]
 [ 1.  0.  0.  0.]]
V: [[ 0.          1.         -0.         -0.         -0.        ]
 [ 0.          0.          1.         -0.         -0.        ]
 [ 0.4472136   0.         -0.         -0.          0.89442719]
 [ 0.          0.         -0.          1.         -0.        ]
 [-0.89442719  0.         -0.         -0.          0.4472136 ]]
s: [ 4.          3.          2.23606798  0.        ]
=== [[ 1.  0.  0.  0.  2.]
 [ 0.  0.  3.  0.  0.]
 [ 0.  0.  0.  0.  0.]
 [ 0.  4.  0.  0.  0.]]

关于什么是奇异值分解,请参考奇异值分解

26、random.rand
例子:

    A=mat([[1,0,0,0,2],[0,0,3,0,0],[0,0,0,0,0],[0,4,0,0,0]])
    print A[:,1]

    #获取3*30-1之间的数字
    rr=random.rand(3,3)
    print rr
    print (rr-0.5)
    print 2.0*(rr-0.5)

结果:

[[0]
 [0]
 [0]
 [4]]
[[ 0.30446153  0.40653841  0.40143809]
 [ 0.77970727  0.57491894  0.85801586]
 [ 0.33509491  0.64652856  0.48276137]]
[[-0.19553847 -0.09346159 -0.09856191]
 [ 0.27970727  0.07491894  0.35801586]
 [-0.16490509  0.14652856 -0.01723863]]
[[-0.39107693 -0.18692318 -0.19712383]
 [ 0.55941453  0.14983789  0.71603172]
 [-0.32981018  0.29305712 -0.03447726]]

random.rand(3,3)随机获取3*3个0-1之间的数字

27、arange
例子:

    delta = 0.25
    x = arange(-3.0, 3.0, delta)
    print x

结果:

[-3.   -2.75 -2.5  -2.25 -2.   -1.75 -1.5  -1.25 -1.   -0.75 -0.5  -0.25
  0.    0.25  0.5   0.75  1.    1.25  1.5   1.75  2.    2.25  2.5   2.75]

arange获取指定起始位置,指定步长的一系列数。

28、nonzero()
例子:

    x =[[1,0,0,0,2],[0,0,3,0,0]]
    print x
    nz=nonzero(x)
    print nz
    print nz[0]

结果:

[[1, 0, 0, 0, 2], [0, 0, 3, 0, 0]]
(array([0, 0, 1]), array([0, 4, 2]))
[0 0 1]

nonzero()函数返回矩阵中非0元素的位置
nz的返回值意义如下:
第一行是所有非零数所在行值
第二行是所有非零值所在列值

29、获取指定位置的元素
例子:

    A=mat([[1,0,0,0,2],[0,0,3,0,0],[0,0,0,0,0],[0,4,0,0,0]])
    sample =A[0,:]
    print sample
    print sample[0]

    ll=mat([3,4,5])
    for i in range(5):
        if sum(ll==i):
            print i

结果:

[[1 0 0 0 2]]
[[1 0 0 0 2]]
3
4

A[0,:]的意义是获取第0行的所有元素
if sum(ll==i):的意义是只要i存在ll矩阵中,if就是True
ll必须是mat转换的矩阵。列表好像不行。并且还是单行的矩阵。

第二个例子:

    A=mat([[1,0,0,0,2],[0,0,3,0,0],[0,0,0,0,0],[0,4,0,0,0]])
    #根据ind序列索引获取矩阵A中的数据
    ind=[2,1,3,0]
    print A[ind,0]

结果:

[[0]
 [0]
 [0]
 [1]]

30、zip
例子:

    ll=[1,2,3,4,5,6]
    #可以互换指定区域的位置
    print ll[3:6]+ll[0:3]
    #成对获取x、y的值
    l1=[1,2,3]
    l2=[4,5,6]
    for x,y in zip(l1,l2):
        print x,y

结果:

[4, 5, 6, 1, 2, 3]
1 4
2 5
3 6

print ll[3:6]+ll[0:3]这行代码实现指定位置的元素进行交换。

31、chr函数,获取指定的字符
例子:

    #获取指定的字符
    for i in range(65,70):
        print str(chr(i))

结果:

A
B
C
D
E

32、random.shuffle
例子:

    ll=range(9)#返回列表
    print ll
    #shuffle函数随机打乱列表中的元素顺序
    print random.shuffle(ll)
    print ll

结果:

[0, 1, 2, 3, 4, 5, 6, 7, 8]
None
[8, 5, 1, 4, 2, 6, 0, 3, 7]

shuffle函数随机打乱列表中的元素顺序

33、vdot 点积
例子:

    #vdot 返回两向量的点积
    l1=[1,2,3]
    l2=[4,5,6]
    ll=[l1,l2]
    print vdot(l1,l2)
    print dot(l1,l2)
    print mat(l1)*mat(l2).T
    print mat(ll)

结果:

32
32
[[32]]
[[1 2 3]
 [4 5 6]]

vdot 和 dot都可以获取向量的点积。

区别:

    ll=[[1,2,3],[4,5,6],[1,2,3]]
    print dot(mat(ll),mat(ll).T)
    print vdot(mat(ll),mat(ll))

结果:

[[14 32 14]
 [32 77 32]
 [14 32 14]]
[[105]]

总结:对于向量来说,vdot dot向量的点积结果相同。
对于矩阵来说,dot是矩阵的点积。
而vdot是对应位置的元素乘积求和。
结果中的105=1^2+2^2+3^3+4^2+5^2+6^2+1^2+2^2+3^3

34、次方计算**
例子:

    print 3*2**2
    print 3*2**0.5
    print (3*2)**2
    print (3*2)**0.5

结果:

12
4.24264068712
36
2.44948974278

可见 **次方计算优先级要高!

35、max函数
例子:

    ll=[3,4,6,2,89,9,3,2]
    print max(ll)
    l2=[[3,4,6,2,89,9,3,2],[3,6,7,8,983,3,5,6]]
    print max(l2[0])
    print max(l2)

结果:

89
89
[3, 6, 7, 8, 983, 3, 5, 6]

对于矩阵的最大值的求法,还不知道怎么搞。稍后更新一下。

36、mgrid函数

    #开始值,结束值,步长。如果步长为虚数,表示产生的个数长度
    print mgrid[-5:5:3j] #结果:[-5.  0.  5.]
    print mgrid[-5:5:3] #结果:[-5 -2  1  4]
    print '*'*20
    print mgrid[-5:5:3j,-5:5:3j]
    print '*'*20
    print mgrid[-5:5:3,-5:5:3]

结果:

[-5.  0.  5.]
[-5 -2  1  4]
********************
[[[-5. -5. -5.]
  [ 0.  0.  0.]
  [ 5.  5.  5.]]

 [[-5.  0.  5.]
  [-5.  0.  5.]
  [-5.  0.  5.]]]
 ********************
 [[[-5 -5 -5 -5]
      [-2 -2 -2 -2]
      [ 1  1  1  1]
      [ 4  4  4  4]]

     [[-5 -2  1  4]
      [-5 -2  1  4]
      [-5 -2  1  4]
      [-5 -2  1  4]]]

mgrid函数说明:开始值,结束值,步长。如果步长为虚数,表示产生的个数长度。
对照结果即可知道函数的功能。

37、ogrid函数
例子:

    print ogrid[-5:5:3j]
    print ogrid[-5:5:3]
    print ogrid[-5:5:3j,-5:5:3j]
    print ogrid[-5:5:3,-5:5:3]

结果:

[-5.  0.  5.]
[-5 -2  1  4]
[array([[-5.],
       [ 0.],
       [ 5.]]), array([[-5.,  0.,  5.]])]
[array([[-5],
       [-2],
       [ 1],
       [ 4]]), array([[-5, -2,  1,  4]])]

ogrid跟mgrid一样。不同的是,一个数矩阵,一个数组形式。第三个擦数如果是实数,说明是步长。如果是虚数,说明是个数。
ogrid[-5:5:3,-5:5:3]第一部分产生多行一列,第二部分产生一行多列。这与mgrid不同。

38、random函数

    print random.seed(1)
    #要每次产生随机数相同就要设置种子,相同种子数的Random对象,相同次数生成的随机数字是完全相同的
    #用于生成一个指定范围内的随机符点数
    print random.uniform(-1,1,5)
    #结果:[ 0.40254497 -0.42350395 -0.67640645 -0.54075394 -0.99584028]
    #均匀分布
    #用于生成一个01的随机符点数: 0 <= n < 1.0
    print random.random()
    print random.random(5)
    print random.rand(2,3)#23列
结果:
    None
0.417022004703
[  7.20324493e-01   1.14374817e-04   3.02332573e-01   1.46755891e-01
   9.23385948e-02]
[[ 0.18626021  0.34556073  0.39676747]
 [ 0.53881673  0.41919451  0.6852195 ]]
    #用于生成一个指定范围内的整数。其中参数a是下限,参数b是上限,生成的随机数n: a <= n <= b
    print random.randint(5, 10)
    print random.randint(5,10,size=(5,5))
结果:
    None
8
[[9 5 6 8 5]
 [5 6 9 9 6]
 [7 9 7 9 8]
 [9 7 9 7 9]
 [6 6 5 6 6]]

更多关于random的函数的解析,请参看这篇文章:http://blog.csdn.net/pipisorry/article/details/39086463

39、tofile、fromfile
例子:

    a =arange(0,12)
    a.shape = 3,4
    print a
    a.tofile("a.bin")
    b = fromfile("a.bin", dtype=float) # 按照float类型读入数据
    print b # 读入的数据是错误的

    print a.dtype # 查看a的dtype
    b = fromfile("a.bin", dtype=int32) # 按照int32类型读入数据
    print b # 数据是一维的
    b.shape = 3, 4 # 按照a的shape修改b的shape
    print b

结果:

None
[[ 0  1  2  3]
 [ 4  5  6  7]
 [ 8  9 10 11]]
[  2.12199579e-314   6.36598737e-314   1.06099790e-313   1.48539705e-313
   1.90979621e-313   2.33419537e-313]
int32
[ 0  1  2  3  4  5  6  7  8  9 10 11]
[[ 0  1  2  3]
 [ 4  5  6  7]
 [ 8  9 10 11]]

使用数组的方法函数tofile可以方便地将数组中数据以二进制的格式写进文件。tofile输出的数据没有格式,因此用numpy.fromfile读回来的时候需要自己格式化数据。
Note:
1. 读入的时候设置正确的dtype和shape才能保证数据一致。
并且tofile函数不管数组的排列顺序是C语言格式的还是Fortran语言格式的,统一使用C语言格式输出。
2. sep关键字参数:此外如果fromfile和tofile函数调用时指定了sep关键字参数的话,
数组将以文本格式输入输出。{这样就可以通过notepad++打开查看, 不过数据是一行显示,不便于查看}
user_item_mat.tofile(user_item_mat_filename, sep=’ ‘)

40、numpy.load和numpy.save
例子:

     a =arange(0,12)
    a.shape = 3,4
    print a
    save('a.npy',a)
    c=load('a.npy')
    print c

结果:

None
[[ 0  1  2  3]
 [ 4  5  6  7]
 [ 8  9 10 11]]
[[ 0  1  2  3]
 [ 4  5  6  7]
 [ 8  9 10 11]]

numpy.load和numpy.save函数(推荐在不需要查看保存数据的情况下使用)
以NumPy专用的二进制类型保存数据,这两个函数会自动处理元素类型和shape等信息,
使用它们读写数组就方便多了,但是numpy.save输出的文件很难和其它语言编写的程序读入。
Note:
1. 文件要保存为.npy文件类型,否则会出错
2. 保存为numpy专用二进制格式后,就不能用notepad++打开(乱码)看了,这是相对tofile内建函数不好的一点

numpy.savez函数
如果你想将多个数组保存到一个文件中的话,可以使用numpy.savez函数。savez函数的第一个参数是文件名,其后的参数都是需要保存的数组,也可以使用关键字参数为数组起一个名字,非关键字参数传递的数组会自动起名为arr_0, arr_1, …。savez函数输出的是一个压缩文件(扩展名为npz),其中每个文件都是一个save函数保存的npy文件,文件名对应于数组名。load函数自动识别npz文件,并且返回一个类似于字典的对象,可以通过数组名作为关键字获取数组的内容:
如果你用解压软件打开result.npz文件的话,会发现其中有三个文件:arr_0.npy, arr_1.npy, sin_array.npy,其中分别保存着数组a, b, c的内容。

    a = array([[1,2,3],[4,5,6]])
    b = arange(0, 1.0, 0.1)
    c = sin(b)
    savez("result.npz", a, b, sin_array = c)
    r =load("result.npz")
    print r["arr_0"] # 数组a
    print r["arr_1"] # 数组b
    print r["sin_array"] # 数组c

结果:

None
[[1 2 3]
 [4 5 6]]
[ 0.   0.1  0.2  0.3  0.4  0.5  0.6  0.7  0.8  0.9]
[ 0.          0.09983342  0.19866933  0.29552021  0.38941834  0.47942554
  0.56464247  0.64421769  0.71735609  0.78332691]

numpy.savetxt和numpy.loadtxt(推荐需要查看保存数据时使用)
Note:savetxt缺省按照’%.18e’格式保存数据, 可以修改保存格式为‘%.8f’(小数点后保留8位的浮点数), ’%d’(整数)等等
总结:
载入txt文件:numpy.loadtxt()/numpy.savetxt()
智能导入文本/csv文件:numpy.genfromtxt()/numpy.recfromcsv()
高速,有效率但numpy特有的二进制格式:numpy.save()/numpy.load()

39、40部分来源:http://blog.csdn.net/pipisorry/article/details/39088003。详细内容请参看这个链接。

41、permutation函数
例子:

    #混淆位置。如果是多维数组,则混淆一维的。例如下面的arr.
    print random.permutation(10)
    print random.permutation([1, 4, 9, 12, 15])

    arr=arange(9).reshape((3,3))
    print arr
    print random.permutation(arr)

结果:

None
[2 9 6 4 0 3 1 7 8 5]
[12  1  4  9 15]
[[0 1 2]
 [3 4 5]
 [6 7 8]]

permutation混淆位置。如果是多维数组,则混淆一维的。例如下面的arr.

42、complex构造复数
例子:

print complex(0,100)#结果:100j

43、norm求范数
例子:

from scipy.linalg.misc import norm
print norm([8,6,10])

结果:

14.1421356237

默认情况下,求二范数。也就是各个元素的平方和再开平方。
当然,也可以求取其他范数。例如p范数、无穷范数等。
我这里import写的是scipy。其实numpy里面也有这个函数。详细的其他范数的求法请参考官网http://docs.scipy.org/doc/numpy-1.10.0/reference/generated/numpy.linalg.norm.html内容的介绍。对于求二范数的结果,我这里介绍的就够用了。

待续。。。。!!

本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系:hwhale#tublm.com(使用前将#替换为@)

python-numpy一些方法总结 的相关文章

随机推荐

  • linux 修改文件用户组和所有者

    目录 1 linux下修改文件用户组 2 linux下修改文件所有者 3 linux下同时修改文件所有者和用户组 1 linux下修改文件用户组 chgrp change group的简写 修改文件所属的用户组 chgrp 用户组名 文件名
  • (转) .net web项目的安装制作

    原 http blog csdn net houlinghouling archive 2005 06 17 396338 aspx 一 创建基本安装部署项目 1 在解决方案资源管理器 右击解决方案 添加 新建项目 安装部署项目 Web安装
  • 在loader中创建GDT,进入保护模式

    回顾 上一节实现了从BIOS中加载MBR MBR从磁盘2扇区读取loader加载到内存0x900处 但loader目前尚未实现任何功能 Q A Q1 loader在OS中主要做什么 答 创建一些系统数据结构 如GDT 页表等 打开进入保护模
  • 在csdn中复制的代码 去掉前面的行号

    在csdn中复制的代码会有行号 如下 1 2 3 4 5 6 解决方法 利用notepad 的替换功能 如下图一个个查找替换便可
  • stable diffusion实践操作-Controlnet

    本文专门开一节写提示词相关的内容 在看之前 可以同步关注 stable diffusion实践操作 文章目录 前言 1 ControlNet是什么 2 常用的模型 3 基本操作 openpose full 1 提示词 2 参数 控制效果参数
  • Thinkpad E580 硬件错误0187、2200、2201解决经历

    我的电脑是Thinkpad E580 最近电脑坏了 以下是具体情况 一天中午 我打开电脑 在屏幕显示完联想的logn之后 它出现了 我从未见过的我的电脑出现这样的情况 它也给我带来了生活上的不便以及精神和金钱上的损失 当然这是后话了 它长这
  • 人工智能(pytorch)搭建模型11-pytorch搭建DCGAN模型,一种生成对抗网络GAN的变体实际应用

    大家好 我是微学AI 今天给大家介绍一下人工智能 pytorch 搭建模型11 pytorch搭建DCGAN模型 一种生成对抗网络GAN的变体实际应用 本文将具体介绍DCGAN模型的原理 并使用PyTorch搭建一个简单的DCGAN模型 我
  • csharp: 用Enterprise Library对象实体绑定数据

    Enterprise Library https msdn microsoft com en us library ff648951 aspx
  • 盒马五年,参悟取舍

    配图来自Canva可画 从盒马举旗首垦算起 新零售已走过五个年头 在众多资本和玩家的集体探索下 整个赛道越来越繁荣 但格局和气象也在不断变化 这两年的变动尤为明显 疫情的打击 少数资本的强势以及头部梯队的强化 为赛道带来了更大的波动 一边是
  • Pandas 使用指南

    Pandas 可以说是在 Python 数据科学领域应用最为广泛的工具之一 Pandas 是 Python 的核心数据分析支持库 提供了快速 灵活 明确的数据结构 旨在简单 直观地处理关系型 标记型数据 Pandas 的目标是成为 Pyth
  • PCB表面处理方式详解

    上一篇文章中介绍了关于PCB设计中的常用基本概念 那么本篇文章就挑选其中的第三点 表面处理 给大家展开说一说 顺便也自我温习一下 什么是表面处理 由于铜长期与空气接触会使得铜氧化 所以我们需要在PCB表面做一些处理 这样才可以保证PCB板的
  • JQuery Dialog(JS模态窗口,可拖拽的DIV)

    代码分析地址 http jingyan baidu com album fb48e8be5347156e622e14ac html 或者http jingyan baidu com article fb48e8be5347156e622e1
  • 升级log4j2遇到的那些坑

    升级log4j2遇到的那些坑
  • 数据分析与挖掘之一:流程

    概念 数据分析 Data Analysis 是以数据为分析对象 以探索数据内的有用信息为主要途径 以解决业务需求为最终目标 包含业务理 解 数据采集 数据清洗 数据探索 数据可视化 数据建模 模型结果可视化 分析结果的业务应用等步骤在内的一
  • vue中读取对象中对象的属性无法读取报错

    问题描述 报错内容 Vue warn Error in render TypeError Cannot read properties of undefined reading cnName 在写博客系统的时候 在文章详情页 需要获取文章中
  • 杂谈随感-6:敢问路在何方?

    不识庐山真面目 只缘生在此山中 人员远虑 必有近忧 敢问路在何方 4G还是5G 是 钻地道 还是 云中飘 3G 已死 没有出路 4G 行将入木 死路一条 5G Hub不死不活 僵尸 RRU与BBU 一个钻地道 一个云中飘 渐行渐远 各走各路
  • 常用Dos命令

    1 dos命令 color all 修改背景字体颜色 cls 清屏 dir 查看当前目录有哪些文件 a 查看隐藏文件 a d 只查看目录不显示文档 r 只读文件 A 准备存档的文件 在内存中写了但是没网硬盘里面写 S 系统文件 c 显示文内
  • Caffe 工程的一些编译错误以及解决方案

    CAFFE深度学习交流群 532629018 整理一下最近遇到caffe工程的一些编译错误以及解决方法 1 cuDNN cuDNN当前最新版本是v5 近两三年的一些caffe工程 使用的版本不尽相同 其中以v2 v3版本的最为常见 所以使用
  • UCenter安装时提示mysql_connect()不支持

    问题描述 安装时 提示mysql connect 不支持 请检查 mysql 模块是否正确加载 如下图 分析原因 查看php官方帮助文档得知 mysql connect是php4有php5中的函数 在php5 5 0已标记为废弃 在php7
  • python-numpy一些方法总结

    1 multiply 例子 x1 1 2 3 x2 4 5 6 print multiply x1 x2 输出 4 10 18 multiply函数得到的结果是对应位置上面的元素进行相乘 2 std 标准方差 var 方差 例子 b 1 3