4. 风场扰动建模与仿真环境

风场扰动建模的目标，是把真实大气中连续、随机、带空间结构的风，转化为控制器设计、仿真验证和鲁棒性评估可用的数学输入。对无人飞行器而言，风不是单一的外力，而是同时作用在气动力、力矩、空速测量和导航解算上的多通道扰动。本节围绕 Dryden 与 von Kármán 湍流模型、离散阵风模型、风场空间相关性、仿真注入方法和典型风场数据集展开，重点说明“模型怎么选、参数怎么定、风怎么注入、结果怎么看”。

4.1 连续湍流：Dryden 与 von Kármán 模型

连续湍流描述的是大气边界层中随机脉动分量。工程上常用功率谱密度（PSD）刻画其频域能量分布。Dryden 谱为有理函数形式，便于生成白噪声驱动的成形滤波器；von Kármán 谱基于更严格的湍流统计理论，在高频段衰减更符合实测，但实现时需处理根号项，通常用近似有理式或频域采样。两者都通过尺度长度 L_u, L_v, L_w 和湍流强度 σ_u, σ_v, σ_w 参数化。

以 Dryden 纵向分量为例，其单边 PSD 可写为：

Φ_u(ω) = σ_u² · (2 L_u / (π V)) · 1 / (1 + (L_u ω / V)²)

其中：
  V   —— 飞行空速（m/s）
  L_u —— 纵向湍流尺度长度（m）
  σ_u —— 纵向湍流强度（m/s）
  ω   —— 角频率（rad/s）

生成时域序列的常用做法是：先产生单位白噪声，再通过传递函数 H_u(s) 成形。Dryden 纵向成形滤波器可写为：

H_u(s) = σ_u · sqrt(2 L_u / (π V)) · 1 / (1 + (L_u / V) s)

横向与垂向分量的滤波器阶次更高，含一对共轭极点。实际仿真中，若采样周期为 Ts，可用双线性变换或零阶保持离散化，得到差分方程后逐拍递推。von Kármán 模型则常写成：

Φ_u(ω) = σ_u² · (2 L_u / (π V)) · 1 / [1 + (1.339 L_u ω / V)²]^(5/6)

由于指数 5/6 不便直接有理化，工程上可用 Dryden 近似或频域采样后 IFFT 生成。选择建议：若控制器带宽较低、关注低频能量，Dryden 足够；若需评估高频载荷或与实测谱对比，优先 von Kármán。湍流强度可按 σ = 0.1 V 量级估算，尺度长度在低空常取 L_u ≈ 200 m、L_v ≈ 200 m、L_w ≈ 50 m，具体应结合高度和风场类别调整。

4.2 离散阵风模型

离散阵风用于描述确定性、短时、大幅值的风突变，常见形式有“1-cos”阵风和锐边阵风。1-cos 阵风在持续时间 T_g 内按余弦包络变化，峰值风速为 V_g，其表达式为：

V_w(t) = 0.5 V_g [1 - cos(2π t / T_g)]，0 ≤ t ≤ T_g
V_w(t) = 0，其他

锐边阵风则近似为阶跃，用于极限情况检查。离散阵风通常沿机体轴或地轴注入，需明确方向定义：本文约定 u_g 为顺风正、v_g 为右风正、w_g 为下洗正。注入时可直接叠加到空速分量或等效迎角、侧滑角上。评估指标包括最大迎角增量、俯仰角速率峰值、高度损失和恢复时间。若阵风峰值 V_g = 10 m/s、T_g = 2 s，则最大等效迎角扰动约为 atan(V_g / V)，在 V = 20 m/s 时约 26.6°，足以触发显著非线性响应。

4.3 风场空间相关性

真实湍流在空间上并非各点独立。空间相关性描述同一时刻不同位置风速的统计关联，常用相关函数或互谱密度表示。对多旋翼或大展弦比固定翼，机翼两端、前后旋翼之间会经历不同风矢量，若忽略相关性，可能高估或低估结构载荷与姿态扰动。工程简化中，可用指数型相关函数：

R(Δx) = exp(-Δx / L_c)

其中 Δx 为两点间距，L_c 为相关长度。仿真实现时，可先生成单点湍流序列，再通过一阶滤波或协方差分解生成多点相关序列。对于小尺寸无人机，若 L_c 远大于机体尺度，可近似认为全机经历同一风矢量；若 L_c 与旋翼间距相当，则必须分点注入。

4.4 仿真注入方法与典型数据集

风场注入通常有三种路径：一是直接叠加到气动模型的速度三角关系，即 V_air = V_ground - V_wind；二是作为外部力与力矩加到刚体动力学；三是注入到传感器模型，模拟空速管、迎角传感器和 GPS 的测量偏差。推荐优先采用第一种，因为它与气动系数计算自然衔接，且便于保持坐标系一致。

典型风场数据集包括：Dryden 标准谱生成的随机序列、MIL-F-8785C 规定的离散阵风、以及基于实测气象塔或激光雷达的湍流样本。仿真前应记录随机种子、采样率、湍流强度、尺度长度和注入方向，确保结果可复现。评估时，可对比无风、常值风、湍流风和阵风四种工况下的姿态误差、位置漂移和控制量饱和情况。

实际使用时，建议先确定飞行包线和评估目标，再选择模型：低频姿态控制验证用 Dryden 加 1-cos 阵风；高频载荷与传感器噪声评估用 von Kármán 加锐边阵风；多旋翼或大尺寸平台需加入空间相关性。注入后应检查空速、迎角和角速率是否出现非物理跳变，必要时限制变化率。最后，用同一随机种子重复运行，比较不同控制器的误差统计量，才能判断抗扰能力是否真正提升。