algebra 0.2.1
Sign up to get free protection for your applications and to get access to all the features.
- checksums.yaml +7 -0
- data/.gitignore +53 -0
- data/.travis.yml +7 -0
- data/CODE_OF_CONDUCT.md +49 -0
- data/Gemfile +7 -0
- data/Gemfile.lock +28 -0
- data/LICENSE.txt +21 -0
- data/README-ja.txt +60 -0
- data/README.md +61 -0
- data/Rakefile +43 -0
- data/algebra.gemspec +33 -0
- data/bin/console +14 -0
- data/bin/setup +8 -0
- data/doc-ja/README-ja.html +69 -0
- data/doc-ja/README-ja.rd +60 -0
- data/doc-ja/README.html +0 -0
- data/doc-ja/README.rd +0 -0
- data/doc-ja/algebra-ja.html +66 -0
- data/doc-ja/algebra-ja.rd +60 -0
- data/doc-ja/algebraic-equation-ja.html +36 -0
- data/doc-ja/algebraic-equation-ja.rd +26 -0
- data/doc-ja/algebraic-extension-field-ja.html +122 -0
- data/doc-ja/algebraic-extension-field-ja.rd +118 -0
- data/doc-ja/algebraic-parser-ja.html +87 -0
- data/doc-ja/algebraic-parser-ja.rd +95 -0
- data/doc-ja/changes.html +161 -0
- data/doc-ja/changes.rd +119 -0
- data/doc-ja/elementary-divisor-ja.html +56 -0
- data/doc-ja/elementary-divisor-ja.rd +45 -0
- data/doc-ja/euclidian-ring-ja.html +78 -0
- data/doc-ja/euclidian-ring-ja.rd +69 -0
- data/doc-ja/finite-group-ja.html +285 -0
- data/doc-ja/finite-group-ja.rd +276 -0
- data/doc-ja/finite-map-ja.html +224 -0
- data/doc-ja/finite-map-ja.rd +227 -0
- data/doc-ja/finite-set-ja.html +393 -0
- data/doc-ja/finite-set-ja.rd +408 -0
- data/doc-ja/index-ja.html +113 -0
- data/doc-ja/index-ja.rd +112 -0
- data/doc-ja/jordan-form-ja.html +106 -0
- data/doc-ja/jordan-form-ja.rd +99 -0
- data/doc-ja/localized-ring-ja.html +121 -0
- data/doc-ja/localized-ring-ja.rd +131 -0
- data/doc-ja/m-polynomial-ja.html +542 -0
- data/doc-ja/m-polynomial-ja.rd +596 -0
- data/doc-ja/matrix-algebra-ja.html +686 -0
- data/doc-ja/matrix-algebra-ja.rd +727 -0
- data/doc-ja/matrix-algebra-triplet-ja.html +124 -0
- data/doc-ja/matrix-algebra-triplet-ja.rd +129 -0
- data/doc-ja/permutation-group-ja.html +159 -0
- data/doc-ja/permutation-group-ja.rd +151 -0
- data/doc-ja/polynomial-converter-ja.html +64 -0
- data/doc-ja/polynomial-converter-ja.rd +54 -0
- data/doc-ja/polynomial-ja.html +446 -0
- data/doc-ja/polynomial-ja.rd +473 -0
- data/doc-ja/residue-class-ring-ja.html +112 -0
- data/doc-ja/residue-class-ring-ja.rd +113 -0
- data/doc-ja/sample-algebraic-equation01.rb.v.rd +11 -0
- data/doc-ja/sample-algebraic-equation02.rb.v.rd +11 -0
- data/doc-ja/sample-algebraic-root01.rb.v.rd +20 -0
- data/doc-ja/sample-algebraicfield01.rb.v.rd +26 -0
- data/doc-ja/sample-algebraicfield02.rb.v.rd +15 -0
- data/doc-ja/sample-cayleyhamilton01.rb.v.rd +14 -0
- data/doc-ja/sample-diagonalization01.rb.v.rd +48 -0
- data/doc-ja/sample-divmod01.rb.v.rd +17 -0
- data/doc-ja/sample-elementary-divisor01.rb.v.rd +48 -0
- data/doc-ja/sample-factorize01.rb.v.rd +11 -0
- data/doc-ja/sample-factorize02.rb.v.rd +13 -0
- data/doc-ja/sample-factorize03.rb.v.rd +14 -0
- data/doc-ja/sample-factorize04.rb.v.rd +14 -0
- data/doc-ja/sample-factorize05.rb.v.rd +53 -0
- data/doc-ja/sample-galois-group01.rb.v.rd +27 -0
- data/doc-ja/sample-gaussian-elimination01.rb.v.rd +19 -0
- data/doc-ja/sample-geometry01.rb.v.rd +33 -0
- data/doc-ja/sample-geometry02.rb.v.rd +34 -0
- data/doc-ja/sample-geometry03.rb.v.rd +45 -0
- data/doc-ja/sample-geometry04.rb.v.rd +33 -0
- data/doc-ja/sample-geometry07.rb.v.rd +62 -0
- data/doc-ja/sample-groebner01.rb.v.rd +13 -0
- data/doc-ja/sample-groebner02.rb.v.rd +18 -0
- data/doc-ja/sample-groebner03.rb.v.rd +19 -0
- data/doc-ja/sample-group01.rb.v.rd +21 -0
- data/doc-ja/sample-jordan-form01.rb.v.rd +56 -0
- data/doc-ja/sample-jordanform01.rb.v.rd +55 -0
- data/doc-ja/sample-lagrange-multiplier01.rb.v.rd +35 -0
- data/doc-ja/sample-m-factorize01.rb.v.rd +15 -0
- data/doc-ja/sample-m-factorize02.rb.v.rd +12 -0
- data/doc-ja/sample-m-polynomial01.rb.v.rd +10 -0
- data/doc-ja/sample-map01.rb.v.rd +9 -0
- data/doc-ja/sample-polynomial01.rb.v.rd +9 -0
- data/doc-ja/sample-polynomial02.rb.v.rd +11 -0
- data/doc-ja/sample-primefield01.rb.v.rd +13 -0
- data/doc-ja/sample-quotientfield01.rb.v.rd +11 -0
- data/doc-ja/sample-quotientfield02.rb.v.rd +18 -0
- data/doc-ja/sample-quotientfield03.rb.v.rd +16 -0
- data/doc-ja/sample-quotientfield04.rb.v.rd +16 -0
- data/doc-ja/sample-set01.rb.v.rd +18 -0
- data/doc-ja/sample-splitting-field01.rb.v.rd +19 -0
- data/doc-ja/samples-ja.html +885 -0
- data/doc-ja/samples-ja.rd +221 -0
- data/doc-ja/style.css +85 -0
- data/doc-ja/todo.html +20 -0
- data/doc-ja/todo.rd +9 -0
- data/lib/algebra.rb +36 -0
- data/lib/algebra/algebraic-equation.rb +56 -0
- data/lib/algebra/algebraic-extension-field.rb +135 -0
- data/lib/algebra/algebraic-parser.rb +160 -0
- data/lib/algebra/algebraic-system.rb +224 -0
- data/lib/algebra/annihilate.rb +52 -0
- data/lib/algebra/array-supplement.rb +53 -0
- data/lib/algebra/auto-require.rb +84 -0
- data/lib/algebra/chinese-rem-th.rb +135 -0
- data/lib/algebra/combinatorial.rb +145 -0
- data/lib/algebra/elementary-divisor.rb +193 -0
- data/lib/algebra/euclidian-ring.rb +161 -0
- data/lib/algebra/factors.rb +305 -0
- data/lib/algebra/finite-group.rb +374 -0
- data/lib/algebra/finite-map.rb +201 -0
- data/lib/algebra/finite-set.rb +456 -0
- data/lib/algebra/galois-group.rb +129 -0
- data/lib/algebra/gaussian-elimination.rb +385 -0
- data/lib/algebra/groebner-basis-coeff.rb +228 -0
- data/lib/algebra/groebner-basis.rb +197 -0
- data/lib/algebra/import-module-single-thread.rb +86 -0
- data/lib/algebra/import-module.rb +491 -0
- data/lib/algebra/jordan-form.rb +114 -0
- data/lib/algebra/linear-algebra.rb +143 -0
- data/lib/algebra/localized-ring.rb +294 -0
- data/lib/algebra/m-index.rb +282 -0
- data/lib/algebra/m-polynomial-factor-int.rb +186 -0
- data/lib/algebra/m-polynomial-factor-zp.rb +114 -0
- data/lib/algebra/m-polynomial-factor.rb +315 -0
- data/lib/algebra/m-polynomial-gcd.rb +40 -0
- data/lib/algebra/m-polynomial.rb +875 -0
- data/lib/algebra/matrix-algebra-triplet.rb +292 -0
- data/lib/algebra/matrix-algebra.rb +929 -0
- data/lib/algebra/numeric-supplement.rb +123 -0
- data/lib/algebra/permutation-group.rb +257 -0
- data/lib/algebra/polynomial-converter.rb +193 -0
- data/lib/algebra/polynomial-factor-alg.rb +148 -0
- data/lib/algebra/polynomial-factor-int.rb +252 -0
- data/lib/algebra/polynomial-factor-zp.rb +165 -0
- data/lib/algebra/polynomial-factor.rb +140 -0
- data/lib/algebra/polynomial.rb +592 -0
- data/lib/algebra/powers.rb +31 -0
- data/lib/algebra/prime-gen.rb +32 -0
- data/lib/algebra/rational.rb +81 -0
- data/lib/algebra/residue-class-ring.rb +219 -0
- data/lib/algebra/sets-system.rb +43 -0
- data/lib/algebra/splitting-field.rb +103 -0
- data/lib/algebra/version.rb +3 -0
- data/sample/Makefile +10 -0
- data/sample/do-testscripts.rb +100 -0
- data/sample/sample-algebraic-equation01.rb +6 -0
- data/sample/sample-algebraic-equation02.rb +6 -0
- data/sample/sample-algebraic-root01.rb +15 -0
- data/sample/sample-algebraicfield01.rb +21 -0
- data/sample/sample-algebraicfield02.rb +10 -0
- data/sample/sample-cayleyhamilton01.rb +9 -0
- data/sample/sample-diagonalization01.rb +43 -0
- data/sample/sample-divmod01.rb +12 -0
- data/sample/sample-elementary-divisor01.rb +43 -0
- data/sample/sample-factorize01.rb +6 -0
- data/sample/sample-factorize02.rb +8 -0
- data/sample/sample-factorize03.rb +9 -0
- data/sample/sample-factorize04.rb +9 -0
- data/sample/sample-factorize05.rb +48 -0
- data/sample/sample-galois-group01.rb +22 -0
- data/sample/sample-gaussian-elimination01.rb +14 -0
- data/sample/sample-geometry01.rb +28 -0
- data/sample/sample-geometry02.rb +29 -0
- data/sample/sample-geometry03.rb +40 -0
- data/sample/sample-geometry04.rb +28 -0
- data/sample/sample-geometry07.rb +58 -0
- data/sample/sample-groebner01.rb +8 -0
- data/sample/sample-groebner02.rb +13 -0
- data/sample/sample-groebner03.rb +14 -0
- data/sample/sample-group01.rb +16 -0
- data/sample/sample-jordan-form01.rb +51 -0
- data/sample/sample-lagrange-multiplier01.rb +30 -0
- data/sample/sample-m-factorize01.rb +10 -0
- data/sample/sample-m-factorize02.rb +7 -0
- data/sample/sample-m-polynomial01.rb +5 -0
- data/sample/sample-map01.rb +4 -0
- data/sample/sample-polynomial01.rb +4 -0
- data/sample/sample-polynomial02.rb +6 -0
- data/sample/sample-primefield01.rb +8 -0
- data/sample/sample-quotientfield01.rb +6 -0
- data/sample/sample-quotientfield02.rb +13 -0
- data/sample/sample-quotientfield03.rb +11 -0
- data/sample/sample-quotientfield04.rb +11 -0
- data/sample/sample-set01.rb +13 -0
- data/sample/sample-splitting-field01.rb +14 -0
- data/sample/test-00-cayley-hamilton.rb +76 -0
- data/sample/test-00-groebner-basis.rb +274 -0
- data/sample/test-00-polynomial-factor-alg.rb +89 -0
- data/sample/test-00.list +4 -0
- data/sample/time-trial.rb +65 -0
- metadata +373 -0
@@ -0,0 +1,686 @@
|
|
1
|
+
<?xml version="1.0" encoding="utf-8" ?>
|
2
|
+
<!DOCTYPE html
|
3
|
+
PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN"
|
4
|
+
"http://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd">
|
5
|
+
<html xmlns="http://www.w3.org/1999/xhtml">
|
6
|
+
<head>
|
7
|
+
<meta http-equiv="Content-type" content="text/html; charset=utf-8" />
|
8
|
+
<title>doc-ja/matrix-algebra-ja.rd</title>
|
9
|
+
<link href="style.css" type="text/css" rel="stylesheet" />
|
10
|
+
</head>
|
11
|
+
<body>
|
12
|
+
<p>[<a href="index-ja.html">index-ja</a>]
|
13
|
+
<a href="#label-0">Algebra::MatrixAlgebra</a>
|
14
|
+
/
|
15
|
+
<a href="#label-60">Algebra::Vector</a>
|
16
|
+
/
|
17
|
+
<a href="#label-78">Algebra::Covector</a>
|
18
|
+
/
|
19
|
+
<a href="#label-96">Algebra::SquareMatrix</a>
|
20
|
+
/
|
21
|
+
<a href="#label-117">Algebra::GaussianElimination</a></p>
|
22
|
+
<h1><a name="label-0" id="label-0">Algebra::MatrixAlgebra</a></h1><!-- RDLabel: "Algebra::MatrixAlgebra" -->
|
23
|
+
<p><em>(行列クラス)</em></p>
|
24
|
+
<p>行列を表現します。実際のクラスを生成するには基底環とサイズを指定して、
|
25
|
+
クラスメソッド <a href="#label-7">::create</a> あるいは関数 <a href="#label-0">Algebra::MatrixAlgebra</a>()
|
26
|
+
を用います。</p>
|
27
|
+
<p>サブクラスとして <a href="#label-60">Algebra::Vector</a>(縦ベクトル),
|
28
|
+
<a href="#label-78">Algebra::Covector</a>(横ベクトル),
|
29
|
+
<a href="#label-96">Algebra::SquareMatrix</a>(正方行列) を持ちます。</p>
|
30
|
+
<h2><a name="label-1" id="label-1">ファイル名:</a></h2><!-- RDLabel: "ファイル名:" -->
|
31
|
+
<ul>
|
32
|
+
<li><var>matrix-algebra.rb</var></li>
|
33
|
+
</ul>
|
34
|
+
<h2><a name="label-2" id="label-2">スーパークラス:</a></h2><!-- RDLabel: "スーパークラス:" -->
|
35
|
+
<ul>
|
36
|
+
<li><var>Object</var></li>
|
37
|
+
</ul>
|
38
|
+
<h2><a name="label-3" id="label-3">インクルードしているモジュール:</a></h2><!-- RDLabel: "インクルードしているモジュール:" -->
|
39
|
+
<ul>
|
40
|
+
<li><a href="#label-117">Algebra::GaussianElimination</a></li>
|
41
|
+
<li>Enumerable</li>
|
42
|
+
</ul>
|
43
|
+
<h2><a name="label-4" id="label-4">関連する関数:</a></h2><!-- RDLabel: "関連する関数:" -->
|
44
|
+
<dl>
|
45
|
+
<dt><a name="label-5" id="label-5"><code>Algebra.MatrixAlgebra(<var>ring</var>, <var>m</var>, <var>n</var>)</code></a></dt><!-- RDLabel: "Algebra.MatrixAlgebra" -->
|
46
|
+
<dd>
|
47
|
+
<a href="#label-7">::create</a>(ring, m, n)と同じです。</dd>
|
48
|
+
</dl>
|
49
|
+
<h2><a name="label-6" id="label-6">クラスメソッド:</a></h2><!-- RDLabel: "クラスメソッド:" -->
|
50
|
+
<dl>
|
51
|
+
<dt><a name="label-7" id="label-7"><code>::create(<var>ring</var>, <var>m</var>, <var>n</var>)</code></a></dt><!-- RDLabel: "::create" -->
|
52
|
+
<dd>
|
53
|
+
<p>環 <var>ring</var> を要素とする, <code> (m, n) </code> 型の行列を
|
54
|
+
表現するクラスを生成します。</p>
|
55
|
+
<p>このメソッドの戻り値は <a href="#label-0">Algebra::MatrixAlgebra</a> クラスのサブクラス
|
56
|
+
です。このサブクラスにはクラスメソッドとして <var>ground</var> と
|
57
|
+
<var>rsize</var>, <var>csize</var>, <var>sizes</var> が定義され、それぞれ、
|
58
|
+
基底となる環 <var>ring</var>、行のサイズ <var>m</var>、列のサイズ <var>n</var>、
|
59
|
+
それらの配列 <code> [m, n] </code> を返します。</p>
|
60
|
+
<p>実際に行列を作るには <a href="#label-8">::new</a>, <a href="#label-9">::matrix</a>,
|
61
|
+
<a href="#label-10">::[]</a> を使います。</p></dd>
|
62
|
+
<dt><a name="label-8" id="label-8"><code>::new(<var>array</var>)</code></a></dt><!-- RDLabel: "::new" -->
|
63
|
+
<dd>
|
64
|
+
<p><var>array</var> を配列の配列とするとき、それを要素とする行列を返します。</p>
|
65
|
+
<p>例:</p>
|
66
|
+
<pre>M = Algebra.MatrixAlgebra(Integer, 2, 3)
|
67
|
+
a = M.new([[1, 2, 3], [4, 5, 6]])
|
68
|
+
a.display
|
69
|
+
#=> [1, 2, 3]
|
70
|
+
#=> [4, 5, 6]</pre></dd>
|
71
|
+
<dt><a name="label-9" id="label-9"><code>::matrix{|<var>i</var>, <var>j</var>| ... }</code></a></dt><!-- RDLabel: "::matrix" -->
|
72
|
+
<dd>
|
73
|
+
<p><var>i</var> と <var>j</var> に行と列のインデックスを与え ... を評価した値を
|
74
|
+
<code> (i, j) </code> 成分にした行列を返します。</p>
|
75
|
+
<p>例:</p>
|
76
|
+
<pre>M = Alebra.MatrixAlgebra(Integer, 2, 3)
|
77
|
+
a = M.matrix{|i, j| 10*(i + 1) + j + 1}
|
78
|
+
a.display
|
79
|
+
#=> [11, 12, 13]
|
80
|
+
#=> [21, 22, 23]</pre></dd>
|
81
|
+
<dt><a name="label-10" id="label-10"><code>::[<var>array1</var>, <var>array2</var>, ..., <var>arrayr</var>]</code></a></dt><!-- RDLabel: "::[]" -->
|
82
|
+
<dd>
|
83
|
+
<p>配列 <code>array1, array2, ..., arrayr</code> をそれぞれ行とする配列を返します。</p>
|
84
|
+
<p>例:</p>
|
85
|
+
<pre>M = Algebra.MatrixAlgebra(Integer, 2, 3)
|
86
|
+
a = M[[1, 2, 3], [4, 5, 6]]
|
87
|
+
a.display
|
88
|
+
#=> [1, 2, 3]
|
89
|
+
#=> [4, 5, 6]</pre></dd>
|
90
|
+
<dt><a name="label-11" id="label-11"><code>::collect_ij{|<var>i</var>, <var>j</var>| ... }</code></a></dt><!-- RDLabel: "::collect_ij" -->
|
91
|
+
<dd>
|
92
|
+
<a href="#label-9">::matrix</a> によく似ていますが、
|
93
|
+
返り値は Algebra::MatrixAlgebra でなく、2重配列(Array の Array)です。</dd>
|
94
|
+
<dt><a name="label-12" id="label-12"><code>::collect_row{|<var>i</var>| ... }</code></a></dt><!-- RDLabel: "::collect_row" -->
|
95
|
+
<dd>
|
96
|
+
<p>第 i 行に ... を評価して得た配列を並べた行列を返します。</p>
|
97
|
+
<p>例:</p>
|
98
|
+
<pre>M = Algebra.MatrixAlgebra(Integer, 2, 3)
|
99
|
+
A = M.collect_row{|i| [i*10 + 11, i*10 + 12, i*10 + 13]}
|
100
|
+
A.display
|
101
|
+
#=> [11, 12, 13]
|
102
|
+
#=> [21, 22, 23]</pre></dd>
|
103
|
+
<dt><a name="label-13" id="label-13"><code>::collect_column{|<var>j</var>| ... }</code></a></dt><!-- RDLabel: "::collect_column" -->
|
104
|
+
<dd>
|
105
|
+
<p>第 j 列に ... を評価して得た配列を並べた行列を返します。</p>
|
106
|
+
<p>例:</p>
|
107
|
+
<pre>M = Algebra.MatrixAlgebra(Integer, 2, 3)
|
108
|
+
A = M.collect_column{|j| [11 + j, 21 + j]}
|
109
|
+
A.display
|
110
|
+
#=> [11, 12, 13]
|
111
|
+
#=> [21, 22, 23]</pre></dd>
|
112
|
+
<dt><a name="label-14" id="label-14"><code>::*(<var>otype</var>)</code></a></dt><!-- RDLabel: "::*" -->
|
113
|
+
<dd>
|
114
|
+
<p>2つの行列の型を掛けた新しいクラスを返します。</p>
|
115
|
+
<p>例:</p>
|
116
|
+
<pre>M = Algebra.MatrixAlgebra(Integer, 2, 3)
|
117
|
+
N = Algebra.MatrixAlgebra(Integer, 3, 4)
|
118
|
+
L = M * N
|
119
|
+
p L.sizes #=> [3, 4]</pre></dd>
|
120
|
+
<dt><a name="label-15" id="label-15"><code>::vector_type</code></a></dt><!-- RDLabel: "::vector_type" -->
|
121
|
+
<dd>
|
122
|
+
<a href="#label-22">rsize</a> と同じサイズの縦ベクトル(Vector)クラスを返します。</dd>
|
123
|
+
<dt><a name="label-16" id="label-16"><code>::covector_type</code></a></dt><!-- RDLabel: "::covector_type" -->
|
124
|
+
<dd>
|
125
|
+
<a href="#label-23">csize</a> と同じサイズの横ベクトル(CoVector)クラスを返します。</dd>
|
126
|
+
<dt><a name="label-17" id="label-17"><code>::transpose</code></a></dt><!-- RDLabel: "::transpose" -->
|
127
|
+
<dd>
|
128
|
+
転置を行った新しい行列クラスを返します。</dd>
|
129
|
+
<dt><a name="label-18" id="label-18"><code>::zero</code></a></dt><!-- RDLabel: "::zero" -->
|
130
|
+
<dd>
|
131
|
+
零行列を返します。</dd>
|
132
|
+
</dl>
|
133
|
+
<h2><a name="label-19" id="label-19">メソッド:</a></h2><!-- RDLabel: "メソッド:" -->
|
134
|
+
<dl>
|
135
|
+
<dt><a name="label-20" id="label-20"><code>[<var>i</var>, <var>j</var>]</code></a></dt><!-- RDLabel: "[]" -->
|
136
|
+
<dd>
|
137
|
+
<code>(i, j)</code> 成分を返します。</dd>
|
138
|
+
<dt><a name="label-21" id="label-21"><code>[<var>i</var>, <var>j</var>] = <var>x</var></code></a></dt><!-- RDLabel: "[]=" -->
|
139
|
+
<dd>
|
140
|
+
<code>(i, j)</code> 成分を x にします。</dd>
|
141
|
+
<dt><a name="label-22" id="label-22"><code>rsize</code></a></dt><!-- RDLabel: "rsize" -->
|
142
|
+
<dd>
|
143
|
+
行サイズを返します。<!-- Reference, RDLabel "::rsize" doesn't exist --><em class="label-not-found">::rsize</em><!-- Reference end --> と同じです。</dd>
|
144
|
+
<dt><a name="label-23" id="label-23"><code>csize</code></a></dt><!-- RDLabel: "csize" -->
|
145
|
+
<dd>
|
146
|
+
列サイズを返します。<!-- Reference, RDLabel "::csize" doesn't exist --><em class="label-not-found">::csize</em><!-- Reference end --> と同じです。</dd>
|
147
|
+
<dt><a name="label-24" id="label-24"><code>sizes</code></a></dt><!-- RDLabel: "sizes" -->
|
148
|
+
<dd>
|
149
|
+
[<a href="#label-22">rsize</a>, <a href="#label-23">csize</a>] の配列を返します。
|
150
|
+
<!-- Reference, RDLabel "::sizes" doesn't exist --><em class="label-not-found">::sizes</em><!-- Reference end --> と同じです。</dd>
|
151
|
+
<dt><a name="label-25" id="label-25"><code>rows</code></a></dt><!-- RDLabel: "rows" -->
|
152
|
+
<dd>
|
153
|
+
<p>各行を要素とする配列を返します。</p>
|
154
|
+
<p>例:</p>
|
155
|
+
<pre>M = Algebra.MatrixAlgebra(Integer, 2, 3)
|
156
|
+
a = M.new([[1, 2, 3], [4, 5, 6]])
|
157
|
+
p a.rows #=> [[1, 2, 3], [4, 5, 6]]
|
158
|
+
p a.row(1) #=> [4, 5, 6]
|
159
|
+
a.set_row(1, [40, 50, 60])
|
160
|
+
a.display #=> [1, 2, 3]
|
161
|
+
#=> [40, 50, 60]</pre></dd>
|
162
|
+
<dt><a name="label-26" id="label-26"><code>row(<var>i</var>)</code></a></dt><!-- RDLabel: "row" -->
|
163
|
+
<dd>
|
164
|
+
i 行目を配列として返します。</dd>
|
165
|
+
<dt><a name="label-27" id="label-27"><code>set_row(<var>i</var>, <var>array</var>)</code></a></dt><!-- RDLabel: "set_row" -->
|
166
|
+
<dd>
|
167
|
+
i 行目を配列 array に入れ換えます。</dd>
|
168
|
+
<dt><a name="label-28" id="label-28"><code>columns</code></a></dt><!-- RDLabel: "columns" -->
|
169
|
+
<dd>
|
170
|
+
<p>各列を要素とする配列を返します。</p>
|
171
|
+
<p>例:</p>
|
172
|
+
<pre>M = Algebra.MatrixAlgebra(Integer, 2, 3)
|
173
|
+
a = M.new([[1, 2, 3], [4, 5, 6]])
|
174
|
+
p a.columns #=> [[1, 4], [2, 5], [3, 6]]
|
175
|
+
p a.column(1) #=> [2, 5]
|
176
|
+
a.set_column(1, [20, 50])
|
177
|
+
a.display #=> [1, 20, 3]
|
178
|
+
#=> [4, 50, 6]</pre></dd>
|
179
|
+
<dt><a name="label-29" id="label-29"><code>column(<var>j</var>)</code></a></dt><!-- RDLabel: "column" -->
|
180
|
+
<dd>
|
181
|
+
j 列目を配列として返します。</dd>
|
182
|
+
<dt><a name="label-30" id="label-30"><code>set_column(<var>j</var>, <var>array</var>)</code></a></dt><!-- RDLabel: "set_column" -->
|
183
|
+
<dd>
|
184
|
+
j 列目を配列 array に入れ換えます。</dd>
|
185
|
+
<dt><a name="label-31" id="label-31"><code>each{|<var>row</var>| ...}</code></a></dt><!-- RDLabel: "each" -->
|
186
|
+
<dd>
|
187
|
+
各行を配列にして <var>row</var> に入れるイテレータです。</dd>
|
188
|
+
<dt><a name="label-32" id="label-32"><code>each_index{|<var>i</var>, <var>j</var>| ...}</code></a></dt><!-- RDLabel: "each_index" -->
|
189
|
+
<dd>
|
190
|
+
各添え字 <code> (i, j) </code> に関するイテレータです。</dd>
|
191
|
+
<dt><a name="label-33" id="label-33"><code>each_i{|<var>i</var>| ...}</code></a></dt><!-- RDLabel: "each_i" -->
|
192
|
+
<dd>
|
193
|
+
各行の添え字 <code>i</code> に関するイテレータです。</dd>
|
194
|
+
<dt><a name="label-34" id="label-34"><code>each_j{|<var>j</var>| ...}</code></a></dt><!-- RDLabel: "each_j" -->
|
195
|
+
<dd>
|
196
|
+
各列の添え字 <code>j</code> に関するイテレータです。</dd>
|
197
|
+
<dt><a name="label-35" id="label-35"><code>each_row{|<var>r</var>| ... }</code></a></dt><!-- RDLabel: "each_row" -->
|
198
|
+
<dd>
|
199
|
+
各行を配列にして <var>r</var> に入れるイテレータです。
|
200
|
+
<a href="#label-31">each</a> と同じです。</dd>
|
201
|
+
<dt><a name="label-36" id="label-36"><code>each_column{|<var>c</var>| ... }</code></a></dt><!-- RDLabel: "each_column" -->
|
202
|
+
<dd>
|
203
|
+
各列を配列にして <var>c</var> に入れるイテレータです。</dd>
|
204
|
+
<dt><a name="label-37" id="label-37"><code>matrix{|<var>i</var>, <var>j</var>| ... }</code></a></dt><!-- RDLabel: "matrix" -->
|
205
|
+
<dd>
|
206
|
+
<a href="#label-9">::matrix</a> と同じです。</dd>
|
207
|
+
<dt><a name="label-38" id="label-38"><code>collect_ij{|<var>i</var>, <var>j</var>| ... }</code></a></dt><!-- RDLabel: "collect_ij" -->
|
208
|
+
<dd>
|
209
|
+
<a href="#label-11">::collect_ij</a> と同じです。</dd>
|
210
|
+
<dt><a name="label-39" id="label-39"><code>collect_row{|<var>i</var>| ... }</code></a></dt><!-- RDLabel: "collect_row" -->
|
211
|
+
<dd>
|
212
|
+
<a href="#label-12">::collect_row</a> と同じです。</dd>
|
213
|
+
<dt><a name="label-40" id="label-40"><code>collect_column{|<var>j</var>| ... }</code></a></dt><!-- RDLabel: "collect_column" -->
|
214
|
+
<dd>
|
215
|
+
<a href="#label-13">::collect_column</a> と同じです。</dd>
|
216
|
+
<dt><a name="label-41" id="label-41"><code>minor(<var>i</var>, <var>j</var>)</code></a></dt><!-- RDLabel: "minor" -->
|
217
|
+
<dd>
|
218
|
+
<var>i</var> 行 <var>j</var> 列を除いた小行列を返します。</dd>
|
219
|
+
<dt><a name="label-42" id="label-42"><code>cofactor(<var>i</var>, <var>j</var>)</code></a></dt><!-- RDLabel: "cofactor" -->
|
220
|
+
<dd>
|
221
|
+
<var>i</var> 行 <var>j</var> 列を除いた小行列式に (-1)**(i+j) を掛けたものを
|
222
|
+
返します。<code>minor(i, j) ** (i + j)</code> と同じです。</dd>
|
223
|
+
<dt><a name="label-43" id="label-43"><code>cofactor_matrix</code></a></dt><!-- RDLabel: "cofactor_matrix" -->
|
224
|
+
<dd>
|
225
|
+
余因子行列を返します。<code>self.class.transpose.matrix{|i, j| cofactor(j, i)}</code> と同じです。</dd>
|
226
|
+
<dt><a name="label-44" id="label-44"><code>adjoint</code></a></dt><!-- RDLabel: "adjoint" -->
|
227
|
+
<dd>
|
228
|
+
<a href="#label-43">cofactor_matrix</a> と同じです。</dd>
|
229
|
+
<dt><a name="label-45" id="label-45"><code>==(<var>other</var>)</code></a></dt><!-- RDLabel: "==" -->
|
230
|
+
<dd>
|
231
|
+
等しいとき真を返します。</dd>
|
232
|
+
<dt><a name="label-46" id="label-46"><code>+(<var>other</var>)</code></a></dt><!-- RDLabel: "+" -->
|
233
|
+
<dd>
|
234
|
+
和を計算します。</dd>
|
235
|
+
<dt><a name="label-47" id="label-47"><code>-(<var>other</var>)</code></a></dt><!-- RDLabel: "-" -->
|
236
|
+
<dd>
|
237
|
+
差を計算します。</dd>
|
238
|
+
<dt><a name="label-48" id="label-48"><code>*(<var>other</var>)</code></a></dt><!-- RDLabel: "*" -->
|
239
|
+
<dd>
|
240
|
+
<p>積を計算します。</p>
|
241
|
+
<p>例:</p>
|
242
|
+
<pre>M = Algebra.MatrixAlgebra(Integer, 2, 3)
|
243
|
+
N = Algebra.MatrixAlgebra(Integer, 3, 4)
|
244
|
+
L = M * N
|
245
|
+
a = M[[1, 2, 3], [4, 5, 6]]
|
246
|
+
b = N[[-3, -2, -1, 0], [1, 2, 3, 4], [5, 6, 7, 8]]
|
247
|
+
c = a * b
|
248
|
+
p c.type #=> L
|
249
|
+
c.display #=> [14, 20, 26, 32]
|
250
|
+
#=> [23, 38, 53, 68]</pre></dd>
|
251
|
+
<dt><a name="label-49" id="label-49"><code>**(<var>n</var>)</code></a></dt><!-- RDLabel: "**" -->
|
252
|
+
<dd>
|
253
|
+
<var>n</var> 乗を計算します。</dd>
|
254
|
+
<dt><a name="label-50" id="label-50"><code>/(<var>other</var>)</code></a></dt><!-- RDLabel: "/" -->
|
255
|
+
<dd>
|
256
|
+
商を計算します。</dd>
|
257
|
+
<dt><a name="label-51" id="label-51"><code>rank</code></a></dt><!-- RDLabel: "rank" -->
|
258
|
+
<dd>
|
259
|
+
階数を返します。</dd>
|
260
|
+
<dt><a name="label-52" id="label-52"><code>dsum(<var>other</var>)</code></a></dt><!-- RDLabel: "dsum" -->
|
261
|
+
<dd>
|
262
|
+
<p>直和を返します。</p>
|
263
|
+
<p>例:</p>
|
264
|
+
<pre>a = Algebra.MatrixAlgebra(Integer, 2, 3)[
|
265
|
+
[1, 2, 3],
|
266
|
+
[4, 5, 6]
|
267
|
+
]
|
268
|
+
b = Algebra.MatrixAlgebra(Integer, 3, 2)[
|
269
|
+
[-1, -2],
|
270
|
+
[-3, -4],
|
271
|
+
[-5, -6]
|
272
|
+
]
|
273
|
+
(a.dsum b).display #=> 1, 2, 3, 0, 0
|
274
|
+
#=> 4, 5, 6, 0, 0
|
275
|
+
#=> 0, 0, 0, -1, -2
|
276
|
+
#=> 0, 0, 0, -3, -4
|
277
|
+
#=> 0, 0, 0, -5, -6</pre></dd>
|
278
|
+
<dt><a name="label-53" id="label-53"><code>convert_to(<var>ring</var>)</code></a></dt><!-- RDLabel: "convert_to" -->
|
279
|
+
<dd>
|
280
|
+
<p><var>self</var> の各成分を行列環 <var>ring</var> にコンバートします。</p>
|
281
|
+
<p>Example:</p>
|
282
|
+
<pre>require "matrix-algebra"
|
283
|
+
require "residue-class-ring"
|
284
|
+
Z3 = Algebra.ResidueClassRing(Integer, 3)
|
285
|
+
a = Algebra.MatrixAlgebra(Integer, 2, 3)[
|
286
|
+
[1, 2, 3],
|
287
|
+
[4, 5, 6]
|
288
|
+
]
|
289
|
+
a.convert_to(Algebra.MatrixAlgebra(Z3, 2, 3)).display
|
290
|
+
#=> 1, 2, 0
|
291
|
+
#=> 1, 2, 0</pre></dd>
|
292
|
+
<dt><a name="label-54" id="label-54"><code>to_ary</code></a></dt><!-- RDLabel: "to_ary" -->
|
293
|
+
<dd>
|
294
|
+
<var>to_a</var> を返します。<var>to_a</var> は <var>Enumerable</var> で定義されています。</dd>
|
295
|
+
<dt><a name="label-55" id="label-55"><code>flatten</code></a></dt><!-- RDLabel: "flatten" -->
|
296
|
+
<dd>
|
297
|
+
<var>to_a.flatten</var> を返します。</dd>
|
298
|
+
<dt><a name="label-56" id="label-56"><code>diag</code></a></dt><!-- RDLabel: "diag" -->
|
299
|
+
<dd>
|
300
|
+
対角成分を配列で返します。</dd>
|
301
|
+
<dt><a name="label-57" id="label-57"><code>transpose</code></a></dt><!-- RDLabel: "transpose" -->
|
302
|
+
<dd>
|
303
|
+
<p>転置行列を返します。</p>
|
304
|
+
<p>例:</p>
|
305
|
+
<pre>M = Algebra.MatrixAlgebra(Integer, 2, 3)
|
306
|
+
a = M.new([[1, 2, 3], [4, 5, 6]])
|
307
|
+
Mt = M.transpose
|
308
|
+
b = a.transpose
|
309
|
+
p b.type #=> Mt
|
310
|
+
b.display #=> [1, 4]
|
311
|
+
#=> [2, 5]
|
312
|
+
#=> [3, 6]</pre></dd>
|
313
|
+
<dt><a name="label-58" id="label-58"><code>dup</code></a></dt><!-- RDLabel: "dup" -->
|
314
|
+
<dd>
|
315
|
+
<p>複製します。</p>
|
316
|
+
<p>例:</p>
|
317
|
+
<pre>M = Algebra.MatrixAlgebra(Integer, 2, 3)
|
318
|
+
a = M.new([[1, 2, 3], [4, 5, 6]])
|
319
|
+
b = a.dup
|
320
|
+
b[1, 1] = 50
|
321
|
+
a.display #=> [1, 2, 3]
|
322
|
+
#=> [4, 5, 6]
|
323
|
+
b.display #=> [1, 2, 3]
|
324
|
+
#=> [4, 50, 6]</pre></dd>
|
325
|
+
<dt><a name="label-59" id="label-59"><code>display([<var>out</var>])</code></a></dt><!-- RDLabel: "display" -->
|
326
|
+
<dd>
|
327
|
+
行列を <var>out</var> に表示します。<var>out</var> が省略されると <var>$stdout</var>
|
328
|
+
に表示します。</dd>
|
329
|
+
</dl>
|
330
|
+
<h1><a name="label-60" id="label-60">Algebra::Vector</a></h1><!-- RDLabel: "Algebra::Vector" -->
|
331
|
+
<p><em>(縦ベクトルクラス)</em></p>
|
332
|
+
<p>ベクトルのクラスです。</p>
|
333
|
+
<h2><a name="label-61" id="label-61">スーパークラス:</a></h2><!-- RDLabel: "スーパークラス:" -->
|
334
|
+
<ul>
|
335
|
+
<li><var>Algebra::MatrixAlgebra</var></li>
|
336
|
+
</ul>
|
337
|
+
<h2><a name="label-62" id="label-62">インクルードしているモジュール:</a></h2><!-- RDLabel: "インクルードしているモジュール:" -->
|
338
|
+
<p>なし</p>
|
339
|
+
<h2><a name="label-63" id="label-63">関連する関数:</a></h2><!-- RDLabel: "関連する関数:" -->
|
340
|
+
<dl>
|
341
|
+
<dt><a name="label-64" id="label-64"><code>Algebra.Vector(<var>ring</var>, <var>n</var>)</code></a></dt><!-- RDLabel: "Algebra.Vector" -->
|
342
|
+
<dd>
|
343
|
+
<a href="#label-66">Algebra::Vector::create</a>(ring, n) と同じです。</dd>
|
344
|
+
</dl>
|
345
|
+
<h2><a name="label-65" id="label-65">クラスメソッド:</a></h2><!-- RDLabel: "クラスメソッド:" -->
|
346
|
+
<dl>
|
347
|
+
<dt><a name="label-66" id="label-66"><code>Algebra::Vector::create(<var>ring</var>, <var>n</var>)</code></a></dt><!-- RDLabel: "Algebra::Vector::create" -->
|
348
|
+
<dd>
|
349
|
+
<p>環 <var>ring</var> を要素とする, <var>n</var> 次元のベクトル(縦ベクトル)
|
350
|
+
表現するクラスを生成します。</p>
|
351
|
+
<p>このメソッドの戻り値は <a href="#label-60">Algebra::Vector</a> クラスのサブクラス
|
352
|
+
です。このサブクラスにはクラスメソッドとして <var>ground</var> と
|
353
|
+
<var>size</var> が定義され、それぞれ、基底となる環 <var>ring</var>、
|
354
|
+
サイズ <var>n</var> を返します。</p>
|
355
|
+
<p>実際にベクトルを作るにはクラスメソッド <var>new</var>,
|
356
|
+
<var>matrix</var>,
|
357
|
+
<var>[]</var> を使います。</p>
|
358
|
+
<p><a href="#label-60">Algebra::Vector</a> は <var>n</var> 行 1 列の
|
359
|
+
<a href="#label-0">Algebra::MatrixAlgebra</a> と同一視されます。</p></dd>
|
360
|
+
<dt><a name="label-67" id="label-67"><code>Algebra::Vector::new(<var>array</var>)</code></a></dt><!-- RDLabel: "Algebra::Vector::new" -->
|
361
|
+
<dd>
|
362
|
+
<p><var>array</var> を配列とするとき、それを要素とす
|
363
|
+
る縦ベクトルを返します。</p>
|
364
|
+
<p>例:</p>
|
365
|
+
<pre>V = Algebra.Vector(Integer, 3)
|
366
|
+
a = V.new([1, 2, 3])
|
367
|
+
a.display
|
368
|
+
#=> [1]
|
369
|
+
#=> [2]
|
370
|
+
#=> [3]</pre></dd>
|
371
|
+
<dt><a name="label-68" id="label-68"><code>Algebra::Vector::vector{|<var>i</var>| ... }</code></a></dt><!-- RDLabel: "Algebra::Vector::vector" -->
|
372
|
+
<dd>
|
373
|
+
<p>第 <var>i</var> 成分を ... にしたベクトルを返します。</p>
|
374
|
+
<p>例:</p>
|
375
|
+
<pre>V = Algebra.Vector(Integer, 3)
|
376
|
+
a = V.vector{|j| j + 1}
|
377
|
+
a.display
|
378
|
+
#=> [1]
|
379
|
+
#=> [2]
|
380
|
+
#=> [3]</pre></dd>
|
381
|
+
<dt><a name="label-69" id="label-69"><code>Algebra::Vector::matrix{|<var>i</var>, <var>j</var>| ... }</code></a></dt><!-- RDLabel: "Algebra::Vector::matrix" -->
|
382
|
+
<dd>
|
383
|
+
第 <var>i</var> 成分を ... にしたベクトルを返します。
|
384
|
+
<var>j</var> には常に 0 が代入されます。</dd>
|
385
|
+
</dl>
|
386
|
+
<h2><a name="label-70" id="label-70">メソッド:</a></h2><!-- RDLabel: "メソッド:" -->
|
387
|
+
<dl>
|
388
|
+
<dt><a name="label-71" id="label-71"><code>size</code></a></dt><!-- RDLabel: "size" -->
|
389
|
+
<dd>
|
390
|
+
次元を返します。</dd>
|
391
|
+
<dt><a name="label-72" id="label-72"><code>to_a</code></a></dt><!-- RDLabel: "to_a" -->
|
392
|
+
<dd>
|
393
|
+
各成分を配列にして返します。</dd>
|
394
|
+
<dt><a name="label-73" id="label-73"><code>transpose</code></a></dt><!-- RDLabel: "transpose" -->
|
395
|
+
<dd>
|
396
|
+
横ベクトル <a href="#label-78">Algebra::Covector</a> に転置します。</dd>
|
397
|
+
<dt><a name="label-74" id="label-74"><code>inner_product(<var>other</var>)</code></a></dt><!-- RDLabel: "inner_product" -->
|
398
|
+
<dd>
|
399
|
+
<var>other</var> との内積を返します。</dd>
|
400
|
+
<dt><a name="label-75" id="label-75"><code>inner_product_complex(<var>other</var>)</code></a></dt><!-- RDLabel: "inner_product_complex" -->
|
401
|
+
<dd>
|
402
|
+
<var>other</var> との内積を返します。
|
403
|
+
<code>inner_product(other.conjugate)</code>と同じです。</dd>
|
404
|
+
<dt><a name="label-76" id="label-76"><code>norm2</code></a></dt><!-- RDLabel: "norm2" -->
|
405
|
+
<dd>
|
406
|
+
ノルムを返します。
|
407
|
+
<code>inner_product(self)</code>と同じです。</dd>
|
408
|
+
<dt><a name="label-77" id="label-77"><code>norm2_complex</code></a></dt><!-- RDLabel: "norm2_complex" -->
|
409
|
+
<dd>
|
410
|
+
ノルムを返します。
|
411
|
+
<code>inner_product(self.conjugate)</code>と同じです。</dd>
|
412
|
+
</dl>
|
413
|
+
<h1><a name="label-78" id="label-78">Algebra::Covector</a></h1><!-- RDLabel: "Algebra::Covector" -->
|
414
|
+
<p><em>(横ベクトルクラス)</em></p>
|
415
|
+
<p>ベクトルのクラスです。</p>
|
416
|
+
<h2><a name="label-79" id="label-79">スーパークラス:</a></h2><!-- RDLabel: "スーパークラス:" -->
|
417
|
+
<ul>
|
418
|
+
<li><var>Algebra::MatrixAlgebra</var></li>
|
419
|
+
</ul>
|
420
|
+
<h2><a name="label-80" id="label-80">インクルードしているモジュール:</a></h2><!-- RDLabel: "インクルードしているモジュール:" -->
|
421
|
+
<p>なし</p>
|
422
|
+
<h2><a name="label-81" id="label-81">関連する関数:</a></h2><!-- RDLabel: "関連する関数:" -->
|
423
|
+
<dl>
|
424
|
+
<dt><a name="label-82" id="label-82"><code>Algebra.Covector(<var>ring</var>, <var>n</var>)</code></a></dt><!-- RDLabel: "Algebra.Covector" -->
|
425
|
+
<dd>
|
426
|
+
<a href="#label-84">Algebra::Covector::create</a> (ring n)と同じです。</dd>
|
427
|
+
</dl>
|
428
|
+
<h2><a name="label-83" id="label-83">クラスメソッド:</a></h2><!-- RDLabel: "クラスメソッド:" -->
|
429
|
+
<dl>
|
430
|
+
<dt><a name="label-84" id="label-84"><code>Algebra::Covector::create(<var>ring</var>, <var>n</var>)</code></a></dt><!-- RDLabel: "Algebra::Covector::create" -->
|
431
|
+
<dd>
|
432
|
+
<p>環 <var>ring</var> を要素とする, <var>n</var> 次元のベクトル(横ベクトル)
|
433
|
+
表現するクラスを生成します。</p>
|
434
|
+
<p>このメソッドの戻り値は <a href="#label-78">Algebra::Covector</a> クラスのサブクラス
|
435
|
+
です。このサブクラスにはクラスメソッドとして <var>ground</var> と
|
436
|
+
<var>size</var> が定義され、それぞれ、基底となる環 <var>ring</var>、
|
437
|
+
サイズ <var>n</var> を返します。</p>
|
438
|
+
<p>実際にベクトルを作るにはクラスメソッド <var>new</var>,
|
439
|
+
<var>matrix</var>,
|
440
|
+
<var>[]</var> を使います。</p>
|
441
|
+
<p><a href="#label-78">Algebra::Covector</a> は 1 行 <var>n</var> 列の <a href="#label-0">Algebra::MatrixAlgebra</a> と
|
442
|
+
同一視されます。</p></dd>
|
443
|
+
<dt><a name="label-85" id="label-85"><code>Algebra::Covector::new(<var>array</var>)</code></a></dt><!-- RDLabel: "Algebra::Covector::new" -->
|
444
|
+
<dd>
|
445
|
+
<p><var>array</var> を配列とするとき、それを要素とす
|
446
|
+
る横ベクトルを返します。</p>
|
447
|
+
<p>例:</p>
|
448
|
+
<pre>V = Algebra.Covector(Integer, 3)
|
449
|
+
a = V.new([1, 2, 3])
|
450
|
+
a.display
|
451
|
+
#=> [1, 2, 3]</pre></dd>
|
452
|
+
<dt><a name="label-86" id="label-86"><code>Algebra::Covector::covector{|<var>j</var>| ... }</code></a></dt><!-- RDLabel: "Algebra::Covector::covector" -->
|
453
|
+
<dd>
|
454
|
+
<p>第 j 成分を ... にしたベクトルを返します。</p>
|
455
|
+
<p>例:</p>
|
456
|
+
<pre>V = Algebra.Covector(Integer, 3)
|
457
|
+
a = V.covector{|j| j + 1}
|
458
|
+
a.display
|
459
|
+
#=> [1, 2, 3]</pre></dd>
|
460
|
+
<dt><a name="label-87" id="label-87"><code>Algebra::Covector::matrix{|<var>i</var>, <var>j</var>| ... }</code></a></dt><!-- RDLabel: "Algebra::Covector::matrix" -->
|
461
|
+
<dd>
|
462
|
+
第 j 成分を ... にしたベクトルを返します。i には常に 0 が代入されます。</dd>
|
463
|
+
</dl>
|
464
|
+
<h2><a name="label-88" id="label-88">メソッド:</a></h2><!-- RDLabel: "メソッド:" -->
|
465
|
+
<dl>
|
466
|
+
<dt><a name="label-89" id="label-89"><code>size</code></a></dt><!-- RDLabel: "size" -->
|
467
|
+
<dd>
|
468
|
+
次元を返します。</dd>
|
469
|
+
<dt><a name="label-90" id="label-90"><code>to_a</code></a></dt><!-- RDLabel: "to_a" -->
|
470
|
+
<dd>
|
471
|
+
各成分を配列にして返します。</dd>
|
472
|
+
<dt><a name="label-91" id="label-91"><code>transpose</code></a></dt><!-- RDLabel: "transpose" -->
|
473
|
+
<dd>
|
474
|
+
横ベクトル <a href="#label-60">Algebra::Vector</a> に転置します。</dd>
|
475
|
+
<dt><a name="label-92" id="label-92"><code>inner_product(<var>other</var>)</code></a></dt><!-- RDLabel: "inner_product" -->
|
476
|
+
<dd>
|
477
|
+
<var>other</var> との内積を返します。</dd>
|
478
|
+
<dt><a name="label-93" id="label-93"><code>inner_product_complex(<var>other</var>)</code></a></dt><!-- RDLabel: "inner_product_complex" -->
|
479
|
+
<dd>
|
480
|
+
<var>other</var> との内積を返します。
|
481
|
+
<code>inner_product(other.conjugate)</code>と同じです。</dd>
|
482
|
+
<dt><a name="label-94" id="label-94"><code>norm2</code></a></dt><!-- RDLabel: "norm2" -->
|
483
|
+
<dd>
|
484
|
+
ノルムを返します。
|
485
|
+
<code>inner_product(self)</code>と同じです。</dd>
|
486
|
+
<dt><a name="label-95" id="label-95"><code>norm2_complex</code></a></dt><!-- RDLabel: "norm2_complex" -->
|
487
|
+
<dd>
|
488
|
+
ノルムを返します。
|
489
|
+
<code>inner_product(self.conjugate)</code>と同じです。</dd>
|
490
|
+
</dl>
|
491
|
+
<h1><a name="label-96" id="label-96">Algebra::SquareMatrix</a></h1><!-- RDLabel: "Algebra::SquareMatrix" -->
|
492
|
+
<p><em>(正方行列環クラス)</em></p>
|
493
|
+
<p>正方行列の作る環を表現するクラスです。</p>
|
494
|
+
<h2><a name="label-97" id="label-97">スーパークラス:</a></h2><!-- RDLabel: "スーパークラス:" -->
|
495
|
+
<ul>
|
496
|
+
<li><var>Algebra::MatrixAlgebra</var></li>
|
497
|
+
</ul>
|
498
|
+
<h2><a name="label-98" id="label-98">インクルードしているモジュール:</a></h2><!-- RDLabel: "インクルードしているモジュール:" -->
|
499
|
+
<p>なし</p>
|
500
|
+
<h2><a name="label-99" id="label-99">関連する関数:</a></h2><!-- RDLabel: "関連する関数:" -->
|
501
|
+
<dl>
|
502
|
+
<dt><a name="label-100" id="label-100"><code>Algebra.SquareMatrix(<var>ring</var>, <var>size</var>)</code></a></dt><!-- RDLabel: "Algebra.SquareMatrix" -->
|
503
|
+
<dd>
|
504
|
+
<a href="#label-102">Algebra::SquareMatrix::create</a>(ring, n) と同じです。</dd>
|
505
|
+
</dl>
|
506
|
+
<h2><a name="label-101" id="label-101">クラスメソッド:</a></h2><!-- RDLabel: "クラスメソッド:" -->
|
507
|
+
<dl>
|
508
|
+
<dt><a name="label-102" id="label-102"><code>Algebra::SquareMatrix::create(<var>ring</var>, <var>n</var>)</code></a></dt><!-- RDLabel: "Algebra::SquareMatrix::create" -->
|
509
|
+
<dd>
|
510
|
+
<p>正方行列表現するクラスを生成します。</p>
|
511
|
+
<p>このメソッドの戻り値は
|
512
|
+
<a href="#label-96">Algebra::SquareMatrix</a> クラスのサブクラス
|
513
|
+
です。このサブクラスにはクラスメソッドとして
|
514
|
+
<var>ground</var> と
|
515
|
+
<var>size</var> が定義され、それぞれ、基底となる環 <var>ring</var>、
|
516
|
+
サイズ <var>n</var> を返します。</p>
|
517
|
+
<p>SquareMatrix は <var>n</var> 行 <var>n</var> 列の Algebra::MatrixAlgebra と同一視されます。</p>
|
518
|
+
<p>実際に行列を作るにはクラスメソッド <var>new</var>,
|
519
|
+
<var>matrix</var>,
|
520
|
+
<var>[]</var> を使います。</p></dd>
|
521
|
+
<dt><a name="label-103" id="label-103"><code>Algebra::SquareMatrix.determinant(<var>aa</var>)</code></a></dt><!-- RDLabel: "Algebra::SquareMatrix.determinant" -->
|
522
|
+
<dd>
|
523
|
+
配列の配列 <var>aa</var> の行列式を返します。</dd>
|
524
|
+
<dt><a name="label-104" id="label-104"><code>Algebra::SquareMatrix.det(<var>aa</var>)</code></a></dt><!-- RDLabel: "Algebra::SquareMatrix.det" -->
|
525
|
+
<dd>
|
526
|
+
<!-- Reference, RDLabel "Algebra::SquareMatrix.determinat" doesn't exist --><em class="label-not-found">Algebra::SquareMatrix.determinat</em><!-- Reference end -->と同じです。</dd>
|
527
|
+
<dt><a name="label-105" id="label-105"><code>Algebra::SquareMatrix::unity</code></a></dt><!-- RDLabel: "Algebra::SquareMatrix::unity" -->
|
528
|
+
<dd>
|
529
|
+
単位行列を返します。</dd>
|
530
|
+
<dt><a name="label-106" id="label-106"><code>Algebra::SquareMatrix::zero</code></a></dt><!-- RDLabel: "Algebra::SquareMatrix::zero" -->
|
531
|
+
<dd>
|
532
|
+
零行列を返します。</dd>
|
533
|
+
<dt><a name="label-107" id="label-107"><code>Algebra::SquareMatrix.const(<var>x</var>)</code></a></dt><!-- RDLabel: "Algebra::SquareMatrix.const" -->
|
534
|
+
<dd>
|
535
|
+
成分が <var>x</var> のスカラー行列を返します。</dd>
|
536
|
+
</dl>
|
537
|
+
<h2><a name="label-108" id="label-108">メソッド</a></h2><!-- RDLabel: "メソッド" -->
|
538
|
+
<dl>
|
539
|
+
<dt><a name="label-109" id="label-109"><code>size</code></a></dt><!-- RDLabel: "size" -->
|
540
|
+
<dd>
|
541
|
+
サイズを返します。</dd>
|
542
|
+
<dt><a name="label-110" id="label-110"><code>const(<var>x</var>)</code></a></dt><!-- RDLabel: "const" -->
|
543
|
+
<dd>
|
544
|
+
成分が <var>x</var> のスカラー行列を返します。</dd>
|
545
|
+
<dt><a name="label-111" id="label-111"><code>determinant</code></a></dt><!-- RDLabel: "determinant" -->
|
546
|
+
<dd>
|
547
|
+
行列式を返します。</dd>
|
548
|
+
<dt><a name="label-112" id="label-112"><code>inverse</code></a></dt><!-- RDLabel: "inverse" -->
|
549
|
+
<dd>
|
550
|
+
逆行列を返します。</dd>
|
551
|
+
<dt><a name="label-113" id="label-113"><code>/(<var>other</var>)</code></a></dt><!-- RDLabel: "/" -->
|
552
|
+
<dd>
|
553
|
+
<code>self * other.inverse</code> を返します。<var>other</var> がスカラーなら
|
554
|
+
各要素を <var>other</var> で割ります。</dd>
|
555
|
+
<dt><a name="label-114" id="label-114"><code>char_polynomial(<var>ring</var>)</code></a></dt><!-- RDLabel: "char_polynomial" -->
|
556
|
+
<dd>
|
557
|
+
<var>ring</var> に多項式環を与えると、特性多項式を返します。</dd>
|
558
|
+
<dt><a name="label-115" id="label-115"><code>char_matrix(<var>ring</var>)</code></a></dt><!-- RDLabel: "char_matrix" -->
|
559
|
+
<dd>
|
560
|
+
<var>ring</var> に多項式環を与えると、特性行列を返します。</dd>
|
561
|
+
<dt><a name="label-116" id="label-116"><code>_char_matrix(<var>poly_ring_matrix</var>)</code></a></dt><!-- RDLabel: "_char_matrix" -->
|
562
|
+
<dd>
|
563
|
+
<var>poly_ring_matrix</var> に多項式成分の行列環を与えると、特性行列を返します。</dd>
|
564
|
+
</dl>
|
565
|
+
<h1><a name="label-117" id="label-117">Algebra::GaussianElimination</a></h1><!-- RDLabel: "Algebra::GaussianElimination" -->
|
566
|
+
<p><em>(ガウスの消去法モジュール)</em></p>
|
567
|
+
<p>ガウスの掃き出し法を実現するモジュールです。</p>
|
568
|
+
<h2><a name="label-118" id="label-118">ファイル名:</a></h2><!-- RDLabel: "ファイル名:" -->
|
569
|
+
<p><var>gaussian-elimination.rb</var></p>
|
570
|
+
<h2><a name="label-119" id="label-119">インクルードしているモジュール:</a></h2><!-- RDLabel: "インクルードしているモジュール:" -->
|
571
|
+
<p>なし</p>
|
572
|
+
<h2><a name="label-120" id="label-120">クラスメソッド:</a></h2><!-- RDLabel: "クラスメソッド:" -->
|
573
|
+
<p>なし</p>
|
574
|
+
<h2><a name="label-121" id="label-121">メソッド</a></h2><!-- RDLabel: "メソッド" -->
|
575
|
+
<dl>
|
576
|
+
<dt><a name="label-122" id="label-122"><code>swap_r!(<var>i</var>, <var>j</var>)</code></a></dt><!-- RDLabel: "swap_r!" -->
|
577
|
+
<dd>
|
578
|
+
<var>i</var> 行と <var>j</var> 行を入れ換えます。</dd>
|
579
|
+
<dt><a name="label-123" id="label-123"><code>swap_r(<var>i</var>, <var>j</var>)</code></a></dt><!-- RDLabel: "swap_r" -->
|
580
|
+
<dd>
|
581
|
+
<var>i</var> 行と <var>j</var> 行を入れ換えたものを返します。</dd>
|
582
|
+
<dt><a name="label-124" id="label-124"><code>swap_c!(<var>i</var>, <var>j</var>)</code></a></dt><!-- RDLabel: "swap_c!" -->
|
583
|
+
<dd>
|
584
|
+
<var>i</var> 列と <var>j</var> 列を入れ換えます。</dd>
|
585
|
+
<dt><a name="label-125" id="label-125"><code>swap_c(<var>i</var>, <var>j</var>)</code></a></dt><!-- RDLabel: "swap_c" -->
|
586
|
+
<dd>
|
587
|
+
<var>i</var> 列と <var>j</var> 列を入れ換えたものを返します。</dd>
|
588
|
+
<dt><a name="label-126" id="label-126"><code>multiply_r!(<var>i</var>, <var>c</var>)</code></a></dt><!-- RDLabel: "multiply_r!" -->
|
589
|
+
<dd>
|
590
|
+
<var>i</var> 行目を <var>c</var> 倍します。</dd>
|
591
|
+
<dt><a name="label-127" id="label-127"><code>multiply_r(<var>i</var>, <var>c</var>)</code></a></dt><!-- RDLabel: "multiply_r" -->
|
592
|
+
<dd>
|
593
|
+
<var>i</var> 行目を <var>c</var> 倍したものを返します。</dd>
|
594
|
+
<dt><a name="label-128" id="label-128"><code>multiply_c!(<var>j</var>, <var>c</var>)</code></a></dt><!-- RDLabel: "multiply_c!" -->
|
595
|
+
<dd>
|
596
|
+
<var>j</var> 列目を <var>c</var> 倍します。</dd>
|
597
|
+
<dt><a name="label-129" id="label-129"><code>multiply_c(<var>j</var>, <var>c</var>)</code></a></dt><!-- RDLabel: "multiply_c" -->
|
598
|
+
<dd>
|
599
|
+
<var>j</var> 列目を <var>c</var> 倍したものを返します。</dd>
|
600
|
+
<dt><a name="label-130" id="label-130"><code>divide_r!(<var>i</var>, <var>c</var>)</code></a></dt><!-- RDLabel: "divide_r!" -->
|
601
|
+
<dd>
|
602
|
+
<var>i</var> 行目を <var>c</var> で割ります。</dd>
|
603
|
+
<dt><a name="label-131" id="label-131"><code>divide_r(<var>i</var>, <var>c</var>)</code></a></dt><!-- RDLabel: "divide_r" -->
|
604
|
+
<dd>
|
605
|
+
<var>i</var> 行目を <var>c</var> 割ったものを返します。</dd>
|
606
|
+
<dt><a name="label-132" id="label-132"><code>divide_c!(<var>j</var>, <var>c</var>)</code></a></dt><!-- RDLabel: "divide_c!" -->
|
607
|
+
<dd>
|
608
|
+
<var>j</var> 列目を <var>c</var> で割ります。</dd>
|
609
|
+
<dt><a name="label-133" id="label-133"><code>divide_c(<var>j</var>, <var>c</var>)</code></a></dt><!-- RDLabel: "divide_c" -->
|
610
|
+
<dd>
|
611
|
+
<var>j</var> 列目を <var>c</var> 割ったものを返します。</dd>
|
612
|
+
<dt><a name="label-134" id="label-134"><code>mix_r!(<var>i</var>, <var>j</var>, <var>c</var>)</code></a></dt><!-- RDLabel: "mix_r!" -->
|
613
|
+
<dd>
|
614
|
+
<var>i</var> 行目に <var>j</var> 行目の <var>c</var> 倍を足します。</dd>
|
615
|
+
<dt><a name="label-135" id="label-135"><code>mix_r(<var>i</var>, <var>j</var>, <var>c</var>)</code></a></dt><!-- RDLabel: "mix_r" -->
|
616
|
+
<dd>
|
617
|
+
<var>i</var> 行目に <var>j</var> 行目の <var>c</var> 倍を足したものを返します。</dd>
|
618
|
+
<dt><a name="label-136" id="label-136"><code>mix_c!(<var>i</var>, <var>j</var>, <var>c</var>)</code></a></dt><!-- RDLabel: "mix_c!" -->
|
619
|
+
<dd>
|
620
|
+
<var>i</var> 列目に <var>j</var> 列目の <var>c</var> 倍を足します。</dd>
|
621
|
+
<dt><a name="label-137" id="label-137"><code>mix_c(<var>i</var>, <var>j</var>, <var>c</var>)</code></a></dt><!-- RDLabel: "mix_c" -->
|
622
|
+
<dd>
|
623
|
+
<var>i</var> 列目に <var>j</var> 列目の <var>c</var> 倍を足したものを返します。</dd>
|
624
|
+
<dt><a name="label-138" id="label-138"><code>left_eliminate!</code></a></dt><!-- RDLabel: "left_eliminate!" -->
|
625
|
+
<dd>
|
626
|
+
<p>左からの基本変形で階段行列に変形します。</p>
|
627
|
+
<p>戻り値は、変形するのに使った正方行列の積とその正方行列の
|
628
|
+
行列式と階数の配列です。</p>
|
629
|
+
<p>例:</p>
|
630
|
+
<pre>require "matrix-algebra"
|
631
|
+
require "mathn"
|
632
|
+
class Rational < Numeric
|
633
|
+
def inspect; to_s; end
|
634
|
+
end
|
635
|
+
M = Algebra.MatrixAlgebra(Rational, 4, 3)
|
636
|
+
a = M.matrix{|i, j| i*10 + j}
|
637
|
+
b = a.dup
|
638
|
+
c, d, e = b.left_eliminate!
|
639
|
+
b.display #=> [1, 0, -1]
|
640
|
+
#=> [0, 1, 2]
|
641
|
+
#=> [0, 0, 0]
|
642
|
+
#=> [0, 0, 0]
|
643
|
+
c.display #=> [-11/10, 1/10, 0, 0]
|
644
|
+
#=> [1, 0, 0, 0]
|
645
|
+
#=> [1, -2, 1, 0]
|
646
|
+
#=> [2, -3, 0, 1]
|
647
|
+
p c*a == b#=> true
|
648
|
+
p d #=> 1/10
|
649
|
+
p e #=> 2</pre></dd>
|
650
|
+
<dt><a name="label-139" id="label-139"><code>left_inverse</code></a></dt><!-- RDLabel: "left_inverse" -->
|
651
|
+
<dd>
|
652
|
+
左からの基本変形による一般逆行列です。</dd>
|
653
|
+
<dt><a name="label-140" id="label-140"><code>left_sweep</code></a></dt><!-- RDLabel: "left_sweep" -->
|
654
|
+
<dd>
|
655
|
+
左からの基本変形で階段行列にして返します。</dd>
|
656
|
+
<dt><a name="label-141" id="label-141"><code>step_matrix?</code></a></dt><!-- RDLabel: "step_matrix?" -->
|
657
|
+
<dd>
|
658
|
+
階段行列であるとき、軸(pivot)の配列を返します。そうでないとき、nil
|
659
|
+
を返します。</dd>
|
660
|
+
<dt><a name="label-142" id="label-142"><code>kernel_basis</code></a></dt><!-- RDLabel: "kernel_basis" -->
|
661
|
+
<dd>
|
662
|
+
<p>右から掛けて零になるベクトルの空間の基底の配列を返します。
|
663
|
+
各基底は <a href="#label-60">Algebra::Vector</a> の要素です。</p>
|
664
|
+
<p>例:</p>
|
665
|
+
<pre>require "matrix-algebra"
|
666
|
+
require "mathn"
|
667
|
+
M = Algebra.MatrixAlgebra(Rational, 5, 4)
|
668
|
+
a = M.matrix{|i, j| i + j}
|
669
|
+
a.display #=>
|
670
|
+
#[0, 1, 2, 3]
|
671
|
+
#[1, 2, 3, 4]
|
672
|
+
#[2, 3, 4, 5]
|
673
|
+
#[3, 4, 5, 6]
|
674
|
+
#[4, 5, 6, 7]
|
675
|
+
a.kernel_basis.each do |v|
|
676
|
+
puts "a * #{v} = #{a * v}"
|
677
|
+
#=> a * [1, -2, 1, 0] = [0, 0, 0, 0, 0]
|
678
|
+
#=> a * [2, -3, 0, 1] = [0, 0, 0, 0, 0]
|
679
|
+
end</pre></dd>
|
680
|
+
<dt><a name="label-143" id="label-143"><code>determinant_by_elimination</code></a></dt><!-- RDLabel: "determinant_by_elimination" -->
|
681
|
+
<dd>
|
682
|
+
体上の正方行列の行列式を掃き出し法で求めます。</dd>
|
683
|
+
</dl>
|
684
|
+
|
685
|
+
</body>
|
686
|
+
</html>
|