DOC: set more suitable parameters for f-v spectra visualization (#160)

Author: Dengda98

docs/source/Advanced/kernel/kernel.rst

@@ -1,6 +1,6 @@
 :author: 朱邓达
 :date: 2025-04-20
-:updated_date: 2025-08-06
+:updated_date: 2025-12-05
 
 核函数 :math:`f-v` 谱图 
 =========================
@@ -21,6 +21,7 @@
 .. note:: 
 
     与之前的示例相比，这里添加了控制波数积分的参数，以显式控制核函数的采样间隔和范围。但这并非是计算格林函数的最佳参数，故此时格林函数计算结果很可能是不收敛的。
+    这里的目的仅是提取特定范围内不同波数下的核函数。
 
 
 .. tabs:: 
@@ -79,4 +80,8 @@ Python端提供了 :py:func:`pygrt.utils.read_kernels_freqs` 函数来完成所
 
 从图上可看到， **核函数** :math:`f-v` **谱图的虚部峰值正是频散曲线的位置** ，尽管其中不同震源的核函数幅值不同，但 :math:`q_m，w_m` 呈现的频散特征一致，这对应的是 **Rayleigh波频散** ，而 :math:`v_m` 对应的是 **Love波频散** 。
 
-不同的震源场点深度、不同的虚频率都会导致幅值变化，积分间隔则会影响插值效果。
+不同的震源场点深度、不同的虚频率都会导致幅值变化，积分间隔则会影响插值效果。
+
+.. note::
+
+    注意其中仍然使用了虚频率，尽管这对于面波频散的分析引入了少许误差，但这使得 :math:`f-v` 谱图中的高亮区域变宽，方便观察。

docs/source/Advanced/kernel/run/run.py

@@ -7,7 +7,7 @@
 
 modarr = np.loadtxt("mod1")
 
-pymod = pygrt.PyModel1D(modarr, depsrc=0.01, deprcv=0.0)
+pymod = pygrt.PyModel1D(modarr, depsrc=0.03, deprcv=0.0)
 
 # 不指定statsidx，默认输出全部频率点的积分过程文件
 # vmin_ref 显式给定参考速度（用于定义波数积分上限），避免使用PTAM

docs/source/Advanced/kernel/run/run.sh

@@ -7,6 +7,6 @@ rm -rf GRN* syn* pygrtstats* *.png
 # -S 后不指定索引表示输出所有频率点的核函数
 # -K+v0.1 显式给定参考速度（用于定义波数积分上限），避免使用PTAM
 # -L20 定义波数积分间隔dk
-grt greenfn -Mmod1 -D0.01/0 -N500/0.02/0.8 -OGRN -R1 -K+v0.1 -S -L20
+grt greenfn -Mmod1 -D0.03/0 -N500/0.02 -OGRN -R1 -K+v0.1 -S -L20
 # END GRN
 # -----------------------------------------------------------------