ベッセル関数計算機

ベッセル関数は、円筒対称性や動径方向の対称性を持つ微分方程式を、より単純な部分に分離すると現れる特殊関数です。円形膜の振動、円筒内の熱流、円形導波管の電磁波、音響モード、回折、流体運動など、多くの数理モデルで使われます。通常のベッセル方程式には独立な二つの解の族があります。n 次の J は第一種ベッセル関数と呼ばれ、非負整数次数では原点で有限です。n 次の Y は第二種の解で、ゼロに特異点があるため、この計算機では正の引数を必要とします。 変形ベッセル関数 I と K は、引数の二乗を含む項の符号を変えた関連方程式を解きます。正の引数に対して I は一般に指数的に増加し、K は指数的に減衰します。この対照的な性質は、一方の動径解を中心付近で有界に保ち、もう一方を源から遠くで減衰させる境界値問題に役立ちます。K もゼロで特異となるため、ここでは正の引数が必要です。種類の選択は計算の都合だけでなく、微分方程式と境界条件に基づいて行います。 整数次数の J と I は収束するべき級数で評価します。通常の浮動小数点精度で新しい項の寄与が無視できるまで加算します。Y と K は安定した近似式で零次と一次を求め、標準的な漸化式で高い整数次数へ進みます。これらの方法は日常的な広い範囲で実用的な工学値を与えます。ただし、極端に大きな次数や引数では、振動、桁落ち、オーバーフロー、アンダーフローにより有効桁が失われ、ブラウザー上の浮動小数点計算が難しくなることがあります。 次数 n は動径モードを決め、このツールでは非負整数でなければなりません。数学関数の引数 x は無次元です。物理への応用では波数と半径の積などを表すことが多く、評価前に単位が相殺されます。結果はモデル全体とその正規化を踏まえて解釈してください。精度に敏感な研究では専門の数値ライブラリと比較し、必要に応じて精度を上げてください。 結果の確認には複数の恒等式が使えます。零次の J と零次の I はゼロでともに一となり、正整数次数の J と I は原点でゼロになります。隣接次数は漸化式を満たし、振動する J と Y の零点は交互に現れます。以下の例では、よく使われる参考値と通常型・変形型の振る舞いの違いを示します。