簡介
克卜勒第三運動定律說明公轉週期的平方與行星到太陽的平均距離三次方成正比。若是行星以橢圓軌道繞行太陽，則行星到太陽的平均距離 ¯d 會是近日點距 r1 加上遠日點距 r2 的平均，也就是橢圓的半長軸，讓我十分好奇該如何證明這件事！
¯d=r1+r22=a
所要做的是就是將橢圓圓周上的每一點，到太陽的距離進行積分後取平均。
數學證明
定義連續函數f(x)在[a,b]區間上的平均值為：
fave=1b−a∫baf(x)dx
在一個以原點為中心的橢圓周上有一點P，其參數式可以表示成P(acosθ,bsinθ)，其中 a,b 為半長軸、半短軸。令太陽的座標為橢圓的其中一個焦點S(−c,0)，如下圖：
則橢圓周上每一點P到太陽S的距離可以表示為：
f(x)=¯PS=√(acos(x)+c)2+(bsin(x))2
如此一來就可以對f(θ)在[0,2π]區間內取平均值：
fave=12π∫2π0f(x)dx=12π∫2π0√(acos(x)+c)2+(bsin(x))2 dx=12π∫2π0√a2cos2(x)+2ac⋅cos(x)+(a2−c2)⋅sin2(x)+c2 dx=12π∫2π0√c2cos2(x)+2accos(x)+a2 dx=12π∫2π0(c⋅cos(x)+a) dx=12π[c∫2π0cos(x) dx+∫2π0a dx]=12π⋅(0+2πa)=a
得證橢圓軌道上地球到太陽的平均距離為半軸長！
程式實作
這次主要是以程式實作積分，透過將角度切成很小塊再疊加起來取平均，得到函數的平均值。以下影片中的橢圓a=13,b=5,c=12，透過將角度切成30度、15度、1度、0.5度觀察疊加的狀況。
角度 = 30度
角度 = 15度
角度 = 1度
角度 = 0.5度
以上四個影片放在原網址，可以看到每段長度加起來的平均值都是半長軸a！
