一种双偏振天气雷达回波信号仿真方法

摘要

本发明涉及雷达仿真技术领域，具体为一种双偏振天气雷达回波信号仿真方法。该仿真方法包括双偏振天气雷达回波信号仿真方法和双偏振雷达体扫模式回波信号仿真方法。其中双偏振天气雷达回波信号仿真方法步骤为：获取双偏振天气雷达水平通道和垂直通道的回波信号功率；进行水平通道和垂直通道回波信号的复频谱建模，建立水平通道和垂直通道回波信号的相位差、强度差、相关性；根据两路正交信号生成时域I/Q回波信号；模拟接收机水平通道和垂直通道回波信号噪声和通道增益。该双偏振天气雷达回波信号仿真方法，建立了双偏振天气雷达系统性能参数与回波信号之间关系，可以仿真雷达性能参数对雷达观测量的影响。

权利要求

1.一种双偏振天气雷达回波信号仿真方法，其特征在于，该方法包括双偏振天气雷达回波信号仿真方法和双偏振雷达体扫模式回波信号仿真方法，其中，所述双偏振天气雷达回波信号仿真方法包括以下步骤：

S1：获取双偏振天气雷达水平通道和垂直通道的回波信号功率；

S2：进行水平通道和垂直通道回波信号的复频谱建模，建立水平通道和垂直通道回波信号的相位差、强度差、相关性；

S3：根据水平通道的两路正交信号生成水平通道的时域I/Q回波信号，根据垂直通道的两路正交信号生成垂直通道的时域I/Q回波信号；

S4：模拟接收机水平通道和垂直通道回波信号噪声和通道增益；

步骤S1中的具体步骤包括：根据天气雷达水平通道反射率因子、雷达波长、天线增益、水平通道发射峰值功率、脉冲宽度、波束宽度、大气损耗、除大气损耗外的总损耗、目标距离雷达的距离、降雨衰减系数、经验常数计算获得所述双偏振天气雷达水平通道回波信号功率；根据天气雷达差分反射率因子、大气损耗、除大气损耗外的总损耗、天气雷达水平通道反射率因子、雷达波长、天线增益、垂直通道发射峰值功率、脉冲宽度、波束宽度、目标距离雷达的距离、降雨衰减系数、经验常数计算获得所述双偏振天气雷达垂直通道回波信号功率；

步骤S2中选用高斯普模型作为回波信号的归一化功率谱，建立水平通道和垂直通道回波信号的相位差、强度差、相关性具体步骤包括：根据频率域的谱宽、速度域的谱宽、多普勒频率、雷达径向速度建立回波信号的功率谱，计算方式为 ，其中 为频率域的谱宽，单位为Hz，；W为速度域的谱宽，即雷达基数据中谱宽数据，单位为m/s； 为多普勒频率 ，单位为Hz， ； 为雷达波长，单位为cm ；V为雷达径向速度，即雷达基数据中谱宽数据，单位为m/s；

使用傅里叶变换获得水平通道和垂直通道功率谱随机化噪声，建立水平通道和垂直通道回波信号的相关性，两者之间的相关性由以下关系式表示 ，，其中 、 分别为水平通道和垂直通道功率谱随机化噪声， 为水平通道零均值的高斯白噪声信号的傅里叶变换， 为垂直通道零均值的高斯白噪声信号的傅里叶变换， 为零滞后相关系数，即双偏振天气雷达的极化参量数据相关系数；

根据回波信号的功率谱、雷达脉冲积累数、脉冲重复频率、差分传播相位、随机相位、水平通道功率谱随机化噪声建立水平通道回波信号的复频率模型，水平通道回波信号的复频谱模型表示方式为 ，其中 为水平通道的回波信号功率，单位为 ； ， N为雷达脉冲积累数， 为脉冲重复频率， 为随机相位，单位弧度；

根据回波信号的功率谱、雷达脉冲积累数、脉冲重复频率、差分传播相位、随机相位、垂直通道功率谱随机化噪声建立垂直通道回波信号的复频率模型，垂直通道回波信号的复频率模型表示方式为 ，其中

为垂直通道的回波信号功率，单位为 ； 为差分传播相位，单位为°；

将水平通道回波信号的不同强度和相位代入水平通道回波信号的复频率模型，将垂直通道回波信号的不同强度和相位代入垂直通道回波信号的复频率模型，获得水平通道和垂直通道回波信号的相位差和强度差。

2.根据权利要求1所述的双偏振天气雷达回波信号仿真方法，其特征在于，步骤S4中，模拟接收机水平通道和垂直通道回波信号噪声和通道增益包括具体步骤：使用接收机噪声系数获得接收机噪声等效功率；

