ndarray – 配列による要素指定

2019-10-27 / tau / コメントする

“Pythonではじめる機械学習”中、27番のコードで気になったので調べた結果。

Scitkit-learnを使ってk-最近傍法の予測結果を保存し、その結果を使って予測されたアヤメの種類を表示させている。

prediction = knn.predict(X_new)
print("Prediction: {}".format(prediction))
print("Predicted target name] {}".format( \
    iris_dataset['target_names'][prediction]))

# Prediction: [0]
# Predicted target name] ['setosa']

prediction = knn.predict(X_new)

print("Prediction: {}".format(prediction))

print("Predicted target name] {}".format( \

iris_dataset['target_names'][prediction]))

# Prediction: [0]

# Predicted target name] ['setosa']

ここでprediction、iris_dataset['target_names']のいずれの型もnumpy.ndarrayであり、ndarrayの引数に整数ではなくてndarrayを使っている。

通常のリストでこれをやると、”リストのインデックスにリストを使ってはダメ”とエラーになる。

a = ['a', 'b', 'c']
print(a)

x = 0
print(a[x])

y = [0]
print(a[y])

# ['a', 'b', 'c']
# a
# ...
# TypeError: list indices must be integers or slices, not list

a = ['a', 'b', 'c']

print(a)

x = 0

print(a[x])

y = [0]

print(a[y])

# ['a', 'b', 'c']

# a

# ...

# TypeError: list indices must be integers or slices, not list

ところがndarrayの場合はその引数にndarrayを渡すことができて、結果はndarrayとして帰ってくる。

import numpy as np

a = np.array(['a', 'b', 'c'])

x = 0
print("{}: {}".format(x, a[x]))

y = np.array([0])
print("{}: {}".format(y, a[y]))

# 0: a
# [0]: ['a']

import numpy as np

a = np.array(['a', 'b', 'c'])

x = 0

print("{}: {}".format(x, a[x]))

y = np.array([0])

print("{}: {}".format(y, a[y]))

# 0: a

# [0]: ['a']

さらに、引数に渡す配列が複数の要素を持つ場合、各要素を引数とした場合に対応する元の配列の要素が配列として返される。

z = np.array([0, 2, 2, 1])
print("{}: {}".format(z, a[z]))

# [0 2 2 1]: ['a' 'c' 'c' 'b']

z = np.array([0, 2, 2, 1])

print("{}: {}".format(z, a[z]))

# [0 2 2 1]: ['a' 'c' 'c' 'b']

なお、元の配列がndarrayであれば、引数に普通のリストを渡しても同じ結果が得られる。

u = [0, 2, 2, 1]
print("{}: {}".format(u, a[u]))

# [0, 2, 2, 1]: ['a' 'c' 'c' 'b']

u = [0, 2, 2, 1]

print("{}: {}".format(u, a[u]))

# [0, 2, 2, 1]: ['a' 'c' 'c' 'b']

Python3 – os/os.path

2019-10-26 / tau / コメントする

主な関数

os.getcwd(): カレントディレクトリー。プログラム実行直後は実行中のディレクトリーを指している。実行ファイルのディレクトリー情報を得るには__file__を使うとよい。
os.path.dirname(dir): dirのディレクトリーのフルパス。
os.listdir(dir): dirにあるファイル・ディレクトリーのリスト。
os.chdir(dir): dirに移動。存在しない場合は’NotADirectoryError’をスローする。

実行例

import sys
import os

if len(sys.argv) > 1:
    target_directory = sys.argv[1]
    try:
        os.chdir(target_directory)
    except NotADirectoryError as e:
        print("Error: Directory not exists")
        sys.exit()

print("Specified directory: {}".format(target_directory))
print()
cwd = os.getcwd()
print('current working directory:', cwd)
print('directory name           :', os.path.dirname(cwd))
print('list of directories      :', os.listdir(cwd))

import sys

import os

if len(sys.argv) > 1:

target_directory = sys.argv[1]

try:

os.chdir(target_directory)

except NotADirectoryError as e:

print("Error: Directory not exists")

sys.exit()

print("Specified directory: {}".format(target_directory))

print()

cwd = os.getcwd()

print('current working directory:', cwd)

print('directory name :', os.path.dirname(cwd))

print('list of directories :', os.listdir(cwd))

実行結果

