January 04, 2004
人工衛星シミュレータ
学科の課題でつくりました。現在のところ、数値の羅列しか出力しないので、つまらないです。時間があったらDirectXやOpenGLを使ってビジュアルに表示できるようにしたいと思います。
ソースはここ
C++で書いたのですが、テンプレート機能に本当に助けられました。ルンゲクッタ積分がそのまま書けたのには涙もんです。
ルンゲクッタ積分[4次](strategy.hより)
template <class Function, class V1, class V2>
V2 nextByRK4(Function f, V1 x, V2 y, V1 h){
V2 k1, k2, k3, k4;
k1 = f(x, y) * h;
k2 = f(x + h/2, y + k1/2) * h;
k3 = f(x + h/2, y + k2/2) * h;
k4 = f(x + h, y + k3) * h;
return y + (k1 + k2*2 + k3*2 + k4)/6;
}
V2 nextByRK4(Function f, V1 x, V2 y, V1 h){
V2 k1, k2, k3, k4;
k1 = f(x, y) * h;
k2 = f(x + h/2, y + k1/2) * h;
k3 = f(x + h/2, y + k2/2) * h;
k4 = f(x + h, y + k3) * h;
return y + (k1 + k2*2 + k3*2 + k4)/6;
}
次の課題はカルマンフィルター(涙)。
コメント
コメントする
スポンサード リンク