3. 参数空间与扫描设计

参数扫描的本质,是在一个由若干参数张成的空间中选取有限个采样点,逐点求解动力系统并记录其行为。扫描设计要回答四个问题:扫哪些参数(维度)、扫到哪里(范围)、隔多远取一点(分辨率)、点怎么排布(采样策略)。这四个问题相互制约:维度增加会使空间体积指数膨胀,分辨率提高会使计算量线性甚至超线性增长,而采样策略决定了在有限预算下能否抓住关键结构。

参数维度:先做敏感性筛选,再决定扫几个

并非所有参数都值得进入扫描网格。一个典型系统可能同时含有激励幅值、频率、阻尼系数、非线性刚度、时滞量等参数,但真正改变定性行为的往往只有两三个。工程上的做法是先做单参数粗扫或局部敏感性分析:固定其余参数在名义值,对每个候选参数单独扫一遍,观察响应指标(如最大李雅普诺夫指数、周期数、振幅、分岔点位置)是否出现显著变化。变化剧烈的参数保留,变化平缓的参数固定。

维度选择还要考虑可视化与解释成本。二维参数平面可以直接用颜色图呈现行为分区,三维以上只能切片或投影。因此常见策略是:主参数二维全扫,次参数取若干离散水平做分层扫描。例如以激励幅值 A 和频率比 Ω 为二维主平面,阻尼比 ζ 取 0.02、0.05、0.10 三个水平,形成三层平面,每层独立成图再对比。

范围:由物理约束与预期分岔位置共同确定

范围定得太窄会漏掉分岔,定得太宽则大量采样点落在无信息区域。确定范围有三条依据:一是物理可行域,如阻尼比必须为正、频率比通常取 0.5~2.0;二是理论或经验预估的分岔位置,如 Duffing 系统的主共振在 Ω≈1 附近,鞍结分岔常在幅值超过某阈值后出现,范围应覆盖阈值两侧各一段;三是先做一维粗扫定位边界,再据此收紧。

范围端点要留出“安全边距”。若预估分岔在 A=1.2 附近,范围取 0.5~2.0 比取 1.0~1.4 更稳妥,因为分岔位置本身可能随其他参数漂移。同时应记录范围选取理由,便于后续复现和调整。

分辨率:由最窄结构尺度决定,而非凭感觉取点数

分辨率的核心判据是:网格步长必须小于待识别结构的最小特征尺度。若周期窗口宽度约为 0.01(以参数单位计),步长取 0.05 就会整体跳过该窗口,扫描图上表现为“干净”的周期区,结论完全错误。反之,若只关心宏观分区,步长 0.02 已足够,盲目加密到 0.001 只会浪费算力。

实际操作中可采用两级分辨率:先用较粗步长(如范围的 1/100)做全局扫描,识别出行为突变带的大致位置;再在突变带附近用细步长(如 1/1000)局部加密。这种“粗扫定位、细扫刻画”的方式,比全程高分辨率更经济,也更符合分岔集通常低维嵌入高维参数空间的几何特征。

采样策略:规则网格、分层与自适应

规则网格是最直接的方式:在二维平面上按等间距取点,形成矩形阵列。优点是实现简单、结果可直接映射为像素图;缺点是对角线方向的结构(如斜向分岔线)会出现阶梯状锯齿,且在高维时点数爆炸。

分层扫描用于处理多参数:对主参数做全网格,对次参数取离散水平,逐层计算并对比层间差异。若层间行为分区发生明显迁移,说明次参数影响显著,应提升其采样密度或转为联合扫描。

自适应采样则根据已算点的局部变化率决定下一步加密位置:在行为指标梯度大的区域插入新点,在平坦区域停止加密。它适合分岔集占比小、预算有限的场景,但实现复杂度较高,且结果依赖初始粗网格的质量。

参数平面扫描设计:粗网格定位 + 局部加密 + 分层对比 ① 粗扫(步长 0.05) A ∈ [0.5, 2.0],Ω ∈ [0.5, 2.0] 红色虚线:预估分岔带 加密 ② 分岔带局部加密(步长 0.005) 仅对虚线框内加密,点数可控 可分辨宽度约 0.01 的周期窗口 分层 ③ 次参数分层(ζ = 0.02 / 0.05 / 0.10) ζ = 0.02:分岔带偏左 ζ = 0.05:分岔带居中 ζ = 0.10:分岔带偏右 层间分区迁移 → 次参数影响显著 应提升 ζ 采样密度或转联合扫描 设计原则:步长 < 最窄结构尺度;范围覆盖分岔两侧并留边距;维度先经敏感性筛选; 粗扫定位 → 局部加密 → 分层对比,在有限算力下兼顾全局覆盖与局部刻画。

一个可执行的扫描配置示例

以受迫 Duffing 振子 ẍ + 2ζω₀ẋ + ω₀²x + βx³ = A·cos(Ωt) 为例,取 ω₀=1、β=0.5,考察 (A, Ω) 平面上的行为分区,ζ 作为分层参数。下面是一份可直接参照的扫描配置(参数为教学示例设定,非实测结果):

主参数:A ∈ [0.5, 2.0],Ω ∈ [0.5, 2.0]
分层参数:ζ ∈ {0.02, 0.05, 0.10}
粗扫步长:ΔA = ΔΩ = 0.05   → 每层 31×31 = 961 点
加密判据:相邻点最大李雅普诺夫指数符号变化,或周期数跳变
加密步长:ΔA = ΔΩ = 0.005  → 仅在突变带 ±0.05 范围内加密
单点积分:四阶龙格-库塔,步长 0.01,舍弃前 200 个周期,统计后 100 个周期
行为标签:周期1 / 周期n / 准周期 / 混沌(由李雅普诺夫指数与庞加莱截面判定)

按此配置,粗扫每层约 961 次积分,三层合计约 2883 次;若突变带面积约占平面的 8%,加密区每层约 400 点,三层再增约 1200 次。总计算量约为全程 0.005 步长(每层 301×301≈9 万点)的 1.5%,却能获得接近细网格的分岔带刻画精度。这正是“粗扫定位、细扫刻画”策略的价值所在。

判断分辨率是否足够的实用检验:将步长减半重扫同一区域,若行为分区图的结构(分区数量、边界走向)基本不变,说明当前分辨率已收敛;若出现新的窄窗口或边界明显移动,则需继续加密。这一“步长收敛检验”比凭经验取点数更可靠。

常见设计失误与规避

  • 范围贴边:分岔恰好落在范围端点,无法判断其两侧行为。规避方法是先做一维粗扫,确认端点外侧仍为同一行为区。
  • 分辨率与结构尺度失配:周期窗口被整体跳过,扫描图呈现虚假的“干净”分区。规避方法是先用小范围高分辨率试扫,估计最窄结构宽度。
  • 维度贪多:四五个参数同时全网格扫描,点数爆炸且结果难以解释。规避方法是先做敏感性筛选,主参数不超过两个,其余分层或固定。
  • 忽略瞬态:积分时间不足,把瞬态误判为稳态行为。规避方法是固定舍弃周期数,并用加倍舍弃周期的方式检验统计量是否稳定。

扫描设计没有唯一正确答案,但有明确的判断依据:维度由敏感性决定,范围由物理约束与预估分岔位置决定,分辨率由最窄结构尺度决定,采样策略由算力预算与结构分布决定。把这四条依据落实到具体数值,并在扫描后做步长收敛检验,才能保证参数扫描的结果既经济又可信。