問題タブ [bessel-functions]

ばかげていたり曖昧に聞こえたら申し訳ありませんが、ベッセル関数に関するコードについて質問があります。私は、ゼロが 1 にスケーリングされる最初の 5 つの 0 次ベッセル関数を使用して、単純な関数 (正確には f(x)=1-x) を表す宿題を課されました。さらに、それらの係数を見つける必要があります。

まず、助けを求める前に、この問題に取り組んだことを示さなければならないことはわかっていますが、どこから始めればよいかわかりません。一次ベッセル関数と二次ベッセル関数については知っていますが、0 次ベッセル関数が何であるかはわかりません。さらに、ベッセル関数で関数を表現できることさえ知りませんでした。テイラー展開またはフーリエ表現を使用して関数を近似できることは知っていますが、ベッセル関数でそれを行う方法がわかりません。このウェブサイトを検索したところ、私の同級生のようで、失礼なことに、私たちの課題をコピーして貼り付けただけで、そのスレッドは閉じられました.

ですから、聖人がここで少なくとも私を正しい方向に向けることができれば、それは素晴らしいことです. これは単純な関数であり、Matlab にはベッセル関数があることを知っているので、コーディングはそれほど難しくありません。微分方程式を解く方法を使用して別の関数を表す方法がわかりません。ああ、そして係数？係数は？私の無知を許して、助けてください！

1 次と 2 次 (iv と kv) の両方の修正ベッセル関数を使用するコードがあります。迷惑なことに、それらには制限があるようです。それらは iv(0,713) と kv(0,697) で、それぞれに 1 を追加すると、それぞれ無限大と 0 になります。これよりも高い値を使用する必要があるため、多くの場合 2000 以上の値を使用する必要があるため、これは私にとって問題です。これらで除算しようとすると、最終的に 0 または無限大で除算されます。つまり、エラーまたはゼロが発生しますが、どちらも必要ありません。

私はscipy bessel functionsを使用しています。はるかに小さい数値とはるかに大きい数値に対応できるより優れた関数、またはこれらの大きな数値を処理するように Python を変更する方法はありますか。なぜPythonが700をはるかに超えてこれらをうまく処理できないのか、ここでの本当の問題が何であるかはわかりません.それは関数ですか、それともPythonですか?

Python が既にそれを行っているかどうかはわかりませんが、たとえば、最初の 5 ～ 10 桁 *10^x だけが必要です。つまり、1000桁すべてを必要とするわけではありません。おそらくこれは、Wolfram Alphaの計算方法と比較して、Pythonの計算方法に問題があるのでしょうか?

この関数をプロットする必要があります

これが私のコードですが、それが正しいかどうかはわかりません。ドットマークがどこにあるべきかわかりません。X軸をデグレスに設定するにはどうすればよいですか？

私のもう1つの質問は、この関数を極座標でプロットする方法です。

私はこのようなものを作りました。

その関数を（y軸で）回転させて3Dプロットを取得することは可能ですか？

この質問の単純さをお詫びします。

Pythonで方程式を実装したいと思います。この式で、K_0 はゼロ次の修正ベッセル関数です。

PythonでK_0を実装する最良の方法は何ですか?

CUDA を使用して Bessel 関数 (例として J0(x) ) を実行しようとしています。式はこちら

結果をいくつかのイプシロン値内に収めようとします。だからここにコードがあります

実行すると黒い画面が表示され、Windows は nVidia ドライバーが応答しなかったと表示し、回復しました。host_resultfまた、コンソール出力では、配列にゼロしかありません。どうしたの？いくつかのイプシロン内で関数を適切に実行するにはどうすればよいですか?

私は自分のクラスの課題に取り組んでおり、選択したプログラム (Matlab を選択しました) を使用してコードを記述し、4 次ルンゲクッタ法を使用してベッセル関数の微分方程式を解くことになっています。参考までに、ベッセル関数 DE は次のとおりです。

x^2*(J_n)''+x*(J_n)'+(x^2-n^2)*J_n=0.

これを 2 つの結合された 1 次 DE に分離するには、次のようにします。

(J_n)'=Z_n および

(Z_n)'+(1/x)*Z_n+[(x^2-n^2)/x^2]*J_n=0.

この課題の前に、Matlab やその他のプログラミング言語の経験はありません。Matlab に 'ode45' コマンドがあることは知っていますが、Matlab のコマンドに依存するのではなく、自分でコードを書くことになっています。これまでのところ、ベッセル関数の n=0 ケースに取り組んできましたが、関数をプロットしようとするとエラーが発生し続けます。私が持っている現在のエラーは、「「double」型の入力引数に対して定義されていない関数またはメソッド「J」です。しかし、このエラーを修正する方法も、コードが正しいかどうかもわかりません。どこが間違っているのか、またはこのコードを書く正しい方法は何か教えてもらえますか?

助けてくれてありがとう

私が取り組んでいる方程式は

$$ E = M_e + \sum_{n = 1}^N\frac{2}{n}\mathcal{J}_n(ne)\sin(nM_e) $$

$\mathcal{J}_n(x)$ は、第 1 種の n 番目のベッセル関数です。

テストとして、最初の 6 つのベッセル関数をプロットしたところ、すべてが正しく機能しました。$n * e$ の引数を入力すると、プロットは期待したものではありません。

プロットは、n = 10 の場合のようになります。

Bessel ローパス フィルターの Java または C 実装を探しています。多くの検索の後、私は多くを見つけることができませんでした。誰かが私を正しい方向に向けることができますか? 前もって感謝します

ブースト ドキュメントを確認しましたが、このエラーは0に 2 番目の引数として渡した結果であることを理解してboost::math::cyl_bessel_kいます。

それを受け取った場合。

私は 98% が C++ 初心者です。プログラムのクラッシュを防ぐために、このエラーをどのように処理できますか?