• ベストアンサー

ベッセル関数って、

大学生です。 最近、勉強していると、振動の分野で「ベッセル関数」なるものが現れました。 振動の分野は全く専門外で、突然ベッセル関数を使われて、こまってます。 また、やさしく説明してくれる本も、探してますが、あまり時間がありません。  この「ベッセル関数」は、    いったいどんな関数で、    この関数をなぜ使うのか?    この関数を使うと、何が出てくるのか、    この関数は何を意味しているのか、  教えてください。 (ちなみに、円筒形の固有振動数を求める式でベッセル関数が出てきました。)

質問者が選んだベストアンサー

  • ベストアンサー
  • Umada
  • ベストアンサー率83% (1169/1405)
回答No.2

chestnutさん、改めましてこんにちは。 さてBessel関数J_n(x)の次数n、これはモードと関係はするのですが関係の仕方がちょっと複雑です。 2次元円形膜の振動はご承知のように、半径方向の関数R(r)と角度方向の関数Θ(θ)に変数分離をして解きますが、R(r)とΘ(θ)のモードはバラバラに決めてよいわけでないのです。 角度方向の関数は 1(定数)・・・0次 角度方向はどこでも同じ値 sin θ・・・1次 角度方向に一周くるりと回ると、節が二つある sin 2θ・・・2次 角度方向に一周くるりと回ると、節が四つある sin 3θ・・・3次 角度方向に一周くるりと回ると、節が六つある ・ ・ となるわけですが、角度方向が0次の振動の場合、半径方向の関数R(r)は、0次のBessel関数J_0(x)と自動的に決まります*。 ただし関数のグラフを見て既にお分かりのように、J_0(x)=0となるようなx(=振動の節となる点)は無数に存在しますから、それらのどれを膜の境界に合わせるか(合うか)によって半径方向にも多くのモードが存在することがお分かり頂けると思います。 同様に角度方向が1次の振動の場合は、半径方向の関数R(r)は、1次のBessel関数J_1(x)と決まります。この場合も同様に、J_1(x)=0となるようなxは無数に存在し、モードも無数にあることになります。 すなわち2次元円形膜の振動のモードは、 (1)角度方向の次数nと、 (2)それにより決まるBessel関数J_n(x)で、J_n(x)=0となるxのどれが境界に当たるかによって決まる次数m の二つで特定されます。 http://www.ep.sci.hokudai.ac.jp/~minobe/mth2ex/ の一番下に「1999講義ノート」というリンクがあり、pdfでダウンロードできるようになっています。モードの図などもあり大変分かりやすいので参考にされると良いと思います。70ページある大作ですが、そのうちの54ページ辺りから読んでみてください。 -------- *角度方向の次数と、Bessel関数の次数を無関係に決められないことは以下のように考えるとすぐ分かります。 2次元円形膜の振動を考え、仮に例えば半径方向の関数がa×J_0(x)、すなわちnに関し0次で、角度方向の関数がsin θ(n=1)だとします(aはある定数)。xは正規化半径に相当します。J_0(x)はx=0で有限のある正の値をとることに注意。 次に0よりわずかに大きい正の値δを考え、θ=0の点とθ=πの点とで位相を含めて振幅を比較します。x=δ、θ=0の点で振幅はa×J_0(δ)、x=δ、θ=πの点では位相がπだけ反転しますから振幅は-a×J_0(δ)になります。δはどんどん0に近づけて構わないわけですが、そうするとθ=0の側から膜の中心に近づいた場合とθ=πの側から近づいた場合では、膜の中心において振幅の不連続が発生することになります。これは直感的にもおかしいことが分かります。 J_0(x)は角度方向の関数が定数である場合(角度方向の次数も0)にのみ当てはまる、というわけです。

chestnut
質問者

お礼

まことにありがとうございました。 勉強いたします。 ほんとにありがとう!!!!!!

その他の回答 (1)

  • Umada
  • ベストアンサー率83% (1169/1405)
回答No.1

1次元の棒や弦の振動に関する微分方程式を解いた時、最後にその振動の様子を表すのにサインやコサインを使いますよね。 2次元の膜(縁が円であるもの)の振動も同様のやり方で解けてその振動の様子も表せるはずなんですが、サインやコサインでは表せないため、サインやコサインの代わりにBessel関数を使います。コサインやサインと同じように変数の変化につれて振動する関数ですが、振動の周期や振幅がだんだんと変化するところが違います。 なぜこんな関数と持ち出すのか? もちろんサインやコサインで用が足りればそちらの方が助かるのですが、2次元の円膜や円筒の振動はこの関数を使わないと表現できないからです。(むしろ、円膜の振動を解いた時に必然的に行き当たるこの一連の関数をBessel関数と名付けたという方が当たっていそうです) この関数のグラフは下記のページなどでごらんください。数式とにらめっこするよりグラフを見た方が早いです。 定義式についてはちょっとここで書くのは大変なので、微分方程式の教科書か理化学辞典などをご覧下さい。 「関数美術館」 http://www.geocities.co.jp/SiliconValley-PaloAlto/4076/scilab/art/ 「gnuplotで描いたサンプル」 http://funada11.denshi.numazu-ct.ac.jp/sei/sotsuken/samples/samples.html

参考URL:
http://www.geocities.co.jp/SiliconValley-PaloAlto/4076/scilab/art/, http://funada11.denshi.numazu-ct.ac.jp/sei/sotsuken/
chestnut
質問者

