標籤:

數學擺(1):周期與角度的關係

數學擺(1):周期與角度的關係

令振幅(即與豎直平面偏轉的角度)為 	heta,對重力進行分解。弦對單擺的拉力計為 T

根據牛頓第二定律,有:

mlfrac{d^2	heta}{dt^2}=-mgsin	heta

約去單擺質量 m 後得:lfrac{d^2	heta}{dt^2}=-gsin	heta

整理得: frac{d^2	heta}{dt^2}+frac{g}{l}sin	heta=0

小角度近似:

	heta 很小時,sin	heta≈	heta ,代入上式,解得: 	heta(t)=Acos(omega t+varphi)A 為振幅,omega 為角頻率, varphi 為相位角。

	heta 求導兩次,得 omega=sqrt{frac{g}{l}} ,根據定義, T=frac{2pi}{omega}=2pisqrt{frac{l}{g}}

非小角度單擺:

現在需要對二階微分方程 frac{d^2	heta}{dt^2}+frac{g}{l}sin	heta=0 進行求解。

為了將二階微分方程轉換成一階微分方程,引入 omega

由於 omega=frac{d	heta}{dt} ,

frac{d^2	heta}{dt^2}=frac{domega}{dt}=omegafrac{domega}{d	heta}

帶入方程, omegafrac{domega}{d	heta}=-frac{g}{l}sin	heta

分離變數得: omega domega=-frac{g}{l}sin	heta d	heta

兩邊同時積分: intomega domega=frac{g}{l}int-sin	heta d	heta

該微分方程的通解為: omega^2=2frac{g}{l}cos	heta+C

t=0 時,令 	heta=alpha ,此時 omega=0

Rightarrow omega^2=2frac{g}{l}(cos	heta-cosalpha)=4frac{g}{l}left[ sin^2 (frac{alpha}{2})-sin^2(frac{	heta}{2}) 
ight]

因為 dt=frac{d	heta}{omega} ,對兩邊同時積分可解出所需時間。考慮 t=frac{1}{4}T ,此時應從0積到最大的位移角,也就是:

frac{1}{4}T=t=int_{0}^{alpha}frac{d	heta}{omega}=frac{1}{2}sqrt{frac{l}{g}}int_{0}^{alpha}frac{d	heta}{sqrt{sin^2(alpha/2)-sin^2(	heta/2)}}

接下來計算下面的積分:

int_{0}^{alpha}frac{d	heta}{sqrt{sin^2(alpha/2)-sin^2(	heta/2)}}

考慮換元:令 sinpsi=frac{sin(	heta/2)}{sin(alpha/2)}

frac{d	heta}{dpsi}=frac{2cos(psi)sin(alpha/2)}{cos(	heta/2)} =frac{2cos(psi)sin(alpha/2)}{sqrt{1-sin^2(psi)sin^2(alpha/2)}}

換元得: int_{0}^{alpha}frac{d	heta}{sqrt{sin^2(alpha/2)-sin^2(	heta/2)}}=int_{0}^{pi/2}frac{frac{d	heta}{dpsi}}{sqrt{sin^2(alpha/2)-sin^2(	heta/2)}}dpsi

則該積分變為:

int_{0}^{pi/2}frac{frac{2cos(psi)sin(alpha/2)}{sqrt{1-sin^2(psi)sin^2(alpha/2)}}}{sqrt{sin^2(alpha/2)-sin^2(alpha/2)sin^2(psi)}}dpsi=2int_{0}^{pi/2}frac{dpsi}{sqrt{1-sin^2(psi)sin^2(alpha/2)}}

查表得該積分為一類橢圓積分。

參考文獻:

The Pendulum?

galileo.phys.virginia.edu圖標
推薦閱讀:

《術數》與傳統數學
關於複數的乘方/開方運算
聽說為什麼數學家會發明"行列式"這種東西呢?會火是因為家庭背景很厲害?
Data Analysis
為什麼有些人的素質這麼差?

TAG:數學 |