• 締切済み

「セカント法」（割線法）についてです.

非線形方程式の解を求めるアルゴリズム「セカント法」（割線法）についてです。 yahoo知恵袋の方でも質問しているのですが、 非線形方程式の解を求めるアルゴリズム「セカント法」（割線法）についてです。 セカント法で解を求める際、その収束次数は黄金比{（1+√5）/2}になると言われていますが、それを確認する方法はないものでしょうか。 セカント法で求めるプログラムはjavaの方で既に組んでおります。 なにかいいアイデアがあれば是非回答よろしくお願いします。 補足 組んだプログラムは与えられた関数の解と反復回数、それぞれの計算における関数値、真の解との誤差をかえすものです。 お願いします。

みんなの回答

  • masa2211
  • ベストアンサー率43% (178/411)
回答No.3

反復回数ごとの真の解との誤差が計算できているなら、あとは次数を計算するだけ。 http://www.cis.twcu.ac.jp/~osada/rikei/rikei2008-5.pdf の、２ページの(5)式 尤も、こちらでは、８倍精度で計算しているようです。

  • 178-tall
  • ベストアンサー率43% (762/1732)
回答No.2

何をお望みか、はかりかねますが、 　secant method convergence でネット検索すると、ゾロゾロ。 手ごろな分量に見える一例だけ、「参考 URL」へ。 　　　　　　　　　　　↓ 　The convergence of the secant method is superlinear 　　

参考URL：
http://homepages.math.uic.edu/~jan/mcs471/Lec5/secant.pdf
  • Tacosan
  • ベストアンサー率23% (3656/15482)
回答No.1

「アイデア」もなにも, 「収束次数」の定義通りに計算していくしかないんじゃないかな.

  1. カテゴリ
  2. 学問・教育
  3. 数学・算数

関連するQ&A

  • ニュートン法

    こんにちは。 aを解とするある非線形方程式に,aの近傍の値x0を初期値としてニュートン法を一反復だけ適用したところ、誤差は |x0-a|=10^-4から|x1-a|＝10^-7に減少した。このとき、ニュートン法をもう一反復適用したら誤差|x2-a|はどの程度になると推定されるのでしょうか？

  • 割線法とはさみうち法の収束性について

    割線法とはさみうち法の収束性について 球根アルゴリズムに割線法とはさみうち法がありますね。 これら二つのアルゴリズムについて質問です。 まず割線法は必ず根に収束するとは限りません。 割線法が根に収束する（またはしない）条件の導出法を教えてください。 次にはさみうち法ですが、はさみうち法にはいくつか亜種があると思います。 まず一つが、初期値a,bによってはさまれたcを求め、bをcに更新してループを回すやり方です。 このやり方ではaは固定されたままです。if f(a)・f(c)<0 then b=c else a=c という条件分岐はやりません。 この場合、収束性はどうなるでしょうか。 二つ目は、上記の条件分岐をやるやり方です。 このやり方は必ず収束するらしいのですがその証明を教えてください。 よろしくおねがいします。

  • ニュートン法でのプログラミング

    ニュートン法を用いてf(x)=X^2-2の正の解を求めなさい。 また、反復回数が増えるにしたがって解が収束していく様子を表示しなさい。 という問題なのですが、いくら反復回数を増やしても誤差がでてしまい、うまくいきません。このプログラミング解いていただきたいです。よろしくお願いします。

  • ごく限られた条件においてニュートン法で収束しません

    未知数が3つの非線形連立方程式をニュートン・ラフソン法で解いています。しかし、ごく限られた条件において、添付図の赤プロットように誤差が減っては増えるを繰り返し、収束しきい値に及びません。すぐ隣の条件では青プロットのように問題なく収束します。 本問題の解決策として、識者の方より、非線形方程式が不連続にならないように弱いばねを与えることを教わり、それによって収束しない範囲は激減しましたが、それでもごく一部残っている状況です。 なお、連立3元1次方程式の解法は掃き出し法で行っています（ピボット操作しています）。 解決策をご存じの方、また同様なご経験をお持ちの方がおられましたら、お知恵をお借りしたく、よろしくお願いいたします。

  • ニュートン法の問題がどうしてもわからないので教えてください。

    非線形方程式の解を求めるアルゴリズムのひとつに、ニュートン法と呼ばれる方式がある。この方法ではf(x)=0の解の近似値を次の漸化式から求める。 x_(i+1) =x_i - f(x_i)/f'(x_i) なぜこのアルゴリズムで解を得られるのかをf(x)=x^2 - e^x -2 =0 を例に説明せよ。 という問題です。 どなたかわかる方がいたらお願いします。

  • はさみうち法とニュートン法について

    はさみうち法とニュートン法のプログラムのついてのことなんですが、ひとつずつ解を求めるプログラムは作れました。例えば三次方程式だったら三つの解をすべて求めるプログラムをつくるにはどうすればいいんでしょうか？

  • 共役勾配法 (Conjugate Gradient)の失敗?

    共役勾配法で線形方程式を解く時に、誤差を10の-8乗 として計算をしています。 たいていの計算はこの誤差で解けるのですが、ある計算 では、誤差が10の-3乗くらいまで減った後に、誤差が増加 していってしまいます。 共役勾配法を使っていると、こういうことが起きるのでしょうか? 数値計算プログラムの作成ミスだろうかと考えてチェックしていますが、 今の所おかしな部分は見つかっておりません。

  • ニュートン法で解が収束しない

    こんにちは。 差分式で表した非線形方程式をニュートン法で解いています。が収束しな解あります。ニュートン法は初期値に依存しているため、初期値を可変的にしてみましたがダメでした。何かいい方法はないでしょうか？ 参考になるか分かりませんが、使っているプログラムのニュートン法の計算の一部は以下のようです。 call ｇ(x,f,df) h=f/df x=x-h if(dabs(h/x)<1.d-14) then 　return endif

  • 微分方程式

    すいません。 微分方程式 dx/dt=y dy/dt=t+0.1-x-(x-t)^{3} を解ける人いたら解いてもらえませんでしょうか？（数値解法でない方法で） 真解が知りたいです。 初期値は適当な値でいいので、t=0,x=4,y=1 でお願いします。 そもそも、これって真解求まるのかすらわかりません。 今、プログラムの課題で上記の微分方程式をおいらー法で解くプログラムを作ってるんですが、真解を求めて、その誤差を測りたいと思っている次第です。 よろしくお願いします。

  • 連立線形方程式の数値解

    Ａ＝｜－２　　０　－１　　８｜　ｂ＝｜８｜ 　　｜　０　　１　　０　　０｜　　　｜９｜ 　　｜　７　－３　－１　　０｜　　　｜２｜ 　　｜　０　－２　　４　　１｜　　　｜９｜ 上のＡ、ｂについて連立線形方程式Ａｘ＝ｂの数値解を解く問題で、この行列式のままだとヤコビやガウス＝ザイデル法では解けない（収束しない）ので、これを収束するように修正したいのですが色々試してみてもできず、ここで聞くに至りました。 収束させるためには具体的にどうしたらいいのか教えて下さい。よろしくお願いします。

注目のQ&A

カテゴリ

専門家に質問してみよう

職業から探して質問する

あなたにピッタリな商品が見つかる！ OKWAVEセレクト

質問する