CN104407329A - 一种用于区域检测的零副瓣二维脉冲压缩方法 - Google Patents
一种用于区域检测的零副瓣二维脉冲压缩方法 Download PDFInfo
- Publication number
- CN104407329A CN104407329A CN201410704941.3A CN201410704941A CN104407329A CN 104407329 A CN104407329 A CN 104407329A CN 201410704941 A CN201410704941 A CN 201410704941A CN 104407329 A CN104407329 A CN 104407329A
- Authority
- CN
- China
- Prior art keywords
- matrix
- row
- dimensional
- code
- echo data
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000007906 compression Methods 0.000 title claims abstract description 25
- 230000006835 compression Effects 0.000 title claims abstract description 25
- 238000000034 method Methods 0.000 title claims abstract description 24
- 238000001514 detection method Methods 0.000 title claims abstract description 20
- 238000001914 filtration Methods 0.000 claims abstract description 36
- 239000011159 matrix material Substances 0.000 claims abstract description 31
- 150000001875 compounds Chemical class 0.000 claims abstract description 4
- 230000005284 excitation Effects 0.000 claims description 10
- 238000004519 manufacturing process Methods 0.000 claims description 3
- 230000015572 biosynthetic process Effects 0.000 claims description 2
- 230000005540 biological transmission Effects 0.000 abstract description 2
- 241000702244 Orthoreovirus Species 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 238000003860 storage Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000005314 correlation function Methods 0.000 description 1
- 238000013144 data compression Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000003306 harvesting Methods 0.000 description 1
- 230000002045 lasting effect Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/28—Details of pulse systems
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/523—Details of pulse systems
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种用于区域检测的零副瓣二维脉冲压缩方法,设置二维编码对,首先找到两对一维的最佳二进序列偶(x1,y1),(x2,y2),使用复合构造法可得二维编码对,其中一组用于编码发射信号,另一组用于编码压缩滤波器;将用于编码发射信号的编码对组成编码发射信号矩阵;取所探测区域之间的回波数据使用B进行滤波进行滤波,从回波数据矩阵中取出探测区域对应的列形成一个矩阵P,将P的第1列至第n列取出,形成一个矩阵Pk1,Pk1与B进行滤波,滑动处理直至将P滤完,经过时间T之后,进行下一次滤波并输出结果,直至完成整个过程;本发明的有益效果是使得二维离散编码信号可以用在雷达的发射和接收中。
Description
技术领域
本发明属于数据压缩技术领域,涉及一种用于区域检测的零副瓣二维脉冲压缩方法。
背景技术
雷达是利用目标对电磁波的反射来进行检测和定位的。声纳是利用目标对机械波的反射来进行检测和定位的。它们通常要用脉冲压缩方法提高探测距离和分辨率。目前脉冲压缩存在较大的副瓣导致干扰相邻目标的检测。降低脉冲压缩副瓣能够提高雷达、声纳的性能。
在目前脉冲压缩的离散编码信号设计中,要求信号具有良好的自相关特性,这样的信号具有能将信号与自身延迟信号区分开来的特性,从而能够从叠加在一起的回波中区分两个相邻的目标。实际中使用的都是一维编码信号,性能最好的是巴克码。
本发明将二维编码使用在脉冲压缩中,用于雷达、声纳对指定区域的目标检测,改善检测性能。
发明内容
本发明的目的在于提供一种用于区域检测的零副瓣二维脉冲压缩方法,首次将二维离散编码信号运用在雷达脉冲压缩中,解决了目前二维离散编码信号无法用在雷达脉冲压缩中的问题。
本发明所采用的技术方案是按照以下步骤进行:
步骤1:设置二维编码对(A,B),其中一组用于编码发射信号,另一组用于编码压缩滤波器;
步骤2:将A用于编码发射信号,组成编码发射信号矩阵;将A复制成[AAA],去除首位两列得
S=[sij]z×(3n-2)=[AC2...ACn AC1 AC2...ACn AC1 AC2...ACn-1],将S按行发射,每一行发射后所得回波数据按行左端对齐组成一个回波数据矩阵;
步骤3:设所探测区域为[z1,z2],z1,z2分别为距离发射机的距离,取出[n·T+t1,n·T+t2]之间的回波数据进行滤波,其中
步骤4:使用B进行滤波;从回波数据矩阵中取出探测区域对应的w列形成一个m×w矩阵P,将P的第1列至第n列取出,形成一个m×n矩阵Pk1,Pk1与B进行滤波,得到F1:接着将P的第2列至第n+1列取出,形成一个m×n矩阵Pk2,Pk2与B进行滤波,得到F2:滑动处理直至将P滤完,经过时间T之后,进行下一次滤波并输出结果,直至完成整个过程;
进一步,所述步骤1中二维编码对(A,B)的产生方法如下:首先找到两对一维的最佳二进序列偶(x1,y1),(x2,y2),使用复合构造法可得:A=[a(i,j)]=[x1(i)·x2(j)],B=[b(i,j)]=[y1(i)·y2(j)];
进一步,所述A表示成列向量的组合,设为:A=[AC1 AC2...ACn],ACi=(a1i,a2i,...,ani)T;
进一步,使用码长较长的二维编码对,能在相同脉宽的情况下降低单个码长持续时间,从而提高分辨率。
本发明的有益效果是使得二维离散编码信号可以用在雷达的发射和接收中,而且此二维离散编码信号的主瓣高度是一维离散零副瓣编码信号的四倍。
附图说明
图1是(S0,B0)完全滤波后的主副瓣示意图;
图2是(S0,B0)限定滤波区域之后的主副瓣;
图3是雷达编码脉冲发射顺序示意图;
图4是发射信号序列及回波存储示意图;
图5是二维脉冲压缩滤波框图;
图6是区域近点滤波结果;
图7是区域远点滤波结果;
图8是区域中点滤波结果。
具体实施方式
下面结合附图和具体实施方式对本发明进行详细说明。
本发明一种用于区域检测的零副瓣二维脉冲压缩方法用于雷达和声纳,为了抑制脉冲压缩的副瓣,本发明首次构建二维脉冲压缩方法,并构造了二维脉冲压缩区域副瓣为零的二维编码对。
本发明的技术方案如下:
步骤1:设置二维编码对,其中一组用于编码发射信号,另一组用于编码压缩滤波器。此类编码对的互相关函数主瓣为16,主瓣周围有一长段副瓣为零的区域,如实施例1中的图1和图2,图1为(S0,B0)完全滤波后的主副瓣,图2为(S0,B0)限定滤波区域之后的主副瓣。
预先选取二维编码对,设所选二维编码对为(A,B),A=(aij)m×n,B=(bij)m×n。二维编码对(A,B)的产生方法如下:首先找到两对一维的最佳二进序列偶(x1,y1),(x2,y2),使用复合构造法可得:A=[a(i,j)]=[x1(i)·x2(j)],B=[b(i,j)]=[y1(i)·y2(j)];将A表示
成列向量的组合,设为:
A=[AC1 AC2...ACn],ACi=(a1i,a2i,...,ani)T;
步骤2:将A用于编码发射信号,组成编码发射信号矩阵;
将A复制成[AAA]:
[A A A]=[AC1 AC2...ACn AC1 AC2...ACn AC1 AC2...ACn]
去除首尾两列,得:
S=[sij]m×(3n-2)=[AC2...ACn AC1 AC2...ACn AC1 AC2...ACn-1]
将S按行发射:
1.从0时刻起发射第一行S1=[s11 s12...s1,3n-2],并在S1发射结束后开始接收和存储回波数据(数据存至RAM-1),接收和存储回波数据至T时刻止。(PRF=1/T为雷达的脉冲重复频率,一般警戒雷达在300Hz左右)。
T时刻即开始发送第二行S2(S1回波接收结束时刻一般为S2开始发送时刻或之前若干时间。S1发射结束时刻到S2开始发射时刻决定了雷达最大不模糊探测距离,通常,S2开始后,S1的回波因为距离远,衰减大,已经低于雷达接收机的灵敏度而收不到了);
2.从T时刻起发送第二行S2=[s21 s22...s2,3n-2],并在S2发射结束后开始接收和存储回波数据(数据存至RAM-2),接收和存储回波数据至2T时刻止。2T时刻起即开始发送第三行;
如是重复,一直将S的m行数据发送完毕。
接着,第m+1次,发送S1,所得回波数据更新存储器RAM-1。第m+2次,发送S2,所得回波数据更新存储器RAM-2,…,第m+m次,发送Sm,所得回波数据更新存储器RAM-m。如是循环。编码发射顺序如图3所示。发射信号序列及回波存储方式如图4所示。待接收完第m行回波数据,即可将这m行数据左端对齐,组成一个矩阵。
步骤3:使用B进行滤波并输出结果。滤波的起始时刻和范围由预先设定的探测区域决定,即想要探测某个区域是否有目标,就对那个区域的回波数据进行滤波。设所探测区域为[z1,z2],z1,z2分别为距离发射机的距离。则应取出[n·T+t1,n·T+t2]之间的回波数据进行滤波(这些数据左端对齐,形成一个矩阵P),其中
步骤4:用B进行滤波的过程:首先,回波数据有m行,将其左端对齐排列,可以形成一个m行的矩阵(此矩阵依情况会有很多列),取出探测区域对应的那些列,例如有w列(一般要求w大于n),形成一个m×w矩阵P,将P的第1列至第n列取出,形成一个m×n矩阵Pk1,Pk1与B进行滤波,得到F1:接着将P的第2列至第n+1列取出,形成一个m×n矩阵Pk2,Pk2与B进行滤波,得到F2:如是滑动处理,直至将P滤完。经过时间T之后,进行下一次滤波并输出结果,如是循环。二维脉冲压缩滤波如图5所示。
为了能充分利用零副瓣区域,使得探测区域的滤波能达到如图6-图8的效果,必须选择合理的码长和单码持续时间:设所探测区域为[z1,z2],称z1为区域近点,z2为区域远点,(z1+z2)/2为区域中点。z1处若存在目标,则滤波结果应如图6所示;z2处若存在目标,则滤波结果应如图7所示,(z1+z2)/2处若存在目标,则滤波结果应如图8所示。由这三幅图可以看出,滤波数据只有13个码位。从区域近点到区域远点移动了10个码位。若z1=15(公里),z2=16(公里),则这10个码位代表的距离即是所探测的区域:16-15=1公里。因此单个码位持续时间为:
通过调整单个码长的持续时间,可以调整分辨率;调整滤波的开始时间和结束时间,可以调整探测区域的定位;使用码长更长的二维编码对,可以在相同脉宽的情况下降低单个码长持续时间,从而提高分辨率。
本发明的优点是使得二维离散编码信号可以用在雷达的发射和接收中,而且此二维离散编码信号的主瓣高度是一维离散零副瓣编码信号的四倍。
下面列举具体实施例对本发明进行说明:
实施例1:
1.为充分利用零相关区,消除副瓣影响,本发明采取区域检测的办法,即预先设定探测范围。举例如下:假设预先需要探测15公里至16公里之间的区域,选取某一编码对A0和B0,先构成新的发射信号S0,其列数为34列,由警戒雷达的PRF为300Hz,得单个码长的持续时长可设为0.6666μs。
2.将A0复制成[A0 A0 A0]形式,并去除掉首尾两列,所得二维编码为:
3.在0时刻,将S0的第一行
[-1 -1 -1 -1 1 1 1 1 1 1 1 -1 -1 -1 -1 -1 1 1 1 1 1 1 1 -1 -1 -1 -1 -1 1 1 1 1 1 1]发送出去,在发射结束时刻接收回波数据并将回波数据按行存储至RAM-1,直至T0时刻。此时发射S0的下一行[-1 -1 -1 -1 1 1 1 1 1 1 1 -1 -1 -1 -1 -1 1 1 1 1 1 1 1 -1 -1 -1 -1 -1 1 1 1 1 1 1]。同样,经过时间T0,接收并存储这一行回波数据至RAM-2。如此进行,一直发射并存储完全部4行数据。
再经过T0,重新发射第一行,并将接收的数据存储于RAM-1,更新RAM-1数据。再经过T0,重新发射第二行,并将接收的数据存储于RAM-2。如是不断循环更新。
4.从第一次完整接收4行数据开始,即使用B0进行二维滤波处理:将4行回波数据左端对齐,形成一个m行的矩阵(由于左端已对齐,此时可以等效的视为每一行的数据都从0时刻开始)。由于要探测的是15-16公里的区域,因此从回波数据中取出从 时刻开始的数据(设为第k1列)(采样间隔与单个码长持续时间相同),顺着取12列,即取至(k1+11)列,得到一个4×12矩阵P1,将P1与B0滤波得F1:接着向右滑动一列取回波数据,即从(k1+1)列开始,取至(k1+12)列,得矩阵P2,将P2与B0滤波得F2:如是滑动处理,直至最后一个取出的回波矩阵的最后一列位置为 止,此时停止滤波,并输出滤波结果(相应的F1,F2,…)。待第一行回波数据更新后,即可进行下一次滤波;第二行回波数据更新后,进行再下一次滤波,如是不断进行循环。
以上所述仅是对本发明的较佳实施方式而已,并非对本发明作任何形式上的限制,凡是依据本发明的技术实质对以上实施方式所做的任何简单修改,等同变化与修饰,均属于本发明技术方案的范围内。
Claims (4)
1.一种用于区域检测的零副瓣二维脉冲压缩方法,其特征在于:按照以下步骤进行:
步骤1:设置二维编码对(A,B),其中一组用于编码发射信号,另一组用于编码压缩滤波器;
步骤2:将A用于编码发射信号,组成编码发射信号矩阵;将A复制成[A A A],去除首位两列得
S=[sij]m×(3n-2)=[AC2 … ACn AC1 AC2 … ACn AC1 AC2 … ACn-1],将S按行发射,每一行发射后所得回波数据按行左端对齐组成一个回波数据矩阵;
步骤3:设所探测区域为[z1,z2],z1,z2分别为距离发射机的距离,取出[n·T+t1,n·T+t2]之间的回波数据进行滤波,其中
τ为单个脉冲的持续时间;
步骤4:使用B进行滤波;从回波数据矩阵中取出探测区域对应的w列形成一个m×w矩阵P,将P的第1列至第n列取出,形成一个m×n矩阵Pk1,Pk1与B进行滤波,得到F1:接着将P的第2列至第n+1列取出,形成一个m×n矩阵Pk2,Pk2与B进行滤波,得到F2:滑动处理直至将P滤完,经过时间T之后,进行下一次滤波并输出结果,直至完成整个过程。
2.按照权利要求1所述一种用于区域检测的零副瓣二维脉冲压缩方法,其特征在于:所述步骤1中二维编码对(A,B)的产生方法如下:首先找到两对一维的最佳二进序列偶(x1,y1),(x2,y2),使用复合构造法得到:A=[a(i,j)]=[x1(i)·x2(j)],B=[b(i,j)]=[y1(i)·y2(j)]。
3.按照权利要求2所述一种用于区域检测的零副瓣二维脉冲压缩方法,其特征在于:所述A表示成列向量的组合,设为:
A=[AC1 AC2 ... ACn],ACi=(a1i,a2i,...,ani)T。
4.按照权利要求1所述一种用于区域检测的零副瓣二维脉冲压缩方法,其特征在于:使用码长较长的二维编码对,能在相同脉宽的情况下降低单个码长持续时间,从而提高分辨率。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410704941.3A CN104407329B (zh) | 2014-11-27 | 2014-11-27 | 一种用于区域检测的零副瓣二维脉冲压缩方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410704941.3A CN104407329B (zh) | 2014-11-27 | 2014-11-27 | 一种用于区域检测的零副瓣二维脉冲压缩方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104407329A true CN104407329A (zh) | 2015-03-11 |
CN104407329B CN104407329B (zh) | 2017-06-16 |
Family
ID=52644974
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410704941.3A Expired - Fee Related CN104407329B (zh) | 2014-11-27 | 2014-11-27 | 一种用于区域检测的零副瓣二维脉冲压缩方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104407329B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113687364A (zh) * | 2021-07-13 | 2021-11-23 | 中国船舶重工集团公司第七一五研究所 | 一种三维成像声纳图像假底的抑制方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1776449A (zh) * | 2005-11-30 | 2006-05-24 | 东南大学 | 构造低峰值旁瓣雷达脉冲压缩波形的方法 |
JP4630767B2 (ja) * | 2005-08-29 | 2011-02-09 | 株式会社東芝 | レーダ装置 |
-
2014
- 2014-11-27 CN CN201410704941.3A patent/CN104407329B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4630767B2 (ja) * | 2005-08-29 | 2011-02-09 | 株式会社東芝 | レーダ装置 |
CN1776449A (zh) * | 2005-11-30 | 2006-05-24 | 东南大学 | 构造低峰值旁瓣雷达脉冲压缩波形的方法 |
Non-Patent Citations (3)
Title |
---|
张兴敢: "逆合成孔径雷达成像及目标识别", 《中国博士学位论文全文数据库 信息科技辑》 * |
毛滔 等: "互补码在高频地波雷达中的应用研究", 《电波科学学报》 * |
龚小年 等: "CORDIC算法和FPGA在线性脉冲数字压缩系统中的应用", 《华东交通大学学报》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113687364A (zh) * | 2021-07-13 | 2021-11-23 | 中国船舶重工集团公司第七一五研究所 | 一种三维成像声纳图像假底的抑制方法 |
CN113687364B (zh) * | 2021-07-13 | 2024-03-19 | 中国船舶重工集团公司第七一五研究所 | 一种三维成像声纳图像假底的抑制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104407329B (zh) | 2017-06-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104160292B (zh) | 雷达装置以及雷达信号处理方法 | |
CN103109203B (zh) | 生成用于距离测量的信号的方法以及用于发射机与接收机之间的距离测量的方法和系统 | |
CN101470202B (zh) | 一种脉冲多普勒雷达系统及其信号处理方法 | |
CN103364770B (zh) | 基于矩阵填充的雷达目标检测系统及其检测方法 | |
CN102193087A (zh) | 发射装置及方法、接收装置及方法、物标探知装置及方法 | |
CN108802718B (zh) | 基于随机发射脉冲的时分mimo雷达空时相位解耦方法 | |
WO2005033728A3 (en) | Ultra-wideband radar system using sub-band coded pulses | |
CN108508438B (zh) | 一种基于单频时变阈值的一比特回波数据采集方法及系统 | |
CN106908773B (zh) | 连续脉冲雷达的信号时域循环编码和恢复方法 | |
EP2778712A8 (en) | Real aperture radar system for use on board a satellite and for maritime surveillance applications | |
US9310483B2 (en) | Device and method for transceiving ultrasonic wave | |
CN102998658A (zh) | 一种脉冲雷达的宽窄正交脉冲补盲方法 | |
CN103792517A (zh) | 一种连续脉冲雷达发射信号编码方法 | |
CN104360346B (zh) | 合成带宽雷达高分辨成像方法 | |
CN104407329A (zh) | 一种用于区域检测的零副瓣二维脉冲压缩方法 | |
CN113759340B (zh) | 回波信号处理方法及装置、激光雷达及存储介质 | |
CN104535970B (zh) | 基于最大值的频率步进雷达信号目标抽取方法 | |
CN106680809A (zh) | 一种穿墙雷达自聚焦压缩感知成像方法 | |
EP1494043A3 (en) | Short-range vehicular radar system | |
CN107329130B (zh) | 一种基于雷达通信一体化系统接收端的信号处理方法 | |
GB2041688A (en) | Surveillance radar systems | |
EP2018578B1 (de) | Verfahren zur ortung von objekten mittels eines abbildenden radarsystems und transponder für die ortung von objekten mittels derartiger radarsysteme | |
JP6641805B2 (ja) | ポラリメトリsar装置および観測方法 | |
CN109490882B (zh) | 同时用于sar成像及运动目标检测的多脉冲组合方法 | |
CN105548994B (zh) | 一种基于压缩感知的水下多亮点目标的检测方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170616 Termination date: 20181127 |
|
CF01 | Termination of patent right due to non-payment of annual fee |