基于时-频域联合分析的非平稳风速反演方法.docxVIP

基于时-频域联合分析的非平稳风速反演方法 1 相干测风激光雷达 正确测量对研究陆军系统的动态模型、数值监测和军事环境的精确报告具有重要意义。它可以用来监测机场的风场变化，监测机场路线前面的风场，为航空航天提供重要的安全保障。 相干测风激光雷达系统采用相干接收体制 利用数据处理算法对仿真的回波信号进行验证、分析和风速反演。目前应用于相干测风激光雷达的数据处理算法主要有：快速傅里叶变换(FFT) 2 目标回波时间模型的特征分析 2.1 胶散射体的回波信号 相干多普勒测风激光雷达基于光外差探测原理 基于端到端的大气分层模型对非稳定大气条件下的相干多普勒测风激光雷达的回波信号进行时域仿真。大气分层模型将探测范围内的气溶胶散射体沿光束传播轴分为M层，每一层的厚度为ΔR，并假设每一层内的气溶胶粒子运动速度大小和方向相同、散射系数恒定(如图1所示)。通过全部分层回波信号的叠加得到总的回波信号。层数M的取值一般为几百到上千，与10 如图1所示，分层厚度ΔR=1/2×c×ΔT，其中c为光速，ΔT为采样时间。层数M=c×P 假设层与层之间的回波信号相互独立，将每一层的回波信号进行非相干累加，可得到总的回波信号为S(t)，其可表示为 式中P K? 式中A是接收望远镜的面积；R 在探测范围内，按仿真所需的风速细节设定v 2.2 时-频域分析 随机信号在理论上可以分为平稳和非平稳两大类。非平稳信号是指统计特性随时间变化的随机信号，其概率密度P(t,x)是时间的函数 式中a 式中f 假设目标运动的径向速度为线性，即v 即在探测时间内(6)式表示的回波信号可看作是啁啾信号，其具有明显的时变特征，是非平稳信号。 傅里叶变换建立了信号从时域到频域的一一对应的映射关系，即 傅里叶变换将一个函数表示成无穷多个最和谐的函数—正弦函数的加权和,对信号的表征要么完全在时域，要么完全在频域，不能揭示某种频率分量出现在什么时候以及随时间的变化情况。因此，常用的快速傅里叶变换方法适于分析信号组成分量的频率不随时间变化的平稳信号。 由以上分析，多普勒回波信号属于非平稳信号的分析领域。为了克服傅里叶变换的局限性，引入时-频域联合分析方法。时频分析方法通过分析时间、频率及能量之间的密切关系，提取信号中所包含的特征信息，以揭示信号中所包含的频率分量及其演化特性 Wigner-Ville分布是常用的时频分析方法，在已提出的各种时频分布中，形式简洁并且具有很好的性质。Wigner-Ville分布具有明确的物理意义，可被看作信号能量在时域和频域中的分布。Wigner在研究量子力学时提出的这种分布可表示为 由于Wigner-Ville分布不是线性的，应用于大气分层模型时域仿真回波信号时会有交叉项产生，一般可通过加入核函数来抑制 3 回波信号仿真 首先对最简单的风速分布情况(常风速)进行分析，常风速是一种特殊的线性风速情况，在探测范围内令各层的风速v 对相干测风激光雷达时域仿真回波信号进行功率谱评估，得到以距离门为空间分辨率的一维风速信息。选取第一个距离门60～120 m为研究对象,由快速傅里叶变换原理的风速反演算法 功率谱峰值对应的频率(52.73 MHz)与声光频移量υ 4 应用期间频繁转换的反演速度 4.1 信号能量的时域分布及优化 通过描述风速细节，可将大气分层模型应用于非平稳的大气环境。下面对v 对比图2与图5可以看到，两图中回波信号都由于散斑的低频调制作用出现了随机起伏，变化趋势基本相同。但v 基于Wigner-Ville时频变换对仿真得到的回波信号进行定量分析，得到频率随时间的变化关系，进而反演出风速值，并与输入的风速值进行比较。系统预置的声光调制频移υ Wigner-Ville时频分布可被看作信号能量在时域和频域中的分布。由每个时间点对应的信号能量的峰值及该能量峰值对应的频率坐标，即信号能量的峰值在时频面上的分布，就可得到一个频率随时间变化的关系，由多普勒测风原理可以反演得到每个时间点对应的风速值。反演得到的风速与线性剪切风速的对比如图7所示。 从图7中可以看到，反演风速在整体上与输入的线性剪切风速有较好的一致性。但是由于系统接收到的大气回波信号非常微弱，加上散斑噪声等因素的影响，反演风速出现了不同程度的波动。因此，采用多个脉冲累加的方式来减少波动，对60组脉冲信号分别作时频变换后累加平均进而反演风速值。图8为60组发射脉冲仿真得到的回波信号作Wigner-Ville变换后累加平均的时频分布。从图8中可以看到，信号能量在时频面上的分布的聚集性得到了明显改善，线性趋势更明显，频率分布更平滑，趋近于线性。噪声对信号及反演风速的影响减小，反演风速的波动得到了明显的改善，如图9所示。分别用单个脉冲和60组脉冲回波信号来反演风速，两次反演风速间的误差随探测距离的变化曲线如图10所示。