使用接收机通道增益来计算接收机通道到数字中频过程中对回波信号的功率增益。

3.根据权利要求1所述的双偏振天气雷达回波信号仿真方法，其特征在于：所述双偏振雷达体扫模式回波信号仿真方法包括连续监测模式、连续多普勒模式、批处理模式、CDX模式、双PRF模式、参差PRF模式、SZ‑2相位编码模式下的双偏振天气雷达回波信号仿真方法。

4.根据权利要求3所述的双偏振天气雷达回波信号仿真方法，其特征在于，所述连续监测模式下的双偏振天气雷达回波信号仿真方法包括具体步骤：根据天气雷达的脉冲重复频率PRF计算雷达最大不模糊速度；

根据当前需要仿真的真实天气雷达回波速度计算回波速度模糊后的速度值；

根据回波速度模糊后速度值模拟仿真产生速度模糊的回波信号特征；

按照S1‑S4所述的双偏振天气雷达回波信号仿真方法，仿真连续监测模式下的双偏振天气雷达回波信号特征。

5.根据权利要求3所述的双偏振天气雷达回波信号仿真方法，其特征在于，所述连续多普勒模式下的双偏振天气雷达回波信号仿真方法包括具体步骤：当回波速度很大并且超过最大模糊速度范围时，按照所述连续监测模式下的双偏振天气雷达回波信号仿真方法，进行连续监测模式下的双偏振天气雷达回波信号特征仿真；

根据天气雷达当前工作的脉冲重复周期PRT获得最大不模糊距离；

根据当前需要仿真的回波位置，判断回波是否会产生距离折叠，若会产生距离折叠则计算会产生多少次距离折叠以及距离折叠的位置；

使用时域混叠方法模拟回波信号的折叠过程。

6.根据权利要求3所述的双偏振天气雷达回波信号仿真方法，其特征在于，所述批处理模式下的双偏振天气雷达回波信号仿真方法包括具体步骤：多次重复发射长脉冲重复周期和短脉冲重复周期的脉冲，仿真产生整个批处理模式下的双偏振天气雷达回波信号。

7.根据权利要求3所述的双偏振天气雷达回波信号仿真方法，其特征在于，所述双PRF模式下的双偏振天气雷达回波信号仿真方法包括具体步骤：（a）对每一个距离库单元，以4：3或3：2的频率多次重复发射长脉冲重复周期和短脉冲重复周期的脉冲产生回波；

（b）在回波信号序列排序和组合输出时，按照先长脉冲重复周期、后短脉冲重复周期的序列组合顺序输出；

其他距离库按照步骤（a）、（b）产生回波，循环所有距离库和所有径向。

8.根据权利要求3所述的双偏振天气雷达回波信号仿真方法，其特征在于，所述SZ‑2相位编码模式下的双偏振天气雷达回波信号仿真方法包括具体步骤：天气雷达发射长脉冲重复周期的发射脉冲信号，扫描一圈，获得长脉冲重复周期的回波信号；

发射具有SZ相位编码的短脉冲重复周期扫描一圈，获取具有相位编码的回波信号；

在信号处理时，结合长脉冲重复周期下不容易距离模糊的回波功率和位置，使用SZ‑2处理算法恢复距离模糊出的反射率值。

说明书

一种双偏振天气雷达回波信号仿真方法

技术领域

[0001] 本发明涉及雷达仿真技术领域，具体涉及一种双偏振天气雷达回波信号仿真方法。

背景技术

[0002] 双偏振天气雷达是当前对灾害性天气进行监测、预警应用最广泛、最有效的工具之一。双偏振天气雷达主要优势在于它除了能够获取降雨过程的强度、速度信息外，还可以
通过探测获取的极化信息，反演降雨的形态、相态等微物理结构，为降雨的定量测量、灾害
性天气的预警提供了丰富的数据基础。

[0003] 随着气象业务朝精细化、精准性的高速发展，对双偏振天气雷的系统性能、处理算法改进、扫描策略也提出了更高的要求。而现在雷达厂家或雷达业务使用部门，在解决这些
问题时，往往需要首先研制雷达设备，再经过长时间观测试验，数据采集和数据分析，才能
得到雷达性能提高、扫描策略改进对雷达观测数据的影响。

[0004] 现有双偏振天气雷达回波仿真，大多数是基于降水回波的物理特性来模拟双偏振天气雷达的谱矩参量和双偏振参量，如根据雨滴大小分布的二维时空随机模型，仿真双偏
振天气雷达反射率、差分反射率、差分相位等，目的是利用已知的雨滴分布，来验证双偏振
天气雷达偏振参量与实际降雨物理特性之间的关系。该方法只能仿真双偏振天气雷达的偏
振参量，不能仿真双偏振天气雷达的回波信号，不能建立与双偏振天气雷达系统参数，体积
扫描工作模式之间关系。