C:\...\python\packages\os>test_os.py "C:\...\python\packages\os"
Specified directory: C:\...\python\packages\os

current working directory: C:\...\python\packages\os
directory name           : C:\...\python\packages
list of directories      : ['dir1', 'dir2', 'test-01.txt', 'test-02.txt', 'test_os.py']

C:\...\python\packages\os>test_os.py "C:\...\python\packages\os"

Specified directory: C:\...\python\packages\os

current working directory: C:\...\python\packages\os

directory name : C:\...\python\packages

list of directories : ['dir1', 'dir2', 'test-01.txt', 'test-02.txt', 'test_os.py']

ファイル／ディレクトリーの判定

ファイルの判定はos.path.isfile(target)。ディレクトリーの判定はos.path.isdir(target)。

print("\nFile or directory")
for file in os.listdir(cwd):
    if os.path.isfile(file):
        print('f ', end='')
    else:
        print('d ', end='')
    print(file)

print("\nFile or directory")

for file in os.listdir(cwd):

if os.path.isfile(file):

print('f ', end='')

else:

print('d ', end='')

print(file)

実行結果

File or directory
d dir1
d dir2
f test-01.txt
f test-02.txt
f test_os.py

File or directory

d dir1

d dir2

f test-01.txt

f test-02.txt

f test_os.py

Python3 – sys

2019-10-26 / tau / コメントする

プログラムの停止

プログラムを停止するのはsys.exit()。

import sys

args = sys.argv

if len(args) == 1:
    print("No arguments, system exits.")
    sys.exit()
else:
    print("Arguments are {}".format(args[1:]))

# C:...\python\test\sys>exit.py
# No arguments, system exits.
# 
# C:...\python\test\sys>exit.py 1 2
# Arguments are ['1', '2']

import sys

args = sys.argv

if len(args) == 1:

print("No arguments, system exits.")

sys.exit()

else:

print("Arguments are {}".format(args[1:]))

# C:...\python\test\sys>exit.py

# No arguments, system exits.

# C:...\python\test\sys>exit.py 1 2

# Arguments are ['1', '2']

sys.argv～コマンドライン引数

sys.argvはコマンドライン引数を配列で返す。0番目の値は実行中のパス付のスクリプト名自体。

args = sys.argv

print('number of afgs {}'.format(len(args)))
print(args)
for arg in args:
    print(arg)

# C:...\python\test\sys>cmd_arg.py one two
# number of afgs 3
# ['C:...\\python\\test\\sys\\cmd_arg.py', 'one', 'two']
# 0 C:\Users\tomo\Google ドライブ\IT_and_Mobile\dev\python\test\sys\cmd_arg.py
# 1 one
# 2 two

args = sys.argv

print('number of afgs {}'.format(len(args)))

print(args)

for arg in args:

print(arg)

# C:...\python\test\sys>cmd_arg.py one two

# number of afgs 3

# ['C:...\\python\\test\\sys\\cmd_arg.py', 'one', 'two']

# 0 C:\Users\tomo\Google ドライブ\IT_and_Mobile\dev\python\test\sys\cmd_arg.py

# 1 one

# 2 two

Python3 – 正規表現 – group()

2019-10-22 / tau / コメントする

たとえば次のようなファイル名のテキストがあるとする。

file_name-10.txt
file_name-1.txt

ファイル名本体の末尾にある番号を2桁に統一したいようなとき、1つ目にはマッチせず、2つ目にはマッチさせたい。

この場合、以下のような正規表現でマッチングできる。

str1 = 'file_name-10.txt'
str2 = 'file_name-1.txt'

pattern = re.compile(r'-\d\.')

match = re.search(pattern, str1)
print(match)
match = re.search(pattern, str2)
print(match)

# None
# <re.Match object; span=(9, 12), match='-1.'>

str1 = 'file_name-10.txt'

str2 = 'file_name-1.txt'

pattern = re.compile(r'-\d\.')

match = re.search(pattern, str1)

print(match)

match = re.search(pattern, str2)

print(match)

# None

# <re.Match object; span=(9, 12), match='-1.'>

ただしこれだけでは、該当するファイル名はわかるが、0の挿入といった部分的な操作ができない。

そのような場合は、正規表現を()で区切り、group()メソッドで部分列を取り出すことができる。

pattern = re.compile(r'(\w+-)(\d\.\w+)')

match = re.search(pattern, str1)
print(match)
match = re.search(pattern, str2)
print(match)

