常微分方程式の数値解法at SIM常微分方程式の数値解法 - 暇つぶし2ch■コピペモード□スレを通常表示□オプションモード□このｽﾚｯﾄﾞのURL■項目テキスト71:名無しさん＠５周年 06/07/31 14:50:12 オイラー法 y[n+1]=y[n]+h*f(x[n], y[n]) hが時間差で ホイン法は、単純な常微分方程式なら k1=h*f(x[n], y[n]) k2=h*f(x[n] + h, y[n] + k1) y[n+1]=y[n] + (k1 + k2)/2 これを連立常微分方程式に拡張するのか f1から出て来た値をぶち込むのがl1 , l2 f2から出て来た値をぶち込むのがm1 , m2 f3は微分方程式じゃないから動かさなくてよろしいｗ z[n]=f3(x[n], y[n]) l1=h*f1(x[n], y[n], z[n]) l2=h*f1(x[n] + l1, y[n] + m1 , z[n]) x[n+1]=f1(x[n], y[n], z[n]) + (l1 + l2)/2 m1=h*f2(x[n], y[n], z[n]) m2=h*f2(x[n] + l1, y[n] + m1 , z[n]) y[n+1]=f2(x[n], y[n], z[n]) + (m1 + m2)/2 で、この問題だと時間差は1みたいだし、2回計算した平均値で 延々回してるだけじゃんかｗ 72:さとし 06/08/22 14:43:50 ブラジウスの公式 http://irws.eng.niigata-u.ac.jp/~chem/itou/ce/diffeq2.html F'''+(1/2)FF''=0 (境界条件:η=0 F=F'=0 および η=∞ F''=1) をルンゲクッタで解くとあるんですけど，実際に解くときにF''，F',Fはどのように計算すれば良いでしょうか？教えて下さいm(_ _)m 次ページ最新レス表示レスジャンプ類似スレ一覧スレッドの検索話題のニュースおまかせリストオプションしおりを挟むスレッドに書込スレッドの一覧暇つぶし2ch