[0005] 现有多普勒天气雷达回波信号仿真方法，建立了基于高斯功率模型的回波信号仿真方法，能够以回波功率、速度和谱宽为基础，仿真单偏振多普勒天气雷达的回波I/Q信号。
该方法不能仿真双偏振天气雷达的双通道回波信号，没有建立回波信号与双偏振天气雷达
系统参数之间的关系，没有考虑体积扫描中工作模式的不同特点，无法仿真天气雷达体扫
过程中各个扫描模式的回波特点。

[0006] 现有方法，不能应用于双偏振天气雷达的雷达系统性能评估、新算法验证与分析、扫描策略分析等领域。

发明内容

[0007] 基于上述问题，本发明提供了一种双偏振天气雷达回波信号仿真方法，能够仿真产生含有偏振信息的垂直通道和水平通道的回波信号，以及仿真产生不同工作模式下的双
偏振天气天气雷达回波信号。

[0008] 本发明的技术方案为：

[0009] 一种双偏振天气雷达回波信号仿真方法，包括双偏振天气雷达回波信号仿真方法和双偏振雷达体扫模式回波信号仿真方法。

[0010] 其中，双偏振天气雷达回波信号仿真的主要目标是从反射率因子、速度、谱宽、差分反射率、差分相位、相关系数六种参数，通过数学建模和雷达性能参数引入，模拟仿真得
到设定雷达性能参数下的双偏振天气雷达水平通道和垂直通道的回波信号。在数学建模
中，引入的双偏振天气雷达的雷达参数包括发射峰值功率、雷达波长、脉冲宽度、波束宽度、
天线增益、接收支路馈线损耗、发射支路馈线损耗、噪声系数、接收机增益、脉冲重复频率、
脉冲积累数等。

[0011] 双偏振天气雷达回波信号仿真方法包括以下步骤：

[0012] S1：获取双偏振天气雷达水平通道和垂直通道的回波信号功率；

[0013] S2：进行水平通道和垂直通道回波信号的复频谱建模，建立水平通道和垂直通道回波信号的相位差、强度差、相关性；

[0014] S3：根据水平通道的两路正交信号生成水平通道的时域I/Q回波信号，根据垂直通道的两路正交信号生成垂直通道的时域I/Q回波信号；

[0015] S4：模拟接收机水平通道和垂直通道回波信号噪声和通道增益。

[0016] 进一步的，步骤S1中的具体步骤包括：

[0017] 根据天气雷达水平通道反射率因子、雷达波长、天线增益、水平通道发射峰值功率、脉冲宽度、波束宽度、大气损耗、除大气损耗外的总损耗、目标距离雷达的距离、降雨衰
减系数、经验常数计算获得所述双偏振天气雷达水平通道回波信号功率，具体计算方式为
，其中
为水平通道的回波信号功率，单位为 ； 为天气雷达水平通道反射率因子，
单位为 ； 为雷达波长，单位为cm；G为天线增益，单位为dB； 为雷达水平通道发
射峰值功率，单位为kw； 为脉冲宽度，单位为 ； 分别为水平方向波束宽度
和垂直方向波束宽度，单位为°； 为大气损耗，单位为dB/km； 为除 外的总损耗，
包括水平通道发射支路的总馈线损耗和水平通道接收馈线支路的损耗总和，单位为dB；R为
目标离雷达的距离，单位为km， 为在雷达距离r处降雨的衰减系数，降雨衰减系数
值由反射率因子计算得到，公式为 ；a和b为经验常数，对于C波段天气
雷达， 和 ，对于X波段天气雷达， 和 。

[0018] 根据天气雷达差分反射率因子、大气损耗、除大气损耗外的总损耗、天气雷达水平通道反射率因子、雷达波长、天线增益、垂直通道发射峰值功率、脉冲宽度、波束宽度、目标
距离雷达的距离、降雨衰减系数、经验常数计算获得所述双偏振天气雷达垂直通道回波信
号功率，具体计算方式为 ，
其中 为垂直通道的回波信号功率，单位为 ； 为差分反射率因子，单位为ｄ
B， 是除 外的垂直通道的总损耗，包括垂直通道发射支路的总馈线损耗和垂直通
道接收馈线支路的损耗总和，单位为ｄB，其它参数参见双偏振天气雷达水平通道回波信号
功率的计算方式。

