CN106019274A - 一种新型多普勒雷达成像装置及方法 - Google Patents
一种新型多普勒雷达成像装置及方法 Download PDFInfo
- Publication number
- CN106019274A CN106019274A CN201610349029.XA CN201610349029A CN106019274A CN 106019274 A CN106019274 A CN 106019274A CN 201610349029 A CN201610349029 A CN 201610349029A CN 106019274 A CN106019274 A CN 106019274A
- Authority
- CN
- China
- Prior art keywords
- signal
- dmd
- frequency
- target
- compressed sensing
- 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
- 238000000034 method Methods 0.000 title claims abstract description 12
- 238000003384 imaging method Methods 0.000 title abstract description 8
- 238000012545 processing Methods 0.000 claims abstract description 38
- 230000003287 optical effect Effects 0.000 claims abstract description 21
- 239000011159 matrix material Substances 0.000 claims description 46
- 229910052782 aluminium Inorganic materials 0.000 claims description 26
- 238000005070 sampling Methods 0.000 claims description 26
- XAGFODPZIPBFFR-UHFFFAOYSA-N aluminium Chemical compound [Al] XAGFODPZIPBFFR-UHFFFAOYSA-N 0.000 claims description 25
- 230000001427 coherent effect Effects 0.000 claims description 15
- 238000001228 spectrum Methods 0.000 claims description 13
- 230000010355 oscillation Effects 0.000 claims description 9
- 230000003321 amplification Effects 0.000 claims description 8
- 238000003199 nucleic acid amplification method Methods 0.000 claims description 8
- 238000005286 illumination Methods 0.000 claims description 7
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 230000011514 reflex Effects 0.000 claims description 3
- 230000015572 biosynthetic process Effects 0.000 claims description 2
- 238000004364 calculation method Methods 0.000 description 10
- 238000011084 recovery Methods 0.000 description 10
- 238000001514 detection method Methods 0.000 description 6
- 238000005259 measurement Methods 0.000 description 6
- 238000005516 engineering process Methods 0.000 description 5
- 230000005540 biological transmission Effects 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000008447 perception Effects 0.000 description 1
- 238000009738 saturating Methods 0.000 description 1
- 238000007493 shaping process Methods 0.000 description 1
- 230000003595 spectral effect Effects 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
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
一种新型多普勒雷达成像装置及方法，该装置包括发射装置，固定在发射装置出光端的目标，固定于目标反射光路上的第一汇聚透镜，固定于第一汇聚透镜汇聚光路上的数字微镜器件DMD，固定于数字微镜器件DMD反射光路上的第二汇聚透镜，固定于第二汇聚透镜汇聚光路上的接收装置；与接收装置连接的混频器，与混频器连接的中频放大器，与中频放大器连接的正交相位检波器，与正交相位检波器连接的信号处理装置；还包括与所述发射装置、混频器和正交相位检波器连接的高稳定信号源，与发射装置、数字微镜器件DMD和信号处理装置连接的同步装置；本发明还公开了新型多普勒雷达成像方法；本发明能够得到目标各个部位的旋转、振动信息，更加有利于目标的识别与跟踪。
Description
技术领域
本发明涉及多普勒测速领域，具体涉及一种新型多普勒雷达成像装置及方法。
背景技术
随着大规模集成电路和数字处理技术的发展，测速雷达广泛应用于机载预警、导航、武器制导、跟踪、战场侦察和武器火控等方面，成为重要的军事设备。装有测速雷达的预警飞机，已成为对付低空轰炸机和巡航导弹的有效军事设备。装有测速雷达的机载火控系统，可以对空中、地面和海上目标边搜索边跟踪，抗干扰性能好，可靠性好。
除了在军事上具有重要应用外，测速雷达也被广泛应用于气象探测、交通管制等和人民生活息息相关的领域。通过对气象回波进行速度分辨，可以获得不同高度大气层中各种空气湍流运动的分布情况，从而进行准确的天气预测。通过对车辆进行速度监测和提醒，可以减少事故的发生，保证出行安全。
目前多普勒测速系统一般都将目标看作是点目标，忽略了目标的不同部位相对于多普勒测速系统的速度变化。而实际情况中，目标不同部位的多普勒频移常常反映了目标的某些振动、旋转特征。这些振动、旋转特征对于目标的检测和识别相当有用。
传统的雷达成像系统对于接收端器件有着很高的要求，大面阵的探测器成为了影响传统雷达成像质量的一个关键因素。但大面阵的探测器往往面临着传感器阵列制造工艺和价格过于昂贵的限制，无法再有更进一步的突破。
压缩感知作为最近几年出现的信号处理技术，一经提出，就引起了学术界的广泛关注，并被美国科技评论评为2007年度十大科技进展。它可以在远小于Nyquist采样率的条件下，每次采样多重传递了若干样本的综合信息，然后通过非线性重建算法完美得重建信号。它巧妙地避开了对原始信号的直接采样，降低了对接收端器件的要求，设计重点由设计昂贵的接收端器件转化为设计新颖的信号恢复算法，极大地降低了雷达成像系统的成本，更加有利于雷达成像系统的小型化与轻型化。
发明内容
为了解决上述问题，本发明的目的在于提供一种新型多普勒雷达成像装置及方法，能够实时地获得目标场景中每一个分辨单元的速度信息，而又不存在测速模糊问题，既适用于单目标探测，又适用于多目标探测，既能得到目标整体的速度信息，又能得到目标不同部位的速度信息，从而得到目标各个部位的旋转、振动信息，更加有利于目标的识别与跟踪。利用单点探测器，有利于装置的小型化与廉价化，同时又能极大地节约传输空间和存储空间。
实现本发明目的的技术思路是：高稳定信号源为发射装置、混频器以及正交相位检波器提供载频一定的正弦信号。同步装置控制发射装置、数字微镜器件DMD以及信号处理装置按周期T同步工作。在每一个周期中，发射装置发射脉冲宽度为t、载频为f0的相参脉冲串信号照在目标上，目标返回来的信号经由第一汇聚透镜传给数字微镜器件DMD，数字微镜器件DMD对照在其上的入射光按照二值随机矩阵Φ的形式进行调制后传给第二汇聚透镜。第二汇聚透镜将数字微镜器件DMD调制过的光信号传给接收装置，接收装置将载频为f0+fdi(i＝1,2,3,…m，m为具有不同运动状态的物体的个数)的光信号转换为电信号，并作为主振信号传给混频器，其中fdi为反应第i个目标运动信息的多普勒频率。在混频器中，主振信号与载频为f0-fI的本振信号进行混频，其中fI为中频信号的频率，并通过选频回路选出频率为fdi+fI的信号输出到中频放大器。中频放大器对其进行功率放大后，进入正交相位检波器，将载频为fdi的信号检出来送入信号处理装置。在信号处理装置中对信号做FFT变换，得到信号的频谱图。信号的频谱图为频率fdi(i＝1,2,3,…m，m为具有不同运动状态的物体的个数)处的sinc函数。在M个周期后，会得到M个信号的频谱图。从中提取M个信号的幅度频谱,并选取M个幅度频谱图中频率为fdi处的幅值作为压缩感知中的测量向量的M个元素，测量矩阵选取调制数字微镜器件DMD使用的二值随机矩阵Φ，通过压缩感知恢复算法计算出目标场景中每一个分辨单元在频率为fdi时的幅值分量，并将所有非零元素提取出来按顺序放入一个行向量中。用同样的方法处理M个幅度频谱图中频率为fdi±n2π/T，n＝1,2,...,T/t处的幅值，得到2T/t+1个行向量。将第w个周期中(w＝1,2,…,M)频率为fdi±n2π/T，n＝1,2,...,T/t处的幅值作为压缩感知中测量向量的元素，将上面得到的2T/t+1个行向量构成矩阵作为压缩感知中的测量矩阵，通过压缩感知恢复算法计算出第w个周期中二值随机矩阵Φ在多普勒频率为fdi的物体处的分量Xwi(w＝1,2,…,M，i＝1,2,3,…m，m为具有不同运动状态的物体的个数)，计算∑i∑jfdiXwij，其中Xwij为Xwi中的第j个元素，将其作为测量向量Y中的第w个元素，测量矩阵选取调制数字微镜器件DMD使用的二值随机矩阵Φ，通过压缩感知恢复算法可以计算得到目标场景的速度图像A。
为达到以上目的，本发明采用如下技术方案：
一种新型多普勒雷达成像装置，包括发射装置1，固定在发射装置1出光端的目标10，固定于目标10反射光路上的第一汇聚透镜3，固定于第一汇聚透镜3汇聚光路上的数字微镜器件DMD2，固定于数字微镜器件DMD2反射光路上的第二汇聚透镜4，固定于第二汇聚透镜4汇聚光路上的接收装置5；与接收装置5连接的混频器6，与混频器6连接的中频放大器7，与中频放大器7连接的正交相位检波器8，与正交相位检波器8连接的信号处理装置9；还包括与所述发射装置1、混频器6和正交相位检波器8连接的高稳定信号源11，与所述发射装置1、数字微镜器件DMD2和信号处理装置9连接的同步装置12；
所述发射装置1，发射周期为T、脉冲宽度为t、初相为0的相参脉冲串光信号照在目标10上；
所述数字微镜器件DMD 2，依照压缩感知二值随机采样矩阵Φ的形式变换其上的铝镜状态，对目标10的反射光进行调制，压缩感知二值随机采样矩阵Φ的元素值即对应数字微镜器件DMD 2的铝镜状态，当矩阵元素值为+1时，数字微镜器件DMD 2对应位置的铝镜偏转+12°，将光信号反射到第二汇聚透镜4并传给接收装置5；当矩阵元素值为0时，数字微镜器件DMD 2对应位置的铝镜偏转-12°，将其上的入射光反射到吸收平面进行入射光吸收；
所述接收装置5，采用单点探测器，接收第二汇聚透镜4传来的光信号并转换为电信号以便后续处理；
所述混频器6，是一种典型的频谱频谱搬移电路，它将接收装置5传来的信号作为主振信号，混频器6的主振信号的频率为f0+fdi，与高稳定信号源11提供的频率为f0-fI的本振信号进行混频，采用选频回路选出频率为fdi+fI的中频信号；
所述中频放大器7，对混频器6传来的中频信号进行功率放大，以便进行后续处理；
所述正交相位检波器8，从中频放大器7传来的放大的中频信号中提取出多普勒频率fdi，得到的信号是频率为多普勒频率fdi的复信号；
所述信号处理装置9，对正交相位检波器8传来的复信号进行FFT变换，并运用压缩感知得到目标10的速度图像；
所述高稳定信号源11，为发射装置1提供载频为f0的信号，为混频器6提供频率为f0-fI的正弦本振信号，为正交相位检波器8提供频率为fI的相参基准信号；
所述同步装置12，保证发射装置1发射光信号、数字微镜器件DMD 2变换其上铝镜状态以及信号处理装置9对这一次采样得到的数据进行处理三者之间的同步进行。
所述正交相位检波器8包括第一乘法器30，与第一乘法器30连接的第一低通滤波器21；第二乘法器31，与第二乘法器31连接的第二低通滤波器22；与第一乘法器30和第二乘法器31连接的移相器20。
上述所述新型多普勒雷达成像装置的成像方法，包括如下步骤：
步骤一，利用信号处理装置9生成一个M×N维的压缩感知二值随机采样矩阵Φ和N×N维的正交稀疏表示矩阵Ψ，其中N＝p×q，p为速度图像的纵向分辨单元的个数，q为速度图像的横向分辨单元的个数，M为压缩感知的采样次数，也即相参脉冲串信号的周期数；
步骤二，信号处理装置9由压缩感知二值随机采样矩阵Φ生成M个测量子矩阵，并编号为1,2,3,…M；其中：M为压缩感知的采样次数；
步骤三，根据速度图像横向分辨单元的个数和速度图像纵向分辨单元的个数对数字微镜器件DMD2上的铝镜进行分块；
步骤四，发射装置1发射载频为f0、周期为T、脉冲宽度为t、初相为0的相参脉冲串光信号照在目标10上，设当前为第i个周期，并令i＝1；
步骤五，目标10的反射光信号经第一汇聚透镜3到达数字微镜器件DMD 2上，同步装置12控制数字微镜器件DMD 2按照编号为i的测量子矩阵设置其上的铝镜状态，完成对目标10的反射光的随机采样；
步骤六，数字微镜器件DMD2上的反射光经第二汇聚透镜4到达接收装置5，接收装置5将接收到的光信号转化为电信号，并作为主振信号传给混频器6；
步骤七，混频器6将主振信号与高稳定信号源11提供的频率为f0-fI的正弦本振信号进行混频，采用选频回路选出频率为fdi+fI的中频信号，并传给中频放大器7；
步骤八，中频放大器7对传来的含有多普勒频率fdi的中频信号进行功率放大，并传给正交相位检波器8；
步骤九，正交相位检波器8将中频放大器7传来的中频信号分别与两个相互正交的相参基准信号相乘后，通过第一低通滤波器21和第二低通滤波器22滤除高频分量后，得到频率为多普勒频率fdi的复信号，并传给信号处理装置9；
步骤十，信号处理装置9对第i个周期得到的复信号进行FFT变换，得到第i个幅度频谱图；
步骤十一，令周期序号数i＝i+1，并判断i是否小于压缩感知的采样次数M，如果i小于等于M，则转至步骤五，如果i大于M，则转至步骤十二；
步骤十二，信号处理装置9根据步骤十一得到的M个幅度频谱图，通过压缩感知计算得出目标10的速度图像。
与现有技术相比，本发明具有以下优点：
第一，能够实时准确地获得目标场景中每一个分辨单元的速度信息而又不存在测速模糊问题，既适用于单目标探测，又适用于多目标探测，更加适合战场等情况复杂多变的场合。
第二，相对于目前大多数的多普勒系统都将目标视为“点目标”，本发明既能得到目标整体的速度信息，又能得到目标不同部位的速度信息，从而得到目标各个部位的旋转、振动信息，更加有利于目标的识别与跟踪。例如直升机的螺旋桨的速度信息，可以作为识别直升机的一个很好的特征。
第三，利用单点探测器作为接收装置，不仅有利于突破传感器阵列制造工艺以及价格对获取高分辨率图像的限制，又有利于装置的小型化与廉价化，同时又能极大地节约传输空间和存储空间。
附图说明
图1为本发明系统的结构示意图。
图2为本发明系统中正交相位检波器的结构示意图。
图3为本发明方法的实现流程图。
具体实施方式
以下结合附图及具体实施例对本发明作进一步的详细描述。
如图1所示，本发明新型多普勒雷达成像装置，包括发射装置1，固定在发射装置1出光端的目标10，固定于目标10反射光路上的第一汇聚透镜3，固定于第一汇聚透镜3汇聚光路上的数字微镜器件DMD2，固定于数字微镜器件DMD2反射光路上的第二汇聚透镜4，固定于第二汇聚透镜4汇聚光路上的接收装置5；与接收装置5连接的混频器6，与混频器6连接的中频放大器7，与中频放大器7连接的正交相位检波器8，与正交相位检波器8连接的信号处理装置9；还包括与所述发射装置1、混频器6和正交相位检波器8连接的高稳定信号源11，与所述发射装置1、数字微镜器件DMD2和信号处理装置9连接的同步装置12。
如图2所示，正交相位检波器8，包括第一乘法器30，与第一乘法器30连接的第一低通滤波器21；第二乘法器31，与第二乘法器31连接的第二低通滤波器22；与第一乘法器30和第二乘法器31连接的移相器20。
本发明新型多普勒雷达成像装置的工作原理如下：
高稳定信号源11为发射装置1、混频器6以及正交相位检波器8提供载频一定的正弦信号。同步装置12控制发射装置1、数字微镜器件DMD2以及信号处理装置9按周期T同步工作。在每一个周期中，发射装置发射脉冲宽度为t、载频为f0的相参脉冲串信号照在目标10上，目标10返回来的信号经由第一汇聚透镜3传给数字微镜器件DMD2，数字微镜器件DMD2对照在其上的入射光按照二值随机矩阵Φ的形式进行调制后传给第二汇聚透镜4。第二汇聚透镜4将数字微镜器件DMD2调制过的光信号传给接收装置5，接收装置5将载频为f0+fdi(i＝1,2,3,…m，m为具有不同运动状态的物体的个数)的光信号转换为电信号，并作为主振信号传给混频器6，其中fdi为反应第i个目标运动信息的多普勒频率。在混频器6中，主振信号与载频为f0-fI的本振信号进行混频，其中fI为中频信号的频率，并通过选频回路选出频率为fdi+fI的信号输出到中频放大器7。中频放大器7对其进行功率放大后，进入正交相位检波器8，将载频为fdi的信号检出来送入信号处理装置9。在信号处理装置9中对信号做FFT变换，得到信号的频谱图。信号的频谱图为频率fdi(i＝1,2,3,…m，m为具有不同运动状态的物体的个数)处的sinc函数。在M个周期后，会得到M个信号的频谱图。从中提取M个信号的幅度频谱,并选取M个幅度频谱图中频率为fdi处的幅值作为压缩感知中的测量向量的M个元素，测量矩阵选取调制数字微镜器件DMD2使用的二值随机矩阵Φ，通过压缩感知恢复算法计算出目标场景中每一个分辨单元在频率为fdi时的幅值分量，并将所有非零元素提取出来按顺序放入一个行向量中。用同样的方法处理M个幅度频谱图中频率为fdi±n2π/T，n＝1,2,...,T/t处的幅值，得到2T/t+1个行向量。将第w个周期中(w＝1,2,…,M)频率为fdi±n2π/T，n＝1,2,...,T/t处的幅值作为压缩感知中测量向量的元素，将上面得到的2T/t+1个行向量构成矩阵作为压缩感知中的测量矩阵，通过压缩感知恢复算法计算出第w个周期中二值随机矩阵Φ在多普勒频率为fdi的物体处的分量Xwi(w＝1,2,…,M，i＝1,2,3,…m，m为具有不同运动状态的物体的个数)，计算∑i∑jfdiXwij，其中Xwij为Xwi中的第j个元素，将其作为测量向量Y中的第w个元素，测量矩阵选取调制数字微镜器件DMD使用的压缩感知二值随机采样矩阵Φ，通过压缩感知恢复算法可以计算得到目标场景的速度图像。
如图3所示，本发明新型多普勒雷达成像方法，包括如下步骤：
步骤一，利用信号处理装置9生成一个M×N维的压缩感知二值随机采样矩阵Φ和N×N维的正交稀疏表示矩阵Ψ，其中N＝p×q，p为速度图像的纵向分辨单元的个数，q为速度图像的横向分辨单元的个数，M为压缩感知的采样次数，也即相参脉冲串信号的周期数；
其中bist为Φ中第i行，第st列的元素，1≤s≤p,1≤t≤q,bist的值取+1或0；
步骤二，信号处理装置9由压缩感知二值随机采样矩阵Φ生成M个测量子矩阵Φi，并令编号i为1,2,3,…M，其中：M为压缩感知的采样次数；每个随机测量子矩阵Φi的形式如下所示：
第i个测量子矩阵Φi即对应原压缩感知二值随机采样矩阵Φ的第i行元素；
步骤三，根据速度图像的横向分辨单元个数p和速度图像的纵向分辨单元个数q对数字微镜器件DMD2上的铝镜进行分块；
设数字微镜器件DMD2上的铝镜个数为m×n，m为数字微镜器件中横向铝镜的个数，n为数字微镜器件中纵向铝镜的个数；将数字微镜器件DMD2上的铝镜分为p×q块，则每块有[min{m,n}/max{p,q}]2个铝镜。以数字微镜器件DMD2上铝镜个数为1024*768和图像分辨率为64*64为例，将数字微镜器件DMD2上的铝镜分为64*64块，每块有144个铝镜，对应于目标场景中的一个分辨单元；
步骤四，发射装置1发射载频为f0、周期为T、脉冲宽度为t、初相为0的相参脉冲串光信号照在目标10上，设当前为第i个周期，并令i＝1；
步骤五，目标10的反射光信号经第一汇聚透镜3到达数字微镜器件DMD 2上，同步装置12控制数字微镜器件DMD2按照编号为i的测量子矩阵设置其上的铝镜状态，完成对目标10的反射光的随机采样，当测量子矩阵相应元素值为+1时，将数字微镜器件DMD2对应位置的铝镜偏转+12°，将其上的入射光反射给第二汇聚透镜4并传给接收装置5；当随机矩阵相应元素值为0时，将数字微镜器件DMD2对应位置的铝镜偏转-12°，将其上的入射光反射到吸收平面进行入射光吸收；
步骤六，数字微镜器件DMD2上的反射光经第二汇聚透镜4到达接收装置5，接收装置5将接收到的光信号转化为电信号，并作为主振信号传给混频器6；
步骤七，混频器6将主振信号与高稳定信号源11提供的频率为f0-fI的本振信号进行混频，采用选频回路选出频率为fdi+fI的中频信号，并传给中频放大器7，具体如下：
接收装置5传来的主振信号可以表示为其中cij表示具有第i种运动状态的物体上的第j个分辨单元的反射率，表示具有第i种运动状态的物体上的第j个分辨单元的距离情况。高稳定信号源11提供的频率为f0-fI的本振信号可以表示为cos(2π(f0-fI)t)，两者混频后可得到如下信号：
选频回路的低通滤波器滤除掉频率为2f0+fdi-fI的信号后，保留频率为fdi+fI的中频信号，传给中频放大器7；
步骤八，中频放大器7对传来的含有多普勒频率fdi的中频信号进行功率放大，增益为K，并传给正交相位检波器8；
步骤九，正交相位检波器8将中频放大器7传来的中频信号分别与两个相互正交的相参基准信号相乘后，通过第一低通滤波器21和第二低通滤波器22滤除高频分量后，得到频率为多普勒频率fdi的复信号，并传给信号处理装置9，具体如下：
设相参基准信号sc(t)＝cos(2πfIt)，其经过移相器20，可得到与其正交的另一个相参基准信号sc'(t)＝cos(2πfIt+π/2)＝-sin(2πfIt)。将中频信号与相参基准信号sc(t)＝cos(2πfIt)相乘，可得信号为：
第一低通滤波器21滤除掉其中的高频分量，可得同相通道信号为
将中频信号与相参基准信号sc'(t)＝-sin(2πfIt)相乘，可得信号为：
第二低通滤波器22滤除掉其中的高频分量，可得正交通道信号为：
由此可获得正交双通道信号X(t)为：
步骤十，信号处理装置9对第i个周期得到的复信号进行FFT变换，得到第i个幅度频谱图；
步骤十一，令周期序号数i＝i+1，并判断i是否小于压缩感知的采样次数M，如果i小于等于M，则转至步骤五，如果i大于M，则转至步骤十二；
步骤十二，信号处理装置9根据步骤十一得到的M个幅度频谱图，通过压缩感知计算得出目标10的速度图像，具体如下。
12.1)选取M个幅度频谱图中频率为fdi处的幅值作为压缩感知中的测量子向量Y的M个元素，测量矩阵选取调制数字微镜器件DMD2使用的压缩感知二值随机采样矩阵Φ，计算压缩感知的全息算子A＝Φ Ψ，其中Ψ为N×N维的正交稀疏表示矩阵，根据公式Y＝AΘ＝ΦΨΘ，通过压缩感知恢复算法计算出Θ，Θ为N×1维的列向量，进而计算出ΨΘ，ΨΘ里面的每个元素为目标10中每一个分辨单元在频率为fdi时的幅值分量，将所有非零元素提取出来按顺序放入一个行向量中，行向量中非零元素的个数为l；
12.2)用同样的方法处理M个幅度频谱图中频率为fdi±n2π/T，n＝1,2,...,T/t处的幅值，得到2T/t+1个行向量。
12.3)将第w个周期中(w＝1,2,…,M)频率为fdi±n2π/T，n＝1,2,...,T/t处的幅值作为压缩感知测量向量y的2T/t+1个元素，将上面得到的2T/t+1个行向量构成矩阵作为压缩感知中的测量矩阵φ，生成l×l维的正交稀疏表示矩阵根据公式通过压缩感知恢复算法计算出θ，进而计算出里面的每个元素为第w个周期中压缩感知二值随机采样矩阵Φ在多普勒频率为fdi的物体处的分量Xwi(w＝1,2,…,M，i＝1,2,3,…m，m为具有不同运动状态的物体的个数)；
12.4)计算∑i∑jfdiXwij，其中Xwij为Xwi中的第j个元素，将其作为测量向量H中的第w个元素，测量矩阵选取调制数字微镜器件DM D使用的压缩感知二值随机采样矩阵Φ，计算压缩感知的全息算子A＝Φ Ψ，其中Ψ为N×N维的正交稀疏表示矩阵，根据公式Y＝Ah＝ΦΨh，通过压缩感知恢复算法计算出h，进而计算出目标10的速度图像C＝Ψh，C的形式如下：
C＝[C11 C12...C1q C21 C22...C2q....Cp1..Cpq]T
变换得到：
C中的每一点都代表目标10中相应分辨单元的速度信息。
Claims (3)
1.一种新型多普勒雷达成像装置，其特征在于：包括发射装置(1)，固定在发射装置(1)出光端的目标(10)，固定于目标(10)反射光路上的第一汇聚透镜(3)，固定于第一汇聚透镜(3)汇聚光路上的数字微镜器件DMD(2)，固定于数字微镜器件DMD(2)反射光路上的第二汇聚透镜(4)，固定于第二汇聚透镜(4)汇聚光路上的接收装置(5)；与接收装置(5)连接的混频器(6)，与混频器(6)连接的中频放大器(7)，与中频放大器(7)连接的正交相位检波器(8)，与正交相位检波器(8)连接的信号处理装置(9)；还包括与所述发射装置(1)、混频器(6)和正交相位检波器(8)连接的高稳定信号源(11)，与所述发射装置(1)、数字微镜器件DMD(2)和信号处理装置(9)连接的同步装置(12)；
所述发射装置(1)，发射周期为T、脉冲宽度为t、初相为0的相参脉冲串光信号照在目标(10)上；
所述数字微镜器件DMD(2)，依照压缩感知二值随机采样矩阵Φ的形式变换其上的铝镜状态，对目标(10)的反射光进行调制，压缩感知二值采样矩阵Φ的元素值即对应数字微镜器件DMD(2)的铝镜状态，当矩阵元素值为+1时，数字微镜器件DMD(2)对应位置的铝镜偏转+12°，将光信号反射到第二汇聚透镜4并传给接收装置(5)；当矩阵元素值为0时，数字微镜器件DMD(2)对应位置的铝镜偏转-12°，将其上的入射光反射到吸收平面进行入射光吸收；
所述接收装置(5)，采用单点探测器，接收第二汇聚透镜(4)传来的光信号并转换为电信号以便后续处理；
所述混频器(6)，是一种典型的频谱搬移电路，它将接收装置(5)传来的信号作为主振信号，混频器(6)的主振信号的频率为f0+fdi，与高稳定信号源(11)提供的频率为f0-fI的本振信号进行混频，采用选频回路选出频率为fdi+fI的中频信号；
所述中频放大器(7)，对混频器(6)传来的中频信号进行功率放大，以便进行后续处理；
所述正交相位检波器(8)，从中频放大器(7)传来的放大的中频信号中提取出多普勒频率fdi，得到的信号是频率为多普勒频率fdi的复信号；
所述信号处理装置(9)，对正交相位检波器(8)传来的复信号进行FFT变换，并运用压缩感知得到目标(10)的速度图像；
所述高稳定信号源(11)，为发射装置(1)提供载频为f0的信号，为混频器(6)提供频率为f0-fI的正弦本振信号，为正交相位检波器8提供频率为fI的相参基准信号；
所述同步装置(12)，保证发射装置(1)发射光信号、数字微镜器件DMD(2)变换其上铝镜状态以及信号处理装置(9)对这一次采样得到的数据进行处理三者之间的同步进行。
2.根据权利要求1所述的新型多普勒雷达成像装置，其特征在于：所述正交相位检波器(8)包括第一乘法器(30)，与第一乘法器(30)连接的第一低通滤波器(21)；第二乘法器(31)，与第二乘法器(31)连接的第二低通滤波器(22)；与第一乘法器(30)和第二乘法器(31)连接的移相器(20)。
3.权利要求2所述新型多普勒雷达成像装置的成像方法，其特征在于：包括如下步骤：
步骤一，利用信号处理装置(9)生成一个M×N维的压缩感知二值随机采样矩阵Φ和N×N维的正交稀疏表示矩阵Ψ，其中N＝p×q，p为速度图像的纵向分辨单元的个数，q为速度图像的横向分辨单元的个数，M为压缩感知的采样次数，也即相参脉冲串信号的周期数；
步骤二，信号处理装置(9)由压缩感知二值随机采样矩阵Φ生成M个测量子矩阵，并编号为1,2,3,…M；其中：M为压缩感知的采样次数；
步骤三，根据速度图像横向分辨单元的个数和速度图像纵向分辨单元的个数对数字微镜器件DMD(2)上的铝镜进行分块；
步骤四，发射装置(1)发射载频为f0、周期为T、脉冲宽度为t、初相为0的相参脉冲串光信号照在目标(10)上，设当前为第i个周期，并令i＝1；
步骤五，目标(10)的反射光信号经第一汇聚透镜(3)到达数字微镜器件DMD(2)上，同步装置(12)控制数字微镜器件DMD(2)按照编号为i的测量子矩阵设置其上的铝镜状态，完成对目标(10)的反射光的随机采样；
步骤六，数字微镜器件DMD(2)上的反射光经第二汇聚透镜(4)到达接收装置(5)，接收装置(5)将接收到的光信号转化为电信号，并作为主振信号传给混频器(6)；
步骤七，混频器(6)将主振信号与高稳定信号源(11)提供的频率为f0-fI的正弦本振信号进行混频，采用选频回路选出频率为fdi+fI的中频信号，并传给中频放大器(7)；
步骤八，中频放大器(7)对传来的含有多普勒频率fdi的中频信号进行功率放大，并传给正交相位检波器(8)；
步骤九，正交相位检波器(8)将中频放大器(7)传来的中频信号分别与两个相互正交的相参基准信号相乘后，通过第一低通滤波器(21)和第二低通滤波器(22)滤除高频分量后，得到频率为多普勒频率fdi的复信号，并传给信号处理装置(9)；
步骤十，信号处理装置(9)对第i个周期得到的复信号进行FFT变换，得到第i个幅度频谱图；
步骤十一，令周期序号数i＝i+1，并判断i是否小于压缩感知的采样次数M，如果i小于等于M，则转至步骤五，如果i大于M，则转至步骤十二；
步骤十二，信号处理装置(9)根据步骤十一得到的M个幅度频谱图，通过压缩感知计算得出目标(10)的速度图像。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610349029.XA CN106019274B (zh) | 2016-05-24 | 2016-05-24 | 一种新型多普勒雷达成像装置及方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610349029.XA CN106019274B (zh) | 2016-05-24 | 2016-05-24 | 一种新型多普勒雷达成像装置及方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106019274A true CN106019274A (zh) | 2016-10-12 |
| CN106019274B CN106019274B (zh) | 2019-01-22 |
Family
ID=57094052
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610349029.XA Active CN106019274B (zh) | 2016-05-24 | 2016-05-24 | 一种新型多普勒雷达成像装置及方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106019274B (zh) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106950554A (zh) * | 2017-03-21 | 2017-07-14 | 中国人民解放军国防科学技术大学 | 实现振动干扰条件下旋转目标参数估计的方法 |
| CN109655821A (zh) * | 2017-10-11 | 2019-04-19 | 西梅奥有限公司 | 用于确定目标的角方位、地点和/或速度的雷达方法和雷达系统 |
| CN113050118A (zh) * | 2021-05-31 | 2021-06-29 | 之江实验室 | 单臂激光外差关联成像雷达空间后端调制方法及装置 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5734347A (en) * | 1996-06-10 | 1998-03-31 | Mceligot; E. Lee | Digital holographic radar |
| CN202485671U (zh) * | 2010-11-26 | 2012-10-10 | 北京空间机电研究所 | 基于压缩感知理论的航天相机成像装置 |
| CN103234479A (zh) * | 2013-04-10 | 2013-08-07 | 浙江大学 | 基于压缩感知脉冲编码调制的三维信息获取装置和方法 |
| CN104142506A (zh) * | 2014-08-15 | 2014-11-12 | 中国科学院上海技术物理研究所 | 一种基于压缩感知的激光雷达成像系统 |
| CN105223582A (zh) * | 2015-09-01 | 2016-01-06 | 西安交通大学 | 一种基于压缩感知的激光雷达成像装置及成像方法 |
| CN105242281A (zh) * | 2015-09-01 | 2016-01-13 | 西安交通大学 | 基于apd阵列的三维激光成像系统及方法 |
-
2016
- 2016-05-24 CN CN201610349029.XA patent/CN106019274B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5734347A (en) * | 1996-06-10 | 1998-03-31 | Mceligot; E. Lee | Digital holographic radar |
| CN202485671U (zh) * | 2010-11-26 | 2012-10-10 | 北京空间机电研究所 | 基于压缩感知理论的航天相机成像装置 |
| CN103234479A (zh) * | 2013-04-10 | 2013-08-07 | 浙江大学 | 基于压缩感知脉冲编码调制的三维信息获取装置和方法 |
| CN104142506A (zh) * | 2014-08-15 | 2014-11-12 | 中国科学院上海技术物理研究所 | 一种基于压缩感知的激光雷达成像系统 |
| CN105223582A (zh) * | 2015-09-01 | 2016-01-06 | 西安交通大学 | 一种基于压缩感知的激光雷达成像装置及成像方法 |
| CN105242281A (zh) * | 2015-09-01 | 2016-01-13 | 西安交通大学 | 基于apd阵列的三维激光成像系统及方法 |
Non-Patent Citations (4)
| Title |
|---|
| 凡守涛: ""单脉冲气象雷达系统分析与设计"", 《中国博士学位论文全文数据库 信息科技辑》 * |
| 刘丽华等: ""连续波多普勒测速雷达射频前端电路设计与仿真"", 《电子科技》 * |
| 迪可新等: ""基于压缩感知的空间稀疏目标成像方法研究"", 《微电子学与计算机》 * |
| 马骏等: ""基于压缩传感理论的单像素成像系统设计"", 《红外技术》 * |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106950554A (zh) * | 2017-03-21 | 2017-07-14 | 中国人民解放军国防科学技术大学 | 实现振动干扰条件下旋转目标参数估计的方法 |
| CN106950554B (zh) * | 2017-03-21 | 2019-04-12 | 中国人民解放军国防科学技术大学 | 实现振动干扰条件下旋转目标参数估计的方法 |
| CN109655821A (zh) * | 2017-10-11 | 2019-04-19 | 西梅奥有限公司 | 用于确定目标的角方位、地点和/或速度的雷达方法和雷达系统 |
| CN109655821B (zh) * | 2017-10-11 | 2024-03-15 | 西梅奥有限公司 | 用于确定目标的矢量速度的雷达方法和雷达系统 |
| CN113050118A (zh) * | 2021-05-31 | 2021-06-29 | 之江实验室 | 单臂激光外差关联成像雷达空间后端调制方法及装置 |
| CN113050118B (zh) * | 2021-05-31 | 2021-08-13 | 之江实验室 | 单臂激光外差关联成像雷达空间后端调制方法及装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106019274B (zh) | 2019-01-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US7986397B1 (en) | FMCW 3-D LADAR imaging systems and methods with reduced Doppler sensitivity | |
| DK2866051T3 (en) | LASER DETECTION AND DISTANCE MEASURING DEVICE FOR DETECTING AN OBJECT UNDER A WATER SURFACE | |
| US6678395B2 (en) | Video search and rescue device | |
| CN107885223A (zh) | 基于激光的无人机回收引导系统 | |
| WO2016106957A1 (zh) | 一种超宽波段图谱关联探测装置与探测方法 | |
| CN109164430A (zh) | 利用激光回波与光斑探测目标位置与姿态的系统及方法 | |
| CN104535996B (zh) | 图像/激光测距/低空调频连续波雷达一体化系统 | |
| US6670920B1 (en) | System and method for single platform, synthetic aperture geo-location of emitters | |
| Fuchs et al. | Conceptual design of the CZMIL data acquisition system (DAS): integrating a new bathymetric lidar with a commercial spectrometer and metric camera for coastal mapping applications | |
| CN106019274A (zh) | 一种新型多普勒雷达成像装置及方法 | |
| CN106443702A (zh) | 一种瑞钠信标结合探测的自适应光学系统 | |
| CN108287350A (zh) | 基于多信息综合的天基空中目标探测关键参数确定方法 | |
| CN104142506A (zh) | 一种基于压缩感知的激光雷达成像系统 | |
| CN106066481B (zh) | 一种基于dmd的四维激光雷达成像装置及方法 | |
| CN107656288B (zh) | 一种新型剪切光束成像系统及目标图像获取方法 | |
| CN111190151A (zh) | 扫描模式下多模态小卫星sar的系统参数设计及发射功率优化方法 | |
| CN204086540U (zh) | 基于压缩感知的激光雷达成像系统 | |
| Zhang et al. | mmHawkeye: Passive UAV Detection with a COTS mmWave Radar | |
| RU2627550C1 (ru) | Стереоскопический когерентный доплеровский локатор | |
| CN104155658A (zh) | 一种基于压缩感知的激光雷达成像系统的图像重构方法 | |
| CN116453277A (zh) | 周界安防入侵事件全覆盖式精准检测方法及系统 | |
| Martin et al. | Design and performance of a distributed aperture millimeter-wave imaging system using optical upconversion | |
| CN106066480B (zh) | 一种基于多维特征的运动目标跟踪装置及方法 | |
| Bounds | The Infrared airborne radar sensor suite | |
| CN114740482A (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 | ||
| GR01 | Patent grant | ||
| CP03 | Change of name, title or address | ||
| CP03 | Change of name, title or address |
Address after: Room 820, 8th floor, Zhongchuang Plaza, No. 385 Aerospace Middle Road, Aerospace Base, Xi'an City, Shaanxi Province, 710100 Patentee after: Xi'an Shenrui Intelligent Technology Co.,Ltd. Country or region after: China Address before: 710199 Shaanxi Xi'an Aerospace base 385 space Road 8, 820 floor, Chong Chuang square. Patentee before: XI'AN SHENQIONG PHOTOELECTRIC TECHNOLOGY CO.,LTD. Country or region before: China |