問題タブ [bessel-functions]

ベッセル関数と積分を含む関数のパラメータを推定したいと思っています。しかし、実行しようとすると、「f(x, ...) のエラー: 関数 "BesselI" が見つかりませんでした」というメッセージが表示されました。私はそれを修正する方法を知りません。関連する提案をいただければ幸いです。

ニュートン ラフソン最大化 反復回数: 0 リターン コード: 100 初期値が範囲外です。

「50 件以上の警告がありました (最初の 50 件を表示するには warnings() を使用してください)」というメッセージが表示され、warnings() を使用すると、次の警告が表示されます。

私のコードでいくつかのガイダンスを得ることを望んでいます。私はPythonにまったく慣れていないので、コード内の変更されたベッセル関数k0に問題があります。実際のコードは 2 つの部分に分かれています

1) 最初のセクションでは、sympy を使用して t を解き、t=0 での x0 の値を返します。

2) 2 番目の部分では、この値 x0 を使用して、さまざまな x の値について t を解きます。

私が困難を抱えているのは2番目のステップです.scipy.special.k0関数は、記載された解決されたx0値に対して機能していないようで、その理由がわかりません

コードを実行するたびに、次のエラーが発生します

どんな助けでも大歓迎です、そして前もって感謝します。

シリンダーが非常に高温に加熱され、周囲温度がゼロ以外の温度、つまり T0 であり、シリンダーの初期温度が Ti であるクエンチング剤に急冷されるとします。

この方程式は、特定の熱伝導率と熱伝達係数に対して、特定の時間 t における表面と中心の間の任意の点の温度を求めるものです。

J0 は第 1 種ベッセル関数で、J1=-d(J0)/dz

固有値 beta(m) は、超越方程式の根から取得されます。

上記の式を使用して冷却曲線を描くにはどうすればよいですか?

私はmatlabが初めてで、これは私のリーグから完全に外れています。他の質問でこのようなものは見たことがありません。助けてくれてありがとう。

式は I ₀ (a*x)=b で、a、b は定数、I ₀ (y) はゼロ次ベッセル関数です。

I ₀ (x)=1/(2*pi)*integral(x*cos(t)) dt [0 ->2*pi から]

a と b が与えられたときの x の値を取得したい。結果が欲しいだけなので、Matlabで解く手順を全て実装する必要はありませんし、おおよその答えでも構いません。

ゼロ次の MATLABベッセル関数 --->を使用しようとしていますbesselj(0,Z)。私が気づいたのは、n乗の実部（たとえば、 ）を処理できますが、複雑な引数でe100（ ）よりも大きい場合、虚部を処理できないことです。e2

例:(2+2*i)処理する複雑な引数。(20000+20000*i)----->しかし、それは「無限大」という答えを返すことはできません。

これは私が得るものです:

誰かがこれを説明できますか？虚数部のべき乗に制限はありますか？

scipy.special次数 n と引数 x のベッセル関数が xjv(n,x)でベクトル化されていることに気付きました。

In [14]: import scipy.special as sp In [16]: sp.jv(1, range(3)) # n=1, [x=0,1,2] Out[16]: array([ 0., 0.44005059, 0.57672481])

しかし、対応する球面ベッセル関数のベクトル化された形式はありませんsp.sph_jn。

さらに、球ベッセル関数は、1 回のパスで N のすべての次数を計算します。したがって、n=5引数 にベッセル関数が必要な場合x=10、n=1,2,3,4,5 が返されます。実際には、jn とその導関数を 1 回のパスで返します。

API にこの非対称性が存在するのはなぜですか。ベクトル化された、または少なくともより迅速に (つまり cython で) 球状ベッセル関数を返すライブラリを知っている人はいますか?

scipy ライブラリを使用して、ベッセル関数のゼロを見つけています。低次のベッセル関数で動作するようです

しかし、カットオフ後、ハングします

上記の行は決して戻りません。これは scipy の予想される動作ですか?