print(match.group(0))
print(match.group(1))
print(match.group(2))

print(match.group(1) + '0' + match.group(2))

# None
# <re.Match object; span=(0, 15), match='file_name-1.txt'>
# file_name-1.txt
# file_name-
# 1.txt
# file_name-01.txt

pattern = re.compile(r'(\w+-)(\d\.\w+)')

match = re.search(pattern, str1)

print(match)

match = re.search(pattern, str2)

print(match)

print(match.group(0))

print(match.group(1))

print(match.group(2))

print(match.group(1) + '0' + match.group(2))

# None

# <re.Match object; span=(0, 15), match='file_name-1.txt'>

# file_name-1.txt

# file_name-

# 1.txt

# file_name-01.txt

パターン設定で()で区切った部分ごとにグルーピングされ、その各部分を後から再利用できる。なお、group(0)はマッチした部分全体になる。

Python3 – enumerateでリストのインデックスが得られる

2019-10-22 / tau / コメントする

通常、リストの要素を順番に操作するときにはforループを使う。

names = ['JP', 'US', 'EU']
for name in names:
    print(name)

# JP
# US
# EU

names = ['JP', 'US', 'EU']

for name in names:

print(name)

# JP

# US

# EU

enumerate()関数を使うと、リストの要素とそのインデックスを同時に得ることができる。

for i, name in enumerate(names):
    print(i, name)

# 0 JP
# 1 US
# 2 EU

for i, name in enumerate(names):

print(i, name)

# 0 JP

# 1 US

# 2 EU

このインデックスにformat()を使うと書式を設定でき、ファイル名のrenameなどに便利。

for i, name in enumerate(names):
    print('{0:02d}'.format(i) + '-' + name)

# 00-JP
# 01-US
# 02-EU

for i, name in enumerate(names):

print('{0:02d}'.format(i) + '-' + name)

# 00-JP

# 01-US

# 02-EU

Python3 – 正規表現の操作

2019-10-22 / tau / コメントする

主なreモジュール関数

まず、reモジュール関数の操作をまとめる。以下の文字列をターゲットにする。

import re
str = 'The rain in spain stays mainly in the plain.'

1 2	import re str = 'The rain in spain stays mainly in the plain.'

re.match()はパターンが先頭でマッチするかどうか。

print(re.match(r'The', str))
print(re.match(r'the', str))

# <re.Match object; span=(0, 3), match='The'>
# None

print(re.match(r'The', str))

print(re.match(r'the', str))

# <re.Match object; span=(0, 3), match='The'>

# None

re.search()はパターンが含まれるかどうか。最初に現れたパターンにのみマッチする。

print(re.search(r'in', str))
print(re.search(r'inn', str))

# <re.Match object; span=(6, 8), match='in'>
# None

print(re.search(r'in', str))

print(re.search(r'inn', str))

# <re.Match object; span=(6, 8), match='in'>

# None

re.findall()は一致するパターン全てのリストを返す。

print(re.findall(r'in', str))
print(re.findall(r'inn', str))

# ['in', 'in', 'in', 'in', 'in', 'in']
# []

print(re.findall(r'in', str))

print(re.findall(r'inn', str))

# ['in', 'in', 'in', 'in', 'in', 'in']

# []

re.finditer()は一致するパターンのイテレータを返す。

for s in re.finditer(r'in', str):
    print(s)

# <re.Match object; span=(6, 8), match='in'>
# <re.Match object; span=(9, 11), match='in'>
# <re.Match object; span=(15, 17), match='in'>
# <re.Match object; span=(26, 28), match='in'>
# <re.Match object; span=(31, 33), match='in'>
# <re.Match object; span=(41, 43), match='in'>

for s in re.finditer(r'in', str):

print(s)

# <re.Match object; span=(6, 8), match='in'>

# <re.Match object; span=(9, 11), match='in'>

# <re.Match object; span=(15, 17), match='in'>

# <re.Match object; span=(26, 28), match='in'>

# <re.Match object; span=(31, 33), match='in'>

# <re.Match object; span=(41, 43), match='in'>

re.sub()は一致するパターンを置き換え。

print(re.sub(r'in', 'IN', str))

# The raIN IN SpaIN stays maINly IN the plaIN.

print(re.sub(r'in', 'IN', str))

# The raIN IN SpaIN stays maINly IN the plaIN.

