問題タブ [differential-equations]

For questions regarding programming in ECMAScript (JavaScript/JS) and its various dialects/implementations (excluding ActionScript). Note JavaScript is NOT the same as Java! Please include all relevant tags on your question; e.g., [node.js], [jquery], [json], [reactjs], [angular], [ember.js], [vue.js], [typescript], [svelte], etc.

0 投票する
1 に答える
289 参照

wolfram-mathematica - キャッシング手数料+プロモーションレートの同等の金利を見つける

私はクレジットカードから「金額」のキャッシングサービスを受け取り、前払いの「料金」(パーセンテージで指定)を支払います。時間「len」のプロモーションレートは「int」です。私は毎月、少なくとも'min'%の未払い額を支払わなければなりません。

「金額」を「p」%の利息を稼ぐ投資口座に入れ、この口座から毎月の支払いも行います。

質問:「p」のどの値に対して、「len」の時間の後でも損益分岐点になりますか?

Mathematicaで設定する方法は次のとおりです。

Mathematicaは「DSolve::bvnul:一般解のいくつかの分岐では、与えられた境界条件は空の解につながる」を返します。これは実際には合理的です。上記の微分方程式の解を生成する「p」の値は1つだけです。 。

Mathematicaにこの値を見つけるように強制するにはどうすればよいですか?

owed [t]を解いてから、owed[t]をhave[t]に置き換えて、owed [len] == have [len]を解いてみましたが、同様のエラーが発生します。「owed[len]== have [len]」でReduceを実行すると、複雑で醜いものが生成されました。

0 投票する
1 に答える
547 参照

c++ - 差分方程式にサンプリング レートの影響を追加するにはどうすればよいですか?

アキュムレータと微分器を設計するために、以下の差分方程式を実装したいと思います。

アキュムレータ:
y[n] = y[n-1] + x[n]、ここy[n]で、 は n 番目の出力、x[n]は n 番目の入力です。

微分子:
y[n] = x[n] - x[n-1]、ここy[n]で、 は n 番目の出力、x[n]は n 番目の入力です。

アキュムレータ:

差別化要因:

入力ベクトルのデータは、連続時間信号から生成されます。つまり、この離散時間信号を取得するためにサンプリングされます。また、サンプリング レートは T dです。

入力信号のサンプリング レートを上げると、アキュムレータの出力の振幅が大きくなります。差別化するサンプルがたまたま多いため、これは予期されることです。一方、サンプリング レートの増加に伴い、微分器の出力振幅は減少します。微分器と積分器の両方の振幅は、サンプリング レートがユニタリ (つまり、1.0 サンプル/秒) の場合、(数学が示す値で) 正確です。

私の質問は;
出力ベクトルでのサンプリング レートの影響を無効にするにはどうすればよいですか?

私の一般的な質問は次のとおりです。
任意の因果離散時間伝達関数 から差分方程式を取得するとしますH(z)。この場合、差分方程式は何でもかまいません。この一般的なケースでサンプリング レートの影響を無効にするにはどうすればよいですか?

(詳細については、このフォーラム スレッドを参照してください: 関連するフォーラム スレッド)

0 投票する
2 に答える
5696 参照

python - 熱方程式とクランクニコルソンの有限差分解における境界条件の適用

以下のコードは、初期条件10 * np.sin(np.pi * x)で両端がゼロ温度に保たれているロッドを表す1D熱方程式を解きます。

この計算では、ディリクレ境界条件(両端の温度がゼロ)はどのように機能しますか?行列Aの上下の行には、ゼロ以外の2つの要素が含まれており、欠落している3番目の要素ディリクレ条件であると言われました。しかし、この条件がどのメカニズムで計算に影響するのかわかりません。Aに要素がない場合、u_{0}またはu_{n}をゼロにするにはどうすればよいですか?

以下の有限差分法では、Crank-Nicholsonを使用しています。

0 投票する
1 に答える
1498 参照

wolfram-mathematica - Mathematicaでドメインの一部で初期条件を指定するにはどうすればよいですか?

単純な熱方程式U_t=U_xxについて、有限体積法の解とMathematicaの解を比較しようとしています。このために、MathematicaでNDSolve関数の初期条件と境界条件を指定する必要があります。境界U=90にしたい。初期条件として、境界を除くすべての正方形ドメインに100を設定したい。これどうやってするの?明らかに私のコードは機能しません。

すべてのドメインで100に等しい初期条件を設定しようとします。

0 投票する
1 に答える
1569 参照

math - 分離可能な変数を使用して IVP を解決する方法

質問は:

x dy/dx = 2y ; y(0)=0

この問題を解くと、積分定数 'c' がゼロになるため、与えられた IVP の解を計算するには、その値を見つける必要があります。

0 投票する
2 に答える
700 参照

matlab - 無限を超える matlab の倍精度積分のサポート


無限サポートに対する二重積分を評価しようとしています。何かのようなもの

quad2d(F, 0, inf, 0, inf)

ただし、quad2d には有限の境界が必要です。無限を扱う代替関数が存在するかどうか知っていますか? 1 つの変数 quadgk でそれを行います。

0 投票する
6 に答える
45101 参照

python - Python の微分方程式

私は DE システムを使用していますが、微分方程式を解くために最も一般的に使用されている Python ライブラリがあればそれを知りたいと思っていました。

私の方程式は非線形一次方程式です。

0 投票する
2 に答える
654 参照

math - How to obtain the numerical solution of these differential equations with matlab

I have differential equations derived from epidemic spreading. I want to obtain the numerical solutions. Here's the equations,

enter image description here

t is a independent variable and ranges from [0,100]. The initial value is

At first, I planned to deal these with ode45 function in matlab, however, I don't know how to express the series and the combination. So I'm asking for help here.

**

The problem is how to express the right side of the equations as the odefun, which is a parameter in the ode45 function.

**

0 投票する
1 に答える
4105 参照

math - MATLABを使用した2次非線形微分方程式

x=0からx=1でmatlabに解くために2階微分方程式を入力しようとしています。方法がわかりません。方程式は次のとおりです。

y'' = 1 + 0.1 \ sqrt {1+(y')^ 2}

初期条件はゼロです。