作者Honor1984 (希望愿望成真)
看板trans_math
标题Re: [积分] beta function
时间Sun Nov 8 13:27:57 2009
※ 引述《midarmyman (midarmyman)》之铭言:
: ∞ x^(p-1) π
: ∫---------------dx=-----------
: 0 1+x sinpπ
这个积分在所有复变课本内容或者习题中是一定会有的
这是基本复变定积分的范围
但是复变最常用在求"定"积分的求值上(优点)
却很难得出相对应的展开式(缺点)
微积分刚好相反
因为这板是讲微积分
我就用微积分的方法做出来
∞ x^(p-1)
∫---------------dx
0 1+x
1 x^(p-1) 1 x^(p-1)
= ∫--------dx + ∫---------dx
0 1+x 0 1+x
∞ (-1)^k (-1)^k
= Σ [--------- + ----------]
k=0 k-p+1 k+p
∞ (-1)^k
= Σ --------
k=-∞ p + k
到这里为止微积分的工作就结束了
这就是我上面说的
这是微积分的缺点
同时也是优点
这个级数同
1
= π----------
sin(pπ)
: 我换到三角函数变这样
: π/2 2
: 2∫(sinθ)^(2p-1)(1-sin θ)^(-p)dsinθ
: 0
: 然後就写不下去了@@
: 帮帮忙罗
你要换算成这个样子的目的是要使用Beta或者Gamma函数
但是除非你已经先把他们某些互相或者各自的关系式记下来
否则你只是写出他们的形式
并不能得到最後的closed form
只能表达定积分可用他们表示
因为会特别立出来Beta函数或者Gamma函数
除了他们的应用外
这些函数的求值一般微积分是没办法直接求出来的
但是既然书上将最初的积分化成这个样子
就表示它要使用Beta或者Gamma已知的关系来"求值"
但是Beta和Gamma函数有很多关系
你只要背其中几点
自然其他等价的关系都应该可以推导出来
其实你给的原积分式本身就是一种Beta函数B(p,1-p)
π/2 2
2∫(sinθ)^(2p-1)(1-sin θ)^(-p)dsinθ
0
t = (sinθ)^2
1
= ∫ t^(p-1) * (1-t)^(-p) dt
0
= B(p,1-p)
Γ(p)Γ(1-p)
= -------------
Γ(1)
π
= -------
sin(pπ)
也可以像你化简成那个三角函数的积分看出是B(p,1-p)
就看你关系式背多少
π
假如你只有背最基本的Beta和Gamma表达式以及Γ(x)Γ(1-x)= --------
sin(πx)
从你变数变换後的式子开始出翻
1 2
2∫(sinθ)^(2p-1)(1-sin θ)^(-p)dsinθ
0
π/2
= 2∫ (sinθ)^(2p-1) * (cosθ)^(-2p+1) dθ
0
∞π/2
2∫∫ r * exp(-r^2) * (sinθ)^(2p-1) * (cosθ)^(-2p+1) drdθ
0 0
= ---------------------------------------------------------------
∞
∫r * exp(-r^2) dr
0
(1/2)Γ(1-p)Γ(p)
= --------------------------
(1/2)
π
= -------
sin(pπ)
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 122.124.106.63
1F:推 G41271:为什麽Γ(1-p)Γ(p)=π/sinpπ 112.104.14.245 11/08 13:41
2F:→ G41271:我的意思是看过Γ(1-p)Γ(p)=π/sin(pπ) 112.104.14.245 11/08 13:41
3F:→ G41271:的证明是由∫x^(p-1)/(1+x)下推出来的 112.104.14.245 11/08 13:43
4F:→ G41271:这题换回Gamma是否有循环论证的毛病 112.104.14.245 11/08 13:44
5F:→ Honor1984:我文中已经说得很明白 你要先知道关系式 122.124.106.63 11/08 13:52
6F:→ Honor1984:用关系式求值 你把文字部分看一遍 122.124.106.63 11/08 13:53
7F:→ Honor1984:Γ(1-p)Γ(p)=π/sinpπ的推导可以由复 122.124.106.63 11/08 13:53
8F:→ Honor1984:变 或者单纯三角函数还有极限概念求得 122.124.106.63 11/08 13:54
9F:→ Honor1984:∫x^(p-1)/(1+x)和从Gamma乘积或者Beta 122.124.106.63 11/08 13:55
10F:→ Honor1984:出发都一样 这些都是等式 得到的答案会 122.124.106.63 11/08 13:56
11F:→ Honor1984:一样是合情合理的 122.124.106.63 11/08 13:56
12F:→ Honor1984:如果G41271看过除了我推文所指π/sinpπ 122.124.106.63 11/08 13:59
13F:→ Honor1984:证明过程的方式以外的方法 特别是光用积 122.124.106.63 11/08 14:00
14F:→ Honor1984:分技巧可以积出来的 麻烦请贴在板上 我 122.124.106.63 11/08 14:00
15F:→ Honor1984:也有兴趣知道^^ 122.124.106.63 11/08 14:00
16F:推 G41271:三角函数法?要怎麽推呢? 谢谢 112.104.14.245 11/08 14:01
17F:→ G41271:哦 我印象中只看过积分x^(p-1)/(1+x)的方法 112.104.14.245 11/08 14:03
18F:→ G41271:偏偏复变没学好,所以一直对此性质不太满意 112.104.14.245 11/08 14:05
19F:→ Honor1984:所以你指的积分是复变积分吗? 122.124.106.63 11/08 14:06
20F:→ Honor1984:另外一个方法是从sinx的展开式得到本文 122.124.106.63 11/08 14:13
21F:→ Honor1984:单纯微积分得到的级数相比差一倍数 但 122.124.106.63 11/08 14:14
22F:→ Honor1984:已经离题很多了 即使用复变 在本板也不 122.124.106.63 11/08 14:14
23F:→ Honor1984:是很恰当 122.124.106.63 11/08 14:15
24F:推 G41271:是的 此积分可用复变求得为π/sinpπ 112.104.14.245 11/08 14:17
25F:→ midarmyman:谢H大 140.117.198.78 11/08 20:37
26F:→ Honor1984:G41271 我在下一篇做出π/sin(pπ)的结140.109.103.151 11/09 14:02
27F:→ Honor1984:果 因为原积分式可通Γ(x)Γ(1-x) 所以140.109.103.151 11/09 14:03
28F:→ Honor1984:意思一样140.109.103.151 11/09 14:03