[0019] 进一步的，步骤S2中选用高斯普模型作为回波信号的归一化功率谱，建立水平通道和垂直通道回波信号的相位差、强度差、相关性具体步骤包括：

[0020] 根据频率域的谱宽、速度域的谱宽、多普勒频率、雷达径向速度建立回波信号的功率谱，计算方式为 ，其中 为频率域的谱宽，单位为Hz，
；W为速度域的谱宽，即雷达基数据中谱宽数据，单位为m/s； 为多普勒频
率 ，单位为Hz， ；V为雷达径向速度，即雷达基数据中谱宽数据，单位为m/s；

[0021] 使用傅里叶变换获得水平通道和垂直通道功率谱随机化噪声，建立水平通道和垂直通道回波信号的相关性，两者之间的相关性由以下关系式表示
，其中
分别为水平通道和垂直通道功率谱随机化噪声， 为水平通道零均值的高斯白噪声
信号的傅里叶变换， 为垂直通道零均值的高斯白噪声信号的傅里叶变换，
为零滞后相关系数，即双偏振天气雷达的极化参量数据相关系数；

[0022] 根据回波信号的功率谱、雷达脉冲积累数、脉冲重复频率、差分传播相位、随机相位、水平通道功率谱随机化噪声建立水平通道回波信号的复频率模型，水平通道回波信号
的复频谱模型表示方式为 ，
其中  N为雷达脉冲积累数， 为脉冲重复频
率， 为随机相位，单位弧度；

[0023] 根据回波信号的功率谱、雷达脉冲积累数、脉冲重复频率、差分传播相位、随机相位、垂直通道功率谱随机化噪声建立垂直通道回波信号的复频率模型，垂直通道回波信号
的复频率模型表示方式为 ，其
中 为差分传播相位，单位为°；

[0024] 将水平通道回波信号的不同强度和相位代入水平通道回波信号的复频率模型，将垂直通道回波信号的不同强度和相位代入垂直通道回波信号的复频率模型，获得水平通道
和垂直通道回波信号的相位差和强度差。

[0025] 进一步的，步骤S３中，根据水平通道的两路正交信号生成水平通道的时域I/Q回波信号，表示方式为 ，其中
和 分别表示水平通道的两路正交信号，ｎ表示每个雷达径向脉冲的序列号；

[0026] 根据垂直通道的两路正交信号生成垂直通道的时域I/Q回波信号，表示方式为，其中 和 分别表
示垂直通道的两路正交信号，ｎ表示每个雷达径向脉冲的序列号。

[0027] 进一步的，步骤S4中，模拟接收机水平通道和垂直通道回波信号噪声和通道增益包括具体步骤：

[0028] 使用接收机噪声系数获得接收机噪声等效功率，水平通道接收机噪声等效功率计算方式为 ，其中 为波尔茨曼常数，
为常用室温温度290K， 为雷达接收机带宽， 为水平通道的接收机噪声系数，垂直通
道接收机噪声等效功率计算方式为 ， 为垂直通道的接收机噪
声系数；

[0029] 使用接收机通道增益来计算接收机通道到数字中频过程中对回波信号的功率增益。

[0030] 考虑接收机噪声和通道增益的水平通道和垂直通道的时域I/Q回波信号表示方式为：

[0031] 其中 和 分别为水平通道和垂直通道的接收机通道增益 和 为高斯随机噪声，用
来表示水平通道和垂直通道的接收机噪声； 和 分别为水平通道和垂直通道的接
收机噪声功率值。

[0032] 通过上述步骤Ｓ１－Ｓ４，实现了一个以反射率因子、速度、谱宽、差分反射率、差分相位、相关系数六种基数据量为输入的双偏振天气雷达I／Q回波信号仿真，同时增加了发射
峰值功率、雷达波长、脉冲宽度、波束宽度、天线增益、接收支路馈线损耗、发射支路馈线损
耗、噪声系数、接收机增益、脉冲重复频率、脉冲积累数等重要雷达参数的定量化计算方法，
使得生成的回波信号，既仿真了天气目标时域和频域特性，又仿真了天气雷达系统的重要
性能。

[0033] 目前新一代多普勒天气雷达业务体扫策略有VCP21,VCP11,VCP31,VCP32四种，广泛采用的是VCP21扫描策略。在四种体扫策略中，一般低仰角采用CS、CD工作模式，中仰角采
用批处理模式，高仰角采用CDX模式。此外，部分业务雷达在CD模式下采用了SZ‑2相位编码，
以减低距离模糊的影响，在批处理模式，采用了双重复频率模式（双PRF）或者参差重复频率
模式（参差PRF），以降低速度模糊的影响。以下为双偏振天气雷达体扫策略中不同工作模式
的回波信号仿真。

