ルンゲクッタ法を用いたバネの単振動を解析する課題をやっています。以下の部分まで作ったのですが、
どうしても最後のほうがわかりません。
関数fとf2のreturn やメイン関数を変更すればいいとおもうのですが、どうやっても満足な結果が得られません。どなたか詳しい方教えてもらえない?
#include<stdio.h>
#include<math.h>
double f(double t,double y1)
{
return y1;
}
double f2(double t2,double y2)
{
double w=5.0,k=0.5;
return -w*w*f(t2,y2)-2.0*k*y2;
}
int main(void)
{
int count,bunkatu,npr;
double t_0,T;
double k1,k2,k3,k4;
double y,t,dt,dt_2,dt_6;
t_0 = 0.0;
T=10.0;
bunkatu=2000;
npr=1000;
y=1.0;
printf("\n# [t_0,T]=[%6.4lf,%6.4lf] N=%d\n",t_0,T,bunkatu);
t=t_0;
dt=(T-t_0)/(double)bunkatu;
dt_2=dt/2.0;
dt_6=dt/6.0;
for(count=0;count <=bunkatu;count++)
{
if(count %(bunkatu/npr)==0)
printf("\n %5.2lf %10.5f",t,y);
k1=f2(t,y);
k2=f(t+dt_2,y+k1*dt_2);
k3=f(t+dt_2,y+k2*dt_2);
k4=f(t+dt,y+k3*dt);
y+=(k1+2.0*k2+2.0*k3+k4)*dt_6;
t+=dt;
}
printf("\n");
return 0;
}
A 回答 (5件)
- 最新から表示
- 回答順に表示
No.4
- 回答日時:
本当に教えてほしいんでしょうかね。
配布されたものらしきプログラムを投げつけ
ルンゲクッタ法を用いたバネの単振動を解析する課題
としか説明せず
みんながよくわからないなりに必死に読んでくれてるにもかかわらず、補足要求に
たいしても
「初期値を決めて解析する方法」
としか言わず。
本当に教えてもらいたいなら
プログラムの注釈をこまかくかいたり、
どこでどのような処理をするプログラムか、説明を書いたり、
物理の様々な式を計算しているなら
最低限計算している式を書くのが常識です。
私も読みましたがルンゲクッタ法がよくわからないのでプログラムをどう改良したらよいかもよくわかりませんでした。
詳しい補足要求お願いします。
No.3
- 回答日時:
物理わからないなりに読んでみましたが、
まず、fとf2の二種類あるのがよくわかりません。
また、"単振動を解析"だけでは何を求めるのかわかりませんが、
仮に求めるのが、時刻tにおける位置yのグラフであるなら、
関数f()は時刻tにおける速度sを返すはずです。
return y1;
が速度でないということはわかりますが、どうでしょうか。
時刻tにおける速度sを求める計算式は、私にはわからないので
ご自身で数学的に求めてください^^。
No.2
- 回答日時:
これはプログラミングするとき全般に言えることですが、
「他人に見せるコードにはわかりやすくコメントしましょう」を守ってください(^^;
私は大学に入って間もないのでルンゲクッタをよく知りません。
それで式と照らし合わそうと努力したのですが、大変な苦労です。
どの変数がルンゲクッタの式の何に対応しているのか教えてください。
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
似たような質問が見つかりました
- 数学 単振り子とルンゲ・タック法 1 2022/07/15 00:05
- C言語・C++・C# c言語でユーザ関数を利用して複素数のべき乗と絶対値の数列を計算するプログラムが作りたいです。 3 2023/01/29 22:13
- C言語・C++・C# Cのdoubleの浮動小数点表示について 3 2023/04/17 13:14
- C言語・C++・C# C言語のエラーについて 2 2022/07/11 13:56
- C言語・C++・C# C 言語の Gauss Jordan 法について 2 2022/12/28 11:16
- C言語・C++・C# 10個の実数に対する降順ソート結果を出力するプログラムを作りたいのですが、以下のプログラムをどう直せ 1 2022/07/09 22:16
- C言語・C++・C# C言語のマクローリン展開ローラン展開のコードについて 3 2022/12/15 14:45
- C言語・C++・C# プログラミングの授業の課題です 1 2023/01/17 22:15
- C言語・C++・C# バイナリファイルをコピーするのにかかる時間を測りたいのですが実行するとFatel error:gli 2 2022/11/03 01:10
- C言語・C++・C# LU分解法のピボッティングについて(C言語/gcc-9) 3 2022/07/11 23:10
このQ&Aを見た人はこんなQ&Aも見ています
関連するカテゴリからQ&Aを探す
おすすめ情報
このQ&Aを見た人がよく見るQ&A
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
c言語で、繰り返し文の中で、0....
-
プログラムでの数字につく”f”の...
-
C言語 関数プロトタイプ宣言の...
-
doubleの変数にintとintの割り...
-
float型とdouble型の変数の違い...
-
2分法で方程式の複数の解を自...
-
C言語の型による処理速度の違い
-
C#イベント中の戻り値の設定の...
-
C 開放してるのにエラー(doubl...
-
difftime()について
-
型について
-
ボール同士の衝突
-
C言語 入力した数値の平均値の...
-
C言語
-
インデックスが配列の境界外です.
-
C言語のマクローリン展開ローラ...
-
C言語を実行すると-infが出てき...
-
関数におけるif文とreturn文に...
-
int とdoubleの比較
-
MATLABで画像のヒストグラムを...
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
プログラムでの数字につく”f”の...
-
doubleの変数にintとintの割り...
-
float型とdouble型の変数の違い...
-
c言語で、繰り返し文の中で、0....
-
C言語を実行すると-infが出てき...
-
C言語 関数プロトタイプ宣言の...
-
C 開放してるのにエラー(doubl...
-
C言語の型による処理速度の違い
-
至急です! マクロ定義で #defi...
-
関数におけるif文とreturn文に...
-
c言語のプログラミングについて...
-
2分法で方程式の複数の解を自...
-
-1.#IND00と出てしまうのですが...
-
doubleは常に%lfとするべきなのか
-
C言語のpow関数の不具合
-
C言語で-23乗を取り扱うには
-
C言語で台形公式を使った二重積...
-
Cで3乗根を求める方法
-
sin(x)の近似について
-
2次方程式の解を求めるプログ...
おすすめ情報