• 締切済み

ベッセル関数について

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のせいで式がベッセル方程式の形にならず、困っています。よろしくお願いします。

みんなの回答

  • f272
  • ベストアンサー率46% (8017/17134)
回答No.2

dz/dx=2x=2√z y'=dy/dx=(dy/dz)(dz/dx)=2√z*Y' y"= (d/dx)(dy/dx) = (d/dx)(2√z*Y') = (2√z)*(d/dz)(2√z*Y') = (2√z)(1/√z*Y')+(2√z)(2√z*Y'') = 2*Y'+4z*Y'' x^2*y''+x*y'+(4*x^4-1/4)*y=0 z*(2*Y'+4z*Y'')+x*(2√z*Y')+(4*z^2-1/4)*y=0 4z^2*Y''+4z*Y'+(4*z^2-1/4)*y=0 z^2*Y''+z*Y'+(z^2-1/16)*y=0

  • f272
  • ベストアンサー率46% (8017/17134)
回答No.1

yをzで1回微分したものをY'とすると、 z^2*Y''+1/2*z*Y'+(z^2-1/16)*y=0 にはならないよ。 z^2*Y''+z*Y'+(z^2-1/16)*y=0 になります。 どこかで計算を間違っている。 y'=2√z*Y' y''=4zY''+2Y' になっていますか?

anonymous_kibou
質問者

補足

すみません y'=2√z*Y'までは求められるのですが どうしてもy''=4zY''+2Y'となる計算がわかりません

