網(wǎng)絡(luò)上似乎沒有公開的內(nèi)潮模態(tài)分解程序，筆者根據(jù)論文?Extraction of Internal Tidal Currents and Reconstruction of Full-Depth Tidal Currents from Mooring Observations （Cao et al., 2015)以及相關(guān)的視頻報告（https://www.bilibili.com/video/BV15L411K7iG?spm_id_from=333.880），編寫了s_modaldecomposition?函數(shù)，可以做內(nèi)潮的動力模態(tài)分解，該函數(shù)添加到了最新的s_tide v1.23 update 4版本里（下載鏈接：https://www.researchgate.net/project/A-non-stationary-tidal-analysis-toolbox-S-TIDE），

下面是這個函數(shù)的一個例子，無論是從時間上還是空間上，模態(tài)分解和重構(gòu)的效果都很不錯。

?load data.mat;?

?% u_k1 means u speed for K1 tide

[N,Fi,dFi,mc,um]=s_modaldecomposition(depth,density,u_k1,4); % 4 baroclinic modes

%um(5,:,:) is the barotropic mode

i=1900; plot(squeeze(um(4,:,i)+um(5,:,i)+um(3,:,i)+um(2,:,i)+um(1,:,i)),depth,'k')

?set(gca,'YDir','reverse');hold on; plot(u_k1(:,i),depth,'r*')

?ylabel('Depth(m)');xlabel('u(m/s)');

?legend('reconstructed','observed')

?

?d=18;

?plot(squeeze(um(4,d,:)+um(5,d,:)+um(3,d,:)+um(2,d,:)+um(1,d,:)),'k*')

?hold on;plot(u_k1(d,:),'r');xlim([1 100])

?ylabel('u(m/s)');xlabel('Time(hour)');

?title('Depth: 85m');legend('reconstructed','observed')