開平法

出典: フリー百科事典『ウィキペディア（Wikipedia）』

開平法（かいへいほう, extraction of square root) あるいは開平算（かいへいざん）とは、正の実数の平方根の小数による近似値を求めるアルゴリズムの一つである。

[編集] 数式による開平法

開平法というと筆算を用いた計算が多いが、ここでは原理的な事も含めて数式で説明する。十進法を用いるが、他の位取り記数法でも同様に計算できる。ここで述べるのと似たような方法で、立方根を求める開立法や、もっと一般に n 乗根を求める事も可能である。

単に筆算による手順を知りたいだけの場合は、筆算による開平法の計算例を参照するだけでよい。

[編集] 帰納的方法

正の実数 x に対し、その平方根 √x が次のように表現されているとする。

$\sqrt{x} = \sum_{k \le n} 10^k a_k = 10^n a_n + 10^{n-1} a_{n-1} + \cdots + 10^2 a_2 + 10 a_1 + a_0 + 10^{-1} a_{-1} + 10^{-2} a_{-2} + \cdots$

ここで、各 a_k は 0 から 9までの整数とする。

n は十分大きな整数で構わないが、大きな k に対して a_k = 0 が続くだけであまり意味は無いので、大抵は、 a_n ≠ 0 となるような最大の n すなわち、 10²ⁿ ≤ x < 10²⁽ⁿ⁺¹⁾ を満たすような n を選ぶ。

ある整数 m があって、 k > m となる a_k の値は分かっているとする。

$p_m = \sum_{m < k \le n} 10^{k-m-1} a_k = 10^{n-m-1} a_n + 10^{n-m-2} a_{n-1} + \cdots + 10^2 a_{m+3} + 10 a_{m+2} + a_{m+1}$

$b_m = \sum_{k \le m-1} 10^{k-m} a_k = 10^{-1} a_{m-1} + 10^{-2} a_{m-2} + 10^{-3} a_{m-3} + \cdots$

と置いて

√x = 10^m (10 p_m + a_m + b_m)

と表記する。 p_m , a_m は整数、b_m は小数である。また、 p_n = 0 である。

x = 10^2m (10 p_m + a_m + b_m)²

10^−2m x − 100 p_m² = (20 p_m + a_m + b_m) (a_m + b_m)

0 ≤ b_m < 1 なので、

(20 p_m + a_m) a_m ≤ 10^−2m x − 100 p_m² < (20 p_m + (a_m + 1)) (a_m + 1)

であり、この式を満たすような a_m を見つければよい。

この段階では、a_m の値を、地道に探す必要があるが、a_m は 0 から 9 までの整数のいずれかなので、多くても 10通りの値を考えればよい。

20 p_m a_m ≤ (20 p_m + a_m) a_m ≤ 10^−2m x − 100 p_m² < (20 p_m + (a_m + 1)) (a_m + 1) ≤ (20 p_m + 10) (a_m + 1)

p_m ≠ 0 であれば

a_m ≤ (10^−2m x − 100 p_m²)/(20 p_m) < (1 + (1/(2 p_m))) (a_m + 1)

なので、 (10^−2m x − 100 p_m²)/(20 p_m) という割り算によって、 a_m がどのくらいかという見当をつけることはできる。

正方形 ABCD の面積は 10^−2m x, 青い正方形の面積は 100 p_m² で、橙色と桃色の部分の面積の合計が (20 p_m + a_m) a_m = r_m である。 p_m の値は既に決まっていて、 a_m をどこまで大きくとれるのかが問題である。

このようにして a_m が求まれば、次は

p_m−1 = 10 p_m + a_m

b_m−1 = 10 b_m − a_m−1

なので

√x = 10^m−1 (10 p_m−1 + a_m−1 + b_m−1)

を考えることにより、a_m と同様に a_m−1 の値が求まる。

[編集] 計算の簡略化

計算を簡略化するために

q_m = 2p_m

r_m = (10 q_m + a_m) a_m

y_m = 10^−2m x − 100 p_m²

とおくと

q_m−1 = 10 q_m + 2a_m

y_m−1 = 10^−2(m−1) x − 100 p_m−1² = 100(10^−2m x − p_m−1²)

=100(10^−2m x − 100 p_m² − (20 p_m + a_m) a_m) =100(y_m − r_m)

となり、 a_m−1 を求めるための不等式は

(10 q_m−1 + a_m−1) a_m−1 ≤ y_m−1 < (10 q_m−1 + (a_m−1 + 1)) (a_m−1 + 1)

となる。

左辺は、r_m−1 に等しいが、 a_m−1 を探しにくくなるので、この不等式では r_m−1 に置き換えない。

この不等式は、両辺が整数なので実際の計算においては、 y_m−1 の小数部分は無視してよい。y_m の小数部分は 10^−2m x の小数部分に等しく、これを β_m とし、y_m の整数部分を z_{m とすれば}

