海洋要素计算实习报告
课程名称：海洋要素计算训练
实验名称：海洋要素计算综合训练
学院：海洋科学与工程学院
专业班级：海洋技术专业
学号：
姓名：
指导教师：
2014年10月
一、实验目的：
根据所给中国近海各个验潮站的水位资料，通过绘制横轴为时间，纵轴为水位的水位过程曲线，初步确定该验潮站的潮汐类型。并使用T-Tide工具箱，进行水位资料的调和分析。通过调和分析再次判断验潮站潮汐类型。
二、验潮站资料信息
验潮站：lianyungang
三、水位过程曲线
由图像可知该验潮站潮汐类型为：半日潮
三、一个月水位资料的调和分析
1、所有显著分潮
tidal amplitude and phase with 95% CI estimates
tide freq amp amp_err pha pha_err snr
MSF 0.0028219 0.0348 0.068 314.12 111.68 0.26
2Q1 0.0357064 0.0263 0.029 228.50 64.44 0.8
*Q1 0.0372185 0.0571 0.029 200.89 29.81 3.8
*O1 0.0387307 0.2035 0.029 214.44 8.44 48
NO1 0.0402686 0.0116 0.029 108.69 107.25 0.16
*K1 0.0417807 0.3078 0.029 270.16 5.51 1.1e+02
J1 0.0432929 0.0086 0.029 185.77 197.89 0.086
OO1 0.0448308 0.0174 0.029 25.32 88.44 0.35
UPS1 0.0463430 0.0092 0.029 81.46 180.54 0.098
*N2 0.0789992 0.3442 0.032 294.05 5.29 1.1e+02
*M2 0.0805114 1.6306 0.032 306.33 1.12 2.6e+03
*S2 0.0833333 0.5894 0.032 4.70 3.13 3.4e+02
ETA2 0.0850736 0.0126 0.032 269.43 147.84 0.15
*MO3 0.1192421 0.0579 0.017 144.39 16.66 12
*M3 0.1207671 0.0294 0.017 73.23 31.98 3.1
*MK3 0.1222921 0.0290 0.017 198.67 32.88 3
*SK3 0.1251141 0.0200 0.017 263.96 47.92 1.5
*MN4 0.1595106 0.0582 0.009 153.86 8.81 41
*M4 0.1610228 0.1333 0.009 168.04 3.86 2.1e+02
*MS4 0.1638447 0.1102 0.009 230.22 4.71 1.5e+02
*S4 0.1666667 0.0209 0.009 305.53 24.97 5.3
2MK5 0.2028035 0.0021 0.004 126.22 97.31 0.34
2SK5 0.2084474 0.0034 0.004 123.35 61.62 0.88
*2MN6 0.2400221 0.0102 0.004 7.42 22.76 6
*M6 0.2415342 0.0126 0.004 31.78 18.57 9.1
*2MS6 0.2443561 0.0186 0.004 97.52 12.67 20
*2SM6 0.2471781 0.0056 0.004 159.12 42.73 1.8
*3MK7 0.2833149 0.0028 0.003 88.52 52.67 1.1
*M8 0.3220456 0.0030 0.002 228.07 32.90 2.8
M10 0.4025570 0.0016 0.002 107.05 63.87 0.74
2、列表写明该验潮站四个主要分潮M2，S2,K1与O1的调和常数tide Amp(振幅) Pha（迟角）
M2 1.6306 306.33
S2 0.5894 4.70
K1 0.3078270.16
O1 0.2035214.44
四、在该月资料中截取25小时（一个太阴日）的资料进行调和分析。
tidal amplitude and phase with 95% CI estimates
tide freq amp amp_err pha pha_err snr
*K1 0.0417807 0.0697 0.069 354.38 56.83 1
*M2 0.0805114 1.9803 0.069 290.75 1.97 8.3e+02
M3 0.1207671 0.0590 0.069 253.18 65.97 0.74
*M4 0.1610228 0.2429 0.069 143.45 15.97 12
2MK5 0.2028035 0.0356 0.069 195.97 109.58 0.27
*3MK7 0.2833149 0.0415 0.019 38.35 26.43 4.5
*M10 0.4025570 0.0286 0.019 50.51 37.59 2.2
1、与题2进行比较，说明资料时间长度对调和分析结果的影响
资料时间长的得到的潮汐数据较为详细准确，能够通过数据判断潮汐类型，而短时间的数据测量得到的数据不够判断验潮站所在区域的潮汐类型。
2、列表写明该验潮站四个主要分潮M2，S2,K1与O1的调和常数（包括振幅和迟角）Tide Amp Pha
M2 1.9083 290.75
K1 0.0697 354.38
五、主要程序
clear all
close all
clc
lusi=load('D:\haiyangyaosu\作业2\sl_lianyungang_19930114_19930214.mat');
lusi= struct2cell(lusi);
lusi= cell2mat(lusi);
lusi_1=double(lusi);
lusi_1=lusi_1/1000.0; %原单位为毫米，先将单位转换成米，lusi1即可作为t_tide 的输入值
%画图，横轴为时间（days），纵轴为水位(m)的水位过程曲线。
A=linspace(1,31,length(lusi_1));
figure(1);
plot(A,lusi_1);
title('lianyungang','FontSize',14);
xlabel('时间/days','FontSize',12);
ylabel('水位/m','FontSize',12);
[tidestruc,pout]=t_tide(lusi_1,'interval',1,'start',[1993,01,14],...
'latitude',31.7,'shallow','M10','error','linear','synthesis',1);
% fieldnames(tidestruc);
% data_2= cell2mat(data);
% madata=double(data_2);
% save result.txt madata -ASCII
lusi_2=lusi_1(1:25);
[tidestruc,pout]=t_tide(lusi_2,'interval',1,'start',[1993,01,14],...
'latitude',31.7,'shallow','M10','error','linear','synthesis',1);
六、根据题2和题3的结果，计算潮型系数A。并根据潮型系数A判断该验潮站的潮汐类型。
由公式：
21
1
M O
K H H
H A +
=
A=0.32<0.5，属于半日潮
A=0.032<0.5，属于半日潮
七、实验感想
通过对潮汐实验数据的计算，得到潮汐水位过程曲线，在这个过程中，对潮汐水位的分布，潮汐类型的特点有了更加深刻的认识，同时，在判断潮汐类型的方法上有了更多的了解，在以后的处理过程中会有更加熟练的操作。