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

用IQ数据重新计算基数据（特定时刻） 2014年7月18日 概述 2014年6、7月份，由北京敏视达公司制造的南京大学C波段双偏振雷达在安徽长丰站进行了持续的降水观测，录取了大量的基数据和IQ数据。 下面，我们针对几个特定时刻的资料（根据“需要定做的数据时刻.rar”中的要求），用Matlab重新对IQ数据进行计算得出基数据，和RVP9信号处理器直接给出的基数据进行初步的对比分析。 RVP9给出的基数据的结果保存在“H:\QXLD\NJU_CPol\BaseData_Select\需要定做的数据时刻目录下。 挑选出的IQ数据保存在“H:\QXLD\NJU_CPol\IQData\需要定做的数据时刻 目录下。 用Matlab重新处理之后的基数据保存在“H:\QXLD\NJU_CPol\BaseData_DSP\需要定做的数据时刻目录下。 IQ计算时的参数设置 Matlab信号处理中的雷达参数，根据2014年5月27日，对雷达进行了仪表测试的结果而设定。测试详见《南大C波段雷达标校情况_2014年5月30日.doc》。 水平通道的dBZ0=-28.31dBZ 垂直通道的dBZ0=-28.58dBZ 水平通道的噪声功率=-75.96dBm（LNA端接噪声源，但不通电时测得） 垂直通道的噪声功率=-76.40 dBm（LNA端接噪声源，但不通电时测得） DSP_Setup.ReceiverPhaseDiff=0; %两个接收通道的相位差 信号处理算法设定的参数如下： DSP_Setup.InterferenceFilterGate=1; %幅度相差多少倍，就认为是异步 DSP_Setup.ClutterFilterSel=2; %对杂波滤除的滤波器的凹口的大小设置 DSP_Setup.USE_R_Calc=true; %采用R0、R1进行计算 DSP_Setup.USE_dBZ_Calc=true; %采用dBZ计算ZDR DSP_Setup.KDP_Average_Distance=2000; %定义KDP计算中，参与平均计算的距离范围 DSP_Setup.ThresholdLOG=3; %信噪比门限，单位:dB。 DSP_Setup.ThresholdCCOR=-25.0; %杂波对消的门限 注意：是负数 相关系数和谱宽的计算，分别采用Multi-Lag1 和常规方法计算。采用Multi-Lag1的计算结果保存在“MultiLag1”子目录下，采用常规方法计算的结果，保存在“normal”子目录下。 计算结果的初步对比 （注意：这里只选取一个文件进行对比，其它的数据对比就由其他研究人员负责完成了。） RVP9给出的结果 从IQ数据用matlab重新计算的结果 Z、V、ZDR的结果 用不同的方法计算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的结果比较可信，采用更多阶数计算的结果存在问题。 用不同的方法计算谱宽 采用标准方法的结果 采用Lag1和2的结果 其它需说明的问题 注意： NJU170145.AR2、NJU170515.AR2、NJU171238.AR2所对应时刻的IQ数据没有保存下来； NJU222923.AR2所对应时刻的IQ数据没有保存下来； 6月11日的IQ数据丢失严重（空调坏，磁盘阵列温度过高），丢失率约50%左右，因此计算结果仅供参考 NJU103618.AR2 基数据丢失严重，无法解出数据； NJU114648.AR2的基数据所对应的IQ数据，以下几个IQ文件缺失：NJU114719_02_RPH.IQNJU115122_11_RPH.IQ NJU115148_12_RPH.IQ NJU115216_13_RPH.IQ NJU115243_14_RPH.IQ NJU115314_15_RPH.IQ 另外，7月11日、12日 的IQ数据还在雷达车上，没有拷贝过来。下次再分析。 下