[0034] 双偏振雷达体扫模式回波信号仿真方法包括连续监测模式、连续多普勒模式、批处理模式、CDX模式、双PRF模式、参差PRF模式、SZ‑2相位编码模式下的双偏振天气雷达回波
信号仿真方法。

[0035] 进一步的，（CS）连续监测模式下的双偏振天气雷达回波信号仿真方法包括具体步骤：

[0036] 根据天气雷达的脉冲重复频率PRF计算雷达最大不模糊速度，计算方式为；

[0037] 根据当前需要仿真的真实天气雷达回波速度V计算回波速度模糊后的速度值 ，计算方式为 ，其中 为雷达回波速
度模糊后的速度值，Ｖ为需要仿真的真实雷达回波速度值；Ｋ为速度模糊次数取值。当V速
度值小于 时，Ｋ取值为0，表示速度不发生折叠， ；当Ｖ速度值为正，大于1倍
，且小于2倍 时，此时速度为1次折叠，Ｋ为取‑1，折叠后的速度值为负速度值；
若Ｖ速度值为负，小于1倍 时且大于2倍 时，此时速度为负方向的1次折叠，
Ｋ为取1次，折叠后的速度值为正速度值。实际过程中，重复频率和待仿真速度值得不同，可
能出现多次速度重叠，速度折叠方法类推；

[0038] 根据回波速度模糊后速度值模拟仿真产生速度模糊的回波信号特征；

[0039] 按照S1‑S4所述的双偏振天气雷达回波信号仿真方法，仿真连续监测模式下的双偏振天气雷达回波信号特征。

[0040] 进一步的，（CD）连续多普勒模式下的双偏振天气雷达回波信号仿真方法包括具体步骤：

[0041] 当回波速度很大并且超过最大模糊速度范围时，按照所述连续监测模式下的双偏振天气雷达回波信号仿真方法，进行连续监测模式下的双偏振天气雷达回波信号特征仿
真。

[0042] 根据天气雷达当前工作的脉冲重复周期PRT获得最大不模糊距离。

[0043] 根据当前需要仿真的回波位置，判断回波信号是否会产生距离折叠，若会产生距离折叠则计算会产生多少次距离折叠以及距离折叠的位置，回波信号的距离折叠位置的计
算方式为 ，其中R为雷达回波信号的实
际位置， 为雷达产生距离折叠后的位置，K为距离模糊次数取值。当R距离值小于
时，K取值为0，表示距离不发生折叠，  =R；当R距离大于1倍 且小于两倍
时，此时距离为1次折叠，K取值为1；当R距离大于2倍 且小于3倍 时，此时
距离为2次折叠，K取值为2；以此类推。

[0044] 使用时域混叠方法模拟回波信号的折叠过程，加入同一方位角度距离位置R1和距离位置R2，且R2大于 ，通过回波信号的距离折叠位置的计算方式折叠后的距离位置
刚好等于R1，则此时雷达收到的回波是距离位置R1和距离位置R2回波信号的叠加，时域
信号叠加计算方式为 ，其中
和 分别表示距离位置R1和距离位置R2垂直通道的原始回波信号，这个信号
是通过考虑接收机噪声和通道增益的水平通道的时域I/Q回波信号计算方式、考虑接收机
噪声和通道增益的垂直通道的时域I/Q回波信号计算方式产生的，n和n‑1表示当前脉冲与
上一个脉冲的远距离回波距离混叠。

[0045] 在天气雷达体扫策略中，通常CS模式和CD模式是配合使用，如先进行CS模式下的方位一圈PPI扫描，然后雷达仰角不变，改变雷达脉冲重复频率和雷达积累数参数，雷达又
进行方位一圈的PPI扫描。对于回波信号仿真来说，产生CS模式和CD模式回波信号，即按照
上述计算过程，基于同一雷达基数据场，分别进行CS模式和CD模式回波信号模拟仿真。

[0046] 进一步的，批处理模式下的双偏振天气雷达回波信号仿真方法包括具体步骤：

[0047] 多次重复发射长脉冲重复周期PRT1和短脉冲重复周期PRT2的脉冲（PRT1< PRT2），仿真产生整个批处理模式下的双偏振天气雷达回波信号。从回波信号产生方面来说，即是
先产生一个径向的CS模式下回波信号，然后又产生一个径向CD模式回波信号，以此类推，仿
真产生整个PPI回波信号。

[0048] 批处理模式下的回波信号可以如下表示：

[0049]

[0050] 其中 和 为长脉冲重复周期PRT1和短脉冲重复周期PRT2下的回波信号序列。因此，批处理模式回波信号的仿真实现可看作是CS模式和CD模
式回波信号排序组合。

