ode

0熱

1答えて

このユーザー定義関数はscipy integrate odeintとどのように機能しますか？

私は、時間の関数としてRLC回路のコンデンサにかかる電圧をプロットする電気工学のプロジェクトに取り組んでいます。次のチュートリアルで私のコードに基づいています：Python-ODE for RLC Circuit。私のコードは完全に機能していますが、ユーザー定義関数とSciPy統合関数の仕組みをよりよく理解したいだけです。私はドキュメントを読んだが、まだいくつかの具体的な質問がある。ここに私のコー

0熱

1答えて

C++ - システムへのodeintの解が大きく変化する

私は、運転中の車両の振動をシミュレートする4つの連立方程式の系を解くためにodeintを使用しています。私の結果はMATLABのものと似ていると思っていましたが、それは残念ながら起こっていません。私は自分の方程式を何度もチェックしているのでエラーはないので、統合中に問題が起こっているはずです。私はC++スクリプトから何を得ているかを検証するためにMATLABのソリューションをコーディングしました

1熱

1答えて

ODE45からのリターンベクトルは、初期条件ベクトルより短い長さを持ち、diffEQを解決しようとしています

ODE45を使用して、垂直平面で回転する2つの回転バーの解決策を見つけようとしていますそれらの間の角度が90度と異なる場合にのみバー上に表示されます。私は、diffEQの定数としてa1-b4を使用しており、その値を関数に送る前にその値を行列に代入しています。私は6つの初期条件を送信しているが、ODE45関数からは5つしか戻っていないというエラーを返す。どのようにこれを修正するための任意のアイデア？

1熱

1答えて

Matlabのode45の結果をリアルタイムで表示するには？

は、我々は一階のODEのシステムのために、以下のソルバーを持っていると仮定してみましょう： % func.m function dydt = func(t,y) dydt = [y(2); (1-y(1)^2)*y(2)-y(1)]; とメインコード： % solver.m tspan=0:1:10; [t,y] = ode45(@func,tspan,[2; 0]);

0熱

1答えて

MATLABオイラー法

オイラー法を使ってフリードマン方程式を近似しようとしていますが、私のコードはうまくいかないようです。私はapでもslopeでも値を得ていません。 N=12; rho0=1; h=0.001; j=N*h; t=0:h:j; a0=1; G=6.67*10^-11; rhop=-3*(sqrt((8*pi*G)/3));

0熱

1答えて

私はこのDopri5実装で

を間違って何をやっている私は、Pythonに完全に新しいですし、次の頌歌を統合しよう： $ \ドット{X} = -2x-y^2 $ $ \ドット{y} = -yx^2 この結果、すべて0の配列になります。私は間違っていますか？それは主にコピーされたコードであり、別のものと結合されたコードではうまく機能しませんでした。 import numpy as np import matplotlib.p

0熱

1答えて

サスペンションの応答のルンゲクッタ積分

ルンゲクッタ積分を使用して、特定の時間間隔の車両サスペンションの加速度応答を解き、プロットする必要があります。応答がプロットされるとき、それは複数の振動を有するべきであるが、減衰のために、加速度応答が消滅するまで振動振幅は小さくなる。私はこの応答を消滅させることはできません。私の最終的な結果を見ると、配列の "z"成分は、それが初期状態の0にとどまり続けます。私はこれが問題かもしれないと思うが、

0熱

1答えて

scipy.integrate.ode、タプルを返す関数の問題

したがって、いくつかのソルバータイプを使用していくつかのODEを統合する関数を作成しました。何らかの理由で、私はあまり扱うことができませんが、ODEのRHSを定義するために書いた関数は、積分関数内で実行されると（またはエラーが示すように）タプルを返していますが、それは別々に返され、配列は...うまくいくように見えますが配列はうまく動作するはずです。 def solve_decay_system(t

1熱

1答えて

この数値メソッドの精度のMATLAB出力は、なぜ正確ではありませんか？

私の問題は構文だと思うのでここに投稿してください：私は2次ODEを分析しようとしています：y '' + 2y '+ 2y = e ^（ - x）* sin（x ）中点メソッドのMATLABコードを使用します。最初にODEを一次方程式の系に変換し、それを下に適用しようとしましたが、離散化[m]が増加すると、出力は.2718で停止します。例えば、M = 11件の収率： ans = 0.2724