数学建模社区-数学中国

标题: Python数据挖掘基础(二):Numpy [打印本页]

作者: 杨利霞    时间: 2020-5-8 17:53
标题: Python数据挖掘基础(二):Numpy
7 s8 Y9 _& K8 P5 @% @5 X! ]

$ r6 ^) G% X, w! Q2 _1 a, m* hPython数据挖掘基础(二):Numpy目录9 e' E# P, w. R- ~! \
$ S) `  r2 L/ G7 C. \; N
1. Numpy的优势: t# n; ^' N. n$ e6 T
2. 数组属性
, K* c& R/ W4 f3. 创建数组+ q8 l9 j9 J5 B( _- i6 g
4. 数组形状与类型变化
3 A) i/ ^4 d/ O5 e1 x. Y4. 数组运算
8 _8 {3 x# z- R& `4.1 逻辑运算& ?  P$ v& \3 b4 V+ X2 b& U' F' ]
4.2 统计运算3 u" h( K0 O1 ^; C
4.3 数组间运算( u& o2 Q, n2 q8 l5 H
4.4 合并分割
- E6 Y- j1 ^: L1 T7 g& ?2 U1. Numpy的优势
' O$ Y$ {9 |, O  W- G/ m4 `- d) r) t) D2 C" }" K7 L
Python已经提供了很多丰富的内置包,我们为什么还要学习NumPy呢?先看一个例子,找寻学习 NumPy 的必要性和重要性。如下:: _* M' o) \: V5 B/ k
, w) g# s+ ~- W3 y% g. u
完成同样的都对元素相加的操作,NumPy比Python快了11倍之多。这就是我们要学好NumPy的一个重要理由,它在处理更大数据量时,处理效率明显快于Python。并且内置的向量化运算和广播机制,使得使用NumPy更加简洁,会少写很多嵌套的for循环,因此代码的可读性大大增强。NumPy计算为什么这么快呢?原因如下:, w) b: q# O! K2 G
2 Q- Y* _4 t# g. z
Python 的 list是一个通用结构。它能包括任意类型的对象,并且是动态类型。
# }2 O3 j' H6 r4 s6 DNumPy 的 ndarray 是 静态、同质 的类型,当ndarray对象被创建时,元素的类型就确定。由于是静态类型,所以ndarray间的加、减、乘、除用 C 和 Fortran 实现才成为可能,所以运行起来就会更快。根据官方介绍,底层代码用 C语言 和 Fortran 语言实现,实现性能无限接近 C 的处理效率。
0 u) X. j  c- B) g" G 2.png
! B* O2 n7 z0 P, V& d从图中我们看出来 NumPy 其实在存储数据的时候,数据与数据的地址都是连续的,这样就给我们操作带来了好处,处理速度快。* {7 N$ {! Q6 ?
支持并行化运算,也叫向量化运算。当然向量是数学当中的概念,我们不过多解释,只需要知道他的优势即可。也就是说 NumPy 底层使用 BLAS 做向量,矩阵运算。2 w# x: d2 S' |6 r+ x4 p0 o5 V! `
numpy的许多函数不仅是用C实现了,还使用了BLAS(一般Windows下link到MKL的,下link到OpenBLAS)  k$ _9 }+ n- {  B" F. |
基本上那些BLAS实现在每种操作上都进行了高度优化
2 `! _! k; P2 Z! V例如使用AVX向量指令集,甚至能比你自己用C实现快上许多,更不要说和用Python实现的比
0 P$ V! ]7 ?% F9 T, a1: w/ s' n6 a; g" V: ^3 e: ]4 l9 h
2( a4 r$ p6 M& _" b  v
3% y8 O1 H( }+ T7 s
由此可见,NumPy 就非常适合做大规模的数值计算和数据分析。( {2 Q; ?' }7 w! W# e4 m, ~8 T& m

# e, u5 I, X( ~' @9 b- w2. 数组属性
7 P5 |) K! ~$ J' @7 L% Z; B: e8 A1 v- b8 f4 y! B- O5 Z5 Q
数组属性反映了数组本身固有的信息。0 E8 O6 p7 U  k# a9 h

  w7 J+ }( l. F" `; I0 F属性名字        属性解释
3 B3 n' ]5 |1 d0 r  _' Sndarray.shape        数组维度的元组  \4 }6 U& W0 k; J
ndarray.flags        有关阵列内存布局的信息
7 \# \9 i  U4 t8 ?/ mndarray.ndim        数组维数! k( X9 S( n2 }3 ?9 ^
ndarray.size        数组中的元素数量
3 y9 r8 v6 w# ^% x% O0 H/ Wndarray.itemsize        一个数组元素的长度(字节)0 O* F* z  P4 N) v/ v2 r
ndarray.nbytes        数组元素消耗的总字节数
* N+ g. v5 ^% Q' f示例代码如下:, h# {* J  t6 n* M
import numpy as np
) P) a& m0 X1 u* X) M& t; l2 I8 C8 R) Y+ \
# 数组的属性0 |1 Y; v5 t& n% P
# 1.创建数组 这里先不用管 后续会详细讲解数组的创建方法
7 N( C5 y& q4 C3 |3 E) U( t% \- T( Qa = np.array([[1,2,3],[4,5,6]])) k% H. L9 L9 Z1 f/ J" h
b = np.array([1,2,3,4]). r: e& j7 P- S
c = np.array([[[1,2,3],[4,5,6]],[[1,2,3],[4,5,6]]])
5 \- W# D- k$ [# 创建数组的时候指定类型) t  G# T- K/ }. F! M
# dtype更多取值: int complex bool object8 c3 k1 M( q7 j* W( s/ Y7 ^3 B5 s
# 还可以显示的定义数据位数的类型,如: int64、int16、float128、complex128。
9 n$ @9 }' B7 k! i1 N, U- Od = np.array([1,2,3,4], dtype=np.float)" H! P( L# O9 ~' F
0 @. z, ]* k% x% F' r( k% b/ N
# 2.测试数组属性
5 Y# w5 O7 @. z2 x0 O! D7 yprint(a.shape)  # 数组形状  (2, 3): 二维数组7 H( ?% v9 s8 }- E9 r
print(b.shape)  # (4,) 一维数组: 有4个元素; e2 [" Y0 ~2 J6 y0 ?% [8 U$ g  ?
print(c.shape)  # (2, 2, 3): 三维数组, }" R  q- p( D( u7 r( E+ Y% R8 k
print(a.ndim)  # 数组维数 2
: e; P: b4 @! {) I3 S% Qprint(a.size)  # 元素的数量 6
8 S$ p$ s; Y. m) jprint(a.itemsize)  # 每一个元素占的位数(字节) 8- F" N! H7 U6 ?' n% B% T
print(a.nbytes)  # 总共占的字节数 6*8 ==> 48
. B2 `- I0 g2 q) o7 ]$ A; _# H" Aprint(a.flags)  # 阵列内存布局
* f7 ?1 g1 a$ N5 N& U* y$ f! Jprint(a.dtype)  # 数组类型 int643 ~5 D7 y' C1 ]& N" I1 c: w8 r
print(d.dtype)  # float642 q, C2 h* I" n

3 V9 S6 B9 L9 C, u) D: V+ v  f
6 F7 q8 T1 P5 A6 i  ~. M2 dimport numpy as np
/ q! L1 A* X  j6 j, e0 `
# E1 `) \) G, j; u- M# u9 b( m$ ?上述代码执行结果如下:  T, q& y" c1 Z+ T
3.png % W5 C/ P+ R# L