関連するQ&A

  • ベッセルの微分方程式について

    ベッセルの微分方程式は標準形で、 x^2(d^2f/dx^2)+x(df/dx)+(x^2-n^2)f = 0 となっています。ある種の物理現象を常微分方程式化してこの形に帰着させると、あとはベッセル関数などを使った定型処理に持ち込めるので、何とかこの形に持ち込むまで頑張るということになると思います。 質問ですが、x^2(d^2f/dx^2)+x(df/dx)+ x c f = 0 という形になった場合、これはベッセルの微分方程式ではない、ということになるでしょうか。左辺最終項の形が標準形と違います。ある文献で”ベッセルの微分方程式になる”と断言されているのですが、標準形にならないので思案しています。例えば、独立変数を変換して式も変換すると標準形になるのでしょうか。よろしくお願いします。 よく、ベッセル、ラゲールなどの微分方程式は級数解を使いますが、いろんな現象から式を変形して標準形に本当にきれいに帰着できるものでしょうか。ちょっと違うとか亜種があるのかなと思いますが。 よろしくお願いします。

  • 変数係数2階線形微分方程式

    変数係数2階線形微分方程式の問題です。 x^2*y(x)''+2x*y(x)'-iαy(x)=0 i:複素数,α:定数 この微分方程式はどのようにして解けばよろしいでしょうか? できるだけ計算過程を詳しくお願いします。 解にはベッセル関数が用いられるみたいです。 自分でベッセルの微分方程式と同様にして解いていっても途中でつまずいてしまいます。 お手数ですがよろしくお願いしたします。

  • ベッセルの微分方程式

    テキストによると、円筒座標系での電磁場のマクスウェル方程式を磁場に関して解いて得られる方程式が f’’+1/x*f’+k^2*f=0 解はベッセル関数 AJ0(kx)+BY0(kx) A,Bは定数 しかしこの方程式は一般的なベッセルの微分方程式と少し違います。 x^2f’’+xf’+x^2f=0 x^2で割り算してるのはともかく、係数kの分だけ違うのです。これでもベッセルの微分方程式であり解はベッセル関数であると言えるのでしょうか?

  • ベッセルの方程式の問題の解き方が分かりません

     次のベッセルの方程式の問題の解き方が分かりません。  数学に詳しい方、よろしければご教示願えないでしょうか。 問題は、  ベッセルの方程式に帰着できるさまざまな方程式がある。示されている置換を 使って、次の微分方程式の一般解を求めよ。 4*x^2*y" + 4*x*y' + (x - ν^2)*y = 0 (√x = z)  このように解いてみました。  ベッセルの微分方程式は、 x^2*y" + x*y' + (x^2 - ν^2)*y = 0 で、  一般解は、 y(x) = A*Jν(x) + B*Yν(x) ここで、A と Bは任意定数、Jν(x)は第1種ベッセル関数、Yν(x)は第2種ベッセル 関数。 √x = z より、 dz/dx = 1 / (2*√x) y'とy"は、 y' = dy/dx = (dy/dz)*(dz/dx) = (dy/dz)/(2*√x) y" = d^2y/dx^2 = (d/dx)*(dy/dx) = (d/dz)/(2*√x)*(dy/dz)/(2*√x) = (d^2y/dz^2)/(4*x) ゆえに、 4*x^2*y" + 4*x*y' + (x - ν^2)*y = 4*x^2*(d^2y/dz^2)/(4*x) + 4*x*(dy/dz)/(2*√x) + (x - ν^2)*y = x*(d^2y/dz^2) + 2*√x*(dy/dz) + (x - ν^2)*y = z^2*(d^2y/dz^2) + 2*z*(dy/dz) + (z^2 - ν^2)*y = 0 となって、第 2項目が z*(dy/dz) にならず、2*z*(dy/dz) になってしまいます。  本の回答をみると、 A*Jν(√x) + B*Yν(√x) となっているので、問題の微分方程式を、 z^2*(d^2y/dz^2) + z*(dy/dz) + (z^2 - ν^2)*y = 0 に変形したのだと思いますが、どのようにすれば良いのでしょうか ?  同様に下記の問題も、 x^2*y" + x*y' + 4*(x^4 - ν^2)*y = 0 (x^2 = z) 同じ解き方をしたため、 z^2*(d^2y/dz^2) + z*(dy/dz) + (z^2 - ν^2)*y = 0 に変形できませんでした。  なにとぞよろしくお願いします。

  • ベッセル関数

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

  • ベッセル関数

    ベッセル関数の問題で困っています。y"+(x^2)*(y)=0で公式y=C1Jr(x)+C2Yr(x)を使いたいのですが、これはx^2y"+xy'+(x^2-y^2)y=0の解に使用できる公式だと思うのですが、どのようにこの形に持っていけばよいのでしょうか?ヒントにはy=u√(x), z=(1/2)*x^2と置けばよいと書かれています。このような問題が他にもあり困っています。どなたかできる方がいましたらこの問題を例にご解法よろしくお願いします。

  • ベッセル関数とノイマン関数の部分で微分方程式の形が分かると、解がすぐに

    ベッセル関数とノイマン関数の部分で微分方程式の形が分かると、解がすぐに出るものがあるそうです。 添付の式1はワイリーの「工業数学 上」p361にある定理1の系1というものらしいです。 これに対して式2が求めたいもので、式1と比べて係数を解決するそうですが、λの二乗の部分から a,b,sの求め方がいまいち理解できません。式3からどう解くのかのヒントをお願いいたします。

  • 変数係数の微分方程式の解き方

    変数係数斉次線形微分方程式y''+P(x)y'+Q(x)y=0の特殊解または、基本解系はどのように求めることができるのでしょうか。

  • ベッセル関数の微分

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

  • 大学院入試の微分方程式の問題がわかりません!

    問題の式を書くとややこしいので画像を添付しました。 【初期条件: y(0)=y0,y'(0)=y1】 画像の微分方程式について (1) 変数変換 u=( x^2 + 2 )y を行って、uに関する微分方程式を導け (2) (1)で導いた微分方程式を解くことで、元の微分方程式の解yを求めよ (3) 【x→∞】lim y(x)を計算せよ また、【x→-∞】lim y(x)が存在するためのy0,y1の条件を求めよ (1)の変数変換を行うときに uを微分してu' u'' を出し それらをy y' y'' の式に直して代入すればできると思うのですが その変形がややこしすぎて何回やっても間違えてしまいます そこで知識ある皆様のお力をお貸しいただければと思い質問しました。 何卒よろしくお願い致します。