y_m = z_m + β_m

z_m−1 = y_m−1 − β_m−1 = 100(z_m − r_m) + 100 β_m − β_{m − 1}

となる。最後の部分の 100 β_m − β_{m − 1} は、β_m の小数第一位、小数第二位を取り出し、それぞれ十の位、一の位とする演算である。

x_i を 0から 9までの整数として

$x = \sum_{i \le j} 10^i x_i$

と表されているとすると

100 β_m − β_{m − 1} = 10 x_2m−1 + x_2m−2

例えば β_m = 0.743 であれば、100 β_m − β_{m − 1} = 74 である。z_m を用いることにより m が大きいうちは、x の上の方の桁だけの計算で済む。

y_m−1 の代わりに z_m−1 を用いた不等式

(10 q_m−1 + a_m−1) a_m−1 ≤ z_m−1 < (10 q_m + (a_m−1 + 1)) (a_m−1 + 1)

を考えれば十分である。

[編集] 計算例

計算例として x = 5630738.132 の平方根を求める。ここで使う式は

$x = \sum_{i \le j} 10^i x_i$
$\sqrt{x} = \sum_{k \le n} 10^k a_k$
$q_m = 2 \sum_{m < k \le n} 10^{k-m-1} a_k = 10 q_{m+1} + 2 a_{m+1}$
r_m = (10 q_m + a_m) a_m
z_m−1 = 100(z_m − r_m) + 10 x_2m−1 + x_2m−2
(10 q_m + a_m) a_m ≤ z_m < (10 q_m + (a_m + 1)) (a_m + 1)

初期値として

q_n = 0
k>n のとき a_k = 0
z_n = 10 x_2n+1 + x_2n

まず、x の値から、 x₆ = 5, x₅ = 6, x₄ = 3, x₃ = 0, x₂ = 7, x₁ = 3, x₀ = 8, x₋₁ = 1, x₋₂ = 2, x₋₃ = 3 となる。これ以外の x_j の値は全て 0 とする。

10⁷ < x < 10⁸ なので、 n = 3 が、a_n ≠ 0 となるような最大の n であり、

$\sqrt{x} = \sum_{k \le 3} 10^k a_k = 10^3 a_3 + 10^2 a_2 + 10 a_1 + a_0 + 10^{-1} a_{-1} + 10^{-2} a_{-2} + \cdots$

の係数 a_k を調べていく。

初期値より、q₃ = 0, a₄ = 0, z₃ = x₇ + x₆ = 5

a₃² ≤ z₃ < (a₃ + 1)²

となるような a₃ を探すと、a₃ = 2 が見つかる。

r₃ = (10 q₃ + a₃) a₃ = 4

q₂ = 2 a₃ = 4

z₂ = 100(z₃ − r₃) + 10 x₅ + x₄ = 163

(10 q₂ + a₂) a₂ ≤ z₂ < (10 q₂ + (a₂ + 1)) (a₂ + 1)

となるような a₂ を探すと、a₂ = 3 が見つかる。

r₂ = (10 q₂ + a₂) a₂ = 129

q₁ = 10 q₂ + 2 a₂ = 46

z₁ = 100(z₂ − r₂) + 10 x₃ + x₂ = 3407

(10 q₁ + a₁) a₁ ≤ z₁ < (10 q₁ + (a₁ + 1)) (a₁ + 1)

となるような a₁ を探すと、a₁ = 7 が見つかる。

このような計算を繰り返すと、各変数の値は次の表のようになる。

m	q_m	a_m	r_m	z_m
3	0	2	4	5
2	4	3	129	163
1	46	7	3269	3407
0	474	2	9484	13838
−1	4744	9	427041	435413
−2	47458	1	474581	837220
−3	474582	7	33220789	36263900
−4	4745834	6	284750076	304311100
−5	47458352	4	1898334096	1956102400
…	…	…	…	…

a_m の列から

√5630738.132 ≒ 2372.91764 …

と分かる。

[編集] 筆算による開平法

開平法は、筆算を用いるとさらに簡単に計算できる。

[編集] 筆算を用いた略記

r_m と q_m−1 の式を比べると、形がとてもよく似ていることがわかる。

r_m = (10 q_m + a_m) a_m
q_m−1 = 10 q_m + 2 a_m = (10 q_m + a_m) + a_m

例えば q₋₁ = 4744, a₋₁ = 9 であれば

47449	47449
×　　9	＋　 9
r₋₁ = 427041	q₋₂ = 47458

というかけ算と足し算である。

この 2つの式を書いていくのは無駄になるため、これを 1つにまとめて

	47449
	9
r₋₁ = 427041	q₋₂ = 47458

