博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
Numpy学习50例
阅读量:5059 次
发布时间:2019-06-12

本文共 3540 字,大约阅读时间需要 11 分钟。

基础部分

1、导入numpy模块

import numpy as np

2、查看numpy版本信息

print(np.__version__)

numpy的主要对象的多维数组Ndarray。Numpy中维度(dimensions)叫做轴(axis),轴的个数叫做秩。

3、通过列表创建一位数组

np.array([1, 2, 3])

4、通过列表创建一个二维数组

np.array([(1, 2, 3),(4, 5, 6)])

5、创建全为0的二维数组

np.zeros((3,3))

6、创建全为1的三维数组

np.ones((2,3,4))

7、创建一维等差数组

np.arange(5)

8、创建二维等差数组

np.arange(6).reshape(2,3)

9、创建二维单位矩阵

np.eye(3)

10、创建等间隔一维数组(1到10之间,共6个数,数字之间间隔相同)

np.linspace(1, 10, num=6)

11、创建二维随机数组

np.random.rand(2,3)

下面说明一下random模块中的一些常用函数的用法

  • numpy.random.rand(d0, d1, ..., dn):生成一个[0,1)之间的随机浮点数或N维浮点数组。
  • numpy.random.randn(d0, d1, ..., dn):生成一个浮点数或N维浮点数组,取数范围:正态分布的随机样本数。
  • numpy.random.randint(low, high=None, size=None, dtype='l'):生成一个整数或N维整数数组,取数范围:若high不为None时,取[low,high)之间随机整数,否则取值[0,low)之间随机整数。
  • numpy.random.seed():使得随机数据可预测。当我们设置相同的seed,每次生成的随机数相同。如果不设置seed,则每次会生成不同的随机数

12、依据自定义函数创建数组

np.fromfunction(lambda i, j: i +j, (5, 6))fromfunction()方法可以根据矩阵的行号列号生成矩阵的元素。 例如创建一个矩阵,矩阵中的每个元素都为行号和列号的和。import numpy as np

def func(i,j):

return i+j
a = np.fromfunction(func,(5,6))
第一个参数为指定函数,第二个参数为列表list或元组tuple,说明矩阵的大小
print(a)
输出:[[ 0. 1. 2. 3. 4. 5.]
[ 1. 2. 3. 4. 5. 6.]
[ 2. 3. 4. 5. 6. 7.]
[ 3. 4. 5. 6. 7. 8.]
[ 4. 5. 6. 7. 8. 9.]]
注意这里行号的列号都是从0开始的,func传进去的实参是矩阵的行列号

13、矩阵乘法运算

np.dot(A, B)

14、使用np.mat将数组定义为矩阵,直接使用进行矩阵乘法运算
np.mat(A)
np.mat(B)

15、矩阵的转置

A.T

16、矩阵求逆

np.linalg.inv(A)

17、常用数学函数计算

np.sin(a)np.exp(a)np.sqrt(a)  #开平方运算np.power(a, 3)  #立方根运算

18、数组的切片和索引

a = np.array([1, 2, 3, 4, 5])a[0], a[-1]a[0:2], a[:-1]

19、更改数组的形状(不改变原始数组)

a.reshape(2,3)

20、更改数组形状(改变原始数组)

a.resize(2, 3)

21、展平数组(将多维变成一维)

a.ravel()       a.flatten()numpy.flatten()返回一份拷贝,对拷贝所做的修改不会影响(reflects)原始矩;而numpy.ravel()返回的是视图(view),会影响(reflects)原始矩阵。

22、垂直拼合数组

np.vstack((a, b))

23、水平拼合数组

np.hstack((a, b))

24、沿轴分隔数组

np.hsplit(a, 3)  #沿横轴切分数组。a是ndarray数组,3表示切分为几份np.vsplit(a,2)  #沿纵轴切分数组。

25、返回数组最值

np.min(a, axis=num)np.max(a, axis=num)np.argmax(a, axis=num)  #返回最大值的索引np.argmin(a, axis=num)