re.subn()は一致するパターンを置き換え、その結果と置換回数のタプルを返す。

print(re.subn(r'in', 'IN', str))

# ('The raIN IN SpaIN stays maINly IN the plaIN.', 6)

print(re.subn(r'in', 'IN', str))

# ('The raIN IN SpaIN stays maINly IN the plaIN.', 6)

コンパイル

パターン文字列をコンパイルしパターンオブジェクト化することで、再利用・高速化が可能。この場合、reモジュール関数と同じ名前のメソッドが使える。

pattern = re.compile(r'in')
print(pattern.search(str))
print(pattern.findall(str))
print(pattern.sub('IN', str))
print(pattern.subn('IN', str))

# <re.Match object; span=(6, 8), match='in'>
# ['in', 'in', 'in', 'in', 'in', 'in']
# The raIN IN SpaIN stays maINly IN the plaIN.
# ('The raIN IN SpaIN stays maINly IN the plaIN.', 6)

pattern = re.compile(r'in')

print(pattern.search(str))

print(pattern.findall(str))

print(pattern.sub('IN', str))

print(pattern.subn('IN', str))

# <re.Match object; span=(6, 8), match='in'>

# ['in', 'in', 'in', 'in', 'in', 'in']

# The raIN IN SpaIN stays maINly IN the plaIN.

# ('The raIN IN SpaIN stays maINly IN the plaIN.', 6)

マッチオブジェクトの利用

上記のうちmatch()、search()、findall()はマッチオブジェクトを返す。マッチオブジェクトの主なメソッドには、マッチした開始点を返すstart()、終了点(+1)を返すend()、それらをタプルで範囲として返すspan()、マッチした文字列を返すgroup()がある。

str = 'Peter piper picked a peck of pickled pepper.'
pattern = re.compile(r'[Pp]\w+r')
m = pattern.search(str)
print(m.start(), m.end(), m.span(), m.group())

# 0 5 (0, 5) Peter

str = 'Peter piper picked a peck of pickled pepper.'

pattern = re.compile(r'[Pp]\w+r')

m = pattern.search(str)

print(m.start(), m.end(), m.span(), m.group())

# 0 5 (0, 5) Peter

複数のマッチした文字列を処理するには、finditer()で順次マッチオブジェクトを取り出すとよい。

for m in pattern.finditer(str):
    print(m.start(), m.end(), m.span(), m.group())

# 0 5 (0, 5) Peter
# 6 11 (6, 11) piper
# 37 43 (37, 43) pepper

for m in pattern.finditer(str):

print(m.start(), m.end(), m.span(), m.group())

# 0 5 (0, 5) Peter

# 6 11 (6, 11) piper

# 37 43 (37, 43) pepper

Python3 – 正規表現～貪欲と非貪欲

2019-10-22 / tau / コメントする

以下のようにタグで囲まれた文字からタグを取り出そうとすると、先頭の'<‘と最後尾の’>’に対して文字列全体がマッチしてしまう。

str = '<> < > <abc> <123>'
print(re.findall('<.*>', str))

# ['<> < > <abc> <123>']

str = '<> < > <abc> <123>'

print(re.findall('<.*>', str))

# ['<> < > <abc> <123>']

これは正規表現のマッチングができるだけ長い文字列でヒットさせようとする貪欲な(greedy)マッチングのため。

これを非貪欲(non-greeding)なマッチングにするためには*→*?とする。この場合、できるだけ短くマッチするようになり、個々のタグが分解される。

print(re.findall('<.*?>', str))

# ['<>', '< >', '<abc>', '<123>']

print(re.findall('<.*?>', str))

# ['<>', '< >', '<abc>', '<123>']

タグの間に1文字以上入っていることを意図して以下のようにしても、やはり貪欲にマッチする。

print(re.findall('<.+>', str))

# ['<> < > <abc> <123>']

print(re.findall('<.+>', str))

# ['<> < > <abc> <123>']

これを非貪欲にマッチさせようとすると、以下のように最初のタグのの'>'ではなく、その次の'>'にマッチする。最初のタグの'<'に対して、その後複数の文字を経て到達するのが2つ目のタグの'>'のため。

print(re.findall('<.+?>', str))

# ['<> < >', '<abc>', '<123>']

print(re.findall('<.+?>', str))

# ['<> < >', '<abc>', '<123>']

