- ベストアンサー
ベッセル関数
円筒座標系での電磁場のマクスウェル方程式を磁場に関して解いて得られる解が複素数を引数とする0次のベッセル関数 AJ0(kr)、kが複素数、Aは実係数、rは実変数 で得られるのですが 引数を実数に変換する方法がわかりません。 純虚数の引数であれば実数の引数の変形ベッセル関数に変換でき、 実数の引数であれば手持ちの本にベッセル関数の値が載っているのですが 複素数の引数の場合の処理方法がわからなくて困っています。 よろしくお願いします。
- pierrot2002
- お礼率66% (2/3)
- 数学・算数
- 回答数3
- ありがとう数7
- みんなの回答 (3)
- 専門家の回答
質問者が選んだベストアンサー
z = x+yi として z^2 = (x^2 - y^2) + 2xyi や r exp(iθ) = r (cosθ+i sinθ) みたいに,実部と虚部をうまく分離したような形で処理できないか, というご質問と思いますが, 一般にはこういうふうにうまくはいきません. 今の例で言えば, J0(z) = u(x,y) + iv(x,y) となるように u,v を定義することはもちろん可能ですが, u,v がよく知られた関数になるわけではありません(と思います). 複素変数のベッセル関数というより仕方がないでしょう. 計算センターなど科学計算用ライブラリにはたいてい 複素変数のベッセル関数のサブルーチンが入っています. あるいは,Cで書かれたソースなどもどこかにあるかも知れませんが ちょっと今は手が回りません.
その他の回答 (2)
- siegmund
- ベストアンサー率64% (701/1090)
siegmund です. http://oshiete1.goo.ne.jp/kotaeru.php3?q=195527 で uni050 さんご紹介のページからちょっとたどった http://www.netlib.org/cgi-bin/search.pl で,bessel と入れてサーチしてみて下さい. 沢山出てきます.
お礼
なるほど、いろいろとありがとうございました。 どうにかやってみます。
- siegmund
- ベストアンサー率64% (701/1090)
siegmund です. 級数展開自体は x を複素数 z にして Jn(z) = (z/2)^n Σ[{(-1)^n (z/2)^2n } / { n! Γ(n + m + 1 )}] でそのまま使えます. n が非整数ですと頭の (z/2)^n のブランチが悩ましいところですが, 整数なら心配は要りませんね. 級数展開は z=0 の一点の周辺の近似になりますから, 実際の数値計算で得かどうかはよくわかりません. 数値計算プログラムとしては,1点の周りで正確であるよりも 対象とする変数の範囲すべてで規定以上の精度を保証しないといけませんから. アルゴリズムはちょっと私の手に余ります. なお,Mathemtica には組み込み関数で BesselJ(n,z) があります. あと,Maple や それに似ている MuPAD にもあるかもしれません(未確認). MuPAD はある条件を満たせばフリーで使えるらしいです (私はまだ使ったことはありませんが). http://home.att.ne.jp/star/mathmodern/link/linkcontents.html などからリンクをたどってみてください.
関連するQ&A
- 複素数を引数とする(?)ベッセル関数
皆様,お忙しい所回答お願いいたします。 一度他の方が似たような質問があったのですが,それだけでは解決できなかったので再度質問させていただきます。 質問はタイトルの通り,複素数を引数zとするベッセル関数Jn(z)がどうなるのか分かりません。 さまざまな参考書では実数もしくは純虚数をを引数とするベッセル関数の理論やプログラムのサブルーチンはあるのですが,複素数に関しては見つかりません。 大きな数学辞典を見ても,載っている数表は引数が実数のものばかりです。 どうしてこのような関数を必要としているのかというと,電磁界の円柱散乱問題の所で円柱媒質が導電率σを持つ損失性媒質の場合,波数kが複素数となり円柱内部電磁界の解析解に含まれるベッセル関数の引数が複素数となってしまうからです。(Jn(kr)という風に) 複素数を引数zとするベッセル関数Jn(z)の理論について,ご教授の方何卒お願いいたします。 また,参考文献等ありましたらそちらもご教授お願いいたします。
- 締切済み
- 数学・算数
- ベッセルの微分方程式
テキストによると、円筒座標系での電磁場のマクスウェル方程式を磁場に関して解いて得られる方程式が f’’+1/x*f’+k^2*f=0 解はベッセル関数 AJ0(kx)+BY0(kx) A,Bは定数 しかしこの方程式は一般的なベッセルの微分方程式と少し違います。 x^2f’’+xf’+x^2f=0 x^2で割り算してるのはともかく、係数kの分だけ違うのです。これでもベッセルの微分方程式であり解はベッセル関数であると言えるのでしょうか?
- ベストアンサー
- 数学・算数
- 平面波exp(-jx)をベッセル関数を用いてあらわすと・・・
円筒波動関数について勉強しています。 今、平面波exp(-jx)を円筒座標系で表そうとしているのですが、その変換式が Σa*Jn(ρ)*exp(jnΦ) (Σはn=-∞~n=∞まで)と表されています。 aは定数、Jn(ρ)はn次の第一種ベッセル関数、ρは円筒座標系の原点から外に伸びていく変数、Φは円筒座標系のxy面上の角度 この式について、わからないことがあります。 なぜこの式がx方向に進む平面波を表すのでしょうか?定性的なことが理解できません。 しかも第一種ベッセル関数は進行波でないのに進行波をあらわしている。 このこともさらに混乱を深めています。 どのように理解すればよいのでしょうか? あまりベッセル関数に関する知識がないのでできれば優しくおねがしします。 勉強している本はR.F.Harringtonのtime-harmonic electromagnetic fields です。
- 締切済み
- 物理学
- ベッセル関数を用いた解析について
円柱状の筒に分布した電界強度を、ベッセル関数で解析しようとしています。円筒の底をXY平面としたXYZ直交座標を設定し、あるZでのXY平面を解析面としています。 しかしそもそも、具体的にベッセル関数で解析するとはどういう事なのでしょうか?? また電界強度分布が径方向により異なる場合、ベッセル関数で解析することが出来るのか、出来るのであればどう解析を進めればよいのでしょうか??? どのような事柄でも構いません。皆さんのご意見をお聞かせ下さい。
- 締切済み
- 物理学
- 変形ベッセル関数の微分について。
変形ベッセル関数の微分について。 添付画像の一行目の式が微分公式です。 そこで質問ですが、二行目の式の左辺のようにベッセル関数の引数に係数が かかった場合は、二行目の式の右辺のようになりますか?
- ベストアンサー
- 数学・算数
- 複素数が引数の円柱関数
引数が複素数の円柱関数の解、及び計算手法がわかりません。|実数部|=|虚数部|のケルビン関数ではなく、 |実数部|≠|虚数部|の場合の解について教えてほしいです。宜しくお願いします。
- 締切済み
- 数学・算数
- フーリエ展開(ベッセル関数)
cos(m_fsinω_st)=J_0(m_f)+2Σ_{n=1}^{∞}J_{2n}cos2nω_st sin(m_fsinω_st)=2Σ_{n=1}^{∞}J_{2n+1}(m_f)sin(2n+1)ω_st J_n(m_f)などは第一種ベッセル関数です。 この2式を導くにはフーリエ余弦・正弦展開すればいいと思うのですが、係数を求める際に、積分がうまくできません。三角関数の引数にさらに三角関数が入っているので、どのように積分したらいいのでしょう。
- ベストアンサー
- 数学・算数
- 偶関数のフーリエ変換
なぜ偶関数の場合はフーリエ変換すると虚数成分がなくなってしまう(複素数計算を必要としない)のですか? 詳しく解説してあるサイトなんかも教えてもらえるとうれしいです。
- ベストアンサー
- 数学・算数
お礼
なるほど、そのまま複素変数のベッセル関数として計算させるしかないわけですね。 ありがとうございます。 サブルーチンを探してみます。 どのような計算手法なのか知りたいので、 第1種ベッセル関数の解が Jn(x) = (x/2)^n Σ[{(-1)^n (x/2)^2n } / { n! Γ(n + m + 1 )}] であるように、複素変数のベッセル関数を数値計算する方法が詳しく載っている本などをもしご存知でしたら教えていただきたいです。