問題タブ [runge-kutta]

このコードから b と c を印刷しようとしていますが、うまくいきません。私が正しければ、このコードは 0.05 のステップ サイズでいくつかのポイントを出力するはずですが、表示されません。このコードから 2 つの値を出力する方法を知っている人はいますか?

車のようなロボットの非ホロノミック モーションに 4 次のルンゲ クッタを実装しようとしています。何が間違っているのかわかりません。基本的に、+-Pi/4 を渡してハードな左折と右折を計算し、異なる軌道を取得しています。しかし、+pi/4 または -pi/4 を渡しても、同じ答えが得られます。私は自分が間違っていることを理解できません。私が使用している制約方程式は次のとおりです。

ここで、s は速度、L は車のようなロボットの長さです。

オイラー法とレンジ クッタ法の両方を使用してばね質量の問題を解決し、プロットを比較しようとしています。Euler と Runge-Kutta の両方の関数を作成しましたが、問題に対して関数を呼び出した後、プロットにデータが表示されないようです。プロットを修正して、コードにエラーがないかどうかを確認してください。ありがとう

ここに私が解決しようとしている問題があります

5 次ルンゲクッタ法と適応ステップサイズ法を使用して、一連の方程式を解きたいと考えています。Taner Akgun によって書かれた便利なコードを見つけました。コードは次のとおりです。

残念ながら、私は Fortran についてまったく詳しくありません。このコードを使用して、次の一連の方程式を解きます。

1. dy/dx=-x
2. dy/dx=-1

コード内では、nvar 変数は方程式の数であり、このコードでは 1 に設定されています。1 以外に変更したい場合は、どのようにコードを変更すればよいですか?

また、すべての x と y の値を出力ファイルに保存したいと考えています。どうすればできますか？

これは、runge-kutta45 を使用した方程式微分のテストです: f(x,y)= (-5*x - y/5)^1/8 + 10

数値結果が異なるのはなぜですか？私が使用した：

と

私は MATLAB の初心者で、ルンゲ クッタ アルゴリズムを使用してフリードマン方程式を解決しようとしていました。

ご存じない場合は、フリードマン方程式は次の形式になります。

膨張宇宙のフリードマン方程式

ここで、曲率kは の値で与えられます-1, 0 or 1。

また、圧力の状態方程式は次のようpに与えられます。

圧力の状態方程式

だから、私はルンゲクッタアルゴリズムを持っています：

、しかし、私は3つの問題があります:

1. 関数の密度を記号的に表す方法がわかりません\rho。
2. g=@(t,x,y)時間に関する a の導関数があるため、2 番目の関数 を定義する方法がわかりません。
3. 最後に、 はどうなりk=sqrt(-1)ますか? その結果は象徴的に必要ですが、結果全体が数値的に必要だからです。

私の質問が基本的なもので申し訳ありませんが、これを行う方法がわからず、アドバイスや助けが必要です.

ありがとうございます ：）

Julia言語で一般的なRunge-Kuttaステップ関数を実装したいと思います. Python では、この RK4 関数が入力として取得する引数の 1 つとして関数を渡すことができます。Juliaでそのようにすると、パフォーマンスが低下しますか?

私の関数は次のようになります。

ここで、Ps はモデルのパラメーターを持つ Julia 型、Vs は ODE の変数の多次元配列、rhs は ODE の右辺 (部分時間導関数) です。

さまざまな方程式との統合のために C でこのコードを使用してきましたが、今日、表示されているものを統合するように変更しました。それは間違ったコーディングの問題ですか、それともこの方程式がこの手順で解決されることを意図していないということですか?

統合ルーチンとメイン コードの両方を次に示します。ありがとう

FIRST THE CODE (各パーツが何をするのかを説明するために切り刻みました。見た目が不快でしたら申し訳ありません)

パラメーター

方程式

主要

出口

初期条件

パラメータの定義 aa.w1=1; aa.mu1=1; aa.w2=6; aa.mu2=1; aa.eps=0;

統合コマンド

書き出す

そして、これが統合ルーチンです（問題ありません）