[0051] 进一步的，在CDX模式下，雷达一般工作在高仰角、高重复频率，雷达回波不容易距离模糊，但高空风速较大，容易速度模糊。所以CDX模式可按照CD模式来仿真产生雷达回波
信号。

[0052] 进一步的，雷达工作在双PRF模式时，从脉冲发射方式来说，与批处理模式类似，雷达先发射一组高脉冲重复频率(PRF1)脉冲，然后再发射一组低脉冲重复频率(PRF2)脉冲，
并以此重复。与批处理模式不同的是，双PRF模式采用的重复频率要高，且PRF1:PRF2通常是
4：3，或3：2等， 然后在信号处理时，双PRF模式与批处理模式在速度计算方法不同。

[0053] 双PRF模式下的双偏振天气雷达回波信号仿真方法包括具体步骤：

[0054] 对每一个距离库单元，以4：3或3：2的频率多次重复发射长脉冲重复周期和短脉冲重复周期的脉冲产生回波；

[0055] 在回波信号序列排序和组合输出时，按照先长脉冲重复周期、后短脉冲重复周期的序列组合顺序输出。

[0056] 双PRF模式下，距离库单位R的水平通道和垂直通道的回波信号排序组合如下表示：

[0057] ；，其中
表示在距离位置R的水平通道、脉冲重复频率PRF1下第n个回波信号。

[0058] 其他距离库按照上述流程产生回波，循环所有距离库和所有径向。

[0059] 进一步的，雷达工作在参差PRF模式时，从脉冲发射方式来说，与批处理模式和双PRF模式稍微不同。参差PRF模式下，雷达是先发射一个高脉冲重复频率(PRF1)脉冲，然后再
发射一个低脉冲重复频率(PRF2)脉冲，并以此重复发射N次。参差PRF模式采用的重复频率
要高，且PRF1:PRF2通常是4：3，或3：2等， 然后在信号处理时，参差PRF模式与双PRF模式在
速度计算方法相同。

[0060] 参差PRF模式下的双偏振天气雷达回波信号仿真方法包括具体步骤：

[0061] 对每一个距离库单元，以4：3或3：2的频率多次重复发射长脉冲重复周期和短脉冲重复周期的脉冲产生回波；

[0062] 在回波信号序列排序和组合输出时，按照先长脉冲重复周期、再短脉冲重复周期、然后长脉冲重复周期、最后短脉冲重复周期的序列组合顺序输出，距离库单位R的水平通道
和垂直通道的回波信号排序组合表示如下

[0063]

[0064]

[0065] 其他距离库按照上述流程产生回波，循环所有距离库和所有径向。

[0066] 进一步的，天气雷达在使用SZ‑2相位编码模式进行退距离模糊处理时，同样采用类似CS、CD模式的分离扫描方式；首先，雷达发射长脉冲重复周期的发射脉冲信号，扫描一
圈，获得长脉冲重复周期的回波信号，然后发射具有SZ（8/64）相位编码的短脉冲重复周期
扫描一圈，获取具有相位编码的回波信号。然后在信号处理时，结合长脉冲重复周期下不容
易距离模糊的回波功率和位置，使用SZ‑2处理算法恢复距离模糊出的反射率值。

[0067] SZ‑2相位编码模式下的回波信号仿真主要包含长脉冲重复周期PRT1、无相位编码的回波信号仿真，和短脉冲重复周期PRT2、SZ(8/64)相位编码的回波信号仿真两个算法过
程。长脉冲重复周期PRT1、无相位编码的回波信号仿真采用CS模式同样的仿真方法和流程，
主要注意速度模糊的仿真处理。

[0068] 短脉冲重复周期PRT2、SZ(8/64)相位编码的回波信号仿真具体过程如下：

[0069] （1）首先按照双偏振天气雷达回波信号仿真方法产生单元距离库的水平和垂直通道的回波信号序列 和 ；

[0070] （2）对 和 进行相位调制，调制公式表示如下：

[0071]

[0072] 其中m和M是相位编码参数。天气雷达常采用SZ(8/64)编码进行调制输出，即m=8,M=64。上述调制公式中的相位以8为周期
变化，8个相位分别为:   

[0073] （3）然后根据CD模式下距离模糊判断方法和距离折叠位置计算方法，对每个距离库进行距离判断和折叠计算，然后，再对编码后回波序列进行距离模糊混叠处理，处理当时
如下：

[0074]  

[0075] （4）然后，其他距离库，按照上述(1)‑(2)流程产生回波，然后对回波序列进行排序组合，循环所有距离库和所有径向，即可仿真得到短脉冲重复周期PRT2、SZ(8/64)相位编码
的回波信号。

