如何构建潮汐信号谐波成分?

Question

我试图重建天文潮时间序列的谐波成分从T_Tide(计算)。我用这个方程,所有知道H (t) =振幅* cos (t * harmonicSpeed +相位滞后),我写在python脚本。输出不匹配与潮流predicction(从T_Tide文件)。正如你所知道的建筑很简单但是我认为我的问题是关于采样间隔,我真的不知道如何解释它´对于这种情况,我写了1/24由于输入时间(在T_TIDE计算,每小时数据),我不确定。

我分享的调和常数文件从T_Tide和潮水predicction值文件(每小时,我想复制)。他们是https://github.com/feraxel/Harmonic-constans/commit/fb11b9ef9909b03e722978ccfa102fc7b9fd1f01 diff - 09639 - ac17104196efd895598afe1af5234cbb5bcac39070e4a3d9133c611a4ea和https://github.com/feraxel/Harmonic-constans/commit/1e46d1d416c49ef131005d8fa68775f9169c0270 diff - 545 b5ec39a1b72b56fa5e37c33afb30a3c3c39ba0e6fe1d583b4a1d77e9fc74a

os.chdir (r 'directory”) #改变文件。con = pd.read_csv (“harmonics_aca txt目录。txt ', 9 = ' ') #我们负载谐波常数

我们创建一个datetime向量连接构造的潮流

日期= pd。date_range (start = ' 01/02/2018 22:00:00”,结束=“05/30/2022 15:00:00”,频率= ' H ') datetime = pd.DataFrame ({Fecha:日期})潮= 0 H =[] #这个向量保存潮t次dt = 1/24 # # # # #采样间隔为np.arange t (0, len (Fecha) * dt, dt):我= con.iloc[2]我#谐波成分振幅Ph = np.deg2rad (con.iloc[4]我)#谐波成分阶段Fr = (con.iloc[1]我* 2 * np.pi) # Frecuency(1 /小时)* 2 *π=我* np。cos (Fr * t - Ph值)潮=潮汐+水平H.append(潮流)潮= 0

Answer 1

我的建议是,你检查提供的潮汐预测代码t_tide (t_predic)。在那里,你会发现这样一段代码:

美联社= tidecon (: 1) / 2。* exp(我* tidecon(:, 3) *π/ 180);我=连词(美联社);

你可以检查你做的那部分正确。

另一个重要部分是天文节点调制修正。这是为您的特定位置的不同阶段。t_vuf中给定的代码。