南京大学C波段双偏振雷达_用IQ数据重新计算基数据_特定时刻_2014年7月18日.doc

用IQ数据重新计算基数据（特定时刻）2014年7月18日1 概述2014年6、7月份，由北京敏视达公司制造的南京大学C波段双偏振雷达在安徽长丰站进行了持续的降水观测，录取了大量的基数据和IQ数据。下面，我们针对几个特定时刻的资料（根据“需要定做的数据时刻.rar”中的要求），用Matlab重新对IQ数据进行计算得出基数据，和RVP9信号处理器直接给出的基数据进行初步的对比分析。RVP9给出的基数据的结果保存在“H:QXLDNJU_CPolBaseData_Select需要定做的数据时刻_20140718”目录下。挑选出的IQ数据保存在“H:QXLDNJU_CPolIQData需要定做的数据时刻_20140718” 目录下。用Matlab重新处理之后的基数据保存在“H:QXLDNJU_CPolBaseData_DSP需要定做的数据时刻_20140718”目录下。2 IQ计算时的参数设置Matlab信号处理中的雷达参数，根据2014年5月27日，对雷达进行了仪表测试的结果而设定。测试详见南大C波段雷达标校情况_2014年5月30日.doc。l 水平通道的dBZ0=-28.31dBZl 垂直通道的dBZ0=-28.58dBZl 水平通道的噪声功率=-75.96dBm（LNA端接噪声源，但不通电时测得）l 垂直通道的噪声功率=-76.40 dBm（LNA端接噪声源，但不通电时测得）l DSP_Setup.ReceiverPhaseDiff=0; %两个接收通道的相位差信号处理算法设定的参数如下：l DSP_Setup.InterferenceFilterGate=1; %幅度相差多少倍，就认为是异步l DSP_Setup.ClutterFilterSel=2; %对杂波滤除的滤波器的凹口的大小设置l DSP_Setup.USE_R_Calc=true; %采用R0、R1进行计算l DSP_Setup.USE_dBZ_Calc=true; %采用dBZ计算ZDRl DSP_Setup.KDP_Average_Distance=2000; %定义KDP计算中，参与平均计算的距离范围l DSP_Setup.ThresholdLOG=3; %信噪比门限，单位:dB。l DSP_Setup.ThresholdCCOR=-25.0; %杂波对消的门限 注意：是负数相关系数和谱宽的计算，分别采用Multi-Lag1 和常规方法计算。采用Multi-Lag1的计算结果保存在“MultiLag1”子目录下，采用常规方法计算的结果，保存在“normal”子目录下。3 计算结果的初步对比（注意：这里只选取一个文件进行对比，其它的数据对比就由其他研究人员负责完成了。）3.1 RVP9给出的结果3.2 从IQ数据用matlab重新计算的结果3.2.1 Z、V、ZDR的结果3.2.2 用不同的方法计算CC CC采用标准方法的结果（不经过SNR的订正）计算的结果如下： CC采用标准方法的结果（经过SNR的订正）对这个结果，根据以下公式进行修正：SNR_Correct_CorrHV=sqrt( ( 1+1./ 10.( H_SNR_Data/10) .* ( 1+1./ 10.( (V_SNR_Data)/10) ); 注意：上面的修正结果中，SNR值就直接利用了Matlab 信号处理给出的SNR值。 CC采用Lag1的结果采用Lag1的结果比较可信，采用更多阶数计算的结果存在问题。3.2.3 用不同的方法计算谱宽 采用标准方法的结果 采用Lag1和2的结果4 其它需说明的问题注意：l NJU.20140531.170145.AR2、NJU.20140531.170515.AR2、NJU.20140531.171238.AR2所对应时刻的IQ数据没有保存下来；l NJU.20140630.222923.AR2所对应时刻的IQ数据没有保存下来；l 6月11日的IQ数据丢失严重（空调坏，磁盘阵列温度过高），丢失率约50%左右，因此计算结果仅供参考l NJU.20140611.103618.AR2 基数据丢失严重，无法解出数据；l NJU.20140611.114648.AR2的基数据所对应的IQ数据，以下几个IQ文件缺失：NJU_20140611_114719_02_RPH.IQNJU_20140611_115122_11_RPH.IQ NJU_20140611_115148_12_RPH.IQ NJU_20140611_115216_13_RPH.IQ NJU_20140611_115243_14_RPH.IQ NJU_20140611_1153