26、数组的统计特性计算

np.median(a, axis=num)  #指定维度的中位数np.mean(a, axis= num)  #算数平均值np.average(a, axis= num)  #加权平均值np.var(a, axis = num)  #方差np.std(a, sxis = num)  #标准差

进阶部分

27、创建一个5*5的二维数组,其中边界值为1,其余值为0
arr = np.ones((5,5)) #方法1
arr[1:-1,1:-1] = 0
z

arr = np.ones((5,5))  #方法二arr = np.pad(arr, pad_width=1, mode='constant', constant_values=0)

28、创建5*5的二维数组,并设置1,2 ,3 ,4落在其对角线下方
arr = np.diag(1+np.arange(4), k = -1)

29、创建一个0-10的一维数组。并将大于等于5的数置为0

arr1 = np.random.randint(0,10,10)arr1[arr1 >= 5] = 0

29、找出两个一维数组中子相同的元素

arr1 = np.random.randint(0,10,10) arr2 = np.random.randint(0,10,10)np.intersect1d(arr1,arr2)

30、打印今天、昨天、明天

yesterday = np.datetime64('today', 'D') - np.timedelta64(1, 'D')today     = np.datetime64('today', 'D')tomorrow  = np.datetime64('today', 'D') + np.timedelta64(1, 'D')print("yesterday: ", yesterday)print("today: ", today)print("tomorrow: ", tomorrow)

31、将数组排序

Z = np.random.random(10)Z.sort()

32、将float32转化为整型

arr = np.arange(10, dtype = np.float32)arr = arr.astype(np.int32, copy = False)

33、将随机二维数组按照第3列从上到下升序排列

arr = np.random.randint(0, 10, (5,5))arr[arr[:,2].argsort()]

34、将二维数组的前两行进行交换

A = np.arange(25).reshape(5,5)A[[0,1]] = A[[1,0]]

35、找出随机一维数组中出现频率最高的值

arr = np.random.randint(0,10,50)np.bincount(arr).argmax()

36、对于给定的55二维数组,在其内部随机放置num个值为1的数
arr = np.zeros((5,5))
np.put(arr, np.random.choice(range(5
5), p, replace=False), 1)

37、获得二维数组点积结果的对角线数组

A = np.random.uniform(0,1,(3,3))A = np.random.uniform(0,1,(3,3))np.sum(A*B.T,axis = 1)  #较快的方法np.diag(np.dot(A, B))  #较慢的方法

38、设置数组中的各元素的值,使其保留两位小数

np.set_printoptions(precision=2)

39、计算两个数组之间的欧式距离

np.linalg.norm(b-a)

转载于:https://www.cnblogs.com/wemo/p/8900565.html

你可能感兴趣的文章
Linux 根文件系统制作
查看>>
IOS--沙盒机制
查看>>
My.Ioc 的性能
查看>>
使用 JointCode.Shuttle 访问任意 AppDomain 的服务
查看>>
hdoj 1846 Brave Game(巴什博弈)
查看>>
Round #345 B. Beautiful Paintings(Div.2)
查看>>
51nod 1018排序
查看>>
sqlite的坑
查看>>
Elastic Search 上市了,市值翻倍,这群人财务自由了!
查看>>
digitalocean --- How To Install Apache Tomcat 8 on Ubuntu 16.04
查看>>
linux swoole
查看>>
An Easy Problem?! - POJ 2826(求面积)
查看>>
【题解】[P4178 Tree]
查看>>
Jquery ui widget开发
查看>>
css3实现循环执行动画,且动画每次都有延迟
查看>>
更改git仓库地址
查看>>
有标号DAG计数 [容斥原理 子集反演 组合数学 fft]
查看>>
Recipe 1.4. Reversing a String by Words or Characters
查看>>
Rule 1: Make Fewer HTTP Requests(Chapter 1 of High performance Web Sites)
查看>>
sql注入
查看>>