[0076] 本发明的有益之处在于：

[0077] （１）该双偏振天气雷达回波信号仿真方法，建立了双偏振天气雷达含有双偏振信息的水平通道和垂直通道回波I/Q信号的仿真方法，解决了水平通道和垂直通道回波信号
之间包含的功率差、相位差、相关性等偏振信息关联问题。可广泛应用于双偏振天气雷达的
算法处理评估中。

[0078] （２）该双偏振天气雷达回波信号仿真方法，建立了双偏振天气雷达系统性能参数与回波信号之间关系，可以仿真雷达性能参数对雷达观测量的影响。在雷达组件或系统硬
件改进之前，可以通过该仿真方法，快速、方便、灵活地评估改进带来的效果，大大地降低了
生产和设计成本，提高试验效率。

[0079] （３）该双偏振天气雷达回波信号仿真方法，考虑了双偏振天气雷达体积扫描策略中各层工作模式的不同特点，如CS模式下，脉冲重复频率低，容易速度模糊，CD模式下，脉冲
重复频率低，容易距离模糊。能够仿真双偏振天气雷达体积扫描中各个工作模式下的回波
仿真信号的不同。可广泛应用于双偏振天气雷达的算法处理、扫描测量验证与分析评估中。

附图说明

[0080] 图1为本发明中双偏振天气雷达回波信号仿真方法流程图。

[0081] 图2为实施例1中H通道回波信号和V通道回波信号的功率谱表示。

[0082] 图3为实施例1中H通道回波信号和V通道回波信号的时域表示，（a）为H通道回波信号的时域表示，（b）为V通道回波信号的时域表示。

[0083] 图4为实施例2中实际天气过程中的反射率因子PPI图。

[0084] 图5为实施例2中实际天气过程中的速度PPI图。

[0085] 图6为实施例2中实际天气过程中的谱宽PPI图。

[0086] 图7为实施例2中实际天气过程中的差分反射率PPI图。

[0087] 图8为实施例2中实际天气过程中的差分相位PPI图。

[0088] 图9为实施例2中仿真的反射率因子PPI图。

[0089] 图10为实施例2中仿真中的速度PPI图。

[0090] 图11为实施例2中仿真中的谱宽PPI图。

[0091] 图12为实施例2中仿真的差分反射率PPI图。

[0092] 图13为实施例2中仿真的差分相位PPI图。

[0093] 图14为实施例3中实际天气过程的反射率因子PPI图。

[0094] 图15为实施例3中实际天气过程的速度PPI图。

[0095] 图16为实施例3中退距离模糊的反射率因子PPI图。

[0096] 图17为实施例3中退距离模糊的速度PPI图。

具体实施方式

[0097] 下面结合附图，对本发明作详细的说明。

[0098] 所举实施例是为了更好地对本发明进行说明，但并不是本发明的内容仅局限于所举实施例。所以熟悉本领域的技术人员根据上述发明内容对实施方案进行非本质的改进和
调整，仍属于本发明的保护范围。

[0099] 实施例1

[0100] 本实施例为双偏振天气雷达回波信号仿真方法，仿真流程参见图1。为实现单个目标点的双偏振回波信号仿真，其具体步骤为：

[0101] （1）读取基数据，读取反射率 、速度 、谱宽 、差分反射率 、差分相位 、相关系数 六种数据，或者自己设置，为本实施例中的公式做数据准备；

[0102] （2）根据仿真需要，设置预定的雷达系统参数 等，为本实施例中的公式做数据准备；

[0103] （3）将上述参数代入水平通道回波信号功率计算方式和垂直通道回波信号功率计算方式，得到回波功率 和 ，再将 和 代入功率谱计算方式和水平通道回波信
号的复频谱模型计算方式中，得到水平通道和垂直通道的回波复频谱 和 ，
最后将所得复频谱，通过考虑接收机噪声和通道增益的水平通道和垂直通道的时域I/Q回
波信号计算方式，即可得到水平通道和垂直通道的回波序列 和 。

[0104] 以某C波段双偏振天气雷达基数据为例，将该天气雷达的系统参数作为仿真算法的输入，如发射功率： ；脉冲宽度： ；水平波束宽度：
，垂直波束宽度： ，脉冲重复频率：
等。然后从其某次实际扫描过程中读取数据，如反射率因子： ；多
普勒速度： ；谱宽： ；差分反射率因子：
；差分相位： ；零滞后相关系数：
。将以上参数按上述双偏振天气雷达回波仿真算法实施，即可得到所需的H通道和V通道的
回波信号。图2即为H通道回波信号和V通道回波信号的功率谱表示，图3即为H通道回波信号
和V通道回波信号的时域表示，对H通道和V通道的回波信号使用信号处理算法计算基数据
可得反射率因子： 多普勒速度： ；
谱宽： ；差分反射率因子： ；差分
相位： 。

