情报科学-lecture.ecc.u.ppt

情報科学 第７回　数値解析(2) 数値解析(2) トピック 連立1次方程式の解 消去法 反復法 乱数とモンテカルロ法 乱数とは？ 疑似乱数 モンテカルロ法 連立1次方程式の数値解法 小規模な連立1次方程式の解法 消去法 Gauss消去法 Gauss-Jordan法 大規模な連立1次方程式の解法 反復法 Jacob法 消去法による解法 小規模な連立1次方程式は消去法によって解くことができる 消去法とは以下の演算を何回か行うことによって未知数を消去し、方程式をより簡単な方程式に変形して解を求める方法（筆算と同様） 1つの方程式に(0でない)ある数を掛ける 1つの方程式にある数を掛けて他の方程式に加える 多くの消去法では、未知数を1回に1個ずつ消去することによって係数行列を三角行列に変形し、後退置換によって未知数の値を逐次求める Gauss消去法 与えられた連立1次方程式を行列演算で表現し、係数行列を定義する Gauss消去法（続き） Gauss消去法（続き） Gauss-Jordan法 Gauss消去法において、係数行列を単位行列に変形する方法 後退置換が不必要である Pivoting 消去法では、係数行列の要素での割算における問題がある 割算の分母が0になる場合（係数行列の対角要素が0） ?行または変数を入れ替える必要がある 割算の分母が0に近い値になる場合 ?割算の結果非常に大きな数字が結果として表れると数値計算の精度が失われる（前回の講義参照） ?この場合も行または変数を入れ替える必要がある 反復法 Gauss消去法は、単純明快という利点があるが、未知数の個数が多い場合演算回数が非常に多くなるという欠点がある 個の連立方程式から、１つの未知数を消去するために に比例する演算が必要となる これを 個の未知数に対して繰り返すためには、合計 　　　に比例する演算が必要となる また、物理現象の解析などの応用においては、希薄な(sparse)係数行列を扱う場合や　　が大きい場合が多く　　　の演算回数は膨大になることがある 多くの未知変数を含む連立方程式の場合は、消去法よりも反復法が用いられることが多い Jacob法などは、行列のかけ算（　　の演算）を繰り返し、解を近似する 　回のかけ算で解が近似できる場合、計算量は となり、　　　　 の場合には反復法の方が有利となる Jacob法 Jacob法は、 適当な初期値ベクトルに、行列の演算を繰り返し、解ベクトルに収束させる方法 Jacob法（続き） 係数行列 を対角行列 と非対角行列 に分解する Jacob法（続き） ノイマン級数 Jacob法では、　　　　　　　　　という方程式を解く際に　　　　　　の逆行列として、無限級数　　を考えた Jacob法の例 Jacob法の例（続き） Jacob法の計算例 Jacob法の計算例（続き） def jacob(a,b) n = b.size x = Array.new(n,0) nx = Array.new(n,0) d = 1.0 while d0.000001 for i in 0..n-1 r = 0 for j in 0..n-1 if j!=i r += a[i][j]*x[j] end end nx[i] = (b[i]-r)/a[i][i] end d = dist(x,nx) print nx=, nx, d=, d, \n for i in 0..n-1 x[i] = nx[i] end end x end 2．乱数とモンテカルロ法 乱数とは？　 一様乱数 正規乱数 疑似乱数 モンテカルロ法 乱数とは？ 乱数とは、ある一つの数が現れたとき、それに続く数が前の数と全く関係なく現れるような列（乱数列）の要素 乱数列とは、ある分布に従う（互いに独立な事象を表す）確率変数の実現値の列をいう 乱数とは、乱数列の各要素である 一様乱数 出現確率が値によらず一定である乱数を一様乱数という（確率分布が一様である場合） 自然現象における乱数列 放射性原子は放射線を出して崩壊する （原子核が崩壊して別の核種になるか、励起状態の原子核 がより低いエネルギー状態へ遷移する） 放射性物質の単位時間あたりの崩壊数を計測する 正規乱数 出現確率が正規分布に従うような乱数を正規乱数という 確率変数　　が 　よりも小さい値をとる確率　　 　が以下の正規分布に従う場合 疑似乱数(pseudorandom number