DOC: add integral and differential in Tutorial

Author: Dengda98

docs/source/Tutorial/dynamic/run/run.py

@@ -165,4 +165,49 @@ def plot_syn(stsyn:Stream, out:str|None=None, sigs:np.ndarray|None=None):
 # 卷积，原地修改
 pygrt.utils.stream_convolve(stsyn, trig)
 # END TIME FUNC
-plot_syn(stsyn, "syn_sf_trig.png", trig)
+plot_syn(stsyn, "syn_sf_trig.png", trig)
+
+
+
+# BEGIN INT DIF
+idx = 2
+stgrn = stgrnLst[idx]   # 选择格林函数
+
+stsyn = pygrt.utils.gen_syn_from_gf_MT(stgrn, M0=1e24, MT=[0.1,-0.2,1.0,0.3,-0.5,-2.0], az=30)
+
+# 使用inplace=False，防止原地修改
+stsyn_int = pygrt.utils.stream_integral(stsyn, inplace=False)
+stsyn_dif = pygrt.utils.stream_diff(stsyn, inplace=False)
+# END INT DIF
+
+def plot_int_dif(stsyn:Stream, stsyn_int:Stream, stsyn_dif:Stream, comp:str, out:str|None=None):
+    nt = stsyn[0].stats.npts
+    dt = stsyn[0].stats.delta
+    t = np.arange(nt)*dt
+
+    travtP = stsyn[0].stats.sac['t0']
+    travtS = stsyn[0].stats.sac['t1']
+
+    fig, axs = plt.subplots(3, 1, figsize=(10, 4), gridspec_kw=dict(hspace=0.0), sharex=True)
+    for i, (st, suffix) in enumerate(zip([stsyn, stsyn_int, stsyn_dif], ["", "_int", "_dif"])):
+        tr = st.select(component=comp)[0]
+
+        ax = axs[i]
+        ax.plot(t, tr.data, c='k', lw=0.5, label=f"{tr.stats.channel}{suffix}")
+        ax.legend(loc='upper left')
+
+        ylims = ax.get_ylim()
+        # 绘制到时
+        ax.vlines(travtP, *ylims, colors='b')
+        ax.text(travtP, ylims[1], "P", ha='left', va='top', color='b')
+        ax.vlines(travtS, *ylims, colors='r')
+        ax.text(travtS, ylims[1], "S", ha='left', va='top', color='r')
+
+        ax.set_xlim([t[0], t[-1]])
+        ax.set_ylim(np.array(ylims)*1.2)
+
+        if out is not None:
+            fig.savefig(out, dpi=100)
+
+for ch in ['Z', 'R', 'T']:
+    plot_int_dif(stsyn, stsyn_int, stsyn_dif, ch, f"syn_mt_intdif_{ch}.png")

docs/source/Tutorial/dynamic/run/run.sh

@@ -39,3 +39,10 @@ grt.syn -GGRN/milrow_2_0_10 -S1e24 -A30 -M33/50/120 -Osyn_dc_zne -N
 grt.syn -GGRN/milrow_2_0_10 -S1e16 -A30 -F1/-0.5/2 -Osyn_sf_trig -Dt/0.3/0.3/0.6
 # END TIME FUNC
 
+
+# BEGIN INT DIF
+# 积分，1表示积分1次
+grt.syn -GGRN/milrow_2_0_10 -S1e24 -A30 -T0.1/-0.2/1.0/0.3/-0.5/-2.0 -I1 -Osyn_mt_int1
+# 微分，1表示微分1次
+grt.syn -GGRN/milrow_2_0_10 -S1e24 -A30 -T0.1/-0.2/1.0/0.3/-0.5/-2.0 -J1 -Osyn_mt_dif1
+# END INT DIF

docs/source/Tutorial/dynamic/run/syn_mt_intdif_R.png

@@ -140,7 +140,7 @@
 
 
 
-时间函数
+卷积时间函数
 ---------------------
 **PyGRT** 内置了一些震源时间函数，例如抛物波、梯形波、雷克子波或自定义，这里以单力源为例。
 
@@ -165,4 +165,37 @@
         其它时间函数以及具体参数用法可在 :py:mod:`pygrt.signals` 模块中查看函数参数。
 
 .. image:: run/syn_sf_trig.png
+   :align: center
+
+
+
+位移对时间积分、微分
+--------------------------------
+这里以矩张量源为例。
+
+.. tabs::  
+
+    .. tab:: C 
+
+        .. literalinclude:: run/run.sh
+            :language: bash
+            :start-after: BEGIN INT DIF
+            :end-before: END INT DIF
+
+
+    .. tab:: Python 
+
+        .. literalinclude:: run/run.py
+            :language: python
+            :start-after: BEGIN INT DIF
+            :end-before: END INT DIF
+
+
+.. image:: run/syn_mt_intdif_Z.png
+   :align: center
+
+.. image:: run/syn_mt_intdif_R.png
+   :align: center
+
+.. image:: run/syn_mt_intdif_T.png
    :align: center