NumPy のバックアップ(No.5)

「.NET 開発基盤部会 Wiki」は、「Open棟梁Project」,「OSSコンソーシアム .NET開発基盤部会」によって運営されています。

戻る > Pythonセカンド・ステップ

目次 †

↑

概要 †

数値計算を効率的に行うための拡張モジュール

効率的な数値計算を行うための型付きの多次元配列のサポートをPythonに加える。
多次元配列を操作するための大規模な高水準の数学関数ライブラリを提供する。

↑

詳細 †

↑

準備 †

↑

インストール †

> python -m pip install numpy

↑

インポート †

>>>import numpy as np

↑

ファースト・ステップ †

↑

NumPyの配列 †

１次元配列はベクトルとも言う。

↑

生成 †

>>>x=np.array([1.0,2.0,3.0])
>>>print(x)
[1. 2. 3.]
>>>type(x)
<class 'numpy.ndarray'>

↑

確認 †

>>> np.ndim(x) # 配列の次元
1
>>> x.shape    # 次元の形状
(3,)
>>> x.shape[0] # 次元の要素数
4

↑

算術計算 †

同じ数の要素を持つ２つのNumPy配列を算術計算で処理する。

>>>x=np.array([1.0,2.0,3.0])
>>>y=np.array([2.0,4.0,6.0])

要素ごとの足し算
```
>>>x+y
array([ 3.,  6.,  9.])
```

要素ごとの引き算
```
>>>x-y
array([-1., -2., -3.])
```

要素ごとの掛け算
```
>>>x*y
array([  2.,   8.,  18.])
```

要素ごとの割り算
```
>>>x/y
array([ 0.5,  0.5,  0.5])
```

↑

NumPyのN次配列 †

↑

行列 †

行列とは？

２次元配列は行列とも言う。
行が、１つの添字で表される１つのべクトル
列は、ベクトルのｎ番目の要素をグループ化したもの。

    列↓
  ┌      ┐
行│○  ○│
→│      │
  │○  ○│
  └      ┘

↑

テンソル †

テンソルとは？

ベクトルや行列を一般化したものをテンソルと呼ぶ。
３次元以上の配列を、テンソル/多次元配列と呼ぶ。

↑

生成 †

>>> a=np.array([[1,2],[3,4],[5,6]])
>>> print(a)
[[1 2]
 [3 4]
 [5 6]]
>>> np.ndim(a) # 配列の次元
2
>>> a.shape    # 行列の形状
(3, 2)         # 行数・列数
>>> a.dtype    # 行列要素のデータ型
dtype('int32')

↑

算術計算 †

要素ごとの足し算

行列数同じ

>>>a=np.array([[1, 2], [3, 4]])
>>>b=np.array([[3, 0], [0, 6]])
>>>a+b
array([[ 4,  2],
       [ 3, 10]])

行列数が異なる

>>>a=np.array([[1, 2], [3, 4]])
>>>b=np.array([10,20])
>>>a+b
array([[11, 22],
       [13, 24]])

要素ごとの掛け算
行列の積でないことに注意

行列数同じ

>>>a=np.array([[1, 2], [3, 4]])
>>>b=np.array([[3, 0], [0, 6]])
>>>a*b
array([[ 3,  0],
       [ 0, 24]])

行列数が異なる

>>>a=np.array([[1, 2], [3, 4]])
>>>b=np.array([10,20])
>>>a*b
array([[10, 40],
       [30, 80]])

↑

要素へのアクセス †

↑

インデックス †

>>> a=np.array([[51, 55], [14, 19], [0, 4]])
>>> print(a)
[[51 55]
 [14 19]
 [ 0  4]]
>>>a[0]
array([51, 55])
>>>a[0][1]
55

↑

ベクトル化 & 抽出 †

ベクトル化

>>> a=np.array([[51, 55], [14, 19], [0, 4]])
>>> a=a.flatten()
>>> print(a)
[51 55 14 19  0  4]

抽出方法１

>>> a[np.array([0,2,4])]
array([51, 14,  0])

抽出方法２

>>>a>15 # 15より大きい値を調査。
array([ True,  True, False,  True, False, False], dtype=bool)
>>>a[a>15] # 15より大きい値を抽出。
array([51, 55, 19])

↑

ブロード・キャスト †

スカラ値や異なる形状の配列と計算を行う。

↑

スカラ値 †

配列とスカラ値との計算では、スカラ値が配列形状に拡張され計算される。

