問題タブ [derivative]

http://www.mit.edu/people/abbe/matlab/ode.htmlのチュートリアルに従い、次のように関数を準備しました。

を使用して呼び出すと

エラーが発生する

私も試しました

XY= ode45(@diffxy,[0 10],[0;1;0])

誰でも何か考えがありますか？

「テイラー級数」を計算するプログラムを作成していて、出力として「inf」を取得しています。プログラムが終了しておらず、すべてがまだ正常に機能していませんが、出力を確認できるようにしたいと思います。出力を短くしたり、表示したりする方法はありますか？助けてくれてありがとう。

GLSLの関数dFdx()と関数を理解しようとしています。dFdy()

私は次のことを理解しています:

1. 導関数は変化率です
2. 2 つのパラメーターを持つ関数の偏導関数は、パラメーターの 1 つを一定に保ちながら関数を微分する場合です。
3. dFdx()dFdy()現在のフラグメントと隣接するフラグメントの間で値が変化する割合を見つけます。

変化率の意味がわかりません。フラグメント座標の変化率ですか？

フラグメント シェーダーの 2 つの呼び出し間の任意の変数の変化率を見つけることができるのでしょうか? シェーダーの呼び出しは、隣接する呼び出しから変数を「読み取っていますか?」（単純な）例：

dそれぞれ -1.0 と 1.0でしょうか?

質量ばねと二重振り子を備えた (やや奇妙な) システムの運動方程式をシミュレートしています。このシステムには、質量行列と関数 f(x) があり、ode45 を呼び出して解決します。

5 つの状態変数 q= [ QDot, phi, phiDot, r, rDot]'; があります。(QDot は現在の QDot に依存するものがないため削除されました。) ここで、いくつかの力を計算するために、ode45 が各積分ステップで計算する rDotDot の計算値も保存したいと思いますが、ode45 はこれを返しません。私は少し調べましたが、私が見つけた唯一の2つの解決策は、a）これを3次問題に変え、状態ベクトルにphiDotDotとrDotDotを追加することです。これは可能な限り回避したいと思います。これはすでに非線形であり、これにより事態がさら​​に悪化し、計算時間が大幅に増加するからです。

b)ここで説明されているように、状態を拡張して関数を直接計算します。ただし、例では、質量行列にゼロの行を追加すると彼は言います。そうしないと導関数を積分し、ある点で評価するだけでなく、質量行列を特異にするので、これは理にかなっています。私にはうまくいかないようです...

これは基本的なことのように思えます (状態ベクトルの微分値を取得するため)。私が考えていない明らかなことはありますか? （または、それほど明白ではないものでも構いません....）

ああ、そしてグローバル変数はそれほど素晴らしいものではありません。なぜなら、ode45 は f() 関数をステップの調整中に何度も呼び出すため、グローバル変数のサイズと返される状態ベクトル q がまったく一致しないからです。

誰かがそれを必要とする場合に備えて、コードは次のとおりです。

FFTW を使用したサインのこの単純な導関数に苦労しています。一見問題ないように見えますが、正確な解と比較するとかなり大きなエラー (5e-6) があります... c2r を取得した後、複雑な入力がすべて台無しになっていることがわかりますが、私には思えます同じ複雑な入力が問題の原因です...何が間違っていますか? ポインタは一切使用せず、すべてをできるだけシンプルにしようとしました。それでも何が悪いのかわかりません。どんな助けでも大歓迎です！ありがとう！！！

与えられた:

1. ハスケル
2. df/dz複素平面上で定義された複素数値関数(を aとしUましょう)。zComplex Double
3. が定義されているz1からのポイント。Udf/dz

質問：

ポイントで導関数f(z)である関数の値を取得する方法は? 私はe。複素平面を想定して、導関数のみを指定して元の関数の値を復元する方法は?df/dzz1

この質問は、複雑な関数の積分の計算に関する以前の質問と多少関連していますが、それらは別のものです。ここで、スカラー値を計算するのではなく、導関数が与えられた元の関数を見つけることに関心があります。基本的に、この導関数の不定積分を計算しています。

私が理解している限り、ピクセル シェーダーはピクセル単位で動作します。しかし、導関数を計算するddxやのような関数があります。ddyしかし、たった 1 つのピクセル座標から導関数を計算するにはどうすればよいでしょうか?? 誰かがこれらについて私を助けることができますか? これらはまた、次のような疑問を提起します

サンプル関数はピクセルごとに計算されるということですか?? サンプラー命令はサブスパンごとに動作すると思いました。

例：

次の 16 ピクセルがあるとします。

そして、私のピクセルシェーダーは次のようになります:

上記のコードは、ピクセル座標の 4 x 4 グリッドで何回呼び出されるでしょうか? ありがとう。

SymPy でシンボルを区別しようとすると、次のようになります。

シンボルをその共役に関して微分すると、結果は次のようになります

ただし、シンボルの共役を微分しようとすると、SymPyはそれを行いません

これでも正しいですが、結果はゼロになるはずです。SimPy に共役の導関数を実行させるにはどうすればよいですか?