No.7ベストアンサー
- 回答日時:
z の値がないのですが、x → y, y → z と考えてよろしいですよね。
ニュートン法で計算してみると確かに収束しないですね。
与えられたパラメータを入れてグラフを書いてみると解はx = B の近傍にありました。これは図を描かなくても f の両辺を C で割ると分かります。
f/C=A/C*x - y*((cos(x)-cos(B))/(x-B))/C-z*((sin(x)-sin(B))/(x-B))/C + 1
右辺の第二項と第三項はいずれも x = B の近くでゼロに近い値となりますし、A/C * B はほとんど -1です。
で、この関数を見てみると x = B を境にして心電図の波形のように一旦跳ね上がってから潜り込み再び跳ね上がるような形をしています。ちょうど解が存在する辺りに変極点を複数持っています。この場合残念ながらニュートン法ではうまく収束しません。
また、(cos(x)-cos(B))/(x-B) や (sin(x)-sin(B))/(x-B) は式の上ではx → B の極限で -sin(B) と cos(B) に収束するのですが、数値計算上はゼロで割ることになるので値が不定となります。
解の近くに変極点が複数あること、解の近くで関数の計算自体が困難であることから、この関数に対してニュートン法は適していないと思います。
これだけではあまり前向きのアドバイスにならないので、一言付け加えるとすると、もし解が x = B の近くにあることが前もって分かっているならば、x=Bの近傍で展開した式
f = A*x + y * sin(B) - z * cos(B) + C
を使うという手があります。この場合は手計算で解けちゃいますね。
ご回答ありがとうございます。
ご指摘のとおり、x→y、y→zです。失礼しました。
今回の多くのアドバイスによりニュートン法の性質を理解することができました。
今後の勉強に役立てて行きたいと思います。
長い間本当にありがとうございました。
No.6
- 回答日時:
関数の情報ありがとうございます。
試しにこちらでも計算してみますから、収束しない場合の A, B, y, z の具体的な値を一組教えてください。この回答への補足
お忙しいところ何度もありがとうございます。
関数を少しだけ訂正させていただきます。
f=A*x-y*((cos(x)-cos(B))/(x-B))-z*((sin(x)-sin(B))/(x-B))+C
で、第1項のxに係数をつけただけです。パラメータが、
A=100.0、B=1.22995195D-08、C=-1.24229252D-06
x=1.0、y=7.87169249D-09
の場合収束しませんでした。ですがパラメータが、
A=100.0、B=1.00699949D-08、C=-1.01710309D-06
x=1.0、y=6.44479672D-09
のときは収束しました。
手間をおかけして申し訳ないですが、よろしくお願いいたします。
No.5
- 回答日時:
>二分法で出した解を初期値としてニュートン法を適応してみましたが、やはり収束しませんでした。
それはなかなか大変ですねぇ。他にやるとしたら単純にx=x-hとせずにx=x-alpha*hとでもして0<alpha<=1の範囲でfの絶対値が小さくなるようなalphaの一次元探索をしてxを決めていく手があります。
>もともとこの関数はニュートン法では収束しない性質があるのでしょうか?もしまだ解決策がありましたらよろしくお願いいたします。
方程式g(x,f,df)を見てみないことには関数の性質がどうであるとか解決策があるのかどうかは言えません。よろしければここに式を書いてください。
この回答への補足
三度ご回答ありがとうございます。
関数g(x、f、df)は一般的に表すと、
f=x-y*((cos(x)-cos(A))/(x-A))-z*((sin(x)-sin(A))/(x-A))+B
で、少し簡略化しました。
A、Bは定数、ニュートン法で求めたいのはxです。yとzも本来は変数ですが、ここでは定数扱いです。
何度も何度も申し訳ありません。
No.4
- 回答日時:
>二分法で確認したところ、解くことができました。
>今後の応用も考えてニュートン法でも解けるようにし>たいと思っているので、何かいい方法ありましたらア>ドバイスお願いいたします。
具体的な関数の形が分からないと原因に関してやはり正確なことは言えないのですが、私ならば二分法を何回か用いて出したラフな精度の解をニュートン法の初期値として用いて精度を上げます。経験上ほぼ間違いなく収束します。
この回答への補足
再度ご回答ありがとうございます。
二分法で出した解を初期値としてニュートン法を適応してみましたが、やはり収束しませんでした。
もともとこの関数はニュートン法では収束しない性質があるのでしょうか?もしまだ解決策がありましたらよろしくお願いいたします。
No.3
- 回答日時:
解いている非線形方程式を書いていただけると一発で答えられます。
まずはNo.2の方がおっしゃってるように解があるかどうかを確認しましょう。そういう時には、収束は遅いですが(最近の計算機は十分速いのでよほど時間のかかる計算を走らせるのではない限り十分だと思いますが)、二分法を使えば便利です。
あるいは解が x=0 だったら収束しないですよ。
その時は、
if ((dabs(x)>=1.d0.AND.dabs(h/x)<1.d-14).OR.
(dabs(x)<1.d0.AND.dabs(h)<1.d-14))
then
return
endif
とでもしましょう。
この回答への補足
ご回答ありがとうございます。
解いている方程式はとても複雑で、実は連立非線形方程式を差分法を使って解いています。しかし今求めようとしている方程式g(x,f,df)の解が求まれば、他の方程式の解も求まってくる連立方程式になっています。(意味が分からなかったらすみません)そのため連立のニュートン法は使っていません。正確には使いましたがやはり収束しない所がありました。
二分法で確認したところ、解くことができました。
今後の応用も考えてニュートン法でも解けるようにしたいと思っているので、何かいい方法ありましたらアドバイスお願いいたします。
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
このQ&Aを見た人はこんなQ&Aも見ています
-
【大喜利】【投稿~9/7】 ロボットの住む世界で流行ってる罰ゲームとは?
ロボットの住む世界で流行ってる罰ゲームとは?
-
フォロワー20万人のアカウントであなたのあるあるを披露してみませんか?
あなたが普段思っている「これまだ誰も言ってなかったけど共感されるだろうな」というあるあるを教えてください
-
映画のエンドロール観る派?観ない派?
映画が終わった後、すぐに席を立って帰る方もちらほら見かけます。皆さんはエンドロールの最後まで観ていきますか?
-
海外旅行から帰ってきたら、まず何を食べる?
帰国して1番食べたくなるもの、食べたくなるだろうなと思うもの、皆さんはありますか?
-
天使と悪魔選手権
悪魔がこんなささやきをしていたら、天使のあなたはなんと言って止めますか?
-
ニュートン法の収束性について
数学
-
sinωTをTで積分。
数学
-
ステップ信号とは何でしょうか?
物理学
-
おすすめ情報
- ・漫画をレンタルでお得に読める!
- ・人生のプチ美学を教えてください!!
- ・10秒目をつむったら…
- ・あなたの習慣について教えてください!!
- ・牛、豚、鶏、どれか一つ食べられなくなるとしたら?
- ・【大喜利】【投稿~9/18】 おとぎ話『桃太郎』の知られざるエピソード
- ・街中で見かけて「グッときた人」の思い出
- ・「一気に最後まで読んだ」本、教えて下さい!
- ・幼稚園時代「何組」でしたか?
- ・激凹みから立ち直る方法
- ・1つだけ過去を変えられるとしたら?
- ・【あるあるbot連動企画】あるあるbotに投稿したけど採用されなかったあるある募集
- ・【あるあるbot連動企画】フォロワー20万人のアカウントであなたのあるあるを披露してみませんか?
- ・映画のエンドロール観る派?観ない派?
- ・海外旅行から帰ってきたら、まず何を食べる?
- ・誕生日にもらった意外なもの
- ・天使と悪魔選手権
- ・ちょっと先の未来クイズ第2問
- ・【大喜利】【投稿~9/7】 ロボットの住む世界で流行ってる罰ゲームとは?
- ・推しミネラルウォーターはありますか?
- ・都道府県穴埋めゲーム
- ・この人頭いいなと思ったエピソード
- ・準・究極の選択
このQ&Aを見た人がよく見るQ&A
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
∞/0って不定形ですか?∞ですか...
-
数列の極限について
-
数学の問題です
-
収束か発散かを示したいです。
-
シグマの問題なのですが。
-
数3の極限です。 0/1の極限は∞...
-
無限大の0乗は、1で正しいですか?
-
1/n^2と1/n^3の無限和の問題を...
-
f : ℝ→ℝ が微分可能で一様連続...
-
離散フーリエ変換での回転子計...
-
Σ_[n=1,∞]1/nは発散?
-
これなぜ最後の不定形が0に収束...
-
無限級数Σ(n=1~∞)(n/n^2+1)の...
-
広義積分の問題です
-
1-1/2+1/3-1/4+1/5-1/6+...
-
多変数関数の極限についてです...
-
定数aのn乗根の極限(n→∞)...
-
数列の収束・発散
-
収束に関する問題
-
極限値lim[n→∞](3^n/(2^n+n^2))...
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
∞/0って不定形ですか?∞ですか...
-
シグマの問題なのですが。
-
数列の極限について
-
数学の問題です
-
数3の極限です。 0/1の極限は∞...
-
極限の問題
-
ラプラス変換後のsの意味って何...
-
ラプラス変換の積分法則 の証明...
-
無限大の0乗は、1で正しいですか?
-
limの問題
-
1/n^2と1/n^3の無限和の問題を...
-
収束
-
ニュートン法で解が収束しない
-
”有界閉区間”という言葉
-
級数の極限
-
極限値lim[n→∞](3^n/(2^n+n^2))...
-
次の条件を満たす数列{an}の...
-
定数aのn乗根の極限(n→∞)...
-
無限級数Σ(n=1~∞)(n/n^2+1)の...
-
高校数学の初歩的な質問ですが(...
おすすめ情報