研究潮汐季节循环需要注意的事项_blue_ocean

http://blog.sina.com.cn/u/2933243165

首页博文目录关于我

个人资料

微博

加好友发纸条

写留言加关注

博客等级：
博客积分：

博客访问：
关注人气：
获赠金笔：0支
赠出金笔：0支
荣誉徽章：

正文字体大小：大中小

研究潮汐季节循环需要注意的事项

(2020-09-14 15:25:47)

标签：

潮汐调和分析

t_tide

s_tide

分类：数学算法

代码所需要的数据kushiro.mat在s_tide工具包里https://www.researchgate.net/project/A-non-stationary-tidal-analysis-toolbox-S-TIDE

clear
load kushiro.mat
for i=1:36

aa=1+720*(i-1);bb=720*i;
[NAME,FREQ,TIDECON,XOUT]=t_tide(kushiro(aa:bb),'interval',1);
K1(i)=TIDECON(6,1);
S2(i)=TIDECON(12,1);
end

plot(K1);hold on;plot(S2,'r')

结果如下图所示，可以看到K1和S2都存在显著的半年变化。需要注意的是逐月的潮汐调和分析是没法把P1和K1区分，以及K2和S2区分。而P1的频率恰好处在K1的半年变化上，而K2的频率恰好处在S2的半年变化上。所以下图实线所示的K1和S2振幅的半年变化是虚假的。

下面的程序使用了差比关系，根据K1和P1以及S2和K2的理论关系将P1和K2分潮强行分辨出来，结果如虚线所示。半年变化强度有所减弱。

infername=['P1';'K2']; inferfrom=['K1';'S2'];
infamp=[.33093;.27215];infphase=[-7.07;-22.40];
for i=1:36

aa=1+720*(i-1);bb=720*i;
[NAME,FREQ,TIDECON,XOUT]=t_tide(kushiro(aa:bb),'interval',1, 'inference',infername,inferfrom,infamp,infphase);
K1(i)=TIDECON(7,1);
S2(i)=TIDECON(13,1);
end
hold on
plot(K1,'--');hold on;plot(S2,'r--')
legend('K1','S2','K1r','S2r');ylabel('Amplitude(mm)');xlabel('Month')

阅读┊ 收藏 ┊ 喜欢 ▼ ┊打印┊举报/Report

前一篇：潮汐分潮名称产生的混乱

后一篇：新型潮汐调和分析工具包S_TIDE教程

新浪BLOG意见反馈留言板　欢迎批评指正