气象杂波.docVIP

气象杂波 基于FPGA的气象杂波图设计与实现,气象杂波,立体杂波图,检测概率,虚警概率,FPGA 1 引 言 无论什么体制的雷达都会受到其工作环境中的噪声和杂波的干扰，从噪声和杂波中发现目标是雷达信号处理的基本任务。在碧空如洗的天空，空中目标检测是最容易的，随着气象变化，会遇到云、雨、雪、冰雹等不同天气，那时目标检测变得相对困难了。各国都有专门的气象雷达用来分析上述气象目标，但对于检测飞机、导弹等飞行目标而言，这些气象目标称之为气象杂波。本文将立体杂波图应用于气象杂波的CFAR处理，根据当前杂波环境的变化实时地产生杂内杂外标志，从而选择不同的信号处理支路处理当前气象杂波，提高了雷达的检测性能，降低了虚警概率。 2 气象杂波 立体杂波包括：气象杂波、箔条杂波、地杂波等杂波，但气象杂波不同于地杂波等其他类型的立体杂波，他有其自身的特殊性，气象杂波一般情况下并非静止不动的，其各处的厚度也不一样，分布不均匀，尤其是稀薄带分布的随机性，气象杂波的这些特点使得他的边缘情况十分复杂，在处理时也增大了难度，因为检测门限稍微偏高的话，稀薄处的气象杂波在处理时就成为了杂外；相反，如果门限稍微偏低，那么稀薄处的气象杂波在处理时就成为了 杂内，不同的判断将导致不同的处理选择。本文将重点讨论这种边缘效应的处理。一般情况下，气象杂波的回波幅度服从瑞利分布，若以x表示气象杂波回波的包络幅度，以δ表示他的平均功率(即标准差)，则x的概率密度函数为： 3建立杂波图的算法 杂波图是雷达威力范围内存储在存储器中的杂波强度分布图。在三坐标雷达中传统的平面杂波图(如图2所示)需要进行扩展，按空间在不同仰角区(如图3所示)分别建立平面杂波图，则组成了立体杂波图(也叫三维杂波图)，本文中主要研究动态杂波图和轮廓杂波图这两种立体杂波图，他们用来共同产生杂内杂外标志。 杂波图的实现是以杂波单元的划分为依据的，对于气象杂波要考虑到气象杂波运动的特点，合理选取杂波单元的大小。本文仅讨论分布相对均匀的大片气象杂波，在均匀杂波背景中，单元平均算法相对其他算法而言，其检测性能是最优的，采用检测单元△V内数据平均，本单元相邻扫描间相关积累。如图1所示，设距离分辨单元的尺寸为r，一个脉冲重复周期内扫描的角度为△α，杂波单元△V内分辨距离单元数为M，脉冲数为N，则△ρ=Mτ，△θ=N△α，任一杂波检测单元的幅度为A(i，j)，(i，j)是杂波单元的坐标，每个脉冲重复周期内，正常信号输入序列为xn(i，j)，那么： 其中，β为叠代因子，合适的叠代因子使每幅杂波图的权重分配均匀。An(i，j)为相邻两个天线扫描周期按式(3)作相关积累运算后存入同一杂波单元△V的数据，所有仰角区的杂波图的单元都以这种方式更新数据。假设噪声服从瑞利分布，则在杂波与噪声同时存在的情况下，递归后幅度分布也为瑞利分布。在设定噪声和杂波的均值和方差的情况下，根据Neyman-pearson准则，利用式(4)和式(5)在Matlab软件下进行仿真，可以得到气象杂波的虚警概率和检测概率与其相应参数的关系曲线图，如图4，图5所示。 其中，M为检测单元内的分辨单元数，图4表明，随着杂噪比SNR和分辨单元数的增大，检测概率越大；图5表明：β取1/8递归时虚警概率相应较小，同时也看到，门限设置越大，虚警概率越小。 4边缘效应的处理 当目标处于大片的气象杂波中时，由于目标周围杂波的平稳性，基本的邻近单元平均CFAR电路是性能是非常好的，但在气象杂波边缘，杂波的变化剧烈，位于各参考单元里的杂波强度有明显的差别。文献对比了CA，GO，OS，MX-CMLD等几种算法抗边缘杂波的性能，通过分析认为，包含GO(选大)逻辑的检测方法的抗边缘杂波性能均较好。为了消除边缘效应，下面将讨论如何建立轮廓杂波图，他不表征杂波的具体信息而是判断杂波相对强弱的。 杂波强弱判断准则采用双门限进行检测，建立轮廓杂波图流程。 首先判断杂波检测单元是否超过噪声门限，然后对超过噪声门限的分辨单元用计数器进行积累，当该检测单元内的计数值超过一定值，即门限2时，认为该检测单元处于强杂波区，输出一位信息1；如果该计数值小于门限2，则认为该检测单元处于弱杂波区，输出一位信息O；为了消除杂波的边缘效应，轮廓杂波图还必须对检测单元进行方位、距离的区域扩展处理，扩展的方法为：对检测单元进行判断杂波强弱的同时，也判断和他在距离、方位上相邻的参考单元的杂波强弱情况，最终将这些单元的判断结果相或，只要其中一个有杂波，就认为该检测单元处于强杂波中，扩展的单元数目一般选为8个，扩展的方法一般采用田字型。 5杂波图的实现 在雷达信号处理中除了上述的轮廓杂波图外还有其他几种类型的杂波图，下面我们只讨论动态杂波图，数字脉压后动态杂波图信号流程图如图9所示。 在正常通道信号处理中