【Python】海洋潮汐与水位控制 回归分析法计算验潮站长期平均海面

一、原理以及算法步骤

1.1回归分析法

  回归分析法，也称线性关系最小二乘法，基本原理是假定两站的短期距平具有比例关系。设比例为k,则有：

 根据之前所学的关于距平的公式：

将（1.2）代入（1.1）中，得：

上式中k未知，需要进一步假设两站的长期平均海面之间为线性比例关系，比例系数仍然为k，则有：

上式中C为位置常数项。将上式带入式（1.3）中，整理得到短期平均海面之间的关系如下 

 对比式（1.4）和（1.5）可知，此时两站的长期平均海面与短期平均海面具有相同的线性关系。

1.2算法步骤

1. 构造观测方程组

  计算两站每天的日平均海面，设为MSLA，MSLB

,其中i=1,2,...,n。式（1.5）是估计k与C的观测方程，将每天的日平均海面代入（1.5）得到观测方程：

 n天的日平均海面序列构建观测方程组如下：

 根据间接平差的原理，观测方程组为：