以下のように、意図する文字構成をできるだけ詳しく記述する方がよい。

print(re.findall('<\w+>', str))

# ['<abc>', '<123>']

print(re.findall('<\w+>', str))

# ['<abc>', '<123>']

Python3 – 文字列の判定～数値か文字か

2019-10-22 / tau / コメントする

文字列を構成する文字がすべて数値か、数値以外の文字かといった判定をするメソッド群。

isdecimal(): 全ての文字が十進数字ならTrue
isdigit(): 全ての文字が数字ならTrue
isnumeric(): 全ての文字が数を表す文字ならTrue
isalpha(): 全ての文字が数字以外の文字ならTrue
isalnum(): 全ての文字が十進数字ならTrue

以下のコードで確認。

strings = ['ABCDE', '12345', 'AB+CD', '①②③', '一二三', '壱弐参', '四＋五']

for str in strings:
    print('for', str)
    print('is decimal ->', str.isdecimal())
    print('is digit   ->', str.isdigit())
    print('is numeric ->', str.isnumeric())
    print('is alpha   ->', str.isalpha())
    print('is alnum   ->', str.isalnum())

strings = ['ABCDE', '12345', 'AB+CD', '①②③', '一二三', '壱弐参', '四＋五']

for str in strings:

print('for', str)

print('is decimal ->', str.isdecimal())

print('is digit ->', str.isdigit())

print('is numeric ->', str.isnumeric())

print('is alpha ->', str.isalpha())

print('is alnum ->', str.isalnum())

半角英字、半角数字の判定結果は想定通り。

for ABCDE
is decimal -> False
is digit   -> False
is numeric -> False
is alpha   -> True
is alnum   -> True
for 12345
is decimal -> True
is digit   -> True
is numeric -> True
is alpha   -> False
is alnum   -> True

for ABCDE

is decimal -> False

is digit -> False

is numeric -> False

is alpha -> True

is alnum -> True

for 12345

is decimal -> True

is digit -> True

is numeric -> True

is alpha -> False

is alnum -> True

半角記号は文字、数字、数を表す文字のいずれでもないと判定。

for AB+CD
is decimal -> False
is digit   -> False
is numeric -> False
is alpha   -> False
is alnum   -> False

for AB+CD

is decimal -> False

is digit -> False

is numeric -> False

is alpha -> False

is alnum -> False

特殊数字は十進数字ではないが数字と判定。

for ①②③
is decimal -> False
is digit   -> True
is numeric -> True
is alpha   -> False
is alnum   -> True

for ①②③

is decimal -> False

is digit -> True

is numeric -> True

is alpha -> False

is alnum -> True

漢数字は十進数字や数字ではないが、数を表す文字と判定。

for 一二三
is decimal -> False
is digit   -> False
is numeric -> True
is alpha   -> True
is alnum   -> True
for 壱弐参
is decimal -> False
is digit   -> False
is numeric -> True
is alpha   -> True
is alnum   -> True

for 一二三

is decimal -> False

is digit -> False

is numeric -> True

is alpha -> True

is alnum -> True

for 壱弐参

is decimal -> False

is digit -> False

is numeric -> True

is alpha -> True

is alnum -> True

全角の記号は文字、数字、数を表す文字のいずれでもないと判定。

for 四＋五
is decimal -> False
is digit   -> False
is numeric -> False
is alpha   -> False
is alnum   -> False

for 四＋五

is decimal -> False

is digit -> False

is numeric -> False

is alpha -> False

is alnum -> False

Python – コマンドライン引数

2019-10-22 / tau / コメントする

コマンドライン引数の取得

以下のような仕様

sysパッケージをインポート
sys.argvに引数のリストが格納されている
引数の1つ目(ゼロ番目)は常にコマンド自身

以下のスクリプトを実行してみる。

import sys

args = sys.argv

print('number of afgs {}'.format(len(args)))
print(args)
for arg in args:
    print(arg)

import sys

args = sys.argv

print('number of afgs {}'.format(len(args)))

print(args)

for arg in args:

print(arg)

引数リストの先頭はコマンド自身(Windowsの\はエスケープされている)。

C:\>cmd_arg.py
number of afgs 1
['C:\\cmd_arg.py']
C:\cmd_arg.py

C:\>cmd_arg.py

number of afgs 1

['C:\\cmd_arg.py']

C:\cmd_arg.py

