这篇文章主要是整理一下我的胡思乱想。
首先假设地球为球体,只讨论水平(也就是沿切面)运动。
原则上非水平运动可以分为地轴(而不是垂直于切面)与水平两个分量,地轴分量对科里奥利力无影响。这可以由叉积的分配律以及平行向量的叉积为
直接得来。
为了推导,建立右手坐标系。以球心为原点,地轴为
轴(球心至北极为正方向),赤道面为
平面。
地理上用纬度
与经度
表示球面上的点。当然,地理上还使用相对高度,这样能够以三个参数表示三维空间上的点。这里不使用相对高度,而使用距离原点的长度
。当然,这个
不一定是地球的半径,也就是可以表示任意高度的点(只不过这时要看成地球的同心球)。
根据坐标系的建立方法,可知北纬为正,南纬为负;东经为正,西经为负。这与 GIS 的习惯也是一致的。
于是可以用向量
作为球上某点的径矢。
由最初假定地球为球体,得该点切面的法向量
。
用
表示地球自转角速度的大小。地球自西向东自转,所以
为正。又地轴为 z 轴,得地球自转的角速度
。
设质点的速度
,沿切面方向。故
。
那么科里奥利力
。
将科里奥利力分为沿切面的水平地转偏向力
与沿法向垂直地转偏向力
,现在我们要推导两者的表达式。
由投影的几何定义得知,
为矩形。故一个向量到某平面的投影向量与到平面法向量的投影向量之和为该向量。
容易知道,
在某个平面的法向量
的投影向量
,在平面的投影向量
。这就解决了上面的问题。
还可以使用代数方法推出。
设
,则
。
由
,得