计算潮流椭圆的源程序，可以用调和常数来计算计算椭圆要素，并画图。-tide ellipse

可以根据给出的不同潮流分潮，进行潮流椭圆参数的计算，输出结果可以直接绘制潮流椭圆图。-it s used to compute tidal ellipse parameters

潮流准调和分析 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（即大 中 小潮）

用matlab自带函数psd对数据进行功率谱分析；画潮流椭圆图(Power spectrum analysis of data using MATLAB self band function PSD and a drawing of a tidal current ellipse.)