CN103389492B - 多通道随机谐波调制采样雷达接收机及其方法 - Google Patents
多通道随机谐波调制采样雷达接收机及其方法 Download PDFInfo
- Publication number
- CN103389492B CN103389492B CN201310316753.9A CN201310316753A CN103389492B CN 103389492 B CN103389492 B CN 103389492B CN 201310316753 A CN201310316753 A CN 201310316753A CN 103389492 B CN103389492 B CN 103389492B
- Authority
- CN
- China
- Prior art keywords
- radar
- target information
- signal
- frequency
- sampling
- 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.)
- Expired - Fee Related
Links
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种多通道随机谐波调制采样雷达接收机及其方法,包括一个混频器、M个模拟乘法器、M个积分器、M个低速AD转换器和一个目标信息重构处理模块,所述的混频器与M个模拟乘法器相连,每个模拟乘法器与一个积分器相连,每个积分器与一个低速AD转换器相连,AD转换器采样得到的数据输入到目标信息重构处理模块,目标信息重构处理模块通过构造观测矩阵、构造雷达稀疏表示基矩阵和重构雷达目标信息三个步骤完成雷达目标信息估计的目的。解决了宽带以超宽带雷达***数字化采样难题;在保证了信号重构效果的同时,有效地降低了调制信号硬件实现的难度与代价。
Description
技术领域
本发明属于雷达技术领域,涉及一种多通道随机谐波调制采样雷达接收机及其方法。本发明可省去传统雷达接收机的匹配滤波处理,简化信号处理流程,并且解决了现有宽带及超宽带雷达信号采样速率过高和信号处理数据量过大的问题。
背景技术
宽带雷达与传统窄带雷达相比具有高距离分辨率、目标识别和成像能力、极强的四抗能力等优势,对宽带雷达***,无论是在军用领域还是在民用领域,提高其距离分辨率一直是雷达技术研究中一个重要方向。但是,对于宽带信号,以Nyquist采样定理为指导的信号采样方法给现有的模拟/数字转换器(ADC)带来了极大的挑战。特别是对于宽带和超宽带雷达回波信号,传统的雷达接收机以Nyquist速率采样的实现成本以及难度都很高。
由于超宽带雷达的信号带宽很宽,其接收机相比于传统的雷达接收机有所区别。与传统的雷达接收机相比,超宽带雷达的接收机结构更加复杂,实现中面临更多的难题。超宽带雷达接收机以数字接收机为主要的研究方向,在该处理机中,输入被下变频成中频或基带信号,然后被高速ADC数字化,再通过数字信号处理器进行信号处理。
为了降低ADC的速度,在目前的超宽带雷达接收机中,常采用的方式是把信号的相关运算即匹配滤波放在A/D转换前执行,这时信号的相关运算是采用模拟形式进行的,而模拟器件的可变性较差,而且容易受到电路失配以及其他的非理想条件的干扰影响。因此,为了达到高性能的超宽带雷达信号接收,需要根据Nyquist采样定理对信号进行采样,由于单个芯片无法达到这样的高采样率和高动态范围,因此多个ADC并行工作是超宽带雷达接收机通常采用的方式。多个ADC构成的接收***中,每个ADC以部分的采样率进行采样。多个ADC接收***可以分为两类:时域交织接收机、频域信道化接收机。对于时域交织接收机,多个ADC采样前要对其输入信号加入准确的时延。每个ADC的输入是占据整个信号频带的信号,因此,每个ADC的采样保持器的速度仍然为原始信号的采样频率,给设计带来了难度;对于频域信道化接收机,各个通道必需是相参的,要求其频率特性一致,并且由于相邻通道滤波器的过渡带不理想会产生通带间干扰,这样的超宽带频域信道化接收机的硬件***结构复杂,实现难度大,成本高。
发明内容
本发明的目的在于克服现有技术的不足,提出了一种多通道随机谐波调制采样雷达接收机,本发明省去传统雷达接收机的匹配滤波处理,简化信号处理流程,并且解决了现有宽带及超宽带雷达信号采样速率过高和信号处理数据量过大的问题。
一种多通道随机谐波调制采样雷达接收机,包括一个混频器、M个模拟乘法器、M个积分器、M个低速AD转换器和一个目标信息重构处理模块,所述的混频器与M个模拟乘法器相连,每个模拟乘法器与一个积分器相连,每个积分器与一个低速AD转换器相连,AD转换器采样得到的数据输入到目标信息重构处理模块,目标信息重构处理模块通过构造观测矩阵、构造雷达稀疏表示基矩阵和重构雷达目标信息三个步骤完成雷达目标信息估计的目的。
所述的多通道随机谐波调制采样雷达接收机,混频器将雷达回波信号与信号载频fc混频,得到带宽为B的雷达基频回波信号。
所述的多通道随机谐波调制采样雷达接收机,模拟乘法器将雷达基频回波信号与随机谐波调制信号相乘,各通道随机谐波调制信号形式如下:
其中,m=1,2…,M,为第m个通道,M为接收机通道个数,km为随机产生的大小在1~2M之间的互不相等的整数,fe为奈奎斯特采样频率,N为本次处理雷达回波信号以fe速率采样所得到的样本个数。
所述的多通道随机谐波调制采样雷达接收机,积分器的积分时间为,其中,fe为奈奎斯特采样频率,N为本次处理雷达回波信号以fe速率采样所得到的样本个数。
所述的多通道随机谐波调制采样雷达接收机, AD转换器采样时钟频率为,其中,fe为奈奎斯特采样频率,N为本次处理雷达回波以fe速率采样所得到的样本个数;得到采样数据具体形式如下:
y=[y1,y2,…,yM]T
其中,ym,m=1,2,…,M为第m个通道的采样数据,M为接收机通道个数。
所述的雷达接收机的目标信息估计方法,包括以下步骤:
A1、混频器将雷达回波信号sR(t)与信号载频fc混频,得到带宽为B的雷达基频回波信号sB(t);
A2、模拟乘法器将雷达基频回波信号sB(t)与随机谐波调制信号相乘,各通道随机谐波调制信号形式如下:
其中,j表示复数的虚部,m=1,2…,M,为第m个通道,M为接收机通道个数,km为随机产生的大小在1~2M之间的互不相等的整数,fe为奈奎斯特采样频率,N为本次处理雷达回波信号以fe速率采样所得到的样本个数;
A3、积分器的积分时间为,其中, fe为奈奎斯特采样频率,N为本次处理雷达回波信号以fe速率采样所得到的样本个数;y=[y1,y2,…,yM]T
A4、AD转换器采样时钟频率为,其中,fe为奈奎斯特采样频率,N为本次处理雷达回波以fe速率采样所得到的样本个数;得到采样数据具体形式如下:
y=[y1,y2,…,yM]T
其中,ym,1,2,…,M为第m个通道的采样数据,M为接收机通道个数;
A5、AD转换器采样得到的数据输入到目标信息重构处理模块,目标信息重构处理模块通过构造观测矩阵、构造雷达稀疏表示基矩阵和重构雷达目标信息三个步骤完成雷达目标信息估计的目的。
所述的目标信息估计方法,目标信息重构处理模块的具体方法为:
1)构造观测矩阵
随机谐波调制信号等效的观测矩阵Φ形式如下:
其中, M为接收机通道个数,k1,k2…,kM为随机产生的大小在1~2M之间的互不相等的整数,fe为等效的采样频率,Ts为等效的采样间隔,N为本次处理雷达回波以fe速率采样所得到的样本个数;
2)构造雷达稀疏表示基矩阵
雷达稀疏基矩阵Ψ的形式如下:
其中s0(t)=rect(t/T0)exp(jπK0t2)为雷达发射信号的包络,t表示时间,c表示光速, rn为待估计距离的第n个距离单元离雷达的距离,在这样的稀疏基矩阵下,基频回波信号表示为:
sB(t)=Ψα+nR(t)
式中,sB(t)为雷达回波基频信号,nR(t)是回波中包含的噪声信号,α为雷达目标稀疏系数;
3)重构雷达目标信息
根据所构造的观测矩阵Φ和雷达稀疏表示基矩阵Ψ,采样数据y和雷达目标信息α的关系如下:
α中包含了雷达目标的距离和多普勒信息,因此,将雷达目标信息提取问题转化为如下的压缩感知理论中的稀疏重构问题:
其中,表示1-范数。
本发明与现有的技术相比具有以下优点:
1)本发明提出的一种多通道随机谐波调制采样雷达接收机可以大大降低对雷达回波信号的采样速率,解决了宽带以超宽带雷达***数字化采样难题。
2)本发明提出的一种多通道随机谐波调制采样雷达接收机省略了传统雷达信号处理中的匹配滤波处理过程,简化信号处理流程。
3)本发明提出的一种多通道随机谐波调制采样雷达接收机采用随机谐波信号调制,在保证了信号重构效果的同时,有效地降低了调制信号硬件实现的难度与代价。
附图说明
图1是本发明的接收机结构图;
图2是本发明的目标信息重构处理流程图;
图3是本发明的目标信息估计仿真图;
具体实施方式
以下结合具体实施例,对本发明进行详细说明。
本发明一种多通道随机谐波调制采样雷达接收机,包括一个混频器、M个模拟乘法器、M个积分器、M个低速AD转换器和一个目标信息重构处理模块,所述的混频器与M个模拟乘法器相连,每个模拟乘法器与一个积分器相连,每个积分器与一个低速AD转换器相连,AD转换器采样得到的数据输入到目标信息重构处理模块,重构得到目标信息估计。本实施例中取M为204。
混频器将雷达回波信号sR(t)与信号载频fc混频,得到带宽为B的雷达基频回波信号sB(t)。本实施例中取fc为3GHz,B为15MHz。
模拟乘法器将雷达基频回波信号sB(t)与随机谐波调制信号相乘,各通道随机谐波调制信号形式如下:
其中,j表示复数的虚部,m=1,2…,M,为第m个通道,M为接收机通道个数,km为随机产生的大小在1~2M之间的互不相等的整数,fe为奈奎斯特采样频率,N为本次处理雷达回波以fe速率采样所得到的样本个数。本实施例中取fe为30MHz,N为1024。
积分器的积分时间为,其中,fe为奈奎斯特采样频率,N为本次处理雷达回波以fe速率采样所得到的样本个数。积分器的积分时间也可以理解为本次处理雷达回波的时间。
AD转换器采样时钟频率为,其中,fe为奈奎斯特采样频率,N为本次处理雷达回波以fe速率采样所得到的样本个数。得到采样数据具体形式如下:
y=[y1,y2,…,yM]T
其中,ym,m=1,2,…,M为第m个通道的采样数据,M为接收机通道个数。本实施例中采样频率降低了1024倍。
结合附图2对目标信息重构处理模块进行进一步说明:
1)构造观测矩阵
随机谐波调制信号等效的观测矩阵Φ形式如下:
其中, M为接收机通道个数,k1,k2…,kM为随机产生的大小在1~2M之间的互不相等的整数,fe为等效的采样频率,Ts为等效的采样间隔,N为本次处理雷达回波以fe速率采样所得到的样本个数。本实施例中Ts为33ns。
2)构造雷达稀疏表示基矩阵
本发明中发射信号选用线性调频信号(LFM)。
雷达稀疏基矩阵Ψ的形式如下:
其中s0(t)=rect(t/T0)exp(jπK0t2)为雷达发射信号的包络,t表示时间。c表示光速, rn为待估计距离的第n个距离单元离雷达的距离。在这样的稀疏基矩阵下,基频回波信号可以表示为:
sB(t)=Ψα+nB(t)
式中,sB(t)为雷达回波基频信号,nR(t)是回波中包含的噪声信号,α为雷达目标稀疏系数。本实施例中T0为5μs,K0为3×1012Hz/s,nR(t)为高斯噪声,r1为4005m,rN为9120m,rn-rn-1为5m,其中n=2,3…,N。
3)重构雷达目标信息
根据所构造的观测矩阵Φ和雷达稀疏表示基矩阵Ψ,采样数据y和雷达目标信息α的关系如下:
α中包含了雷达目标的距离和多普勒信息,因此,将雷达目标信息提取问题转化为如下的压缩感知理论中的稀疏重构问题:
其中,表示1-范数。本实施例中采用的稀疏重构算法为稀疏贝叶斯学习(SBL)算法。
本发明的效果通过以下仿真试验进一步说明:
1、仿真条件
雷达参数如下:雷达发射线性调频信号,载频f0=3GHz,脉冲重复频率PRF=10KHz,脉冲宽度T=5μs,带宽B=15MHz,Nyquist采样率30MHz,信噪比SNR定义为后向散射系数最小处的散射点信噪比。
图3中,3个静止目标分别位于[6000,6500,7000]m,后向散射系数分别为[1,0.8,0.6],信噪比为10dB。
2. 仿真内容
验证本发明目标信息估计的正确性。
3. 仿真结果分析
从图3可以看出,本发明的雷达接收机,在采样率为30KHz,N/M=5的情况下,能够有效的重构出目标稀疏系数,证明了本发明的有效性和正确性。
应当理解的是,对本领域普通技术人员来说,可以根据上述说明加以改进或变换,而所有这些改进和变换都应属于本发明所附权利要求的保护范围。
Claims (5)
1.一种多通道随机谐波调制采样雷达接收机,其特征在于,包括一个混频器、M个模拟乘法器、M个积分器、M个低速AD转换器和一个目标信息重构处理模块,所述的混频器与M个模拟乘法器相连,每个模拟乘法器与一个积分器相连,每个积分器与一个低速AD转换器相连,AD转换器采样得到的数据输入到目标信息重构处理模块,目标信息重构处理模块通过构造观测矩阵、构造雷达稀疏表示基矩阵和重构雷达目标信息三个步骤完成雷达目标信息估计的目的,模拟乘法器将雷达基频回波信号与随机谐波调制信号相乘,各通道随机谐波调制信号形式如下:
其中,m=1,2...,M,为第m个通道,M为接收机通道个数,km为随机产生的大小在1~2M之间的互不相等的整数,fe为奈奎斯特采样频率,N为本次处理雷达回波信号以fe速率采样所得到的样本个数。
2.根据权利要求1所述的多通道随机谐波调制采样雷达接收机,其特征在于,混频器将雷达回波信号与信号载频fc混频,得到带宽为B的雷达基频回波信号。
3.根据权利要求1所述的多通道随机谐波调制采样雷达接收机,其特征在于,积分器的积分时间为其中,fe为奈奎斯特采样频率,N为本次处理雷达回波信号以fe速率采样所得到的样本个数。
4.根据权利要求1所述的多通道随机谐波调制采样雷达接收机,其特征在于,AD转换器采样时钟频率为其中,fe为奈奎斯特采样频率,N为本次处理雷达回波以fe速率采样所得到的样本个数;得到采样数据具体形式如下:
y=[y1,y2,...,yM]T
其中,ym,m=1,2,...,M为第m个通道的采样数据,M为接收机通道个数。
5.根据权利要求1-4任一所述的雷达接收机的目标信息估计方法,其特征在于,包括以下步骤:
A1、混频器将雷达回波信号sR(t)与信号载频fc混频,得到带宽为B的雷达基频回波信号sB(t);
A2、模拟乘法器将雷达基频回波信号sB(t)与随机谐波调制信号相乘,各通道随机谐波调制信号形式如下:
其中,j表示复数的虚部,m=1,2...,M,为第m个通道,M为接收机通道个数,km为随机产生的大小在1~2M之间的互不相等的整数,fe为奈奎斯特采样频率,N为本次处理雷达回波信号以fe速率采样所得到的样本个数;
A3、积分器的积分时间为其中,fe为奈奎斯特采样频率,N为本次处理雷达回波信号以fe速率采样所得到的样本个数;
A4、AD转换器采样时钟频率为其中,fe为奈奎斯特采样频率,N为本次处理雷达回波以fe速率采样所得到的样本个数;得到采样数据具体形式如下:
y=[y1,y2,...,yM]T
其中,ym,m=1,2,...,M为第m个通道的采样数据,M为接收机通道个数;
A5、AD转换器采样得到的数据输入到目标信息重构处理模块,目标信息重构处理模块通过构造观测矩阵、构造雷达稀疏表示基矩阵和重构雷达目标信息三个步骤完成雷达目标信息估计的目的。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310316753.9A CN103389492B (zh) | 2013-07-25 | 2013-07-25 | 多通道随机谐波调制采样雷达接收机及其方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310316753.9A CN103389492B (zh) | 2013-07-25 | 2013-07-25 | 多通道随机谐波调制采样雷达接收机及其方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103389492A CN103389492A (zh) | 2013-11-13 |
CN103389492B true CN103389492B (zh) | 2015-09-02 |
Family
ID=49533816
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310316753.9A Expired - Fee Related CN103389492B (zh) | 2013-07-25 | 2013-07-25 | 多通道随机谐波调制采样雷达接收机及其方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103389492B (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103763233B (zh) * | 2014-01-08 | 2017-02-08 | 西安电子科技大学 | 基于数字调制解调的gsm外辐射雷达信号处理方法 |
CN104020450B (zh) * | 2014-03-26 | 2019-01-08 | 北京跟踪与通信技术研究所 | 一种基于邻波相关和相位分层的雷达回波重构方法 |
CN103986478B (zh) * | 2014-05-13 | 2017-06-09 | 天津大学 | 一种适用于微网谐波监测的压缩感知重构方法 |
CN105471525B (zh) * | 2015-11-12 | 2019-03-08 | 中国电子科技集团公司第四十一研究所 | 一种矢量网络分析仪的四通道压缩感知数字接收机信号处理方法 |
TWI676043B (zh) * | 2018-11-08 | 2019-11-01 | 立積電子股份有限公司 | 超寬頻雷達收發器及其運作方法 |
CN110058219B (zh) * | 2018-12-30 | 2021-11-09 | 上海微波技术研究所(中国电子科技集团公司第五十研究所) | Mimo谐波雷达 |
CN111338252A (zh) * | 2020-03-09 | 2020-06-26 | 中国电子科技集团公司第十四研究所 | 一种多通道兼容雷达信号处理的fpga架构设计 |
CN111355493B (zh) * | 2020-04-03 | 2023-05-23 | 哈尔滨工业大学 | 一种面向调制宽带转换器的支撑集筛选重构方法 |
CN112014639B (zh) * | 2020-09-02 | 2022-07-05 | 安徽一天电能质量技术有限公司 | 一种交流电力谐波向量测量方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101951270A (zh) * | 2010-08-24 | 2011-01-19 | 哈尔滨工业大学深圳研究生院 | 一种脉冲超宽带信号的压缩采样接收***及方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2007050680A2 (en) * | 2005-10-25 | 2007-05-03 | William Marsh Rice University | Method and apparatus for on-line compressed sensing |
-
2013
- 2013-07-25 CN CN201310316753.9A patent/CN103389492B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101951270A (zh) * | 2010-08-24 | 2011-01-19 | 哈尔滨工业大学深圳研究生院 | 一种脉冲超宽带信号的压缩采样接收***及方法 |
Non-Patent Citations (1)
Title |
---|
"Beyond Nyquist: Efficient Sampling of Sparse Bandlimited Signals";Joel A. Tropp et al.;《IEEE TRANSACTIONS ON INFORMATION THEORY》;20100131;第56卷(第1期);第520-544页 * |
Also Published As
Publication number | Publication date |
---|---|
CN103389492A (zh) | 2013-11-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103389492B (zh) | 多通道随机谐波调制采样雷达接收机及其方法 | |
CN107064883A (zh) | 一种基于压缩感知的线性调频雷达信号处理方法 | |
Chen et al. | A novel radar waveform compatible with communication | |
Chen et al. | Time-frequency transforms for radar imaging and signal analysis | |
CN104062640B (zh) | 一种外辐射源雷达距离徙动补偿的快速实现方法 | |
CN103412295B (zh) | 基于回波精确模型的高速机动弱目标检测方法 | |
CN102830409B (zh) | 一种基于压缩感知的导航信号采集方法 | |
CN109782250A (zh) | 基于有限新息率采样的雷达目标参数提取方法 | |
CN103235295B (zh) | 基于压缩卡尔曼滤波的小场景雷达目标距离像估计方法 | |
CN109738887B (zh) | 一种适用于微动干扰场景的目标人体运动状态识别方法 | |
CN102798856B (zh) | 一种基于小波包的24GHz LFMCW雷达测距方法 | |
CN108008369A (zh) | 一种非线性调频信号欠采样处理方法 | |
Hao et al. | Research on the method of smart noise jamming on pulse radar | |
CN104007421A (zh) | 基于全变差和压缩感知的罗兰c被动雷达toa估计方法 | |
CN104698456A (zh) | 基于预调制aic的太赫兹雷达成像方法 | |
Cui et al. | Reconstruction filter design for stepped-frequency continuous wave | |
CN103399301B (zh) | 一种宽带sar信号的接收装置及接收方法 | |
CN103630899B (zh) | 地面运动目标高分辨雷达压缩感知成像的方法 | |
CN114167423A (zh) | 基于深度回归网络的雷达海浪参数测量方法 | |
CN102841337A (zh) | 一种天波超视距雷达回波信号去非线性相位污染方法 | |
CN103412289A (zh) | 基于模拟sinc核采样的雷达高分辨距离成像方法 | |
Yan et al. | Through-the-wall human respiration detection using impulse ultra-wide-band radar | |
CN103727960A (zh) | 一种基于drfm的无线电高度表干扰信号产生方法 | |
CN106443630A (zh) | 一种高效的组合式顺序等效采样方法 | |
Jihong et al. | High-speed target ISAR imaging via compressed sensing based on sparsity in fractional Fourier domain |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150902 Termination date: 20160725 |
|
CF01 | Termination of patent right due to non-payment of annual fee |