numpy 中的一些函数用法

  • 对于一个一维数组Z,其以下表达式的结果如下:
import numpy as np
Z = np.arange(10)
Z = Z**Z

用以上代码生成的Z数组为[0,1,2,3,4,5,6,7,8,9],而Z**Z的内部计算方式为 [0^0,1^1,2^2,……9^9].

  • 对于一个一维数组,以下表达式的计算结果如下:
Z =  np.arange(10)
print(2<<Z>>2)

>>>[ 0 1 2 4 8 16 32 64 128 256]

其计算过程可以简单描述为,[2^(0+1)/4,2^(1+1)/4 …… 2^(9+1)/4]
即 2<<Z 表示为[2^(0+1),2^(1+1) …… 2^(9+1)],而Z>>2 就是简单的每个元素除以4了。

  • 复数域和实数域中sqrt(-1)的差别
    有如下表达式,其计算结果分别如下所示:
print(np.sqrt(-1))  #输出结果为Nan
print(np.emath.sqrt(-1)) #输出结果为1j

>>> nan
>>>1j

  • np.add.reduce()的详细计算过程可以由以下代码实现,reduce的作用,是把add函数作用于矩阵a的每一维的对应位置上的元素相加,计算这些元素。
import numpy as np
a = np.arange(45).reshape(3,3,5)
print(a)
b = np.add.reduce(a)
print(b,'\n')
c = []
for j in range(3):
    d = []
    for i in range(5):
        print(a[:,j,i],'\t',sum(a[:,j,i]))
        d.append(sum(a[:,j,i]))
    c.append(d)
print(np.array(c))

详细的计算过程,如下输出结果所示:

>>>[[[ 0 1 2 3 4]
[ 5 6 7 8 9]
[10 11 12 13 14]],
[[15 16 17 18 19]
[20 21 22 23 24]
[25 26 27 28 29]],
[[30 31 32 33 34]
[35 36 37 38 39]
[40 41 42 43 44]]],
计算结果:
[[45 48 51 54 57]
[60 63 66 69 72]
[75 78 81 84 87]]
计算过程:
[ 0 15 30] 45
[ 1 16 31] 48
[ 2 17 32] 51
[ 3 18 33] 54
[ 4 19 34] 57
[ 5 20 35] 60
[ 6 21 36] 63
[ 7 22 37] 66
[ 8 23 38] 69
[ 9 24 39] 72
[10 25 40] 75
[11 26 41] 78
[12 27 42] 81
[13 28 43] 84
[14 29 44] 87
最后输出:
[[45 48 51 54 57]
[60 63 66 69 72]
[75 78 81 84 87]]

  • 关于求取最大值的说明:
Z = np.random.random(10)
print(Z.max())  #向量的最大值
print(Z.argmax())  #向量最大值的索引
  • 关于np.bincount(X,W,minlength)函数,有如下用法:
  1. 不指定权重参数W,np.bincount(X)输出是一个其索引值出现的次数,如下代码所示:
X = [0,1,3,2,5,6,3,3]
print(np.bincount(X))

运行如上代码,输出结果为:[1,1,1,3,0,1,1],需要注意的是,输出数组的长度比数组中的最大值+1

  1. minlength参数被指定之后,输出数组的长度为minlength,如下代码所示:
X = [0,2,4,4,2,4,5,6,6]
print(np.bincount(X,minlength = 10))

运行以上代码,输出结果为:[1,0,1,0,3,1,2,0,0,0]

  1. 当weigths参数被指定时,运行结果可以用以下代码详解:
# 关于np.bincount指定权重向量的用法
import numpy as np
X = [1,2,3,4,5,6]
I = [1,3,9,3,4,1]
#最大值是9,所以输出F包含10个元素
#F[0] = 0
#F[1] = X[I[i=1]] = X[0] + X[6] = 1 + 6 = 7
#f[2] = 0 ,I中没有出现2
#F[3] = X[I[i = 3]] = X[1] + X[4]  = 6
#F[4] = X[I[i = 4]] = X[5] = 5
#F[5] = 0 I中没有5
#F[6] = 0 I中没有出现6
#F[7] = 0 I中没有出现7
#F[8] = 0 I中没有出现8
#F[9] = X[I[i=9]] = X[4] = 3  
F = np.bincount(I,X)
print(F)

>>>[0. 7. 0. 6. 5. 0. 0. 0. 0. 3.]

©著作权归作者所有,转载或内容合作请联系作者
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 229,732评论 6 539
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 99,214评论 3 426
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 177,781评论 0 382
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 63,588评论 1 316
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 72,315评论 6 410
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 55,699评论 1 327
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 43,698评论 3 446
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 42,882评论 0 289
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 49,441评论 1 335
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 41,189评论 3 356
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 43,388评论 1 372
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 38,933评论 5 363
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 44,613评论 3 348
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 35,023评论 0 28
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 36,310评论 1 293
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 52,112评论 3 398
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 48,334评论 2 377

推荐阅读更多精彩内容