/ D; y; Q3 h9 k! Y5 M3. 创建数组) \0 ]5 ?2 z( l( |7 H: B
3 S  w. l, e+ F! y; Z
1. 创建0和1的数组 示例代码如下:+ B, }5 p- Z% F8 P3 k( V
4.png
! a( ^$ c9 V( l, F1 n* S2. 从现有的数据中创建 示例代码如下:
( ^+ h/ @( f, Y; |5 _2 a! L6 n. X
7 S  Q* k# m2 V4 ^5 rnum_list = [[1,2,3], [4,5,6]]
9 U: \+ \6 Z3 [0 D4 E* [a = np.array(num_list)
3 p) W5 K( N: @( Sa1 = np.array(a)  # 创建了一个新的数组' I0 B) o! M3 }3 F( W& ]# U
a2 = np.asarray(a)  # 还是引用原来的数组, _5 O! A/ \" x, `  I
print(a)5 b/ {$ o7 I% S2 A2 |; Y( ^. ^
print(a1)
+ n: J" b% R; E8 e6 iprint(a2)
' D% a/ @; y/ K% D0 E' `a[0] = 10
5 n$ w1 c: c, pa, a1, a2
4 w0 N# M! K* P. S! z1 ^& Z
) _  M4 v! W  l5 D3 @上述代码执行结果如下:- W2 D6 K8 i* U" c. r
5.png
( Z: |" m) t; K, E2 ?+ L3. 创建固定范围的数组,语法格式如下:
1 C( \+ K; E  R4 L1 \- U' [$ Y0 g" v- C
np.linspace(start, stop, num, endpoint, retstep, dtype) 生成等间隔的序列- Y3 }( T: L3 k8 E' k
start:序列的起始值; {& V" |" v, Z% R% k
stop:序列的终止值
5 {) W# H. P6 Jnum:要生成的等间隔样例数量,默认为50( H7 B- v8 J7 ]* T2 @
endpoint:序列中是否包含stop值,默认为True2 V  @# I$ B: E  ]
retstep:如果为True,返回样例,以及连续数字之间的步长
% _! I& b! H& L; tdtype:输出ndarray的数据类型6 N' y9 H$ i% _; \0 {4 x
17 C  S/ t9 r0 \* n; D
23 s+ s( c' e" d; j- p+ X! F/ i
3% J, D" C+ v8 L8 \6 J% k, S
43 s! C8 J7 [  h0 [  m5 w
57 X6 o( J! G2 W5 W
6
9 H" j) _+ z  R1 H" b- O( r7
. d& O* F) \$ s& t6 E& l4 C3 \& w示例代码如下:0 R1 [& ~/ e9 g9 J+ k' u
' v& E* L5 N" x8 G
arr = np.linspace(0, 10, 10)
8 t+ R, q, V: q$ d7 C3 iarr/ Q* Z5 I" ?* F% _6 _. M
16 A, W4 }! O  x7 D# H) N+ k& |6 W
2
1 V0 U* c  X7 J# E执行结果如图所示:6 }, _5 t. Q# X

- y5 s& P/ n+ X5 {# O2 f7 w其它的还有:
, C7 d' X9 O+ Q4 p* P  q: Z( f1 r. K! {$ ]
numpy.arange(start,stop, step, dtype) 示例代码如下:* {$ A7 W- X( W$ h
np.arange(1, 10, 2)
1 N* a& ~# @; V; l$ R& N5 V( q1
9 z% O& C  f. d6 o& e; D3 v" U运行结果如图所示:
: y% [. |* f# C7 @4 c- @ 6.png
% O7 Q4 O* a$ Y. o9 Rnumpy.logspace(start,stop, num, endpoint, base, dtype) 构造一个从10的-2次方 到 10的2次方 的等比数列,这个等比数列的长度是 10 个元素,示例代码如下:$ ~% D  v& j/ A# n
np.logspace(-2,2,10)
# H/ C+ ]# j, B1
4 Q- p8 U( Q$ [4 o/ z6 |运行结果如图所示:2 [5 C- b% [  R) n6 B& Q
7.png
$ L; J1 Q, A: o' Y3 U如果不想是10的次方,也就是想改变基数,那么可以这么写,代码如下:
, N- i4 z- K% S. e. ?# T+ pnp.logspace(-2,2,10,base=2)9 T0 {- L  O2 \8 p" v* S" ~
1
2 f3 J  k4 b3 Q$ J运行结果如图所示:
) D0 w. G& a3 c# p  i 8.png 7 e& n1 ]2 e! J
4. 创建随机数组5 [$ M: J0 \' r8 g( j* e# j

0 |8 O4 S* g; y0 J4 B4 I% W  snp.random 模块生成随机数组,更加方便,示例代码如下:
* x% C% W* L+ H( R' k* {
  T) Z% J) l9 K1 ^上面产生的数据是属于一个均匀分布。那么什么是均匀分布呢?在概率论和统计学中,均匀分布也叫矩形分布,它是对称概率分布,在相同长度间隔的分布概率是等可能的。 均匀分布由两个参数a和b定义,它们是数轴上的最小值和最大值,通常缩写为U(a,b)。! W- H8 }" Q/ P% c* X% V
% x4 ~, t( {; W0 j$ F; A7 d9 K
正态分布?给定均值/标准差/维度的正态分布,示例代码如下:
) D9 ?: _/ c! D* [; K9 r% f9 }+ [% o! ~) T  z, `0 _6 I
正态分布是一种概率分布。正态分布是具有两个参数μ和σ的连续型随机变量的分布,第一参数μ是服从正态分布的随机变量的均值,第二个参数σ是此随机变量的方差,所以正态分布记作N(μ,σ )。$ u/ L6 [2 f2 p; A* a

: c- A: }* n1 ^1 \$ Q* H3 c) O生活、生产与科学实验中很多随机变量的概率分布都可以近似地用正态分布来描述。μ决定了其位置,其标准差σ。决定了分布的幅度。当μ = 0,σ = 1时的正态分布是标准正态分布。
! c' F3 s$ L. N; |5 j8 H3 N6 }
( p6 o  J. l: i* h3 M$ r5 ?标准差如何来?方差是在概率论和统计方差衡量一组数据时离散程度的度量。
' r6 R( {6 d" L
3 b9 |9 |  g: e4 f& N其中M为平均值,n为数据总个数,s为标准差,s^2可以理解一个整体为方差。. `* G; L+ v/ e, n: h

6 B" ~* \: j! W' b/ R2 N9 Z  ]4 ?通过索引切片等获取数组中的值,一维数组示例代码如下:
1 s: Q$ C5 R" n. j( T( l& f! F/ V( H
二维数组示例代码如下:! p# ]: b2 J6 X6 l4 ^) Q
( J; f6 X/ O+ e) X3 E( ]
三维数组示例代码如下:
  ?3 O2 [+ `$ o5 @" g3 {9 w
7 u, K, h9 a( ]1 [5 U5 A0 @1 I# L1 l( j0 ~! o5 S8 K
4. 数组形状与类型变化
5 W  \% w* W  h
- h$ O+ j3 w. |6 }/ x& D1. ndarray.reshape(shape[, order]) Returns an array containing the same data with a new shape. 示例代码如下:
% U8 S# Q, f& ^) E( F) Y
  v) w, h: u9 a' ?/ l+ ?2. ndarray.resize(new_shape[, refcheck]) Change shape and size of array in-place. 示例代码如下:/ @0 _5 F- F: G+ T6 G
. e0 B- @( l; u
3. 修改类型 ndarray.astype(type) 示例代码如下:
& d. A' Q3 R  N9 |8 `' p' b& T! W( G/ z4 J
4. 修改小数位数 ndarray.round(arr, out) Return a with each element rounded to the given number of decimals. 示例代码如下:/ P+ ~; q6 D% k, n8 ?6 H& |, `7 @. I6 l

3 Z5 K3 q. F% f8 b0 ]5. ndarray.flatten([order]) Return a copy of the array collapsed into one dimension. 示例代码如下:6 M* E% M( H* k6 W; M5 u/ L( N
+ s; g/ s% \$ \7 }# F" ]/ B. A0 D3 `
6. ndarray.T 数组的转置 将数组的行、列进行互换 示例代码如下:
! q+ w8 y' `7 F  c9 I. C. }: o1 K( L; u0 C
7. ndarray.tostring([order])或者ndarray.tobytes([order]) Construct Python bytes containing the raw data bytes in the array. 转换成bytes3 T; y$ Z8 Y5 n9 Y' `' ^
6 H7 e9 f6 R# O
8. ndarray.copy([order]) Return a copy of the array. 当我们不想修改某个数据的时候,就可以去进行拷贝操作。在拷贝的数据上进行操作,示例代码如下:
5 _1 M- f5 \) c3 e* p( f& F1 ?4 j1 q7 _% s8 M- s: b. j
9 d, \' @( j! K# ?: d" ]
4. 数组运算
% [3 U4 d' y1 U* v/ G' c, n6 s
6 b. ~% a6 I7 {2 `4.1 逻辑运算3 F9 R8 s9 n" f0 I+ d) @0 z( A+ h- Q

3 x% j& k% c5 Y; ?
3 Z5 }" d4 G, z通用判断函数,np.all(),示例代码如下:7 E6 Y2 r8 \" ^  S( n) p

+ {8 U' v6 q8 q* E9 Znp.unique():返回新的数组的数值,不存在重复的值,示例代码如下:
2 @$ N! b1 ]7 @, W( T( P  A& I. W; I/ C; M% o" U; p0 [; w! k) P
np.where (三元运算符):通过使用np.where能够进行更加复杂的运算,示例代码如下:
% N# T0 ]6 w2 Z6 u( I$ a* k2 @9 O. v( ?& N" b

5 ], r9 [# @2 k: q4.2 统计运算& @5 V$ }9 E' H: ?* J! N
9 W7 c, ~% u1 ]! w: e0 ?/ E
在数据挖掘/机器学习领域,统计指标的值也是我们分析问题的一种方式。注意:进行统计的时候,axis轴 的取值并不一定, NumPy中不同的API轴的值都不一样,在这里,axis 0代表列, axis 1代表行 去进行统计。常用的指标如下:
- ]" y+ c$ u2 m, ?% E% ]; _
9 B: n4 z! N7 K/ P5 u" Y+ Y) {; Vmin(a[, axis, out, keepdims]) Return the minimum of an array or minimum along an axis. 示例代码如下:- c5 Y4 H6 K1 i0 p! K# B

) N; p  B. p4 ~4 A( @max(a[, axis, out, keepdims]) Return the maximum of an array or maximum along an axis. 示例代码如下:
4 E4 e9 r8 N6 A6 |2 b' [% g5 a" a: y* q
median(a[, axis, out, overwrite_input, keepdims]) Compute the median along the specified axis. 示例代码如下:) ~+ O% ]6 \. |2 a2 D

. [1 J/ t( d8 `$ l# j2 h' H; L1 emean(a[, axis, dtype, out, keepdims]) Compute the arithmetic mean along the specified axis. 示例代码如下:7 \9 V  A/ a1 Y6 l$ B7 A8 o" {4 ^

- O6 Y8 z) f0 g  E, }- ostd(a[, axis, dtype, out, ddof, keepdims]) Compute the standard deviation along the specified axis. 示例代码如下:
- J) [- n. R& K1 H% q
- r+ X" `3 y3 x1 U0 bvar(a[, axis, dtype, out, ddof, keepdims]) Compute the variance along the specified axis. 示例代码如下:  h3 k! @# h5 n2 D% B8 p3 l4 D: a

0 ?/ U; Q- [2 w2 U9 z/ W: }, ?4 Knp.argmax(temp, axis=) 示例代码如下:
; O6 J# i; z! a, k5 H* f
6 l5 ], F# ~0 [6 S. a! _np.argmin(temp, axis=) 示例代码如下:
) D: D6 N) v  H1 }
* S: d* X0 o& _. _4.3 数组间运算
) u0 f( O2 I# H, U
3 I3 Y" O! A) J& H数组与数的运算,示例代码如下:
1 K6 T3 o5 Z+ A. J( J2 s/ {& ]$ {, Y. E/ K9 P
矩阵运算,什么是矩阵?矩阵,英文matrix,和array的区别矩阵必须是2维的,但是array可以是多维的。示例代码如下:
* Y, n' {$ m3 w! O) s) b8 Z' V
9 l, r& ^; T5 ~1 A# Q* M1 s& A1 I
- n0 Q2 @4 F9 I" i4.4 合并分割
( \) g0 B. E- |$ X1 \6 w
9 D; I( |1 X8 P( Bnumpy.concatenate((a1, a2, …), axis=0) 示例代码如下:
3 e* x+ c* G. ]3 g' E5 F6 ~- q. A6 k# @  b4 }/ A
numpy.hstack(tup) Stack arrays in sequence horizontally (column wise). 示例代码如下:1 p4 ?3 r1 u$ |* }1 h8 D+ i
5 q- \8 F8 b6 ~: p6 f+ G2 {
numpy.vstack(tup) Stack arrays in sequence vertically (row wise). 示例代码如下:
' h& P; S3 b* a6 @# [# m* Q$ Q
+ ^' U8 x# g/ m2 |7 v: j1 fnumpy.split(ary, indices_or_sections, axis=0) Split an array into multiple sub-arrays. 示例代码如下:/ d! }) L* T+ U4 ^+ w/ R3 t

' W8 r4 D$ [+ a$ H————————————————/ j7 X& x, P# W. B9 O' U1 Y) X& `
版权声明:本文为CSDN博主「Amo Xiang」的原创文章,遵循CC 4.0 BY-SA版权协议,转载请附上原文出处链接及本声明。
, G. _: s0 n7 E6 a0 G/ q% n原文链接:https://blog.csdn.net/xw1680/article/details/1059313131 r: D8 I5 L- v

9 |  G+ i8 }7 r: _7 Z" I: \4 B* h

1.png (142.46 KB, 下载次数: 450)

1.png

2.png (164.52 KB, 下载次数: 430)

2.png


作者: 建模小白12138    时间: 2020-5-9 14:30
谢楼主!!!!!!!!!!!!
* u) P. f* k. x8 j) R. j, `1 J& m4 h




欢迎光临 数学建模社区-数学中国 (http://www.madio.net/) Powered by Discuz! X2.5