引数は入力したままの形。

C:\>cmd_arg.py 123 abc 'abc'
number of afgs 4
['C:\\cmd_arg.py', '123', 'abc', "'abc'"]
C:\cmd_arg.py
123
abc
'abc'

C:\>cmd_arg.py 123 abc 'abc'

number of afgs 4

['C:\\cmd_arg.py', '123', 'abc', "'abc'"]

C:\cmd_arg.py

123

abc

'abc'

引数がない場合の処理と、数値・文字別の処理を組み込んだ例。

import sys

args = sys.argv

print('number of afgs {}'.format(len(args)))
args = args[1:]
if len(args) == 0:
    print('no arguments')
    sys.exit()

print(args)
for arg in args:
    if arg.isdecimal():
        print('{} is decimal -> {}'.format(arg, int(arg) * 2))
    else:
        print('{} is not decimal'.format(arg))

import sys

args = sys.argv

print('number of afgs {}'.format(len(args)))

args = args[1:]

if len(args) == 0:

print('no arguments')

sys.exit()

print(args)

for arg in args:

if arg.isdecimal():

print('{} is decimal -> {}'.format(arg, int(arg) * 2))

else:

print('{} is not decimal'.format(arg))

実行結果。

C:\>cmd_arg_dight_check.py
number of afgs 1
no arguments

C:\>cmd_arg_dight_check.py ABC 123
number of afgs 3
['ABC', '123']
ABC is not decimal
123 is decimal -> 246

C:\>cmd_arg_dight_check.py

number of afgs 1

no arguments

C:\>cmd_arg_dight_check.py ABC 123

number of afgs 3

['ABC', '123']

ABC is not decimal

123 is decimal -> 246

numpy – 連立方程式

2019-10-21 / tau / コメントする

導入

以下のように行列表示された連立方程式を考える。

$\begin{equation*} \left( \begin{array}{ccc} 1 & 2 & -3 \\ 2 & -1 & 3 \\ -3 & 2 & 1 \end{array} \right) \left( \begin{array}{c} x \\ y \\ z \end{array} \right) = \left( \begin{array}{c} 5 \\ 4 \\ 1 \end{array} \right) \end{equation*}$

この方程式の解は、 $(x, y, z) = (2, 3, 1)$ 。

逆行列による方法

係数行列、未知数ベクトル、定数ベクトルをそれぞれ $\boldsymbol{A},\boldsymbol{x}, \boldsymbol{b}$ と表す。

$\begin{equation*} \boldsymbol{A} \boldsymbol{x} = \boldsymbol{b} \end{equation*}$

このとき係数行列に逆行列が存在するなら、未知数ベクトルは以下で解ける。

$\begin{gather*} \boldsymbol{A}^{-1} \boldsymbol{A} \boldsymbol{x} = \boldsymbol{A}^{-1} \boldsymbol{b} \\ \boldsymbol{x} = \boldsymbol{A}^{-1} \boldsymbol{b} \\ \end{gather*}$

これをnumpy.linalgパッケージの行列操作で解いてみる。

行列にベクトルを掛けるのに、reshape()で行ベクトルから列ベクトルに変換している点に注意。

A = np.array([[1, 2, -3],
              [2, -1, 3],
              [-3, 2, 1]])

b = np.array([5, 4, 1])
print(LA.inv(A) @ b.reshape(3, 1))

# [[2.]
#  [3.]
#  [1.]]

A = np.array([[1, 2, -3],

[2, -1, 3],

[-3, 2, 1]])

b = np.array([5, 4, 1])

print(LA.inv(A) @ b.reshape(3, 1))

# [[2.]

# [3.]

# [1.]]

numpy.linalg.solve()による方法

solve(A, b)関数は、第1引数に係数行列、第2引数に定数ベクトルを与えて、連立方程式の解のベクトルを得ることができる。

A = np.array([[1, 2, -3],
              [2, -1, 3],
              [-3, 2, 1]])

b = np.array([5, 4, 1])
print(LA.solve(A, b))

# [2. 3. 1.]

A = np.array([[1, 2, -3],

[2, -1, 3],

[-3, 2, 1]])

b = np.array([5, 4, 1])

print(LA.solve(A, b))

# [2. 3. 1.]

非正則行列の場合

以下のような、明らかな非正則行列の場合（連立方程式が不定の場合）、逆行列を計算しようとする時点でsingular matrixのエラーになる。