[0105] 实施例2

[0106] 本实施例为一个仰角所有目标点的双偏振天气雷达回波信号仿真方法。其实际天气过程中的反射率因子PPI图参见图4，实际天气过程中的速度PPI图参见图5，实际天气过
程中的谱宽PPI图参见图6，实际天气过程中的差分反射率PPI图参见图7，实际天气过程中
的差分相位PPI图参见图8。为实现一个仰角所有目标点的双偏振回波信号仿真，首先对该
仰角内的所有目标循环执行实例1的步骤，即可完成一个仰角所有目标点的双偏振回波信
号仿真。

[0107] 以某C波段双偏振天气雷达基数据为例，对该次扫描过程所有目标进行回波仿真并对得到的信号进行使用信号处理算法计算基数据，可得图9‑13所示仿真结果的PPI图。仿
真中的反射率因子PPI图参见图9，仿真中的速度PPI图参见图10，仿真中的谱宽PPI图参见
图11，仿真中的差分反射率PPI图参见图12，仿真中的差分相位PPI图参见图13。

[0108] 实施例3

[0109] 本实施例为不同体扫模式下的双偏振雷达回波信号仿真，首先仿真CS模式下的双偏振天气雷达回波信号，其他模式下以此为基础进行仿真。

[0110] 在CS模式下，取仿真时的脉冲重复频率： ，此时的最大不模糊速度 ，其它系统参数与具体实施例（1）中的参数一致。对于实
施例（1）中所用的回波数据，此时的脉冲重复频率： ，此时的
最大不模糊速度 。故在CS模式下需对大于 的散射体进
行速度模糊仿真。如在实际天气雷达数据中选择反射率因子： ；
多普勒速度： ；谱宽： ；差分反射率因子： ；
差分相位： ；零滞后相关系数： 。由于此时的速度
，故需对仿真的数据做速度模糊处理。具体步骤如下：

[0111] （1）首先按照水平通道回波信号功率计算方式和垂直通道回波信号功率计算方式，求回波功率 和 。

[0112] （2）根据 大小判断为速度一次模糊， ，然后根据回波速度模糊后的速度值，计算的模糊后的速度值 。

[0113] （3）再将 和 代入功率谱计算方式和水平通道回波信号的复频谱模型计算方式，得到水平通道和垂直通道的回波复频谱 和

[0114] （4）最后将所得复频谱，通过考虑接收机噪声和通道增益的水平通道和垂直通道的时域I/Q回波信号计算方式，即可得到水平通道和垂直通道的回波序列 和
。

[0115] （5）在一个仰角对所有目标执行，循环以上步（1）,（2）,（3）,（4）步骤，即可完成CS模式下的所有目标的回波信号仿真。

[0116] 在CD模式下，取仿真时的脉冲重复频率： ，此时的最大不模糊距离 。 此时若选择探测距离大于Rmax的数据，则会发生
距离模糊，因此需对仿真数据做距离模糊处理。具体步骤如下：

[0117] （1）首先按照CS模式计算方式，得到水平通道回波信号和垂直通道回波信号，

[0118] 和 。

[0119] （2）根据 与 大小判断判断回波信号是否会产生距离折叠，若会产生距离折叠，根据回波信号距离折叠位置的计算公式，得到发生折叠的位置R1和R2。

[0120] （3）根据距离叠加位置，按照时域信号叠加计算方式，得到叠加后的回波信号和 。

[0121] （4）在一个仰角对所有目标执行，循环以上步（1）,（2）,（3）步骤，即可完成CD模式下的所有目标的回波信号仿真。

[0122] 从某C波段天气雷达实际扫描中读取数据，分别进行CS模式和CD模式下的回波信号仿真，并对仿真结果做退距离模糊处理。本实施例中的实际天气过程的反射率因子PPI图
参见图14，实际天气过程的速度PPI图参见图15，退距离模糊的反射率因子PPI图参见图16，
退距离模糊的速度PPI图参见图17。

[0123] 最后说明的是，以上实施例仅用以说明本发明的技术方案而非限制，尽管参照较佳实施例对本发明进行了详细说明，本领域的普通技术人员应当理解，可以对本发明的技
术方案进行修改或者等同替换，而不脱离本发明技术方案的宗旨和范围，其均应涵盖在本
发明的权利要求范围当中。