• ベストアンサー

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

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

  • erara
  • お礼率37% (189/500)

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

  • ベストアンサー
回答No.3

1.問題点が絡み合っています。まず、  BESSELJ(x,n)は普通の第一種ベッセル関数です。ベッセル関数の根y'01=2.4 に対しては  BESSELJ(2.4,0)が0になればよいのですが、当然、誤差のため0に近い値となるだけです。  というわけで、まず、引数の与え方が間違っています。 2.つぎに、円筒形の共振器(導波管でも同じ)で出てくるTE波の境界条件はベッセル関数を微分  したものになります。だから慣習として、根y'01に、「'」がついているのです。  この時の根は、「3.832」です。根「2.405」はTM波の境界条件となります。  どこかで、何かが混乱したのでしょう。 3.なお、ベッセル関数の微分はあまり出っ張っていないので、公式 Jm’=(Jm-1 - Jm+1)/2  を使うとよいでしょう。また、根はamn のようにあらわされており、BESSELJ(amn,m)=0  となるので、混乱しないように注意してください。 4.さらに、計算にはフリーソフトmaximaを使うとよいです。マセマテイカ相当の超強力な  ツールです。マニュアルもあります。  日本語訳  http://maxima.sourceforge.jp/maxima.html

erara
質問者

お礼

皆様、ご回答ありがとうございます。 勉強不足、知識不足が甚だしく、質問自体にもミスがあったりと申し訳ありません。 とりあえず、目の前の問題を解くの必要なだけのベッセル関数の根だけは 計算することが出来ました。 (自分でではなく機械がしてくれたわけなので、理解は不十分ですが…) 教えて頂いたソフトについても順に勉強して行こうと思います。 本当に助かりました。

その他の回答 (2)

  • tadys
  • ベストアンサー率40% (856/2135)
回答No.2

2.4÷0.76=3.158… ⇒ π です。 どこかにπが出てくるのでしょう。

  • spring135
  • ベストアンサー率44% (1487/3332)
回答No.1

ベッセル関数は多くの種類があって、そのうちのどれかを確認する必要があります。さらに指数が整数のものはなじみ安いですが非整数次のものは目が回ります。 空洞共振器の解を与えるベッセル関数はJとは異なります。 電磁気の基礎方程式から導かれるベッセル方程式がどれなのかを四季を追いながら注意深く確認してください。

関連するQ&A

  • ベッセル関数って、

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

  • ベッセル関数

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

  • FDTDを用いた共振モードの検出

    円筒共振器の共振モードを3次元のFDTD(Meep)を使って求めようとしています. いくつかの共振モードは見つけることができたのですが(TE030,040,410,420),解析的にはTE510,320,610・・・ともっと多くの共振モードがあるはずなのですが, 見つけることができません. 解析的に求まる共振モードをすべてFDTDで見つけるのは不可能なのでしょうか? それとも,なにか工夫をしてFDTDを行えば,共振モードはみつかるものなのでしょうか? 漠然とした質問ですいません. 要は,FDTDで解析的にでてくる共振モードをすべて見つけられるか,ということです. (補足) 解析的というのは教科書のようなものにのっている式から手計算で求めたということです.

  • ベッセル関数

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

  • 空洞共振器と温度の関係について

    空洞共振器は温度の影響を受けてQ値が変動すると聞きました。 これはなぜでしょうか? 共振器の温度係数τfと関係がありますでしょうか? 温度係数に関する式は τf=1/f0・⊿f(φD, L)/⊿T f、f0は共振周波数、⊿fは共振周波数の変動量、φDは共振器の径、Lは共振器の軸長です。 温度係数についての式は下記サイトから引用しました。 http://www.spc.co.jp/spc/pdf/giho19_06.pdf ご回答宜しくお願い致します。

  • ベッセル関数の近似式

    VBAを使っていて、ベッセル関数を使いたいのですが、 worksheetfunctionを使わずに計算したいので、近似式を探しています。 第一種0次ベッセル関数に関してはこのページで見つけることができました。 http://soudan1.biglobe.ne.jp/qa5670519.html 他のベッセル関数や特殊関数の近似式が書かれてある webページがあれば教えていただけないでしょうか?

  • ベッセル関数の計算

    ベッセル関数の計算 νについてnへの極限をとった際に、添付の式のように求まるというらしいのですが、 どういう計算で下の式になっているのか分かりません。 漸化式などを使っているのでしょうか? よろしくお願いいたします。

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

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

  • ベッセル関数の微分

    ある膜の振動を数値計算で解こうとしているのですが分からないところがあるので質問します。 計算の途中ベッセル関数J_ν、Y_νをνで微分しなければならないのですが、微分した後導関数がどんな形になるか分からないです。νは実数です。 どなたか分かる方がいらっしゃれば宜しくお願いします。

  • ベッセル関数について

    x^2*y''+x*y'+(4*x^4-1/4)*y=0の一般解をz=x^2の変数変換を行って求めてください。 yをzで1回微分したものをY'とすると、 z^2*Y''+1/2*z*Y'+(z^2-1/16)*y=0 ここまでは計算できるのですが、1/2*z*Y'の項の係数1/2のせいで式がベッセル方程式の形にならず、困っています。よろしくお願いします。