数值计算入门.pdfVIP

1 数値計算入門 武尾英哉 1. 離散数学と数値計算 数学的解法の中には理論計算では求められないものもある．例えば，定積分は，まず は積分（被積分関数の原始関数をみつけること）できなければ値を得ることはできない． また，ある関数の所定の値における微分値を得るには，まずその関数の微分ができなけ ればならない．さらに代数方程式の解を得るためには，解析的に代数方程式を解く必要 がある．ところが，これらは必ずしも解析的に導けるとは限らない．そこで本章では， 計算機を使って，離散的な数値（とびとびの数値を扱う）を扱って，近似的に解を導く 手法について説明する．このことを離散数学 という． 本稿では，連続的な数値による数学（これまでの解析的数学）を，離散的な数値によ る数学として近似的に扱うことで解を求めることができる手法について説明する．この ように，離散数学を用いて種々の手法による近似計算により解を求めることを数値 計算 とよぶ． 2. 数値代数方程式（ニュートン法） 数値計算を使って代数方程式を解く代表的な手法であるニュートン法について説明す る． x 代数方程式を解く（f (x ) 0 を満足する を求める）ということは，グラフにおいて x x a 関数y f (x ) の 軸切片を求めることと等価である．すなわち，図 1 における を 求めることである．以下，ニュートン法による手順について説明する． ①まず初めに，予想される真の解に近いと思われる値をひとつとる（x の初期値）． x の初期値をx0 とする． ②次に，グラフの接線を考え，そのx 軸切片を計算する．図 1 において，x x0 にお y f (x )(x x ) f (x ) x ける接線の方程式は 0 0 0 であるので，この接線が 軸と交差す f (x ) 0 る点は，x x x がx0  のときである．この を新たに 1 とする． f (x0 ) f (x ) x1 x 0  0 ただし， f (x 0 ) 0 とする． f (x0 ) x この 軸切片の値は，予想される真の解により近いものとなるのが一般である． ③②の手法を繰り返すと，限りなくx a に近い値となる．すなわち，一般的にx x n 2 から次のステップのx xn1 を求める式は次のようになる． f (x n ) xn1 xn  ただし， f (x n ) 0 とする． (1) f (x n ) ④この繰り返し計算を希望する精度の解を得るまで行う．すなわち，実際には，次の