发明名称 一种船用雷达海杂波参数获取方法
摘要 本发明公开了一种船用雷达海杂波参数获取办法,通过对船用雷达回波信号进行采样量化、抽取或内插、快累计、去均值、累加和求平均值处理,来获取海杂波参数。本发明所涉及的海杂波参数获取办法,应用于船用雷达信号处理中的海杂波抑制,能够简单、快速的获取海杂波参数,避免了传统海杂波参数获取办法在实际应用中计算复杂和不能实时反映实际海域海杂波分布的缺点。
申请公布号 CN102540161B 申请公布日期 2013.09.04
申请号 CN201110392059.6 申请日期 2011.12.01
申请人 宁波成电泰克电子信息技术发展有限公司 发明人 田丹;李建;李浩;陈曦;陶吉怀;李宏波;冉元进;陈闹
分类号 G01S7/41(2006.01)I 主分类号 G01S7/41(2006.01)I
代理机构 宁波奥圣专利代理事务所(普通合伙) 33226 代理人 邱积权
主权项 1.一种船用雷达海杂波参数获取方法,其特征在于具体步骤如下:S1.定义如下各参数:自然数n<sub>0</sub>为当前雷达量程下单条雷达回波信号的采样点数,常数N为当前雷达量程下雷达显示屏上一条雷达扫描线上的点数,常数M为雷达扫描一周的扫描线条数,实数k为当前雷达量程下对采样数据进行抽取或内插处理的倍数,<img file="FDA00003383098900011.GIF" wi="225" he="134" />实数s为当前雷达量程的直径距离,实数c为电磁波在真空中的传播速度,自然数f为雷达回波信号的采样速率,实数数列<img file="FDA00003383098900012.GIF" wi="499" he="64" />为AD采样器对M条雷达回波信号进行采样得到的数据,实数数列b<sub>00</sub>,b<sub>01</sub>,b<sub>02</sub>……b<sub>(M-1)(N-1)</sub>为对<img file="FDA00003383098900013.GIF" wi="497" he="58" />进行抽取或内插处理后得到的数据,实数数列c<sub>00</sub>,c<sub>01</sub>,c<sub>02</sub>……c<sub>(M-1)(N-1)</sub>为对b<sub>00</sub>,b<sub>01</sub>,b<sub>02</sub>……b<sub>(M-1)(N-1)</sub>进行快累计处理后得到的数据,实数数列d<sub>00</sub>,d<sub>01</sub>,d<sub>02</sub>……d<sub>(M-1)(N-1)</sub>为对c<sub>00</sub>,c<sub>01</sub>,c<sub>02</sub>……c<sub>(M-1)(N-1)</sub>进行去均值处理得到的数据,实数数列c<sub>i0</sub>,c<sub>i1</sub>,c<sub>i2</sub>……c<sub>i(N-1)</sub>为第i条雷达扫描线回波信号经过了快累计处理后得到的数据,average<sub>i</sub>为数列c<sub>i0</sub>,c<sub>i1</sub>,c<sub>i2</sub>……c<sub>i(N-1)</sub>的平均值,实数数列d<sub>0j</sub>,d<sub>1j</sub>,d<sub>2j</sub>……d<sub>(M-1)j</sub>为M条雷达扫描线去均值处理后第二下标相同的数据,average<sub>j</sub>为数列d<sub>0j</sub>,d<sub>1j</sub>,d<sub>2j</sub>……d<sub>(M-1)j</sub>的平均值,实数数列para<sub>0</sub>,para<sub>1</sub>……para<sub>N-1</sub>为当前雷达量程下的海杂波参数值,实数Sum<sub>ij</sub>为对数列d<sub>0j</sub>,d<sub>1j</sub>,d<sub>2j</sub>……d<sub>(M-1)j</sub>中第一下标从0到i、第二下标为j的数的和,i、j分别为数列的两个下标,且0≤i&lt;M、0≤j&lt;N;S2.初始化i=0、j=0,Sum<sub>ij</sub>=0;S3.对第i条雷达回波信号进行采样,获得n<sub>0</sub>个采样值<img file="FDA00003383098900014.GIF" wi="438" he="58" />S4.对采样值<img file="FDA00003383098900015.GIF" wi="412" he="59" />进行k倍的内插或抽取处理,得到b<sub>i0</sub>,b<sub>i1</sub>,b<sub>i2</sub>……b<sub>i(N-1)</sub>,内插或抽取过程如下:当k&gt;1,对<img file="FDA00003383098900023.GIF" wi="411" he="63" />进行k倍的抽取,当k&lt;1时,对<img file="FDA00003383098900024.GIF" wi="413" he="64" />进行k倍的内插,当k=1时,<img file="FDA00003383098900025.GIF" wi="411" he="61" />保持不变;本步骤处理后得到的数据为b<sub>i0</sub>,b<sub>i1</sub>,b<sub>i2</sub>……b<sub>i(N-1)</sub>;S5.对b<sub>i0</sub>,b<sub>i1</sub>,b<sub>i2</sub>……b<sub>i(N-1)</sub>进行快累计处理,获得实数数列c<sub>i0</sub>,c<sub>i1</sub>,c<sub>i2</sub>……c<sub>i(N-1)</sub>,快累计处理过程如下:当j&lt;L时,c<sub>i0</sub>=b<sub>i0</sub>,c<sub>i1</sub>=b<sub>i1</sub>,……,c<sub>ij</sub>=b<sub>i(L-1)</sub>,当j≥L时,<img file="FDA00003383098900021.GIF" wi="320" he="209" />j=0,1,2,……767;S6.对c<sub>i0</sub>,c<sub>i1</sub>,c<sub>i2</sub>……c<sub>i(N-1)</sub>进行去均值处理,得到数列d<sub>i0</sub>,d<sub>i1</sub>,d<sub>i2</sub>……d<sub>i(N-1)</sub>,去均值处理过程如下:计算平均值<img file="FDA00003383098900026.GIF" wi="363" he="197" />当c<sub>ij</sub>≥average<sub>i</sub>,d<sub>ij</sub>=c<sub>ij</sub>-average<sub>i</sub>,当c<sub>ij</sub>&lt;average<sub>i</sub>,d<sub>ij</sub>=0;S7.第i条扫描线上的回波信号经过步骤S6处理后,进行一次累加处理,累加处理计算公式如下:Sum<sub>ij</sub>=Sum<sub>(i-1)j</sub>+d<sub>ij</sub>,j=0,1,2……(N-1);S8.如果i=M-1,计算平均值<img file="FDA00003383098900022.GIF" wi="393" he="139" />j=0,1,2……(N-1),否则回到步骤S3;S9.当前量程下的海杂波参数为para<sub>j</sub>=average<sub>j</sub>,j=0,1,2……(N-1)。
地址 315040 浙江省宁波市高新区扬帆路999弄5号7-4室