潮汐调和分析方法的探讨张凤烨.pdfVIP

潮汐调和分析方法的探讨 1,2 1,2 1,2 1,2 1,2 张凤烨 , 魏泽勋 , 王新怡 , 王永刚 , 方国洪 (1. 海洋环境科学与数值模拟国家海洋局重点实验室, 山东 青岛 266061; 2. 国家海洋局 第一海洋研究所, 山东 青岛266061) 摘要: 为实现对不等时距潮汐资料的分析, 基于 Matlab 内部函数功能, 提出了一种调和分析方法。基 于这种方法, 分别对大连、北海两个站位 1985 年的全年等时间间距取样的资料和非等时间间距取样的 资料进行了调和分析, 结果显示, 由等时间距资料和非等时间距资料计算的调和常数基本吻合。对大 连、北海两个站位的全年资料进行多个不同时间间距取样分析, 发现当分潮频率大于取样频率的二分 之一时, 分潮发生频率混淆。若分潮周期明显大于样品长度, 该分潮的分析结果产生很大误差。最后 得出结论为: 此调和分析方法, 适合对非等时间间距、非连续潮汐潮流资料进行调和分析, 并且能够获 得与传统方法精度相当的结果。 关键词: Matlab; 潮汐; 调和分析; 取样 中图分类号: P731.23 文献标识码: A 文章编号: 1000-3096(2011)06-0068-08 在实际的海洋调查中, 尤其是长时间的观测, 本文采用基于 Matlab 内部函数功能的最小二乘 由于受到仪器故障、恶劣天气、地理位置制约, 观测 法实现潮汐资料的调和分析, 并应用于大连、北海两 方式等因素的影响, 很难得到从观测初始时刻到结 站位的水位观测资料分析, 通过对两站位全年资料 束时刻这段时间内完整的高质量数据资料, 也可以 的调和分析, 从非理论推导的角度, 分析基于 说获得的海洋资料是不完整的、不连续的。对于潮 Matlab 内部函数功能的调和分析方法是否可应用于 汐资料, 调和分析方法作为一种主要的分析方法, 等时间间隔和不等时间间隔的观测资料的调和分析, 发展至今, 理论基础已经很成熟, 但是传统的调和 以及观测资料的取样间隔对于调和分析精度的影响, 分析大多是对观测时间间隔为等距的资料进行调和 然后通过对两站位 7 月份资料的调和分析, 来分析 分析, 这是由于非等时序列数据的法方程系数矩阵 采样长度对潮汐调和分析精度的影响。 很大, 所需的计算量也非常大, 通常要求观测数据 1 改进的潮汐调和分析方法介绍 等时间间隔以减少计算量, 而且常用的观测资料也 [5] 是以观测时间间隔为等距离为主。若资料中出现数 调和分析是建立在分潮的概念上 , 将潮汐看 据缺测漏测, 或者因为异常天气和其他原因使得部 成是以不同频率传播的各种潮波叠加产生的现象, 分实测数据相对正常潮汐存在大的误差, 就首先要 调和分析的目的就是求出各个分潮的振幅和迟角, 对数据进行前期处理[1], 然后应用最小二乘法对潮 迟角是指: 某时刻、某一地点实际的分潮的相角与理 汐进行调和分析, 这种方法在计算时, 首先将原始 论上该时刻的分潮相角的差值, 迟角分为3 种: 有区 矛盾方程进行求导处理, 从而得到矛盾方程的法方 时迟角