广义二项式与 y‘‘ = −y 的求解

旧事重提，温故而知新。

公式背得熟，$y^{\prime \prime }=-y$ 一看就是 $y=\sin x$ （或 $y=\cos x$），而 $y^{\prime }=y$ 则不由分说就是 $y=e^x$，照这个思路作级数展开，解方程求 $\pi$，代入式求 $e$ 就显得自然而无聊了。

说无聊是因为这只是在描述一个流程，而不是一个想法，我们无法从这个流程中获得更多的东西。这一切都是已知的，如果没有级数概念和相关运算法则，我们便一筹莫展。

但二项式展开不一样，它是一个朴素概念，在朴素二项式基础上，牛顿 “简单猜” 出了广义二项式展开式的形式，事后看，这不需求多少已知知识，只需要释放并有能力验算的大胆的想法。我前期写过一篇随笔 二项式的哲学和秩序 来描述牛顿的这个创举以及背后的形而上学含义。

广义二项式几乎可以拟合一切连续的，光滑的行为，这真是一个伟大的思想。所以只要我们形成诸如 $y^{\prime }=y，y^{\prime \prime }=-y$ 等动力学方程，就可以用二项式拟合 $y$，然后得到系数的通项，这本质上也是一种级数。

通过具象化的 $y^{\prime \prime }=-y$，我们知道有多种方法可以解出 $y$：

• 通过特征方程 $\alpha y^{\prime \prime }+\beta y^{\prime }+\gamma y=0$ 获得通解，但这属于事后论证，没特征方程认知的时代不好使；
• 数值方法，属于古代方法，通过大量数据探索 $y^{\prime \prime }$ 和 $y$ 的关系式，在当代计算机算力下，威力巨大；
• 二项式拟合，四海皆准的现代算术方法，无需任何前置知识和算力基础，本质上它就是一个级数方法；
• …

简单的二项式，插值是通用技法，所以还是得推一下万能的二项式拟合，代入二项式展开求系数，然后算出 $\pi$，这个在 动力学方程中的 π 和 e 中由于时间关系被省略。

设 $y$ 的二项式拟合式为：

$y(x)=a_0+a_{1}x+a_2{x}^2+a_3{x}^3+...{\sum }_{n=0}^{\infty }{a}_n{x}^n$

经 $y^{\prime }(x)$ 进而求得 $y^{\prime \prime }(x)$：

$y^{\prime \prime }(x)=2a_2+6a_{3}x+12a_4{x}^2+20a_5{x}^3+...$

代入 $y^{\prime \prime }=-y$ 并以初始值 $y(0)=a_0=0$ 启动递推：

$a_0=0，a_1=1，a_2=0，a_3=-\frac{1}{6}，a_4=0，a_5=\frac{1}{120}，...$

很容易确定阶乘关系，从而获得级数：

$\begin{array}{rl}y(x)& =x-\frac{x^3}{3!}+\frac{x^5}{5!}-\frac{x^7}{7!}+...\\ & =\sum _{n=0}^{\infty }(-1{)}^n\frac{x^{2n+1}}{(2n+1)!}\end{array}$

这是从 $y^{\prime \prime }=-y$ 这个具有具象化动力学含义的方程中推导出的结果，和 $y^{\prime }(x)=y(x)$ 自然增长的极限 $e$ 直接蕴含在 $x=0$ 的 $y(x)$ 中不同，很难从以上的 $y(x)$ 表达式获得自然常量，因为 $y^{\prime \prime }=-y$ 是个对称结构，控制变量的正是 $x$ 而不是 $y$，因此常量需要反解 $0,1,\cdots =y(x)$ 中的 $x$.

在另一条与动力学完全无关的路径上，数学家们早就知道了 $y=\sin x$ 的展开式就是上式：

$\sin x={\sum }_{n=0}^{\infty }(-1{)}^n\frac{x^{2n+1}}{(2n+1)!}$

由此可知，$y^{\prime \prime }=-y$ 的解就是 $y=\sin x$。

为了锁定对称常量，就需要求 $y=\arcsin x$ 在 $x$ 等于循环边界时 $y$ 的值。

将 $y=\sin x$ 的反函数 $y=\arcsin x$ 变换一下写法：

$x=\sin y$

同时对 $x$ 求导，根据链式法则，有：

$1=y’\cos y $

变形就是：

$y^{\prime }=\frac{1}{\cos y}$

接着处理右边，使之成为 $x$ 的表达式。根据 $\cos y=\sqrt{1-\sin y}，x=\sin y$，代入得：

$y^{\prime }=\frac{1}{\sqrt{1-x^2}}$

于是：

$(\arcsin x{)}^{\prime }=y^{\prime }=(1-x^2{)}^{-\frac{1}{2}}$

利用广义二项式定理展开：

$\begin{array}{rl}(\arcsin x{)}^{\prime }& =(1-x^2{)}^{-\frac{1}{2}}\\ & =1+\frac{1}{2}{x}^2+\frac{1}{2}\cdot \frac{3}{4}{x}^3+\frac{1}{2}\cdot \frac{3}{4}\cdot \frac{5}{6}{x}^4+...\\ & =\sum _{n=0}^{\infty }\frac{(2n-1)!!}{(2n)!!}{x}^{2n}\end{array}$

将其逐项积分就是 $\arcsin x$：

$\begin{array}{rl}\arcsin x& ={\int }_0^x(1+\frac{1}{2}{t}^2+\frac{1}{2}\cdot \frac{3}{4}{t}^3+\frac{1}{2}\cdot \frac{3}{4}\cdot \frac{5}{6}{t}^4+...)\text{d}t\\ & =x+\frac{1}{3}\cdot \frac{x^3}{3}+\frac{1}{3}\cdot \frac{3}{4}\cdot \frac{5}{6}\cdot \frac{x^7}{7}+...\\ & =\sum _{n=0}^{\infty }\frac{(2n-1)!!}{(2n)!!(2n+1)}{x}^{2n+1}\end{array}$

若不是古人就知道 $\pi$ 的存在，人们并不知道上面的级数蕴含着一个常量，可即使古人不知道 $\pi$，有了上面的级数和 $y=\sin x$ 的知识，人们也会很容易发现 $y$ 的循环节，卡住每个周期，$x$ 便是周期的整数倍，而这个周期就是 $\pi$。

代入 $x=1$ ，级数就收敛到一个常量，即 $\frac{\pi }{2}$，即使代入任何 $x$，均是常量，这就是 $y^{\prime \prime }=-y$ 刻画的行为的本质。然而若不是人们真的在古代就知道了 $\pi$，且确信三角函数一定蕴含 $\pi$，从 $y=\sin x$ 反求 $y=\arcsin x$ 的二项式展开的动机便不会那般强烈。

还是那句话，发现 $\pi$ 的方式有一万种，因为 $\pi$ 就在那里。

对比 莱布尼兹 “凑出” 的级数，我们发现它们何其相似，但由于 $x$ 的幂次衰减不够快，与莱布尼兹 $\pi$ 级数同样，$\pi =2{\sum }_{n=0}^{\infty }\frac{(2n-1)!!}{(2n)!!(2n+1)}$ 的收敛速度也是非常缓慢。

浙江温州皮鞋湿，下雨进水不会胖。