No.6ベストアンサー
- 回答日時:
No.5のお礼に対する解答。
私の言う特殊状態というのは、計算を続けていくと定常状態に達し、さらにそのまま計算を続行すると特殊な状態に到達する、という意味でした。
今回のご返事によって、計算の続行ではなく、初期条件を変えて再計算するものと考え直すことにしました。
ご質問のプログラムは、「相性のいい」初期データを与えれば収束して定常状態に達し、「相性の悪い」初期データを与えると定常状態にならず、振動してしまう性質があるようですね。
予想通りの値で定常状態に達する事例があるということは、計算式のプログラム化については問題ないないでしょうが、この実験には他に隠れた問題…"観測点"という言葉で推測できるのですが、「境界値問題」が含まれていませんか?
これは、式を適用する際のパラメタ(初期値)に制限があり、それを超えるような条件を与えるときは結果が保証されないというものです。微分方程式を解く場合などでは、境界値問題に目を向ける必要があります。
いま少しプログラミング作業から離れて、「境界値問題」のキーワードで実験を見直してみるか、発散,収束に関係なくもっと観測点を増やし、解の傾向を分析されてみてはいかがでしょうか?
お礼が遅れてしまいすみません.
度々,回答していただきありがとうございます.
プログラムのミスとかそのような問題ではないようですので,境界値問題をキーワードに発散の原因を調べていきたい思います.
見ず知らずの私に対してここまで検討していただき本当にありがとうございました.
No.5
- 回答日時:
No.3です。
「定常状態を過ぎると発散」ということは、初期条件を変えるからですか?
「初期条件から出発→定常状態→特殊状態」という遷移の仕方であれば計算の方針は合っていると思います。
特殊状態に達するまでに何回の加減算,乗除算をやっているか見当がつきますか?
その回数は回答者側にとって参考になるので、できれば開示してください。
その他、余りにも大きな数と小さな数同士の加減算は、桁落ちどころか情報落ちという現象が発生して、加減算をしていないのと同じ事態になりますので、プログラムを工夫する必要がでてきます。私の推測ですが、目的の値が振動したりしていませんか?
発散なら、きれいに発散しているのですか?
2度も回答ありがとうございます.
媒質の値の設定により発散が起こらなかったりします.特殊状態というのは発散が起こる状態のことでしょうか.特殊状態に達するまでに何回の加減算,乗除算をやっているかは,式やその式に使う係数の数が多く,また係数にも加減算,乗除算を行っているので回数を数えるのがちょっと大変で今すぐはわかりません.すみません.発散の状況はある観測点において,値が収束せずに他の観測点よりも大きな値で振動してしまうという状況です.一応定常状態になるということは,理論的な考えはあっているけれども,プログラムの組み方が悪いということでしょうか.毎日プログラムをいじっていますが未だ解決せずです・・・(泣)
No.4
- 回答日時:
発散の意味をはっきりさせた方が良いと思いますが。
データ解析を行うプログラムだとしますと、解析的に答えが出るような擬似的なデータを作り、そのプログラムによる解析結果と理論計算とを比較すれば、プログラムが間違っていればすぐ判ります。
プログラムが正しくて、実測データの解析結果が予測と違うなら、理論(モデル)が間違っているのでしょう。
回答ありがとうございます.
電磁界解析プログラムで,簡単なモデルを考え,その理論値と比較したりしています.定常状態直後でその理論値と比較した所,一致しました.ただ,定常状態を過ぎてしばらく経ったあたりから発散が起きてしまうのです.発散の状況はある観測点において,値が収束せずに他の観測点よりも大きな値で振動してしまうという状況です.毎日プログラムをいじっていますが未だ解決せずです・・・(泣)
No.3
- 回答日時:
> 発散がおき,
・べき級数計算のプログラムで、収斂するはずが発散した。
・装置を用いてデータ収集したところ、理論的には収斂するはずが、
発散(または振動)した。
どちらとも取れます。
前者なら計算誤差の累積、後者なら雑音(ノイズ)や分解性能を疑ってかかりましょう。
回答ありがとうございます.
電磁界解析プログラムで,簡単なモデルを考え,その理論的と比較したりしています.定常状態直後でその理論値と比較した所,一致しました.ただ,定常状態を過ぎてしばらく経ったあたりから発散が起きてしまうのです.モデルの値を変えると発散したりしなかったりします.やはり計算誤差の累積の為なのでしょうか・・・
No.2
- 回答日時:
> 発散がおき,
・数値計算のプログラムでアルゴリズム的には収束して欲しいのに発散した。
・考えが発散し、集中できず、研究がはかどらない。
どちらとも取れます。
前者っぽいですが、どちらでしょう?
前者なら、
・モデルを単純化して実験し、理論値との比較から問題を切り分け、プログラムのバグを絞り込む。
・手計算やExcelなんかによる計算結果と比較、検討。
とか。
後者なら、大学の保健センターなんかが身近ですし、経験豊富で頼りになると思います。
回答ありがとうございます.
「・数値計算のプログラムでアルゴリズム的には収束して欲しいのに発散した。
・考えが発散し、集中できず、研究がはかどらない。
どちらとも取れます。
前者っぽいですが、どちらでしょう?」
・数値計算のプログラムでアルゴリズム的には収束して欲しいのに発散した。の方です.
電磁界解析プログラムで,簡単なモデルを考え,その理論的と比較したりしています.定常状態直後でその理論値と比較した所,一致しました.ただ,定常状態を過ぎてしばらく経ったあたりから発散が起きてしまうのです.今もプログラムとにらめっこしています.
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
似たような質問が見つかりました
- その他(自然科学) 科学技術計算の仕事について 2 2023/02/04 18:09
- 大学院 大学研究室選び 失敗 5 2022/04/29 18:50
- 教育・学術・研究 情報系の学部で卒業研究中の大学四回生です。 現在、卒業研究中なのですが不安に感じることがあるので質問 2 2022/09/09 00:19
- 大学院 社会人大学院へ行かれた方教えてください 3 2023/07/15 20:38
- その他(教育・科学・学問) 大学理系研究室における器材の貸し出しや取り決めについて 1 2022/12/01 22:04
- 大学院 共同研究相手にアイデアを盗まれたかもしれない件 1 2022/06/19 22:51
- 大学院 やっぱり大学院進学したいです。 工学部4年で研究にはあまり興味がなかったため推薦を断って、就活をして 2 2022/09/30 10:41
- 大学院 修士課程(社会人コース)行くべきか 6 2023/07/15 16:18
- その他(ニュース・社会制度・災害) 内閣府目標 人が身体、脳、空間、時間の制約から解放された社会を実現って・・出来るのでしょうか? 1 2022/08/13 23:43
- 大学受験 AIの研究者になるための進路 4 2023/01/30 00:14
関連するカテゴリからQ&Aを探す
おすすめ情報
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
65536は2の何乗なのでしょうか?
-
VBAで一時的にオーバーフローを...
-
変化させるセルが変化しない
-
再帰呼び出しの計算量
-
EXCELなどで「返す」という表現
-
VBでReplace
-
エクセル以外で麻雀の成績を管...
-
Delphiでバッチ処理を直ぐに書...
-
数値計算の高速化 (cos, sin, exp)
-
スレッド処理からダイアログを...
-
schemeです
-
2つの住所から直線距離を計算...
-
スライムがつぶれていく様子を...
-
prolog計算
-
for文である数の倍数になるまで...
-
Visual C++でdebugとreleaseで...
-
連立一次方程式を解くプログラム
-
プログラムの発散
-
fortranのプログラムについて
-
骨折リスク評価のFRAXについて...
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
65536は2の何乗なのでしょうか?
-
変化させるセルが変化しない
-
排他的論理和 BCC(水平パリテ...
-
VBAの再計算が反映されない件に...
-
VBAで関数をつくる
-
モジュラス103の計算とは何でし...
-
バッチファイルでウインドウを...
-
EXCELなどで「返す」という表現
-
数値計算の高速化 (cos, sin, exp)
-
傾いた四角形内の範囲の条件式
-
骨折リスク評価のFRAXについて...
-
Excel VBAにてFFT
-
matlab計算での進捗状況を知りたい
-
C言語についてです。 再帰を使...
-
C言語について 下の画像は do-w...
-
アドオン利率を実質年率に変換
-
エクセルで特定のセルのみを任...
-
電卓でmodの計算
-
y=(x^2 +3x+1)^4を微分の定義を...
-
「評価」「評価値」という言葉...
おすすめ情報