NumPyの配列

>>>x=np.array([1.0,2.0,3.0])
>>>x/2.0
array([ 0.5,  1. ,  1.5])

NumPyのN次配列

>>>a=np.array([[1, 2], [3, 4]])
>>>a*10
array([[10, 20],
       [30, 40]])

↑

異なる形状 †

→ 行列の積算

↑

セカンド・ステップ †

↑

行列の積算 †

Ａ行列 * Ｂ行列 =Ｃ行列

Ａ行列の列数とＢ行列の行数が一致ししている必要がある。
計算後のＣ行列は、Ａ行列の行数 * Ｂ行列の列数の行列になる。

↑

2行2列の場合 †

┌      ┐┌      ┐   ┌                    ┐
│a1  b1││a2  b2│   │a1a2+b1c2  a1b2+b1d2│
│      ││      │ = │                    │
│c1  d1││c2  d2│   │c1a2+d1c2  c1b2+d1d2│
└      ┘└      ┘   └                    ┘
  Ａ行列    Ｂ行列              Ｃ行列
  2行2列    2行2列              2行2列

↑

x行y列の場合 †

ベクトルと行列との計算では、
ベクトルが行列形状に拡張され計算される。

┌     ┐┌  ┐   ┌     ┐┌    ┐  ┌         ┐
│a1 b1││a2│   │a1 b1││a2 0│  │a1a2+b1b2│
│     ││  │ = │     ││    │ =│         │
│c1 d1││b2│   │c1 d1││b2 0│  │c1a2+d1b2│ 
└     ┘└  ┘   └     ┘└    ┘  └         ┘ 
  Ａ行列   Ｂ行列                       Ｃ行列
  2行2列   2行1列                       2行1列

Ａ行列の列数とＢ行列の行数が一致ししている必要がある。
計算後のＣ行列は、Ａ行列の行数 * Ｂ行列の列数の行列になる。

┌      ┐┌              ┐   ┌                                          ┐
│a1  b1││a2  b2  c2  d2│   │a1a2+b1e2  a1b2+b1f2  a1c2+b1g2  a1d2+b1h2│
│      ││              │   │                                          │
│c1  d1││e2  f2  g2  h2│ = │c1a2+d1e2  c1b2+d1f2  c1c2+d1g2  c1d2+d1h2│
│      │└              ┘   │                                          │
│e1  f1│                     │e1a2+f1e2  e1b2+f1f2  e1c2+f1g2  e1d2+f1h2│
└      ┘                     └                                          ┘
  Ａ行列        Ｂ行列                             Ｃ行列
  3行2列        2行4列                             3行4列

↑

メソッド †

NumPyのdot（ドット積）メソッドを使用する。

>>> a=np.array([[1,2,3],[4,5,6]])
>>> b=np.array([[1,2],[3,4],[5,6]])
>>> c=np.array([[1,2],[3,4]])
>>> d=np.array([7,8])
>>> a
array([[1, 2, 3],
       [4, 5, 6]])
>>> b
array([[1, 2],
       [3, 4],
       [5, 6]])
>>> c
array([[1, 2],
       [3, 4]])
>>> d
array([7, 8])
>>> np.dot(a, b)
array([[22, 28],
       [49, 64]])
>>> np.dot(b, c)
array([[ 7, 10],
       [15, 22],
       [23, 34]])
>>> np.dot(b, d)
array([23, 53, 83])
>>> np.dot(a, c) # Ａ行列の列数とＢ行列の行数が一致しない。
Traceback (most recent call last):
  File "<stdin>", line 1, in <module>
ValueError: shapes (2,3) and (2,2) not aligned: 3 (dim 1) != 2 (dim 0)

↑

参考 †

NumPy - Wikipedia
https://ja.wikipedia.org/wiki/NumPy

目次 †

概要 †

詳細 †

準備 †

インストール †

インポート †

関数 †

数列の生成 †

統計関数 †

行列演算 †

信号処理 †

ファースト・ステップ †

NumPyの配列 †

生成 †

確認 †

算術計算 †

NumPyのN次配列 †

行列 †

テンソル †

生成 †

算術計算 †

要素へのアクセス †

インデックス †

ベクトル化 & 抽出 †

ブロード・キャスト †

スカラ値 †

異なる形状 †

セカンド・ステップ †

行列の積算 †

2行2列の場合 †

x行y列の場合 †

メソッド †

参考 †