以下の積分を不完全ベータ関数を使用し超幾何関数に変換してで表すことができるようです。
導き出す手順を教えてください。
∫{ (x-b)^(a-1)}/x dx
a, b は実数です。
(参考)
http://oshiete.goo.ne.jp/qa/8067505.html
(解)
自力で計算できなかったので、Wolfram Matheaticaのサイトで積分したところ答えが
{ (x-b)^a/ab }{ 1 - (1-b/x)^-a 2F1(-a, -a; 1-a; b/x) }
とガウスの超幾何級数を用いて算出されました。
導き方がわかりません。
A 回答 (4件)
- 最新から表示
- 回答順に表示
No.4
- 回答日時:
クンマー変換は、ひとつの変換ではなくて、
x → 1-x と x → 1/x が生成する 6 個の変換からなる
変換群 { x, 1/x, 1-x, 1/(1-x), (x-1)/x, x/(x-1) }
のこと。これが何をしているかというと、
2F1 は 0, 1, ∞ に特異点を持つ複素関数なので、
一次分数変換で、三個の特異点を移し換えている。
それぞれの変換を施した微分方程式を解くと、
x(1-x)(d^2y/dx^2)+{γ-(α+β+1)x}(dy/dx)+αβy=0
に対して、解の 6 個の表示
2F1(α,β,γ; x),
{x^(1-γ)}・2F1(α-γ+1,β-γ+1,2-γ; x),
2F1(α,β,α+β-γ+1; 1-x),
{x^(γ-α-β)}・2F1(γ-α,γ-β,γ-α-β+1; 1-x),
{x^(-α)}・2F1(α,α-γ+1,α-β+1; 1/x),
{x^(-β)}・2F1(β,β-γ+1,β-α+1; 1/x)
が出るが、ここから 2F1 をマクローリン展開して得られる
6 個の級数は、上から 2 個ずつ
x=0, x=1, x=∞ の近傍でのみ収束し、
そのペアで、各領域での解空間の基底となる。
(解析接続すると、一次独立なものは、2 個。)
…受け売り。私の計算は、正直、追いついていません。
"クンマーの変換" よりも "クンマーの関数等式" で
検索したほうがよかったかも。
No.3
- 回答日時:
"超幾何微分方程式 クンマーの変換" で google.
ご回答ありがとうございます。
クンマーの変換というものがあるのですね。また一つ勉強になりました。
ご教示いただいたもののクンマーの変換は
F(α,β,γ:x)=(1-x)^(γ-αーβ)F(γ-α,γ-β,γ:x)
F(α,β,γ:x)=(1-x)^(-α)F(α,γ-β,γ:x/(1-x))
でしょうか?
この方程式を使用して頑張って変換しようとしましたが、
2F1(1-a, 0, 1; x/b)→2F1(-a, -a; 1-a; b/x)
とはいきませんね(^^;
少なくとも最後の引数 x/b→b/x と逆数になるような変換はできそうにありません。
引き続きアドバイス頂ければ幸いです。
No.2
- 回答日時:
x(x-b)(d^2y/dx^2)+{(-a+2)x-b}(dy/dx)=0
を x=bz で置換すれば、
z(1-z)(d^2y/dz^2)+{1+(a-2)z}(dy/dz)=0
となる。これって、
z(1-z)(d^2y/dz^2)+{γ+(a-2)z}(dy/dz)+αβy=0,
α=1-a, β=0, γ=1
ってことじゃない?
y = (c_1) 2F1(1-a, 0, 1; x/b)
でいいような気がするけど…
ご回答ありがとうございます。
すいません。自分の思慮が足りなく、助かります。
質問ばかりで申し訳ありませんが本結果から、
{ (x-b)^a/ab }{ 1 - (1-b/x)^-a 2F1(-a, -a; 1-a; b/x) }
の解になるには、少なくとも 2F1(1-a, 0, 1; x/b)→2F1(-a, -a; 1-a; b/x)
の変換が必要になると思いますが、式をいじくってもできそうにありません。
どのようにすれば解を導けるでしょうか。
ご教示いただければ幸いです。
No.1
- 回答日時:
y = ∫{ (x-b)^(a-1)}/x }dx
なら、
y' = (x-b)^(a-1)/x,
y'' = (x-b)^(a-2){(a-1)x-(x-b)}/x^2.
両方式から (x-b)^(a-2) を消去すると、
x(x-b)y'' + {(-a+2)x-b}y' + 0y = 0.
これって、超幾何微分方程式でしょ。
変数変換でうまく係数を調整したら、
ガウス型にならないかな。
ご回答ありがとうございます。
ガウスの微分方程式
x(1-x)y"+{γ-(α+β+1)x}y'+αβy=0
ですね。
この場合の一般解はy=c_1 2F1(α,β,γ;x)+c_2 x^(1-γ) 2F1(α+1-γ,β+1-γ,2-γ;x)
ですね。
bが1なら良いのですが、、、がんばってみましたが変数変換できそうにありません。
試しにy=Σ[k=0-∞] c_k x^(r+k)とおいて解を探そうとしましたが、複雑になりすぎて
私には導けませんでした。
うまい変数変換の方法がありましたらご教示ください。
不完全β積分を使って超幾何関数へ一般化する方法を前回ご提案頂きましたが、
定積分範囲が0→1で定義されているβ関数からΓ関数へ変換するところで躓いてしまいました。
もしこちらから導ける場合は、その方法でも良いです。
引き続きご教示いただければ幸いです。
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
似たような質問が見つかりました
- 物理学 ベクトルと座標系につきまして 1 2022/04/03 06:23
- 物理学 『数か物か』 4 2022/06/13 06:54
- 小学校 算数の問題で悩んでいます。 2つの数A,Bを四捨五入して整数の概数にすると、順に25と3になりました 5 2023/08/21 15:05
- 数学 ある方から 「一応 「①の被積分関数の 1/(z'-z) の部分を以下のように、等比級数 の公式を使 3 2023/02/16 05:30
- 数学 「f(z)=1/(z^2-1)に関して ローラン展開を使う場合、マクローリン展開を使う場合、テイラー 3 2022/08/27 19:56
- 数学 複素関数にロピタルの定理を使おうとしている回答者は、複素関数論はおろか微積分学もよく分かっていない、 5 2022/12/28 18:02
- 数学 「違います 質問11 n≦-2ではz≠π/2で g(z)=tan(z)/(z-π/2)^(n+1) 3 2022/07/16 18:12
- 物理学 力学の運動方程式につきまして 4 2023/07/17 14:43
- Excel(エクセル) エクセル・スプレッドシートで、一定数を超えたらゼロから再累計する方法 8 2022/05/28 03:52
- C言語・C++・C# C言語 3 2022/10/04 15:07
おすすめ情報
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
5進法を10進法への直し方
-
50以下は“50”も入るのですか?
-
デジタル信号処理の問題です。...
-
10進数の50を2進数で表すといく...
-
偏微分の記号をタイプするため...
-
二進数の循環小数を十進数に変...
-
Excel 16進数
-
HEX2BIN関数の使い方。
-
デシベルから加速度の変換について
-
EXCELで10進数表記をB...
-
dBm→dBμV/mの換算について
-
8進数から16進数 16進数から8進数
-
dBm/HzからdBm/MHzへの単位変換
-
相似変換とユニタリ変換
-
16進小数0.Cを10進数小数に変換...
-
インチをセンチに換算
-
2進数
-
マイナスの進数変換
-
10進小数→2進小数、16進小数が...
-
フーリエ変換・逆変換の虚数成...
マンスリーランキングこのカテゴリの人気マンスリー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アダプターの消費電力の件
-
幾何と代数は同じ数学でしょうか
-
フーリエ変換・逆変換の虚数成...
-
二進数の循環小数を十進数に変...
おすすめ情報