计算机编程/物理/加速物体的位移函数

维基文库:源代码

加速物体的位移可以用数学函数 $\mathbf {s} (t)$ 描述。可以使用泰勒级数 获得广义函数

$\mathbf {s} (t+t_{0})=\sum {\frac {t^{n}}{n!}}\mathbf {s} ^{(n)}(t_{0})$,

其中 $\mathbf {s} ^{(n)}$ 是 $nth$ 阶导数 $\mathbf {s}$

• $\mathbf {s} ^{(0)}={\frac {d^{0}\mathbf {s} }{dt^{0}}}=\mathbf {s}$
• $\mathbf {s} ^{(1)}={\frac {d^{1}\mathbf {s} }{dt^{1}}}=\mathbf {v}$
• $\mathbf {s} ^{(2)}={\frac {d^{2}\mathbf {s} }{dt^{2}}}=\mathbf {a}$
• 等等。

此函数的精度取决于所用项数，因为 ${\frac {1}{n!}}$ 迅速减小。此外，时间 $t$ 可以同步，使得 $t_{0}=0$ （麦克劳林级数）。

注意，对于恒定加速度，大多数项将变为零，只剩下

$s(t)={\frac {1}{0!}}\mathbf {s} _{0}+{\frac {t}{1!}}\mathbf {s} _{0}^{(1)}+{\frac {t^{2}}{2!}}\mathbf {s} _{0}^{(2)}$

或者

$\mathbf {s} (t)=\mathbf {s} _{0}+\mathbf {v} _{0}t+{\frac {1}{2}}\mathbf {a} t^{2}$

C++

template<class Vector,class Number>
Vector PositionAcceleratingBody(Vector *s0,Number t,size_t Accuracy)
{
     Vector s(0);     //set to zero if int, float, etc. or invoke the
                      //     "set to zero" constructor for a class
     Number factor(1);//0!==1 and t^0==1
     for(size_t n(0);n<Accuracy;n++)
     {
          if(n)factor*=(t/n);//0!==1 and t^0==1
          s+=(factor*s0[n]); //s0 is the array of nth derivatives of s
                             //     at t=t0=0
     }
     return s;
}

使用泰勒级数的理由

泰勒级数可以通过系统地选择哪些变量是常数，然后将其外推到无限极限来推导。

• 恒定位置

$s(t)=s_{0}$

或者

$s(t)={\frac {t^{0}}{0!}}s^{(0)}(t_{0}=0)$

• 恒定速度

$v(t)=v_{0}$

$s(t)=\int vdt=\int v_{0}dt$

$s(t)=v_{0}t+s_{0}$

或者

$s(t)={\frac {t^{0}}{0!}}s^{(0)}(t_{0}=0)+{\frac {t^{1}}{1!}}s^{(1)}(t_{0}=0)$

• 恒定加速度

$a(t)=a_{0}$

$v(t)=\int adt=\int a_{0}dt$

$v(t)=a_{0}t+v_{0}$

$s(t)=\int vdt=\int (a_{0}t+v_{0})dt$

$s(t)={\frac {1}{2}}a_{0}t^{2}+v_{0}t+s_{0}$

或者

$s(t)={\frac {t^{0}}{0!}}s^{(0)}(t_{0}=0)+{\frac {t^{1}}{1!}}s^{(1)}(t_{0}=0)+{\frac {t^{2}}{2!}}s^{(2)}(t_{0}=0)$

• 加速度的变化率恒定

$a^{(1)}(t)=a_{0}^{(1)}$

$a(t)=\int a^{(1)}dt=\int a_{0}^{(1)}dt$

$a(t)=a_{0}^{(1)}t+a_{0}$

$v(t)=\int adt=\int (a_{0}^{(1)}t+a_{0})dt$

$v(t)={\frac {1}{2}}a_{0}^{(1)}t^{2}+a_{0}t+v_{0}$

$s(t)=\int vdt=\int ({\frac {1}{2}}a_{0}^{(1)}t^{2}+a_{0}t+v_{0})dt$

$s(t)={\frac {1}{6}}a_{0}^{(1)}t^{3}+{\frac {1}{2}}a_{0}t^{2}+v_{0}t+s_{0}$

或者

$s(t)={\frac {t^{0}}{0!}}s^{(0)}(t_{0}=0)+{\frac {t^{1}}{1!}}s^{(1)}(t_{0}=0)+{\frac {t^{2}}{2!}}s^{(2)}(t_{0}=0)+{\frac {t^{3}}{3!}}s^{(3)}(t_{0}=0)$

• 等等。