ガウスの微分方程式
No.1の回答は途中で間違えたので、書き直して、決定方程式まで書いた。式(3)以外の変換の計算が完了してないが一応投稿する。
f(x)=x (1−x)y ″ +(γ−(α + β +1)x)y ′ − α β y = 0__(1)
最初にあなたの計算した変換が最も簡単なので、(2)(3)の変数変換を行う。
x=ξ^(−1) __(2)
y=η__(3)
(2)をξで微分すると(4)になるので、No.1で書いたように、xで微分する演算子d/dxを(5)とする。
dx/dξ=−1/(ξ^2) __(4)
d/dx=(1/(dx /dξ))(d/dξ)= (−ξ^2)d/dξ__(5)
(5)の演算子を使って、(3)を微分してy ′とy ″を計算すると、(6)(7)となる。
y ′=dη/dx= (−ξ^2)dη/dξ__(6)
y ″=(−ξ^2)(d/dξ)(−ξ^2)dη/dξ=(ξ^2)(d/dξ)((ξ^2)dη/dξ)
=(ξ^2)((2ξ)dη/dξ+(ξ^2)d^2η/dξ^2)__(7)
(1)に(2)(3)(6)(7)を入れて、ξとηの微分方程式(8)にする。
f(x)=x (1−x)y ″ +(γ−(α + β +1)x)y ′ − α β y = 0
=(1/ξ)(1−1/ξ)(ξ^2)((2ξ)dη/dξ+(ξ^2)d^2η/dξ^2)+(γ−(α + β +1)(1/ξ))(−ξ^2)dη/dξ− α βη
=(ξ−1)((2ξ)dη/dξ+(ξ^2)d^2η/dξ^2+(γξ−(α + β +1))(−ξ)dη/dξ− α βη= 0__(8)
式(8)に(−1)をかけて2階微分,1階,0階の順に並べると− f(x)の方程式(9)になる。
− f(x)= (1−ξ)((2ξ)dη/dξ+(ξ^2)d^2η/dξ^2)+(γξ−(α + β +1))ξdη/dξ+ α βη__
= (1−ξ)((ξ^2)d^2η/dξ^2)+(1−ξ)(2ξ)dη/dξ+(γξ−(α + β +1))ξdη/dξ+ α βη
= (1−ξ)((ξ^2)d^2η/dξ^2)+(2−2ξ)ξdη/dξ+(γξ−(α + β +1))ξdη/dξ+ α βη
= (1−ξ)((ξ^2)d^2η/dξ^2)+{(γ−2)ξ−(α + β−1)}ξdη/dξ+ α βη=− f(x)__(9)
(9)はあなたの計算と一致するが、分数は級数解を出すのに邪魔なので、1−ξをかけて分母をはらう。
ηの級数解Σ(λ=0→∞)anξ^(κ+λ)の最低次の項はa0ξ^κで、最低次だけ計算するために
η= a0ξ^κ__(10)__とすると(11)が成立する。
dη/dξ= a0κξ^(κ−1),d^2η/dξ^2)= a0κ(κ−1)ξ^(κ−2)__(11)__(9)に(10)(11)を入れると
a0(1−ξ)((ξ^2)κ(κ−1)ξ^(κ−2))+a0{(γ−2)ξ−(α + β−1)}ξκξ^(κ−1)+a0α βξ^κ
= a0[(1−ξ)κ(κ−1)(ξ^κ)+{(γ−2)ξ−(α + β−1)}κ(ξ^κ)+ α β(ξ^κ)]=0__(12)
式(12)の括弧を外すと、(ξ^κ)と(ξ^(κ+1))の項が出るが、最低次数だけを計算するために
(ξ^(κ+1))の項は無視すると、(13)になる。a0の係数が0なら、解の初期条件a0を自由にきめても
式(13)が成り立つ。無視した項はa1の項を加えて消せる。a1,a2,・・・の係数を漸化式から決めて(9)を成立させる。(ξ^κ)で割って(14)を得る。
a0[κ(κ−1)(ξ^κ)−(α + β−1)}κ(ξ^κ)+ α β(ξ^κ)]=0__(13)
κ(κ−1)−(α + β−1)}κ+ α β=κ^ 2−κ−(α + β−1)}κ+ α β=(κ−α)(κ− β)=0__(14)
(14)が決定方程式で、根はαとβである。
No.2ベストアンサー
- 回答日時:
回答が付かないのと、教えてグーからの締め切りの催促で
質問を締め切ってしまい、ご迷惑をおかけしたことをお詫びいたします。
丁寧に書いていただき有難うございます。
x=ξ^(−1) __(2)
y=η__(3)
この変換だと、
(κ−α)(κ− β)=0__(14)
となって、
α
β
が根になります。
その通りだと思います。
私の計算では、
吉田耕作さんの変換
x=ξ^(−1) __(2)
y=(ξ^(−α))η__(3)
だと
0
ーα+β
が根になります。
また
寺沢さんの変換
x=ξ^(−1) __(2)
y=(ξ^(α))η__(3)
だと
2α
α+β
が根になります。
これらの根の間の関係は
αーα=0
βーα=-α+β
となり、
α+α=2α
β+α=α+β
となっている。
わたしは、吉田さんの変換も寺沢さんの変換も
根が α、β だと書いてあるのに、そうはならないので間違いだと思っています。
とりあえず、積分方程式論の48ページについては勝手に修正して
今は、96ページを読んでいます。
他のミスプリント20個(?)も勝手に修正していますが、
48ページの記述はミスプリントとはいえないと思ったので
質問させていただきました。
丁寧に教えていただき、有難うございました。
また、質問する機会があろうかと思いますので、
これからもご指導よろしくお願いします。
間違った回答のままでは、他の人に見られても恥ずかしい。修正が届いてうれしいです。
御研鑽のご様子を伺いました。Good Luck!
No.1
- 回答日時:
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
似たような質問が見つかりました
- 物理学 ポテンシャルが有限で不連続の時、右側の波動関数をφ1(x)、左側をφ2(x)とする。境界条件の「波動 2 2023/06/04 13:53
- 数学 x=r・cosθの2回微分 θ=ωtとすると? 5 2022/05/10 23:53
- 数学 連立微分方程式の解き方について 7 2022/12/16 13:39
- 数学 (1+x^2)y'=1 の微分で教えて下さい 2 2022/08/30 10:23
- 数学 【完全微分方程式⠀】 分数で分母が0になり定義できない場合、分母を仮にtと置いてそれを極限t→0とし 1 2022/05/06 14:43
- 数学 微分について教えてください 放物線y=x^2のx=1における微分係数を定義に従って求め、その点におけ 5 2023/04/16 15:38
- 高校 数学III 積分 数学IIIの積分でf(ax+b)の積分公式がありますが b=0の時どのように考えれ 4 2022/09/30 02:06
- 数学 多変数関数の微分とテイラー展開について 5 2022/04/24 16:55
- 数学 『Cの微分.2』 3 2023/02/15 19:47
- 数学 【全微分について】 z=f(x,y) の全微分は df=(∂f/∂x)dx+(∂f/∂y)dy と表 1 2023/02/25 05:49
おすすめ情報
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
50以下は“50”も入るのですか?
-
5進法を10進法への直し方
-
デジタル信号処理の問題です。...
-
Excel 16進数
-
10進数の50を2進数で表すといく...
-
HEX2BIN関数の使い方。
-
二進数の循環小数を十進数に変...
-
16進小数0.Cを10進数小数に変換...
-
デシベルから加速度の変換について
-
8進数から16進数 16進数から8進数
-
2進数の0.101101101101・・・...
-
進数の小数点をどう変換するのか
-
偏微分の記号をタイプするため...
-
dBm→dBμV/mの換算について
-
相似変換とユニタリ変換
-
ラプラス逆変換
-
相関関係 クロス集計 散布図
-
2進数
-
ヤコビアンが0になってしまう場...
-
EXCELで10進数表記をB...
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
50以下は“50”も入るのですか?
-
5進法を10進法への直し方
-
16進小数0.Cを10進数小数に変換...
-
Excel 16進数
-
HEX2BIN関数の使い方。
-
8進数から16進数 16進数から8進数
-
1分45秒75で289,995円稼ぐA君が...
-
「じじょう」が正しい読み方?
-
偏微分の記号をタイプするため...
-
dBm/HzからdBm/MHzへの単位変換
-
n進法→m進法への変換
-
EXCELで10進数表記をB...
-
dBm→dBμV/mの換算について
-
10進数の50を2進数で表すといく...
-
ヤコビアンが0になってしまう場...
-
算数計算 大至急お願いします
-
ACアダプターの消費電力の件
-
幾何と代数は同じ数学でしょうか
-
フーリエ変換・逆変換の虚数成...
-
二進数の循環小数を十進数に変...
おすすめ情報
質問者:uyama33
質問日時:2018/02/17 08:27
からの質問が間違えた回答のまま締め切られてしまった。
上記のように修正したい。