我在这篇《潮汐的科普与计算（更新）》中重新计算了潮汐高度随角度的变化。计算思路是，先计算潮汐力，也就是外天体引力在地心外一点的引力与在地心处引力的差，然后计算这个潮汐的势。由于水平面是等势的（否则水会从高势流向低势），整个水面的形状也就计算出来了。

$\frac{\vec{f}(\vec{r})}{Gm} =\frac{r}{D^3} \left(2\cos\theta \hat{x} - \sin\theta\hat{y}\right)$

单位径向矢量为

$\hat{r} = \cos\theta \hspace{1mm} \hat{ x} + \sin\theta \hspace{1mm} \hat{y}$

单位切向矢量为

$\hat{\theta} = - \sin\theta \hspace{1mm} \hat{ x} + \cos\theta \hspace{1mm} \hat{y}$

潮汐势：

$V(\vec{r}) =- \int_{0}^{r} \vec{f} \cdot d\vec{r} = - Gm\frac{2 \cos^2\theta - \sin^2\theta}{D^3}\int_0^r r \hspace{1mm}dr \\ = -Gm\frac{r^2}{4D^3} \left(3\cos 2\theta+1\right)$

$H_{tide}(\vec{R}) =\frac{m R^4}{4M_e D^3} (3\cos 2\theta+1)$

其中是造成潮汐的天体质量，为地球质量，为地球半径，为地球到天体的距离，是该地点与地心的连线与地心--天体连线的夹角。由此公式计算出月球产生的最大潮汐高度为40厘米，太阳引发的潮汐高度为16厘米。

其随位置（角度）变化图示如下

2. 现在我们用徐教授的概念进行计算

由于潮汐力的切向分量（也就是海面的水平分量），在力的方向水要高出一点，这样高出的水的部分产生的压力才能抵消水平潮汐力。考虑一块面积为 A、厚度 为 $R d\theta$ 的水 （因为地球半径远大于水深、扇梯形可近似为长方形），令水涨 dh 压力抵消潮汐力，我们有（左边是水块受到 潮汐力，右边是初中物理的水压强公式乘以面积）：

$\rho A R d\theta f_{\theta} = A \rho g dh$

两边的水的密度rho 与面积A消去了。因此，

$dh = R f_{\theta} d\theta /g$

因此，从90度到 theta 度，水总共涨高量为：

$\Delta H= \int dh = R/g \int_{\pi/2}^\theta f_{\theta} d\theta \\ =\frac{Gm R^2}{gD^3} \int_{\pi/2}^\theta (-2\cos\theta\sin\theta - \sin\theta\cos\theta) d\theta\\ = \frac{3mR^4}{4M D^3} (\cos 2\theta+1)$

显然，我的计算与用徐教授的概念进行计算，两者得出的水面高度差完全一致：比如0度处与90度处高度差，两个结果都是  3/2 。不同的是，用势的计算中得出了水的绝对高度变化，而不只是高度差。

为什么两个计算会得到同样的结果呢？ 稍微想想发现道理非常简单。在我的潮汐势计算中，是直接从地心到地表一点计算，一条直路。但我也可以走一条弯路：先从地心垂直往上到地球表面，然后沿着表面转到角度的一点，两条路计算出来的势应该是一样的。这就是所谓保守力，势只与位置有关与路径无关。但在后一条路中，第一部分垂直向上潮汐力做功是固定的，与角度无关，全部与角度有关的贡献来自后半部分切向力的贡献。换言之，以90度处为零点沿着地球表面计算切向潮汐力的势，与以地心为零点用径向潮汐力计算势，两条路是一样的，只是势能零点不同。

殊途同归！