不等时距潮汐资料调和分析.docxVIP

不等时距潮汐资料调和分析 在实际的海洋研究中，尤其是长期观测，由于仪器故障、恶劣天气、地理位置限制、观测方法等因素的影响，很难获得完整的数据集。在观测开始和结束之间的时间内获得的数据集的质量分析是不完整和不连续的。对于潮汐资料,调和分析方法作为一种主要的分析方法,发展至今,理论基础已经很成熟,但是传统的调和分析大多是对观测时间间隔为等距的资料进行调和分析,这是由于非等时序列数据的法方程系数矩阵很大,所需的计算量也非常大,通常要求观测数据等时间间隔以减少计算量,而且常用的观测资料也是以观测时间间隔为等距离为主。若资料中出现数据缺测漏测,或者因为异常天气和其他原因使得部分实测数据相对正常潮汐存在大的误差,就首先要对数据进行前期处理,然后应用最小二乘法对潮汐进行调和分析,这种方法在计算时,首先将原始矛盾方程进行求导处理,从而得到矛盾方程的法方程,这样大大增加了计算的难度。如果我们从线性代数这门学科的角度来考虑,矛盾方程组只是一个线性代数中的超定方程组,只需要直接解这个超定方程组即可,基于Matlab内部函数功能求解线性方程组能很好地解决这一问题,并且方程组所得解即为最小二乘解,从而简化了调和分析的算法。有多少样本,就有多少个方程,更能体现最小二乘法的灵活性,完全不受等时或非等时取样的限制。 本文采用基于Matlab内部函数功能的最小二乘法实现潮汐资料的调和分析,并应用于大连、北海两站位