频率的单位转换:
1 hz = 1/s; 3hz = 3/s
1cph = 1/(1h) =1/3600s, 因此 3600cph = 1/s = 1hz ; 1cph = 1/3600 hz; 1hz = 3600cph
1cpd = 1/(24h) = 1/24cph = 1/(3600×24)s ; 1cpd = 1/24 cph = 1/(3600×24) hz; 1hz = 3600*24cpd
1cph = 24cpd
s2分潮的周期为12小时,该周期对应的频率为 24/12 = 2cpd,
f=1/(12×3600) hz = 3600*24/(12×3600) = 2cpd= 3600×1/(12×3600)=1/12cph
在功率谱的计算结果中:功率谱密度的单位 为 [power/freq],如水位序列(m)单位为 [m2/hz],
此时需要将功率谱密度的转换为[m2/cph]或cpd
1[m2/hz] = 1/3600 [m2/cph]
1[m2/cph] = 1/24[m2/cpd]
1[m2/hz] = 1/3600 [m2/cph] = 1/(3600 ×24)[m2/cpd]
因此功率谱的输出结果 单位的转换
p = p/(3600*24) ; 功率谱密度 由 [m2/hz]→[ m2/cpd]
freq = freq×(3600×24); 频率 由 [hz] → [cpd]
rdn = rdn/(3600*24) ; 红噪声检验 由 [m2/hz]→[ m2/cpd]
==================================================================================
若对时长为8766hour的水位数据作调和分析,可得到68个分潮,长周期6个,全日分潮21个,半日分潮19个,浅水分潮22个,他们在显示的时候标注原有分潮的名称,并用 以下语句标识:
text(tx,ty,‘
Q
1
Q_{1}
Q1’,‘interpreter’,‘latex’)
text(tx,ty,naD1{4},‘interpreter’,‘latex’)
长周期分潮 (6个)
分潮名称简写
naDL ={'$S_{a}$';...
'$S_{sa}$';...
'$M_{sm}$';...
'$M_{m}$';...
'$M_{sf}$';...
'$M_{f}$'};
全日分潮 (21个)
naD1 ={'$\alpha_{1}$';...
'$2Q_{1}$';...
'$\sigma_{1}$';...
'$Q_{1}$';...
'$\rho_{1}$';...
'$O_{1}$';...
'$\tau_{1}$';...
'$\beta_{1}$';...
'$NO_{1}$';...
'$\chi_{1}$';...
'$\pi_{1}$';...
'$P_{1}$';...
'$S_{1}$';...
'$K_{1}$';...
'$\psi_{1}$';...
'$\varphi_{1}$';...
'$\theta_{1}$';...
'$J_{1}$';...
'$SO_{1}$';...
'$OO_{1}$';...
'$\upsilon_{1}$'};
半日分潮D2(19个)
naD2 = {'$OQ_{2}$';...
'$\varepsilon_{2}$';...
'$2N_{2}$';...
'$\mu_{2}$';...
'$N_{2}$';...
'$\nu_{2}$';...
'$\Gamma_{2}$';...
'$H_{1}$';...
'$M_{2}$';...
'$H_{2}$';...
'$MKS_{2}$';...
'$\lambda_{2}$';...
'$L_{2}$';...
'$T_{2}$';...
'$S_{2}$';...
'$R_{2}$';...
'$K_{2}$';...
'$MSN_{2}$';...
'$\eta_{2}$'};
浅水分潮(22个)
naDH = {'$MO_{3}$';...
'$M_{3}$';...
'$SO_{3}$';...
'$MK_{3}$';...
'$SK_{3}$';...
'$MN_{4}$';...
'$M_{4}$';...
'$SN_{4}$';...
'$MS_{4}$';...
'$MK_{4}$';...
'$S_{4}$';...
'$SK_{4}$';...
'$2MK_{5}$';...
'$2SK_{5}$';...
'$2MN_{6}$';...
'$M_{6}$';...
'$2MS_{6}$';...
'$2MK_{6}$';...
'$2SM_{6}$';...
'$MSK_{6}$';...
'$3MK_{7}$';...
'$M_{8}$'};