潮流准调和分析



Dd.m [DO1 dO1 DK1 dK1 DM2 dM2 DS2 dS2 DM4 dM4 DMS4 dMS4]=Dd(y,D,Y,tm,A1) 计算6个准调和分潮的振幅系数和初相角公式



aM.m [N n aa liu_u liu_v]=aM() 基础系数矩阵，目前D和d值是取用的各段数据的中间时间值，t是以中间为0左右对称取；

视差潮龄A1需要预先得知

同时本程序负责读取原始测流数据，格式为

站位名称

观测数据段数N（即大 中 小潮）

第一段起始时间YYYY MM DD HH mm ss

第一段个数n1

speed direction

speed direction

...

第二段起始时间 YYYY MM DD HH mm ss

第二段个数n2

speed direction

speed direction

...

第三段起始时间 YYYY MM DD HH mm ss

第三段个数n3

speed direction

speed direction

...



UVg.m [U V Ug Vg]=UVg() 求调和常数U V和迟角Ug Vg



Tuoyuan.m 计算各准调和分潮的椭圆要素W theta K k

-The quasi harmonic analysis of tidal current

Dd.m [DO1 dO1 DK1 dK1 DM2 dM2 DS2 dS2 DM4 dM4 DMS4 dMS4]=Dd (y, D, Y, TM, A1)　calculation of 6 quasi harmonic division of amplitude and initial phase angle formula of the coefficient of tide

AM.m [N n AA liu_u liu_v]=aM ()　foundation coefficient matrix, the current D and D values are the middle time of each section of data taken　the value of T is in the middle for symmetric about 0

Age of parallax A1 need to know in advance

At the same time, this program is responsible for reading the original flow measurement data, format for

Site name

The observation data segment number N (that is, in the low tide)

The first section of the YYYY MM DD HH mm the starting time of SS

The first section number N1

Speed direction

Speed direction

...

The second section starting time YYYY MM DD HH mm SS

The second section number N2

Speed direction

Speed direction

...

The third section starting time YYYY MM DD HH mm SS

The third section number