|
|
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ 阳光男孩羊: 金币+8, ★有帮助 2015-08-15 10:48:21
要计算数值,就要先把不定积分改为定积分:Integral{dφ/sqrt[Sin(φ)] , α0, α}=C'- β*t
其中为直线的初始位置的角度;α为t时刻直线的位置的角度。当α=π/2时直线就转到了X轴的位置。
由 t=0时α=α0,故代入方程中后得:C'=0
故:Integral{dφ/sqrt[Sin(φ)] , α0, α}=- β*t
当=π/2时,t=T
T=-1/β*Integral{dφ/sqrt[Sin(φ)] , α0, π/2}
现在用高斯勒让德数值计分方法计算。
令φ=(α0+π/2)/2+(π/2-α0)/2*ξ ,dφ=(π/2-α0)/2*dξ
原积分=-1/β*Integral{(π/2-α0)/2*dξ/Sin[(α0+π/2)/2+(π/2-α0)/2*ξ] ,-1 , 1}
=-(π/2-α0)/(2*β)*Integral{dξ/Sin[(α0+π/2)/2+(π/2-α0)/2*ξ] ,-1 , 1}
=-(π/2-α0)/(2*β)*Integral{dξ/f(ξ) ,-1 , 1}
=-(π/2-α0)/(2*β)*Sum{δk*1/f(ηk), k=1~n}
其中,n为高斯点数量,根据需要选取。原则是计算精度要求高时可多取,反之可少取。δk为插值点ηk处对应的被积函数的值的系数。
比如,n=1时,η=0,δ=2;若n=2,则η1=-1/sqrt(3),δ1=1,η2=1/sqrt(3),δ2=1 ;n=3时,......。具体可参考下面链接:
http://wenku.baidu.com/link?url= ... RZmLWeZFFALlBgoKTvC |
|