differential-equations

2熱

1答えて

コードを実行すると、サブルーチンにエラーメッセージが表示されます。このコードは、Kincaid & Gauss Seidelの楕円偏微分方程式を解く方法に関するCheneyの本の練習からのものです。エラーメッセージがある：（1）における変数定義コンテキスト（割当）でINTENT（IN）と仮引数 'U'。私は以下のコードの（1）を参照します。エラーメッセージが表示されないようにサブルーチンを

0熱

2答えて

MATLABでネストされた非線形ODEを定義して解く

にはa(t) = k*y'^2を代入します。ここで、yはt（時間）の関数です。私の初期条件はy(0) = 600;とMATLABで y'(0) = 0;私は ydd = diff(y,t,2) == -g + a(t)/m; でy''を定義する方法を知っているが、私は、これは「ネストされた」非線形微分方程式であり、私は事実で迷ってしまいましたMATLABでそれを解決することはもちろん、それをど

3熱

1答えて

積分微分方程式

Iは、フォームの微分方程式のセットを有する： x1dot = x3; x2dot = x2; x3dot = x1; x4dot = x2 + integral(x1,t,tend) Iはx1、tstartでx2とための境界条件を有する、tendでx4。積分項がなければ、それはBVP4Cを使った単純な実装です。積分に使うことができるBVPソルバの状態に対して以前の解を持つことが可能かど

-1熱

1答えて

pythonのケプラー軌道数値解の微分方程式は？

私はケプラーポテンシャルV（r）= - 1/rの軌道の微分方程式を解くコードを書こうとしています。V（r）= - 1/r M/M^2 = 0 ここで、u = 1/R と我々は - D^2U/D（FI）^ 2 + U：このようになり微分方程式最終的にr（fi）を探しますと今は def func(y,fi): m=4 M=5 return [y[1],m/M^2-y[0

0熱

1答えて

ループによって生成された方程式を受け入れるためのodeINTの取得に問題がある

ループによって生成された複数の方程式の集合を統合するためにPythonのodeINTを使いたいとします。方程式はすべて結合されているので、odeINTを1回呼び出すことで同時に積分する必要があります。問題は、初期条件（ "y0"）はリストのリストまたは行列（単にリストではない）でなければならないということです。 odeINTは、「初期条件y0は1次元でなければならない」というエラーを出します。私は

0熱

1答えて

2次微分方程式は、私が一次方程式に、次の2次方程式を変換する方法を理解していない

1次への変換：のx「」（T）= 686 -0.5 * X '（tは^ 2）+ 10（X（t）は+ 150）私は、置換法を使用しますがx（t）を変換する方法を確認していない： X '（T）= Y（t）は、 x ''（t）= y '（t）、 x（t）= y^2（t）/ 2, Y '（T）= 686 -0.5 * Y（T^2）+ 10（y^2（T）/ 2 + 150）

0熱

1答えて

fminconの目的関数と制約関数の単一計算

私のプログラムでは、線形制約もない最適化ルーチンで境界値問題を解決する必要があります。この問題にはfminconを使用し、目的関数と非線形制約関数を評価するにはBVPの解が必要です。現在、私は目的関数と制約関数でBVPを解いています。目的関数内のBVPを1回評価した後、もう1つのBVP評価を減らすために、その関数を制約関数に渡すことができる、より効率的な方法がありますか？任意のアイデア

0熱

1答えて

MATLAB：私は方程式</p> <p><a href="https://i.stack.imgur.com/pOTax.gif" rel="nofollow noreferrer"><img src="https://i.stack.imgur.com/pOTax.gif" alt="enter image description here"></a></p> <p>のシステムの下に持ってMathWorks社のMATLAB

と非線形二次ODEのシステムを解くことは、MATLABの関数ode45でこれを解決することは可能ですか？私は2次方程式を2次の1次方程式に変換する必要があることを知っていますが、私の混乱は、のとthetaの派生品であるという用語に由来します。 P.S. Yは、Yが（3）、（2）シータの誘導体である、（1）シータことがYを考慮 function dy = pend(t,y) beta = 1;

1熱

1答えて

MathWorks社のMATLABプロットのベクトルが同じ長さ

でなければならない私は、撮影-二分法を使用して、次の境界値問題の解決に関するプログラムを書いています： y''-y+x=0, y(0)=y(1)=0. 私は最初の一次方程式のシステムにこれを変換し、 y'=z を設定し、その後、私はdydtはベクトル（Y 'Z'）を表し、およびスクリプトファイルを思い付くしてみましょう： function dydt=shoot(t,y) dydt=[y(

2熱

2答えて

Pythonユーザ入力方程式

私はEulers法を使って微分方程式を解く単純なPythonスクリプトを書いていますが、今は新しい方程式を解くたびにソースコードを変更する必要があります。（安全性が、これは個人的なプロジェクトであるとして、現時点での任意の関心事ではありません。）私の質問：すると、ユーザーの入力として数式を入力し、それを作ることも可能です変数？たとえば、を入力し、eを後で使用するために入力する場合は、e=in