ニュートン法による非線型方程式の解.pptVIP

ニュートン法による非線型方程式の解 ニュートン法 方程式f(x)=0の解を求める 処理手順 1. 初期値 x の決定 x = x0　を決める 2. 接線と x 軸の交点の計算 x = x0 における y = f(x) の接線を引き、今度はこの接線と y=0 (x軸)の交点を x1 とする すなわち、 xn+1 = xn - (f(x)/f ’(x)) を計算する 3. 2. を繰り返して値が収束したらそれを解とする ニュートン法 ニュートン法の弱点 ニュートン法の弱点 ニュートン法での収束条件 ニュートン法では現在の xn がどれだけ真の値に近いかは，一般には分からない 収束条件として、ある小さな正の数εに対して 　　　　　　　　　　　xn+1-xn 　　　　　　　　　　　　xn+1　 　　となった時点で計算を終了し xn+1 を解とする ニュートン法のプログラム 入力 初期値　x0 計算精度　ε 方程式　f(x) f(x)の導関数　f ’(x) 繰り返し回数の上限　number 出力 解(計算過程) ニュートン法の注意点 初期値をいくつにするか? 初期値の設定の際、あまりに解と掛け離れた値を与えると、収束するのに時間がかかったり、収束しなかったりする 収束条件をどうするか? どの程度の精度で計算するかを決定していないと，繰り返しをいつ終えるか決まらない 収束しない場合はどうするか? 虚数解は求まらない 例題　ニュートン法のプログラム 初期値，計算精度、繰り返し上限回数を読み込んで，f(x)=x2 –2 をニュートン法で解くプログラムを作成する f(x) = x2 - 2 g(x) = f ’(x) = 2x x : 現在のx new_x : 次のx g(x) が 10 以下なら重解とする #include stdio.h #include math.h void function(void); double f(double x); double g(double x); /* f(x)の出力 */ void function() { printf(f(x) = pow(x,2) - 2\n); } /* f(x) */ double f(double x) { return pow(x,2) - 2; } /* f(x)の導関数 */ double g(double x) { return 2 * x; } int main(void) { double x; double new_x; double eps; int number, i; char buf[100]; function(); /* f(x)を表示 */ printf(初期値 : ); fgets(buf, 100, stdin); sscanf(buf,%lf, x); printf(計算精度 : ); fgets(buf, 100, stdin); sscanf(buf,%lf, eps); printf(繰り返し上限回数 : ); fgets(buf,100,stdin); sscanf(buf,%d, number); printf(繰り返し回数\tnew_x\t\tf(x)\t\tg(x)\n); for(i = 0 ; i number ; i++) { new_x = x - f(x) / g(x); if(fabs(new_x - x) eps * fabs(new_x)) { printf(x = %lf\n,new_x); break; } printf(%2d\t\t%lf\t%lf\t%lf\n,i,new_x,f(x),g(x)); if(fabs(g(x)) 1.0e-4){ printf(x = %lf(重解)\n,new_x); break; } x = new_x; } if(i == number) printf(繰り返し上限\n); } 課題 * 牧之内研「C/C++プログラミング入門」Webページ http://www.db.is.kyushu-u.ac.jp/adp/ ?初期値x=a をとる グラフに接線１を引く 接線１が横軸をきる点x=b aよりも真の値に