专利名称：一种基于希尔伯特黄变换的风廓线雷达去噪方法
技术领域：
本发明涉及信号处理与检测领域，尤其涉及一种基于希尔伯特黄变换的风廓线雷达去噪方法。
背景技术：
风廓线雷达是利用大气湍流对电磁波的散射作用对大气进行探测的一种遥感设备，诞生于20世纪80年代，多年来的运行结果表明，风廓线雷达雷达满足观测精度要求，其时空分辨率超过任何高空风测量系统。但由于雷达探测的目标大气风场中的湍流对电磁波的散射非常微弱，其易受噪声、地物杂波、降水干扰和间歇性杂波干扰等影响。地物杂波是 由周围树木、建筑物等地物造成的干扰，这种干扰使雷达基础数据受到污染，进而影响到基本产品以及导出产品，其抑制主要依赖于降低天线副瓣电平和加装天线屏蔽网措施。降水杂波是在降水天气情况下，风廓线雷达接收到的信号主要是由降水粒子散射的回波组成，其功率谱密度、垂直径向速度、信噪比、谱宽等比晴空天气下都有明显的增大。对噪声的干扰抑制，Strauch等人在Colorado风廓线雷达组网时采用一种经典的处理方法，由中频接收机输出的I/Q数据经过相干积累、WFFT和谱平均等步骤，提高了风谱的可检测性，后续很多改进算法也是基于这种方法。改进算法鉴于风廓线回波窄带性质，努力提高其谱分辨率的措施，Boyer等人比较了复倒谱、参数建模谱估计算法和MUSIC算法在提取风谱中的应用，MUSIC算法有较高的谱分辨率，MUSIC算法仅仅完成对风谱的检测，风速等参量的估计需要另行处理，该算法的运算量很大，硬件实现困难。在干扰抑制，尤其是杂波对消的研究方面，在杂波功率谱服从高斯包络的前提下，Siggia和Passarelli提出了基于模型的杂波处理方法，提出用小波变换算法去除地杂波和间歇性杂波；May和Strauch、Jordan等人提出用多项式拟合对消杂波。对于降水杂波的干扰抑制，国内的卜祥元等人把分数阶傅里叶变换应用到WPR杂波抑制，选定最优的分数阶数，在相干积累的基础上进行分数阶FRFT，利用各阶矩的统计特性确定杂波干扰检测阈值，最后进行FRFT逆运算。王勇、卜祥元等利用小波变换对风廓线雷达间歇性杂波抑制方法进行了研究。张雯雯等利用提升的小波静态变换对风廓线雷达的地杂波抑制进行了研究。以上方法对噪声和杂波干扰都有一定的抑制，很多算法具有优良的性能，某些算法对硬件的要求很高，目前还不能实际采用，在噪声消除和杂波抑制处理上自适应性和边界效应改善不大。希尔伯特-黄变换(Hilbert-HuangTransform,简称 HHT)是 Huang 等人在 1998年提出的一种处理非平稳信号的新方法，其分析非稳态资料方法独特，基于信号的局部时间尺度进行信号分解，相对于其他信号处理方法的最大优势在于自适应性，因此很适合对高噪声背景下的非线性、非平稳信号的分析。
发明内容
本发明的目的在于克服现有技术的不足，提出一种基于希尔伯特黄变换的风廓线雷达去噪方法，其具体技术方案如下所述基于希尔伯特黄变换的风廓线雷达去噪方法，基于风廓线雷达去噪处理机，所述处理机包括改进型经验模式分解�？椋６乇浠荒？橐约昂笮藕糯砟？椋龇椒òㄈ缦虏街�
1)风廓线雷达去噪处理机接收回波信号，利用改进型经验模式分解模块获得回波信号的一系列的固有模式函数分量；
2)每一个固有模式函数分量通过希尔伯特变换�？榻邢６乇浠唬玫皆夹藕诺哪芰科字担�
3)通过后续信号处理�？槎栽夹藕诺哪芰科字到斜浠坏玫饺コ肷蟮姆缋呃状锸�。所述的基于希尔伯特黄变换的风廓线雷达去噪方法进一步设计在于，所述改进型经验模式分解�？榈牟僮鞑街璋ǎ�
a)计算原始回波信号的所有极值点的局部平均值使用公式(I)计算出两相邻极值点间的局部均值序列；
b) |在原始数据介于:和之间，利用公式(2)计算Λ|对应的时间；
|漏厂JC |χ( - 》
= ,J- +II* %i e Uj- J/frl ]⑵
c)利用公式(3)加权平均计算1^!处极值点的局部均值爾(Im);
也)二执卜r丨·
其中，=k(tM)=^ (4) iO- ~V2 ' V
d)使用公式(5)将局部均值Htlill)从原始信号JTii中分离出来，得到剩余部分!nil,
補=塌4Q⑶
利用公式(6)对第一个固有模式函数分量^进行第二次过滤处理，可消除叠加波从而
确保波形对称，^jfj _______ O 二 (6)
重复i次，直到满足公式(7)停止筛�。�
、Z着.)I 31^2^3^***···· (7)
I I5^lV
其中If为离散信号序列的总长度，i为重复次数，f取O. 2-0. 3之间的一个值，%是一个基本模式分量公式(8)；
权利要求
1.一种基于希尔伯特黄变换的风廓线雷达去噪方法，基于风廓线雷达去噪处理机，所述处理机包括改进型经验模式分解�？椋６乇浠荒？橐约昂笮藕糯砟？椋涮卣髟谟冢龇椒òㄈ缦虏街� 1)风廓线雷达去噪处理机接收回波信号，利用改进型经验模式分解模块获得回波信号的一系列的固有模式函数分量； 2)每一个固有模式函数分量通过希尔伯特变换�？榻邢６乇浠唬玫皆夹藕诺哪芰科字担� 3)通过后续信号处理模块对原始信号的能量谱值进行变换得到去除噪声后的风廓线雷达数据。
2.根据权利要求I所述的基于希尔伯特黄变换的风廓线雷达去噪方法，其特征在于，所述改进型经验模式分解�？榈牟僮鞑街璋� a)计算原始回波信号的所有极值点的局部平均值 m(ti+1)使用公式⑴计算出两相邻极值点间的局部均值序列；
3.根据权利要求I所述的基于希尔伯特黄变换的风廓线雷达去噪方法，其特征在干，所述希尔伯特变换�？榈牟僮鞑街璋� A)获取原始信号的时频參数，对每个固有模式函数分量分别作希尔伯特变换后得到信号的瞬时频率，
全文摘要
本发明涉及一种基于希尔伯特黄变换的风廓线雷达去噪方法，基于风廓线雷达去噪处理机，所述处理机包括改进型经验模式分解�？椋６乇浠荒？橐约昂笮藕糯砟？椋龇椒òㄈ缦虏街�1)风廓线雷达去噪处理机接收回波信号，利用改进型经验模式分解�？榛竦没夭ㄐ藕诺囊幌盗械墓逃心Ｊ胶至浚�2)每一个固有模式函数分量通过希尔伯特变换�？榻邢６乇浠唬玫皆夹藕诺哪芰科字担�3)通过后续信号处理�？槎栽夹藕诺哪芰科字到斜浠坏玫饺コ肷蟮姆缋呃状锸荨Ｆ溆幸嫘Ч哂辛己玫氖视π裕芄幻飨缘娜コ肷确瞧轿雀扇牛绕淠芄挥行诮德实幕斓跋旌捅呓缧в�。
文档编号G01S13/95GK102736070SQ20121023960
公开日2012年10月17日 申请日期2012年7月12日 优先权日2012年7月12日
发明者单慧琳, 周杰, 张银胜, 李家强 申请人:南京信息工程大学
专业数控铣床产品供应商
20年生产经验，实力工厂，精工品质，质量有保障
