differential-equations

0熱

1答えて

この割り当てのタスクは、レイヤの入力に対する損失の偏微分を計算することです。チェーンルールを実装する必要があります。私は概念的にどのように機能を設定するのが難しいのですか。アドバイスやヒントをいただければ幸いです！関数変数のサンプルデータは、最下部にあります。 def dense_grad_input(x_input, grad_output, W, b): """Calcula

-1熱

1答えて

Microsoft Oslo Libraryを使用して剛性方程式を解く方法はありますか？

Microsoft Oslo Libraryのどの句は、スティッフ常微分方程式を解くことができますか？には、それらをまったく解決する別の方法がありますか？

2熱

1答えて

math.jsでLotka-Volterra方程式を解く

math.jsのホームページにある微分方程式を解く例がありますが、それはかなり複雑で、数学を適用するには個人的に十分な情報を提供していません。他の同様の問題についてはjs。だから私がやろうとしているのは、Lotka–Volterra equations for predator-prey simulationを解決することです。システム内の2つの方程式があります、しかしhttp://mathjs.

1熱

1答えて

intまたは分数を期待し、sympyで％と％を取得

私は微分方程式y '= 6e ^（2x-y）を解く必要があります。私はdsolve（）でsympyでそれをしようとしています。 sol = dsolve(Derivative(f(x), x) - 6 *(e**(2*x-f(x))), f(x)) しかし、常にエラー期待int型または画分を取得するには、7.38905609893065022723042746058と6 問題は何であるので

0熱

1答えて

ポアンカレ区をプロットする方法は？（Duffing Oscillator）

Duffing方程式の単純なリミットサイクルを成功裏に示すプログラムを作成しました。しかし、私は今、このケースのためにポアンカレのセクションをプロットする必要があります。フェーズスペースダイアグラムのスナップショットを一定の時間間隔（例えば、t*omega = 2*pi*n）で取得することで、これを行う必要があります。この場合、オメガが1に設定されているので、これはちょうどt = 2*pi*nの

0熱

1答えて

numpy.ndarrayオブジェクトを取得するときに呼び出し可能なエラーが発生しない

エラーが発生しました。なぜそれが起こっているのか100％確信していません。私は、オイラーとヒュンの方法の両方を実行し、互いに対してプロットしようとしています。どんな助けや提案もすばらしいでしょう！私はPythonでかなり初心者ですので、これは私には少し新しくなっています。 import numpy #eulers yi = 0 h = 0.2 te = [0.0,2.0,0.2] def

0熱

1答えて

ScipyでODEの集合を解く

システムに到達するまでの時間の経過とともに変化する細胞、基質および生成物（すなわち、）の濃度を予測するアルゴリズム（scipy.integrate.odeint（）を使用）を開発しようとしています。定常状態（約100時間または200時間）。バイオリアクター中の細胞の初期濃度は0.1 /であり、最初は反応器中にグルコースまたは生成物が存在しない。 0.01 /ℎと0.25/betweenの間のさまざ

-1熱

2答えて

この2つのコードの違いは何ですか？

差分問題を解くための2つのコードを書きました。どちらも予測子補正ツールを使っています。最初は差分系の解法でもありますので、基本的には系で解く方程式の数に等しい大きさのベクトルyを使いますので、基本的には各工程を上書きしています前のステップを保存します。代わりに、ステップ数に等しいサイズのベクトルを使用するので、使用している要素を簡単に制御できます。私はデバッグの多くをやったが、私は私が間違って

1熱

1答えて

メソッドを書く別の方法

私は、微分方程式の右辺の項を定義するクラスを持っています：このクラスは、rhs関数とその派生物を計算する方法を提供します関数はこのようにベクトルコンテナに格納されますこのクラスは微分方程式の系にも適しています。ここ私は変更するだろうメソッドが定義されているインタフェース template <typename Type = double> class rhsODEProblem {