お礼

ありがとうございました。 本当に助かりました。 もうひとつだけ、いいでしょうか。 このいくつかの条件があるようですが、 これは虚数か何か関係しているのでしょうか。 振動の「モード」に当たるようなものでしょうか? よければ、教えていただけませんか? 簡単に書いてある本が見つけられず、本当に助かりました。 理解して勉強を進められそうです。ありがとうございました。

関連するQ&A

  • 平面波exp(-jx)をベッセル関数を用いてあらわすと・・・

    円筒波動関数について勉強しています。 今、平面波exp(-jx)を円筒座標系で表そうとしているのですが、その変換式が Σa*Jn(ρ)*exp(jnΦ) (Σはn=-∞~n=∞まで)と表されています。 aは定数、Jn(ρ)はn次の第一種ベッセル関数、ρは円筒座標系の原点から外に伸びていく変数、Φは円筒座標系のxy面上の角度 この式について、わからないことがあります。 なぜこの式がx方向に進む平面波を表すのでしょうか?定性的なことが理解できません。 しかも第一種ベッセル関数は進行波でないのに進行波をあらわしている。 このこともさらに混乱を深めています。 どのように理解すればよいのでしょうか? あまりベッセル関数に関する知識がないのでできれば優しくおねがしします。 勉強している本はR.F.Harringtonのtime-harmonic electromagnetic fields です。

  • ベッセル関数が難しく

    電磁気の勉強でベッセル関数が出てきました。はじめて見た関数だったため、この関数についていろいろと調べたのですが、どのような時に用いる関数なのか、この関数を用いるとどう便利なのか、そういったことが理解・イメージができません。 どうか、このベッセル関数に関しての、わかりやすい説明、あるいは参考URL等、ご教授いただけないでしょうか?

  • ベッセル関数

    周波数変調波の式を第一種ベッセル関数を用いて表わすにはどうすればいいのか分かりません。何かいいサイトなどあればよろしくお願いします。

  • ベッセル関数

    円筒座標系での電磁場のマクスウェル方程式を磁場に関して解いて得られる解が複素数を引数とする0次のベッセル関数 AJ0(kr)、kが複素数、Aは実係数、rは実変数 で得られるのですが 引数を実数に変換する方法がわかりません。 純虚数の引数であれば実数の引数の変形ベッセル関数に変換でき、 実数の引数であれば手持ちの本にベッセル関数の値が載っているのですが 複素数の引数の場合の処理方法がわからなくて困っています。 よろしくお願いします。

  • 円筒関数とは?

    ずっと円筒関数とベッセル関数は全く同じだと思っていたのですが、 手元にある特殊関数の本には ベッセル関数は円筒関数の一つの特殊な形である と書かれていますが、では一体円筒関数は一体どういう定義を持っているのでしょうか? ベッセル関数との違いなどについて教えて下さい。

  • 空洞共振器とベッセル関数の根

    空洞共振器の共振周波数を求める問題で躓いています。 直方体の共振器はようやく理解出来たつもりですが、 円筒形の共振器の計算で出てくるベッセル関数の根というのが全く解りません。 例えばTE011モードの場合、式の途中にベッセル関数の根y'01というのが出てきますが、 テキストを読み進みますとこの値は「2.4」となってます。 エクセルの関数にある「BESSELJ」で計算すると「BESSELJ(0.1)」は「0.76…」となるなど 全く異なります。 これとは別の関数なのでしょうか? どうやって求めたらよいか解らず、他のモードでの計算が出来ず困っています。 よろしくお願いします。

  • ベッセル関数と環状の膜の振動のモード

    ベッセル関数を調べています。 検索するとウィキペディアに説明文がありました。 その中の「応用」のところに 「環状の膜の振動のモード」と書かれてありました。 この「環状の膜の振動のモード」をベッセル関数で表すと どのような表現になりますでしょうか。 書籍、またはホームページを教えていただけたら助かります。 よろしくお願いします。

  • 解析学(ベッセル関数について)おねがいします!

    大学の解析学の授業でベッセル関数について課題が出され困っているところです。その内容ですがまずどんな関数か定義を書くのですがこれは本で調べたのでなんとなくわかったのですが次の直交関数系の証明、完全系かどうかの証明がまったくわかりません。解き方あるいは参考となるようなアドバイスを頂ければありがたいのですが・・・どなたかお願い致します!

  • ベッセル関数とcosの積分

    \int_{0}^{∞} J0(nx) cos(mx) dx という式、ここでJ0(nx)は0次のベッセル関数です。 山口勝也さんという方の「詳細微分方程式・特殊関数演習」という本のp244に上記の解が2通りかかれています。 1つは 1 / √(n^2 - m^2) もう1つは 0 しかしながら、その本にはそれぞれの解になる「条件」が明記されているようではないのですが、どういう条件でそれぞれの解になるかご存知でしょうか?

  • 固有値問題について・・・

    建築の構造解析をしているものです。最近この分野を勉強しはじめました。そこで、6次の行列の固有値を算出する段階です。そこで出てきた固有振動数、固有ベクトルのことなんですが、固有振動数が6個、固有ベクトルが6個ずつ出てきました。 そのような場合の固有振動数の意味合いはどういったものですか?できましたら詳しく教えていただけませんか?よろしくお願いいたします。