简介TMD 2.05 是运行于 MATLAB 环境的专业潮汐分析与预测工具箱面向海洋科学研究者、物理海洋学专业学生以及从事港口规划、海岸工程和海洋灾害预警的工程师。该工具箱通过图形用户界面与谐波分析、潮流成分解析等模型支持潮汐成分识别、长期与短期水位预测以及环境影响评估。压缩包内共 82 个文件主体为 71 个 MATLAB 脚本涵盖主程序入口、潮汐系数提取、插值、椭圆参数计算及各环节工具函数另含模型配置数据、坐标文件与全球潮汐模型数据整体结构清晰、开箱即用压缩包仅 379KB轻量易部署。已有 1012 人下载学习。解压并配置 MATLAB 路径后用户可直接调用工具箱完成从数据导入、参数设置到结果可视化与预测的全流程操作也可参考内部脚本学习潮汐调和分析与模型建立的实现思路。1. TMD2.05一套被低估的潮汐预测工具箱我在处理中国近海某站点的逐时潮位预报时最初想自己实现FES2014的调用结果被数据格式和球谐展开折腾了两周。后来换成TMD2.05这个zip包配合TPXO8模型数据当天就出了可用的时间序列。这个压缩包里没有花哨的界面核心是一个个独立的MATLAB函数从文件名就能看出整个潮汐计算的流程读取模型、提取调和常数、预测水位、计算潮流椭圆。它解决的主要问题是把全球潮汐模型高效地插值到任意经纬度并输出工程上需要的潮位或流速序列。适合物理海洋、港口航道、海洋环境评估方向的研究生和工程师也适合想快速验证潮汐模型效果的MATLAB用户。2. 函数族与数据文件TMD2.05的文件结构决定了调用顺序解压TMD2.05.zip后你会看到一组.m文件、两个数据目录Model_load7.2、tpxo8_atlas30和若干经纬度文本。这个结构不是随便堆的它对应了潮汐预测的三个阶段读取模型、空间插值、时间预测。理解文件之间的调用关系比直接运行TMD.m重要得多因为很多报错都源于调用顺序和参数格式不对。2.1 核心函数族从tmd_tide_pred到TMD.mtmd_tide_pred.m是使用频率最高的预测入口它的输入是模型句柄、经纬度、时间数组输出是该点的天文潮位时间序列。TMD.m更像是主控脚本负责组织整个流程包括加载模型、设置坐标、调用预测函数并绘图。tmd_tide_pred_mapts.m则是在tmd_tide_pred基础上扩展出的多测点版本适合批量计算。我通常不会直接改tmd_tide_pred而是把它当作黑盒重点看它依赖的函数tmd_get_coeff.m负责从模型文件中读取分潮系数BLinterp.m完成双线性空间插值nodal.m和astrol.m计算天文幅角和节点因子这四个环节是潮汐预测的数学核心。再看配合文件TMD_check_lat_lon.m和TMD_check_date.m会在传入非法经纬度或日期时直接报错拦掉大部分低级错误。虽然这些函数看起来只是工具但少了它们实际预测结果会出现不可预期的NaN或跳变。% 典型调用链子模型 - 时间轴 - 逐时潮位 model tmd_mk_submodel(DATA/Model_tpxo8_atlas30, ...); % 构建子模型 tlist datenum(2024,1,1):1/24:datenum(2024,1,31); % 小时分辨率时间轴 Z tmd_tide_pred(model, 30.0, 122.0, tlist); % 预测潮位这里的model是包含网格坐标、水深和分潮复振幅的结构体lat/lon用十进制度数tlist使用MATLAB的datenum序列。返回的Z单位是米。要注意的是tmd_tide_pred内部可能期望的时间基准是MJDModified Julian Date所以直接传datenum有时会报错或产生偏移这时候就需要date_mjd.m作为中转。2.2 模型数据目录Model_load7.2与tpxo8_atlas30解压后出现Model_load7.2和tpxo8_atlas30两个文件夹它们存的是不同来源的网格数据。Model_load7.2可能是某个区域模型的二进制文件而tpxo8_atlas30是TPXO8的30弧分全球模型。tmd_mk_submodel.m用来把这些原始网格裁剪成局部区域裁剪后生成的文件会小很多加载速度快也便于分发。tmd_submodel_atlas.m和tmd_submodel_atlas_compact.m是不同压缩策略的实现compact版本会减少内存占用但牺牲部分精度。lat_lon.dat以及lat_lon、lat_lon_time、lat_lon_0、lat_lon_1、lat_lon_2这些文件我理解是测绘坐标点文件用于测试或者定义兴趣区。tmp.mat、ltmp这类带临时字样的文件往往是脚本运行过程中生成的缓存变量不影响主流程可以直接忽略。真正要关心的是OUT目录因为子模型默认会输出到这里。2.3 支持函数的功能速查表很多人在论坛上卡住就是因为不知道某些函数的存在。下表是从TMD2.05的函数堆里筛出来的高频实用项。函数/文件作用使用场景tmd_extract_HC.m提取某点的调和常数振幅/迟角站点潮汐特征分析tmd_get_bathy.m提取水深校核模型网格tmd_ellipse.m由M2/K1等分潮流计算椭圆参数潮流分析与余流场剥离tmd_mk_submodel.m制作子模型局部区域快速计算tmd_timeOnOff.m控制时间开关时变模型处理rdModFile.m读取OTIS格式模型文件自定义数据源时使用grd_in.m / xy_ll.m网格坐标与经纬度互换插值和绘图前坐标系转换tmd_changeCaxis.m调整绘图色标结果出图InferMinor.m推算次要分潮浅水分潮补偿这个速查表对应了完整工作流先用tmd_get_coeff读取系数再用tmd_extract_HC提取某站点的调和常数最后用tmd_tide_pred做时间预测。每个函数前都有检查逻辑比如TMD_check_lat_lon.m会在经纬度超出网格时提示避免你拿着一个空矩阵空算半小时。3. 构建子模型并完成第一次潮汐预测这一步是整个过程里最需要谨慎的。直接把全球模型加载进内存不是不行但每次插值都访问大文件会很慢所以习惯先对目标海域做一个子模型。子模型既能保留主要分潮信息又能把计算区域限制在兴趣范围内调试时也更容易排查问题。3.1 解压、路径设置与模型选择解压zip后把TMD主目录和DATA目录加入MATLAB路径。比较稳妥的是直接在TMD根目录运行addpath(genpath(pwd));这会把当前目录及所有子目录加入搜索路径避免因为嵌套目录导致的函数找不到问题。然后检查是否存在Model_tpxo8_atlas30文件夹这是TPXO8的全球网格分辨率为30弧分覆盖范围广适合作为子模型的源数据。如果手头还有更高分辨率的区域网格比如Model_load7.2也可以替换但要注意tmd_get_coeff.m是否能识别对应格式。3.2 用tmd_mk_submodel裁剪局部网格假设要研究东海某区域经度122125°E纬度2832°N。常见做法是src DATA/Model_tpxo8_atlas30; % 源模型文件 dst OUT/east_china_model.mat; % 输出子模型 lon_range [122, 125]; lat_range [28, 32]; model tmd_mk_submodel(src, lon_range, lat_range, dst);src指向模型目录或NetCDF文件dst是输出的.mat文件后两个参数控制裁剪范围。执行后生成的model结构体里包含网格坐标、水深和分潮复振幅矩阵。注意裁剪范围不要跨越±180°经线否则需要先做坐标偏移TMD_check_lat_lon.m会拦掉这种输入。裁剪时还有一个隐藏坑边界网格点的插值依赖外围一个网格所以裁出来的区域实际比最外层略小plot后你会发现边缘有空白带这是正常的。3.3 生成预测时间序列有了model就可以预测任意点的潮位。我最常用的是逐小时序列便于后续调和分析对比lat 30.5; lon 123.5; t0 datenum(2024,3,1); tf datenum(2024,3,31); tlist t0:1/24:tf; [time, tide] tmd_tide_pred(model, lat, lon, tlist); plot(time, tide);注意tmd_tide_pred返回的时间可能是UTC或本地时区需要结合TMD_cor_date.m做一个时间校正否则和实测站数据比对时会差几个小时。这一点我多次在项目里栽过。TMD_cor_date.m的输入参数通常包含原始时间数组和时区偏移量具体以函数注释为准。另一个容易忽略的是水深点与潮位点的区别海洋模型中有些网格点是干出点潮间带水深为负值tmd_tide_pred在这些点上的插值结果可能异常需要用到h_in_p.m和h_in.m的掩膜逻辑。3.4 从子模型到全模型哪种预测更可靠当研究区跨多个时区或近岸地形复杂时子模型的分辨率可能不够。这时直接使用原始模型文件而不裁剪tmd_tide_pred在每次调用时会从模型网格上做双线性插值精度更好但耗时明显增加。我通常在方案论证阶段用子模型快速出图正式计算时改用完整模型。如果想兼顾速度与精度可以尝试tmd_get_coeff.m只提取少数几个主要分潮比如只用M2、S2、K1、O1这四个忽略浅水分潮这样预测序列呈现的误差会稍大但趋势和峰值相位仍然可信。浅水区需要额外加入InferMinor.m推算M4、MS4等分潮否则高潮位会偏低。4. 从预测结果到调和常数验证、可视化与排错预测只是第一步我更关心某个站点的调和常数因为M2、S2、K1、O1等分潮的振幅和相位决定了这个站的潮汐性质也是后续工程报审需要的数据。TMD2.05提供了直接从模型网格提取调和常数的接口不需要自己做谱分析。4.1 用tmd_extract_HC提取分潮复振幅TMD内部把分潮的振幅和相位存储成复数值实部为余弦分量虚部为正弦分量。提取某点调和常数的常见写法[HC, names] tmd_extract_HC(model, 30.5, 123.5); % HC是一个复值向量对应constit.m定义的分潮顺序 for i 1:length(names) A abs(HC(i)); ph angle(HC(i)) * 180/pi; fprintf(%s: 振幅 %.3f m, 相角 %.2f deg\n, names{i}, A, ph); end这里abs得到振幅angle得到格林尼治迟角。names列表来自constit.m它定义了工具箱内使用的分潮列表常见的M2、S2、K1、O1都在其中。注意相位的参考时间是模型使用的历元和实测调和分析得到的迟角不一定直接可比需要做天文幅角归算。如果要输出给水动力模型做开边界最好保留原始的复数形式而不是振幅相位格式这样对方模型读取时可以避免相位单位转换产生的误差。4.2 用实测数据校验模型预测拿到站点实测逐时潮位后把TMD预测值减去实测值差值中高位部分可能是非潮汐成分。我常用下面的方式obs load(obs_30.5_123.5.txt); % [mjd, height] pre interp1(time, tide, obs(:,1)); % 对齐时间 err pre - obs(:,2); std(err) % 看均方根误差如果std超过0.3 m就需要检查是否接近河口或浅水区那里的非线性分潮如M4往往没被全球模型覆盖。另外还要检查时间和水位基准面是否一致。TMD默认输出的是平均海面MSL上的潮位而中国沿海的工程基准多为理论最低潮面两者可能相差几厘米到几十厘米不能直接拿来做工程高程换算。4.3 常见报错与处理下表是我实际调试中遇到的高频问题以及对应的排查路径。报错关键词原因处理办法lon/lat out of range经纬度超出网格范围用TMD_check_lat_lon.m检查边界并裁剪模型model not loaded未读取模型句柄改用tmd_get_coeff先加载模型结构体date format error时间格式非MATLAB日期统一用datenum必要时转MJDNaN in bathymetry水深数据缺测用smooth.m和getridofNaNss.m预处理subscript indices error子模型边界内缩扩大裁剪范围1个网格最隐蔽的是时间格式问题。tmd_tide_pred内部可能要求MJD直接传datenum有时得到错误时间轴。TMD_InFileFormat.m里有详细说明但很多人不会打开看。我建议先读这个文件里面有输入参数的格式说明比在GUI里乱点高效得多。5. 进阶潮流椭圆、批量预测与节点修正做完潮位预测后很多工程场景还要求潮流流速和流向比如海底管线冲刷评估、通航条件论证。TMD2.05自带的tmd_ellipse.m可以基于调和常数构造潮流椭圆这是它最容易被忽略的功能。5.1 用tmd_ellipse计算潮流椭圆参数当模型包含分潮流的uv分量时tmd_ellipse.m可以根据M2分潮的U、V复振幅计算椭圆长轴、短轴、方位角和旋转方向。调用前需要先通过tmd_extract_HC获得东西向和南北向的复振幅再组合成椭圆参数[u, names] tmd_extract_HC(model, lat, lon, u); [v, ~] tmd_extract_HC(model, lat, lon, v); % 假设M2是第一个分潮 ell tmd_ellipse(u(1), v(1));ell各字段分别是长轴长度、短轴长度、椭圆倾角和旋转方向。注意浅水区M2椭圆接近直线这时方位角对噪声非常敏感需要结合K1等分潮做综合判断不要只靠单一分潮。5.2 批量多点预测与内存优化如果要画一张区域最大天文潮分布图需要遍历网格所有节点。逐个调用tmd_tide_pred会很慢我一般用tmd_tide_pred_mapts.m它一次传入多个经纬度点内部共享模型插值结果明显提速。也可以配合parfor并行但要注意每个worker都能访问TMD函数路径否则会报“Undefined function”错误。5.3 别忘了天文学修正最后说nodal.m和nodal_arg.m。潮汐预报不是简单的常数叠加月球的18.6年交点周期会造成振幅变化。TMD2.05内置了这两个函数在tmd_tide_pred中默认生效。如果你从TMD提取调和常数然后自己合成潮位一定要保留节点因子修正否则长序列预测在十几年后会漂移。TMD_timeOnOff.m可以控制这些修正是否开启对比开关前后的结果你能直观看到节点修正对M2分潮振幅的影响这也是检验你理解tmd_tide_pred内部逻辑最快的方式。本文还有配套的精品资源点击获取