下のURLで出している円板の固有振動数の式のKmnの出し方を教えてください
http://www.math.ryukoku.ac.jp/~iida/lecture/gr/g …
No.3ベストアンサー
- 回答日時:
掲題のサイトの議論は以下のような疑問があり、ここでの質問に意味がないように
思います。
1.境界条件
詳しくないのですが、R(a)=0はありません。周囲はフリーです。
中心でとめる場合も固定していないようなので境界条件は|R(r)|<∞
しかない。したがって、R'(a)=0 だけになります(これも詳しくないが)。
2.まず、方程式はパラメータの正負によって次の2つの式がでます。
R''+R'/x+(1-n^2/x^2)R=0 (1)
R''+R'/x-(1+n^2/x^2)R=0 (2)
R''+R'/x-(n^2/x^2)R=0 (3)
3.(1)の解は R=AJn+BYn
(2)の解は R=AIn+BKn
(3)の解は R=Ar^n+Br^(-n)
掲題のサイトでは違うパラメータの方程式の解(1)(2)をごっちゃにしています。
4.r=0で有限という条件から
(1)の解は R=AJn
(2)の解は R=AIn
(3)の解は R=Ar^n
となり、R'(a)=0を使うと(2)(3)はA=0 となります(In'>0,r^(n-1)>0)。結局、解は
R=AJn
だけとなる。
5.以上の解でよければ求めるknmは(R'(a)=0を使うと)
Jn'(knm)=0
です。これは、#1のようにして求めることもでききますが、円形導波管の
電磁界の教科書に載っています。
参考URL:http://tpweb2.phys.konan-u.ac.jp/~susa/keisan_bu …
No.1
- 回答日時:
数値計算するしかありません。
フリーソフトMaximaの結果を示します。基のサイトはJm, kmn としていますが、以下ではJn, knmと考えています。
まず、g(n,x)=Jn'(x)/Jn(x)-In'(x)/In(x)
としてg(n,x)=0の解をニュートン法で求めます。
関数を変形したのは、(解が求めやすいと思って?)素直なグラフにするためです。
ベッセル関数の形は次のコマンドを実行すれば表示します。
wxplot2d([bessel_j(1,x),bessel_j(2,x)], [x,0,50]);
wxplot2d([bessel_i(1,x),bessel_i(2,x)], [x,0,5]);
1.まず、JnとInの微分は公式を使って、関数bj(n,x)とbi(n,x)を定義する。
bj(n,x):=(bessel_j(n-1,x)-bessel_j(n+1,x))/2;
bi(n,x):=(bessel_i(n+1,x)+bessel_i(n-1,x))/2;
次に計算関数として
g(n,x):=bj(n,x)/bessel_j(n,x)-bi(n,x)/bessel_i(n,x);
を定義する。つぎのコマンドで、グラフを描いて大体の零点をもとめ、
ニュートン法の開始値とする。
[wxplot2d([g(0,x),g(1,x),g(2,x)],[x,0,20],[y,-10,10],
[legend, "g(0,x)", "g(1,x)","g(2,x)"])]$
2.ニュートン法のパッケージをロードする。
load(mnewton);
3.変数に値をセット
[n:0, x1:3, x2:7, x3:10]$
ニュートン法の計算
[mnewton(g(n,x),x,x1),mnewton(g(n,x),x,x2),mnewton(g(n,x),x,x3)];
結果は次のように表示されます。
(%o39) [[[x=3.196220616582541]],[[x=6.306437047688424]],[[x=9.439499137876405]]]
すなわち、k01=3.196220616582541
k02=6.306437047688424
k03=9.439499137876405
となります。
同様に、次を実行すれば各値が求まります。
[n:1, x1:4, x2:7.5, x3:11]$
[mnewton(g(n,x),x,x1),mnewton(g(n,x),x,x2),mnewton(g(n,x),x,x3)];
(%o41) [[[x=4.610899879049056]],[[x=7.799273800811232]],[[x=10.9580671919195]]]
[n:2, x1:6, x2:9, x3:12.5]$
[mnewton(g(n,x),x,x1),mnewton(g(n,x),x,x2),mnewton(g(n,x),x,x3)];
(%o43) [[[x=5.905678235420522]],[[x=9.196882599635321]],[[x=12.40222096686439]]]
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
このQ&Aを見た人はこんなQ&Aも見ています
-
天使と悪魔選手権
悪魔がこんなささやきをしていたら、天使のあなたはなんと言って止めますか?
-
フォロワー20万人のアカウントであなたのあるあるを披露してみませんか?
あなたが普段思っている「これまだ誰も言ってなかったけど共感されるだろうな」というあるあるを教えてください
-
映画のエンドロール観る派?観ない派?
映画が終わった後、すぐに席を立って帰る方もちらほら見かけます。皆さんはエンドロールの最後まで観ていきますか?
-
海外旅行から帰ってきたら、まず何を食べる?
帰国して1番食べたくなるもの、食べたくなるだろうなと思うもの、皆さんはありますか?
-
天使と悪魔選手権
悪魔がこんなささやきをしていたら、天使のあなたはなんと言って止めますか?
-
円板の振動モードについて
その他(自然科学)
おすすめ情報
- ・漫画をレンタルでお得に読める!
- ・人生のプチ美学を教えてください!!
- ・10秒目をつむったら…
- ・あなたの習慣について教えてください!!
- ・牛、豚、鶏、どれか一つ食べられなくなるとしたら?
- ・【大喜利】【投稿~9/18】 おとぎ話『桃太郎』の知られざるエピソード
- ・街中で見かけて「グッときた人」の思い出
- ・「一気に最後まで読んだ」本、教えて下さい!
- ・幼稚園時代「何組」でしたか?
- ・激凹みから立ち直る方法
- ・1つだけ過去を変えられるとしたら?
- ・【あるあるbot連動企画】あるあるbotに投稿したけど採用されなかったあるある募集
- ・【あるあるbot連動企画】フォロワー20万人のアカウントであなたのあるあるを披露してみませんか?
- ・映画のエンドロール観る派?観ない派?
- ・海外旅行から帰ってきたら、まず何を食べる?
- ・誕生日にもらった意外なもの
- ・天使と悪魔選手権
- ・ちょっと先の未来クイズ第2問
- ・【大喜利】【投稿~9/7】 ロボットの住む世界で流行ってる罰ゲームとは?
- ・推しミネラルウォーターはありますか?
- ・都道府県穴埋めゲーム
- ・この人頭いいなと思ったエピソード
- ・準・究極の選択
このQ&Aを見た人がよく見るQ&A
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
物理の有効数字2桁
-
整数で答えよ、のように書かれ...
-
この問題分かる人解き方教えて...
-
この問題の正の向きと負飲む木...
-
FEMのモード法と直接法の違...
-
単振動の微分方程式の解き方に...
-
高さHのビルの屋上から初速v₀で...
-
ダイアモンドの結晶構造因子に...
-
大学4年生です。 卒業に124単位...
-
未だに卒論の準備してないのは...
-
大学生です。7000~10000字のレ...
-
卒論が最低20000字何ですが、引...
-
大学の卒論の不合格認定について
-
【大学のレポートの書き方】「...
-
レポートの字数に参考文献は含む?
-
卒論提出せずに卒業してしまい...
-
ゼミ必修ではないのでゼミを辞...
-
大学のゼミの1回目のの講義で5...
-
論文などで怒られるのは皆さん...
-
【卒論】大学生が書いた卒業論...
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
物理の有効数字2桁
-
bcc(面心立方体)の四面体位置で...
-
整数で答えよ、のように書かれ...
-
FEMのモード法と直接法の違...
-
pHについてです
-
慣性モーメント
-
音波の振動数の求め方を教えて...
-
水素原子の波動関数の直交性に...
-
1リットルは1辺が1cmの立方体が...
-
斜面を上る運動 摩擦なし
-
保存力→ポテンシャル→力 の計...
-
たたみ込み積分って名前つける...
-
変数分離が成功したからといっ...
-
この問題の有効数字の合わせ方...
-
高校の理科の答えの書き方
-
円板の固有振動数
-
てこの計算問題 輪軸の応用問題
-
全微分のdx,dyの意味
-
写真の問題の解き方を調べたと...
-
なぜ写真①の場合と②の場合で自...
おすすめ情報