《南京大学c波段双偏振雷达_对双通道接收机的噪声功率_分析报告_2014年10月21日》由会员分享，可在线阅读，更多相关《南京大学c波段双偏振雷达_对双通道接收机的噪声功率_分析报告_2014年10月21日（33页珍藏版）》请在金锄头文库上搜索。
1、南京大学C波段双偏振雷达对双通道接收机的噪声功率分析报告南京大学2014年10月22日目 录1 概述12 分析方法33 接收机茅草的分析结果43.1 对每个IQ数据文件的单独分析43.2 对每个IOP的综合分析93.2.1 IOP193.2.2 IOP2113.2.3 IOP3123.2.4 IOP4143.2.5 IOP5153.2.6 IOP6173.2.7 IOP7183.2.8 IOP8193.2.9 IOP9203.2.10 IOP10233.2.11 IOP11244 初步结论255 下一步工作265.1 把雷达车上的空调制冷修好265.2 更换更强大制冷能力的恒温箱，实时记录温度
2、265.3 在关于低SNR下，LDR结果出现系统偏差的问题_分析报告中指出，研究更高精度的接收机噪声的实时估计方法，并应用到信号处理之中去265.4 把分析的噪声功率起伏的数据代入信号处理的参数中，重新计算基数据266 附录1：恒温箱的性能（根据该产品的说明书）277 附录2：在溧水山上观测时，该雷达的接收机噪声稳定性281 概述2014年67月份，由北京敏视达公司制造的南京大学C波段双偏振雷达在安徽长丰站进行了持续的降水观测，总计有11个IOP观测。在每次体扫观测中，都进行了垂直90度扫描。我们对垂直90度观测的数据进行了系统的分析，绘制了各种统计分析的图片和曲线。以下是垂直90度的部分分析
3、结果（全部分析结果见南京大学C波段双偏振雷达_垂直90度观测_分析总结报告_第1部分_2014年8月5日_7日_21日补充.doc和南京大学C波段双偏振雷达_垂直90度观测_分析总结报告_第2部分_2014年9月22日.doc）雷达自身的RVP9信号处理器的输出，其IOP6、IOP9、IOP10的ZDR变化曲线如下：ZDR存在“小时”量级的变化此处ZDR变化，是因为时间上相隔了12个小时此处ZDR变化，是因为时间上相隔了8个小时9月份，我们利用数字中频接收机输出的IQ数据，重新对各个IOP的观测结果进行了计算。其IOP6、IOP9、IOP10的ZDR变化曲线如下：备注：从IQ数据重新计算的结果
4、，其ZDR的系统偏差和RVP9给出的有一点差别。这个原因是因为：基于IQ数据重新计算，其计算中的雷达常数、接收机灵敏度等值，是严格按照用仪表测定的参数设定的。而RVP9计算中的参数，其具体的值不清楚。从这些图中（包括RVP9给出的结果，和用IQ数据重新计算的结果）可以看出，存在两个主要的问题：l ZDR的偏差存在“小时”量级的波动的，波动范围约为0.30.4dB。l 第9个IOP，有11个VCP的ZDR系统偏差为0.3dB；第10个IOP，前7个VCP的ZDR的系统偏差为0.2dB。这些偏差结果和其它IOP下、其它VCP的结果有0.30.4dB的差异。由于在IQ数据计算出基数据的过程中，雷达常
5、数、接收机灵敏度等值是固定的。因此也就说明此时雷达的确发生了雷达常数、增益、灵敏度等硬件上的变化。最有可能跟该雷达双通道接收机的温度没有稳定有关。这个可以通过对IQ数据进行茅草的统计来分析是否是这个原因。本文就是对接收机茅草进行了分析，结果表明：l 从IOP1IOP11的分析结果可以发现，该雷达接收机两个通道的噪声功率和功率之差存在一定的起伏的。 噪声功率的起伏，会直接造成回波的dBZ也会有同样的起伏（由于在信号处理器的计算中，采用了仪表测定的噪声功率值）。 噪声功率之差的起伏，会直接造成ZDR的系统偏差出现同样的起伏。l 特别是对于IOP9、IOP10和IOP11，噪声功率起伏达到了12dB
6、，功率之差也有较大的起伏（0.30.5dB），这应该就是在垂直90度观测中发现ZDR的系统偏差存在波动的原因。l 我们认为，由于雷达硬件上存在的两个问题，共同造成了这个现象： 雷达发射机开机之后，整个舱室的温度会升高，会达到35度左右，甚至更高（由于当时雷达车上的空调有问题，无法达到很好的制冷效果）。 恒温箱在环境温度比较高的条件下，难以做到恒温，只能将内部的温度和周围的环境温度有一个差值，如果周围的环境温度高了几度，它内部也会高几度，无法实现恒定在25度。l 下一步工作： 把雷达车上的空调制冷修好 更换更强大制冷能力的恒温箱，并对温度进行监测记录 研究更高精度的接收机噪声的实时估计方法，并应
7、用到信号处理之中去 把分析的噪声功率起伏的数据代入信号处理的参数中，重新计算基数据同时，本文的分析结果也说明了雷达的硬件指标的重要性、精妙性，以及各种问题的防不胜防。这些经验教训，希望能用于今后布网雷达的双偏振升级改造之中。2 分析方法Matlab程序见“Read_NJU_CPol_Radar_Receiver_NoisePower_Analysis.m”。分析的方法是利用雷达发射的重复周期的远区部分的IQ数据，来计算噪声的功率。所谓远区的部分，是指在这个区域，所对应的海拔高度很高，不会有气象回波。我们在一个VCP扫描中，选取两个仰角层面的数据进行分析：l 当天线处于仰角为14度扫描时（即VC
8、P扫描的第12个仰角层面），距离是120150km的数据，相当于海拔高度是2936km。l 当天线处于垂直90度扫描时（即VCP扫描的第15个仰角层面），距离（即相当于海拔高度）是1620km的数据；也就是说，我们在每个径向上，选取最后的1/5的数据进行计算。但是，对于更低的低仰角而言，由于在远距离处，仍有可能会存在气象回波，因此对接收机噪声功率的统计就存在一定的困难。计算步骤为：1. 遍历所有的IOP，去掉RHI模式，去掉仰角层数没有达到15层的数据； 2. 找出仰角为第12层和第15层所对应的IQ文件3. 然后，打开对应的IQ数据文件，取最后的1/5 距离处，来进行茅草的统计，计算远处的噪
9、声的功率值（反映在数字中频接收机输入端的功率值），单位:dBm4. 在对全方位总共360根径向上的噪声功率值，进行平均，作为该体扫的噪声功率值。并求出方差，方差可以用于判断数据是否可信；5. 最后，将计算结果保存至txt文件中，并绘制各种曲线图。3 接收机茅草的分析结果3.1 对每个IQ数据文件的单独分析所有图片在“E:CJJWRP900AZK敏视达雷达对双通道接收机茅草的分析图片对每个IQ文件单独分析”目录下。这些图片，将每根径向上的噪声功率值绘制出来，以方便分析噪声功率随方位的变化。其中，NJU.20140531.003208.AR2.bz2所对应的第12和第15仰角的噪声随方位的变化图如
10、下：注意：上图中，所谓的“基数据中，水平通道噪声”和“基数据中，垂直通道噪声”是指RVP9记录在数据文件结构中的两个参数，这两个参数和真正的噪声功率并不一样，不知道RVP9为何如此记录。注意：上图中，可以发现接收机的噪声功率水平通道为-77.5dBm，垂直通道为-77.9dBm。但是，我们在2014年5月27日，对雷达用仪表测试得出的噪声功率为：l 水平通道的噪声功率=-75.96dBm（LNA端接噪声源，但不通电时测得）l 垂直通道的噪声功率=-76.40 dBm（LNA端接噪声源，但不通电时测得）这两者之间有约1.5dB的差异。但是这个差异是可以解释的：由于天线仰角比较高，波束对准了太空，
11、而太空的微波辐射是非常微弱的，因此其接收机输出的噪声功率自然要比用仪表测试时（此时接收机输入端接的是没有通电的噪声源，可以认为是一个室温下的黑体）要低。在观察生成的这些图片的时候，我们发现有时在第12个仰角上会出现噪声功率突然变大的现象，如下图：此处噪声强，是因为天线正好对准了太阳经过分析，认为这个噪声应该是此时天线正好对准了太阳，接收到了太阳的辐射能量所致。下面是计算出的此时此刻的太阳位置（雷达的位置是用Google Earth得到的）：北京时间可见，此时太阳位于仰角为14.35度、方位是287.32度，和上面图中所显示的方位是一致的。备注：国外有篇文章“Determining Weathe
12、r Radar Antenna Pointing Using Signals Detected from the Sun at low antenna elevations.pdf”，就是利用雷达在正常体扫工作中，雷达天线每次扫过太阳时的噪声功率和方位、仰角的关系，来对天线方位、仰角的定位精度进行在线的在线高精度监测。另外，我们还发现，有时还会出现噪声功率随方位变化，有较大起伏的现象，如下图：垂直90度的数据也有同样的现象：通过观察这个时刻的PPI的图，如下：从上面的PPI图可以发现，此时在雷达四周有很强的降水，而强大的降水雨滴自然也会发出微波辐射（就像地物一样），这些微波辐射被雷达的天线和高
13、灵敏度的接收机所收到，就会抬高了接收机的噪声了。从上面的噪声变化曲线中可以发现，这个抬高的程度竟然达到了1dB！今后在进行精确的CC、LDR、ZDR等的计算中，还必须要将其考虑进去。3.2 对每个IOP的综合分析所有图片在“E:CJJWRP900AZK敏视达雷达对双通道接收机茅草的分析图片对每个IOP综合分析”目录下。这些图片，将每个体扫下，接收机两个通道的噪声功率和功率之差，都绘制了出来。同时，还将噪声功率随方位变化的方差也绘制出来，可以用于判断数据是否可信，是否存在外界的太阳辐射，或者是降水雨滴的微波辐射。各个IOP的接收机噪声功率的曲线如下：3.2.1 IOP1天线没有达到指定的仰角就开
14、始观测记录了上图中，方差非常大的数据，是因为此时天线没有达到指定的仰角就开始观测记录了，自然其数据就完全不正常了（有一个专门的报告谈到这个雷达伺服的问题）。此时的PPI如下图：天线没有达到指定的仰角就开始观测记录了3.2.2 IOP2从上图可以看出，当天线位于第12层时，接收机噪声功率的方差为0.035dB左右。此时，噪声统计是对(150e3-120e3)/75*64=25600个采样点进行平均的结果（75是一个距离库的大小，64是积累数）。而当天线位于第15层（垂直90度）时，由于采样点比较少，因此接收机噪声功率的方差增加到0.065dB左右。此时，噪声统计是对(20e3-16e3)/75*
15、128=6800个采样点进行平均的结果，采样点降低到6800/256001/4，因此方差则增加到约2倍。3.2.3 IOP3从上图可以看出，接收机噪声功率随时间变化（从6月14日16日）有约0.5dB的起伏，水平和垂直的噪声功率之差有约0.2dB 的起伏。另外，还可以发现当天线位于第12层（仰角为14度）时，噪声功率比起天线位于第15层（仰角为90度），要高0.1dB。这也是很容易解释的：当天线位于仰角14度时，由于受到外界低层大气自身辐射的影响，噪声自然要高一点点。仰角升高，噪声慢慢降低下面是该雷达曾经测得的接收机噪声功率随仰角的变化图，可以发现，随着仰角的抬高，噪声会慢慢降低。这是太阳的辐射能量图 31 噪声随仰角的变化天线扫描正好经过了太阳注意：上图中，第12层的噪声的方差有时会达到0.3dB