A = np.array([[1, 1], [2, 2]])
b = np.array([2, 4])
print(LA.inv(A))

# numpy.linalg.LinAlgError: Singular matrix

A = np.array([[1, 1], [2, 2]])

b = np.array([2, 4])

print(LA.inv(A))

# numpy.linalg.LinAlgError: Singular matrix

このようなケースでは、linalg.solve()関数でも同様のエラーとなる。

A = np.array([[1, 1], [2, 2]])
b = np.array([2, 4])
print(LA.solve(A, b))

# numpy.linalg.LinAlgError: Singular matrix

A = np.array([[1, 1], [2, 2]])

b = np.array([2, 4])

print(LA.solve(A, b))

# numpy.linalg.LinAlgError: Singular matrix

以下のようなケースはややこしい。同じ係数行列と定数ベクトルに対して、逆行列による解とsolve()による解の値が異なっている。

A = np.array([[1, 2, 3],
              [4, 5, 6],
              [7, 8, 9]])
b = np.array([1, 2, 3])
print(LA.det(A))
print(LA.inv(A) @ b.reshape(3, 1))
print(LA.solve(A, b))

# 6.66133814775094e-16
# [[ 2.]
#  [-4.]
#  [ 0.]]
# [-0.33333333  0.66666667  0.        ]

A = np.array([[1, 2, 3],

[4, 5, 6],

[7, 8, 9]])

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

print(LA.det(A))

print(LA.inv(A) @ b.reshape(3, 1))

print(LA.solve(A, b))

# 6.66133814775094e-16

# [[ 2.]

# [-4.]

# [ 0.]]

# [-0.33333333 0.66666667 0. ]

行列Aの行列式は理論上はゼロであることが確認できる。

$\begin{eqnarray*} \left| \begin{array}{ccc} 1 & 2 &3 \\ 4 & 5 & 6 \\ 7 & 8 & 9 \end{array} \right| &=& 1 \cdot (5 \cdot 9 - 6 \cdot 8) - 2 \cdot (4 \cdot 9 - 6 \cdot 7) + 3 \cdot(4 \cdot 8 - 5 \cdot 7) \\ &=& 45 - 48 -2(36 - 42) + 3(32 - 35) \\ &=& -3 +12 + 9 = 0 \end{eqnaray*}$

この方程式を掃き出し法で解いていくと以下の通り。

$\begin{equation*} \left( \begin{array}{ccc} 1 & 2 & 3 \\ 4 & 5 & 6 \\ 7 & 8 & 9 \end{array} \right) \left( \begin{array}{c} x \\ y \\ z \end{array} \right) = \left( \begin{array}{c} 1 \\ 2 \\ 3 \end{array} \right) \end{equation*}$

$\begin{equation*} \left( \begin{array}{ccc} 1 & 2 & 3 \\0 & -3 & -6 \\ 0 & -6 & -12 \end{array} \right) \left( \begin{array}{c} x \\ y \\ z \end{array} \right) = \left( \begin{array}{c} 1 & -2 & -4 \end{array} \right) \end{equation*}$

$\begin{equation*} \left( \begin{array}{ccc} 1 & 2 & 3 \\0 & -3 & -6 \\ 0 & 0 & 0 \end{array} \right) \left( \begin{array}{c} x \\ y \\ z \end{array} \right) = \left( \begin{array}{c} 1 & -2 & 0 \end{array} \right) \end{equation*}$

途中省略するが、ここでzを消去すると、 $y = -2x$ となり、先の計算結果と符合するが、zは0である必要はない。

不定連立方程式において、逆行列やsolve()関数を使って解く場合には注意が必要。

TauStation

タグ: Python

ndarray – 配列による要素指定

Python3 – os/os.path

主な関数

ファイル／ディレクトリーの判定

Python3 – sys

プログラムの停止

sys.argv～コマンドライン引数

Python3 – 正規表現 – group()

Python3 – enumerateでリストのインデックスが得られる

Python3 – 正規表現の操作

主なreモジュール関数

コンパイル

マッチオブジェクトの利用

Python3 – 正規表現～貪欲と非貪欲

Python3 – 文字列の判定～数値か文字か

Python – コマンドライン引数

コマンドライン引数の取得

numpy – 連立方程式

導入

逆行列による方法

numpy.linalg.solve()による方法

非正則行列の場合