問題タブ [bessel-functions]

次数 0 の球状ベッセル関数に関する問題があります。独自の球状ベッセル関数を作成しました。

これは、私のすべてのテストケースで Mathematicas に組み込まれているものと一致しているようです。問題はnuまたはにありx =0ます。Mathematica は正しく 1 を返しますが、私の MATLAB スクリプトは NaN を返します。コードを修正するにはどうすればよいですかx = 0:1:5?

私のカスタム関数はこれを行うための最良の方法ですか?

ありがとう

この Matlab コードのベクトル化を手伝ってくれる人はいますか? 具体的な問題は、ベクトル入力の合計関数とベッセル関数です。ありがとうございました！

一連のリングに囲まれた中心スポットを持つ円形の回折パターンを作成しようとしています。それを行うには、コードで定義されているベッセル積分が必要です。

私の問題は、コードが実行されるのを 10 分待ったのに何も表示されなかったように、時間がかかりすぎることです。私のベッセル積分にはポイントごとに1000回の反復があるため、これを手伝ってくれる人はいますか?

私は正しい軌道に乗っていますか？

Mark Newmans の本 Computational Physics で Python と計算物理学を独学しようとしています。演習は計算物理学の 5.4 です。この章へのリンクは次のとおりです。9ページにあります。 http://www-personal.umich.edu/~mejn/cp/chapters/int.pdf

これが私が作ろうとしているイメージです。

.

私のコード:

Unity3d で第 1 種のベッセル関数を使用する必要があります (C# を使用)。Unity で mathnet.numerics .dll を使用しようとした後、明らかに unity3d が .NET 4 をサポートしていないために多くのエラーが見つかりました。その後、jetbrains dotPeek デコンパイラーを使用してベッセル関数のソースを見つけることにしました。

しかし、ベッセル関数は間違った出力を出し、正直なところ、数学は私の給与水準を超えているため、デバッグは非常に困難でした. 誰でも私が使用したもののエラーを見ることができますか?

たとえば、x = 20 の出力は 43558282.5571362 であり、0.167025 であると想定されています

私は試しboost::math::cyl_bessel_k(x,y) * exp(y)ました。ほとんどの場合、これは Matlab の scaled と同じbesselk(x,y,1)です。しかし、場合によっては (例: x=1、y=2000) と の両方besselk(x,y)=0でboost::math::cyl_bessel_k(x,y)=0、Matlab のスケーリングされbesselk(x,y,1)たバージョンでは異なる値が得られ、 が変化し10^-3ます。しかし、boost::math::cyl_bessel_k(x,y) * exp(y)戻ります-nan。Matlab のbesselk(x,y,1). どうすればこれを処理できますか?

Matlab では、ベッセル関数 j1(x) を微分すると、x= 1:10 で 10 ではなく 9 つの値が得られるのはなぜですか?