非常接近的情况下,可以通过对其进行级数展开来近似计算相互作用力。而且不止是经典n体问题,同样可以引申到相对论性n提问题。”
听到乔泽的话,彼得舒尔茨也忍不住站了起来,凑到了乔泽身边开始看起了他在稿纸上的演算。
“这是椭圆模型?”
“对,先预设一个三体问题,将三体系统的每个物体的位置表示为椭圆函数的解。”
说完,乔泽又在手稿上写下了三个公式。
[x_ita_icosomega_it+phi_i,]
[y_itb_isinomega_it+phi_i,]
[z_itc_icosomega_it+phi_i,]
然后开口解释道:“其中a_i,b_i,c_i分别是椭圆的半长轴、半短轴和半高轴,omega_i是椭圆的角频率,phi_i是初始相位。”
洛特杜根先是露出恍然的神色,随后又皱着眉头问道:“但这如何影响相互作用力的计算?”
“通过级数展开来逼近相互作用力,比如我们先考虑物体i和j之间的引力,那么定义相互作用力为”
说话间,乔泽又在手稿上写出一串公式。
[mathbff_ij-gfracm_im_j|mathbfr_i-mathbfrj|^2hatmathbfrij,]。
然后说道:“g是引力常数无需解释,m_i,m_j分别是物体i和j的质量,hatmathbfr_ijmathbfr_j-mathbfr_i|mathbfr_j-mathbfr_i|是单位矢量。”
“这怎么做级数展开?超越几何学还涉及到力学的计算?”
“超越几何学中允许使用逐项逼近技术,可以用于做级数展开,具体可得”
[frac1|mathbfr_i-mathbfrj|sumk0^inftyfracpsi_kr^k+1,]
“对了,psi_k就是系数。”
看着手稿上开始逐渐变得丰满的公式,彼得舒尔茨突然感觉人有些不好了,拧着眉头说道:“不对,这样会有级数展开的截断误差,这个误差是不可控的吧?”
“有办法的,你忘了今天我讲解论文的时候是怎么解决的吗?设置一个截断参数n,仅考虑级数展开的前n项。只要n的值够大,