と書き、r₋₁ の値は、z₋₂ の計算に用いるので他の所へ書くことになる。

さらに、q_m−1 を求める式は、 q_m を 1桁ずらして、 2 a_m を加えるという演算を繰り返し行うものである。

例えば q₃ = 0, a₃ = 2, a₂ = 3, a₁ = 7, a₀ = 2 であれば

2
2	← a₃ を縦に 2つ並べる。
43
3	← a₂ を右端に縦に 2つ並べる。
467
7	← a₁ を右端に縦に 2つ並べる。
4742
2	← a₀ を右端に縦に 2つ並べる。
4744

のように、足し算と a_m を並べる作業を、一連の筆算として繰り返すことで、記述を簡略化できる。

[編集] 計算例

この節では、以上の節で定義された変数も用いて開平法の筆算の手順を説明をするが、変数を用いた説明は飛ばしても差し支えないように書いたので、手順を知るためだけであれば、これまでの節を読む必要は無い。

計算例として x = 5630738.132 の平方根を求める。

まず、 x の値を書き、小数点を基準にして 2つずつの桁に区切る。

√5 63 07 38.13 2

説明の便宜上、分けた一つ一つのかたまりをブロックと呼ぶことにする。

これは、漸化式

z_m−1 = 100(z_m − r_m) + 10 x_2m−1 + x_2m−2

の最後の 2項において行われる x から 2桁ずつ取り出して最後尾に付加する操作を見やすくするために、分けている。

二乗しても（左端のブロックの） 5を超えない最大の整数を探すと 2であるので、右側の方に 2を縦に 2つ並べて書く。それを足し算の筆算と思って計算し、すぐ下に結果の 4を書く。この筆算を足し算ではなく、かけ算と思って計算した場合の結果（この場合は足し算と同じ 4）を 5の下に書く。左の筆算で引き算を行い、その結果の 1をすぐ下に記入する。

2	2
√5 63 07 38.13 2	2←かけ算が5を超えない最大の整数2
4←右の筆算をかけ算として計算	4 ←足し算の筆算として計算
1 ← 引き算の筆算として計算

z₃ = 5で、それを元に a₃ = 2 を探した。それから r₃ = 4 を計算し、左側の筆算に書き、引き算を行った。さらに q₂ = 4 を計算し、右側の筆算に記述した。

次のブロックをおろし、右の筆算のかけ算が 163を超えないように右端の数を探すと 3である。63のブロックの上に 3と書き、右の筆算に 3を縦に 2つ並べる。さらに、右の筆算のかけ算と足し算の結果を、左右の筆算に書き加え、左の筆算では引き算を行う。

2 3	2
√5 63 07 38.13 2	2
4 ↓次のブロックをおろす	43
1 63	3←かけ算が163を超えない最大の整数3
1 29←右の筆算をかけ算として計算	46 ←足し算の筆算と思って計算
34 ← 引き算の筆算として計算

z₂ = 163で、それを元に a₂ = 3 を探した。それから r₂ =129 を計算し、左側の筆算に書き、引き算を行った。さらに q₁ = 46 を計算し、右側の筆算に記述した。

このような計算を繰り返すことによって

2 3 7 2. 9 1 7 6 4	2
√5 63 07 38.13 2	2
4	43
1 63	3
1 29	467
34 07	7
32 69	4742
1 38 38	2
94 84	47449
43 54 13	9
42 70 41	474581
83 72 20	1
47 45 81	4745827
36 26 39 00	7
33 22 07 89	47458346
3 04 31 11 00	6
2 84 75 00 76	474583524
19 56 10 24 00	4
18 98 33 40 96	474583528
57 76 83 04

などと計算する。

左側の筆算において、おろしてくる数が無くなった場合は 0 をおろしてくる。小数点の位置は、 x の小数点の位置にそろえる。

以上により

√5630738.132 ≒ 2372.91764 …

と分かる。

得られた値を平方してみると

2372.91764² = 5630738.1262231696

2372.91765² = 5630738.1736815225

なので

2372.91764 < √5630738.132 < 2372.91765

であることが分かる。

[編集] 関連項目

ネイピアの骨：ネイピアの骨を使った場合の開平法

"http://ja.wikipedia.org../../../%E9%96%8B/%E5%B9%B3/%E6%B3%95/%E9%96%8B%E5%B9%B3%E6%B3%95.html" より作成

カテゴリ: 算術 | 初等数学 | 数学に関する記事

開平法

出典: フリー百科事典『ウィキペディア（Wikipedia）』

目次

[編集] 数式による開平法

[編集] 帰納的方法

[編集] 計算の簡略化

[編集] 計算例

[編集] 筆算による開平法

[編集] 筆算を用いた略記

[編集] 計算例

[編集] 関連項目

Views

ナビゲーション

ヘルプ

検索