CN103630878B - 一种基于迭代自适应方法的频带融合方法及装置 - Google Patents
一种基于迭代自适应方法的频带融合方法及装置 Download PDFInfo
- Publication number
- CN103630878B CN103630878B CN201310038953.2A CN201310038953A CN103630878B CN 103630878 B CN103630878 B CN 103630878B CN 201310038953 A CN201310038953 A CN 201310038953A CN 103630878 B CN103630878 B CN 103630878B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- msup
- target echo
- radar
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 53
- 238000007500 overflow downdraw method Methods 0.000 title claims description 9
- 230000004927 fusion Effects 0.000 claims abstract description 42
- 230000001427 coherent effect Effects 0.000 claims abstract description 28
- 238000012545 processing Methods 0.000 claims abstract description 15
- 230000003044 adaptive effect Effects 0.000 claims abstract description 14
- 239000011159 matrix material Substances 0.000 claims description 17
- 238000005070 sampling Methods 0.000 claims description 14
- 238000001228 spectrum Methods 0.000 claims description 12
- 230000003595 spectral effect Effects 0.000 claims description 11
- 230000004044 response Effects 0.000 claims description 8
- 239000013598 vector Substances 0.000 claims description 7
- 230000006978 adaptation Effects 0.000 claims description 2
- 230000017105 transposition Effects 0.000 claims description 2
- 238000003384 imaging method Methods 0.000 abstract description 10
- 238000002156 mixing Methods 0.000 abstract 2
- 238000010586 diagram Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000007499 fusion processing Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000001629 suppression 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
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
-
- 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
- G01S7/285—Receivers
- G01S7/288—Coherent receivers
- G01S7/2883—Coherent receivers using FFT processing
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
本发明公开了一种基于迭代自适应方法的多频带融合方法,其包括:(1)建立多频带雷达信号融合模型;(2)利用全相位快速傅里叶变换方法估计各频带雷达信号的初始相位差及线性相位差,并对各频带雷达信号进行相干补偿;(3)对相干处理后的各频带雷达信号应用迭代自适应方法实现多频带信号的融合。本发明采用全相位快速傅里叶变换方法估计各频带雷达信号的相位差,提高了相位差的估计精度,避免了常规算法中由噪声引起的雷达信号模型参数估计误差对后续相干补偿精度和最终融合成像质量的影响。同时对相干处理后的各频带雷达信号直接应用迭代自适应方法进行融合,避免了常规频带融合算法中需要进行空白频域数据预测的难题。
Description
技术领域
本发明属于目标检测领域,尤其涉及一种基于迭代自适应方法的频带融合方法及装置。
背景技术
在雷达成像中,距离分辨率由雷达发射信号的带宽决定,高距离分辨率的获取需要增加雷达系统发射信号的带宽,大带宽的发射信号则会增加雷达系统的复杂度,同时给接收端的数据采样和后期的数据传输及处理带来了压力。为了克服大带宽带来的难度,在单雷达观测基础上,研究人员提出了许多超分辨率成像算法来实现超分辨成像,但是由于单雷达观测包含的目标信息有限,其分辨率改善能力因此受到了限制,因此更多的研究集中在多频带雷达成像,通过融合多频带的信息等效实现大带宽的高分辨率成像。传统的多频带融合算法大多基于谱估计方法进行,尽管这些方法的估计精度较高,但是估计精度易受噪声的影响,稳健性较差,而且此类算法均需要散射中心的数目作为先验知识。
全相位傅立叶变换(apFFT)的原理如下,假设离散信号x[n]形式为:
其中,N为x[n]点数,频率ω0=nβ2π/N,β可以为小数,则x[n]的傅立叶变换(DFT)谱X[k]为:
apFFT则是从考虑包含某样点所有循环移位后的数据的DFT谱衍生的。对于时间序列中的样点x[0],存在且只存在N个包含该点的N维向量:
x0=[x(0),x(1),…,x(N-1)]T
x1=[x(-1),x(0),…,x(N-2)]T (3)
……
xN-1=[x(-N+1),x(-N+2),…,x(0)]T
将上述每个向量进行循环移位,直至样点x(0)位于首位,然后将移位后的N个向量相加并取平均,即可得到全相位数据向量:
对式(4)进行离散傅立叶变换,即得到apFFT的输出:
对比式(2)与(5)可以发现,apFFT谱幅值为传统FFT幅值的平方,这就意味着旁谱线相对主谱线的比值也按照这种平方关系发生衰减,从而使主谱线更加突出,因而具有很好的谱泄漏抑制性能。更为重要的是,相对于传统FFT各谱线的相位值与频率偏离值(β-k)相关的特点,apFFT谱的相位值为即为中心样点的理论相位值,也就是说apFFT具有相位不变性。而传统FFT在非整周期采样时则需要进行校正才能知道信号的初始相位,而且其精度还不是很高。
发明内容
本发明要解决的技术问题是:针对目前常规谱估计方法估计精度易受噪声的影响,稳健性较差,同时需要目标散射中心个数作为先验知识的问题,提供一种基于迭代自适应方法的多频段融合方法。该方法采用了全相位FFT算法估计各频带雷达信号的相位差,提高了相位差的估计精度,避免了常规算法中由噪声引起的雷达信号模型参数估计误差对后续相干补偿精度和最终融合成像质量的影响,并消除了常规算法在估计雷达信号模型阶数时引入的估计误差。同时对相干处理后的各频带雷达信号直接应用迭代自适应方法进行融合,有效避免了常规算法中需要进行空白频域数据预测的难题。
本发明公开了一种基于迭代自适应方法的多频带融合方法,其包括:
步骤S1:接收多频带雷达目标回波信号,利用全相位快速傅里叶变换方法估计所述多频带雷达目标回波信号间的初始相位差及线性相位差,进而对所述多频带雷达目标回波信号间的初始相位差和线性相位差进行相干补偿;
步骤S2:根据相干补偿后的多频带雷达目标回波信号,建立多频带雷达信号融合模型;
步骤S3:根据所述多频带雷达信号融合模型,应用迭代自适应方法获得融合后的多频带雷达融合信号。
本发明还公开了一种基于迭代自适应方法的多频带融合装置,其包括:
相干补偿装置,其接收多频带雷达目标回波信号,利用全相位快速傅里叶变换方法估计所述多频带雷达目标回波信号间的初始相位差及线性相位差,进而对所述多频带雷达目标回波信号间的初始相位差和线性相位差进行相干补偿;
融合模型建立装置,其根据相干补偿后的多频带雷达目标回波信号,建立多频带雷达信号融合模型;
回波信号融合装置,其根据所述多频带雷达信号融合模型,应用迭代自适应方法获得融合后的多频带雷达融合信号。
本发明基于迭代自适应方法实现了多频带雷达信号的融合,与现有技术相比优点在于:
(1)本发明不需要对各频带雷达信号构造匹配信号模型,避免了信号模型参数估计过程引入的误差,同时不需要目标散射中心个数作为先验知识;
(2)本发明采用全相位FFT算法估计各频段回波信号的初始相位差,较常规算法估计精度高、抗噪声能力强;
(3)本发明对各频段相干数据直接应用迭代自适应方法实现频带融合,避免了常规算法中为了填补频带空白进行的频域数据预测引入的误差。
附图说明
图1为本发明中基于迭代自适应方法的频带融合方法的流程图;
图2为无相位误差时频带融合后的一维距离像示意图;
图3为有相位误差时频带融合后的一维距离像示意图;
图4为初始相位差估计精度示意图;
图5为相位误差随信噪比的变化曲线图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,以下结合具体实施例,并参照附图,对本发明进一步详细说明。
本发明提出了一种基于迭代自适应方法的频带融合方法,如图1所示吗,其具体包括:
步骤S1:接收多频带雷达的回波信号,利用全相位快速傅里叶变换方法估计所述多频带雷达目标回波信号间的初始相位差及线性相位差,进而对所述多频带雷达目标回波信号间的初始相位差和线性相位差进行相干补偿;
步骤S2:根据补偿后的多频带雷达目标回波信号,建立多频带雷达信号融合模型;
步骤S3:根据所述多频带雷达信号融合模型,应用迭代自适应方法获得融合后的多频带雷达目标回波融合信号。
在步骤S1中,具体包括:
步骤S101:接收多频带雷达目标回波信号,进行去调频处理;
假设多频带雷达的发射信号为线性调频信号,雷达i的发射信号载频为fci,脉冲宽度为Ti,发射信号带宽为Bi,调频率为Kri,则雷达i照射目标的发射信号可表示为:
其中,t表示时间,rect(·)为矩形窗函数,i=1,2,…,L,L为雷达数目。
当多频带雷达满足邻近配置时,认为观测目标的散射特性一致,则雷达i接收到的由目标反射回来的目标回波信号为:
其中,σq为目标观测区域内目标等效散射中心q的后向散射系数,Q为目标等效散射中心数目,τiq为目标等效散射中心q到雷达i的时间延迟。
采用去调频(DeChirp)处理后的目标回波信号可表示为:
其中,τis为雷达i接收到的目标回波信号采用DeChirp处理的参考延时,Δτiq=τiq-τis为第q个目标等效散射中心的相对延时。令则目标回波信号可表示为:
步骤S102:对去调频处理后目标回波信号采用全相位傅里叶变换方法处理,获得估计多频带雷达目标回波信号间的初始相位差及线性相位差。
由式(4)的指数项可知,不同频段雷达接收的回波信号存在相位差,需要以其中一个雷达接收的回波信号为参考基准,对其他雷达接收的回波信号进行相位补偿。由于剩余视频相位通常很小,故可以从式(4)中忽略。同时,假定雷达i存在初始相位αi,则由式(4)可得,雷达i(i≠1)经DeChirp后的目标回波信号与雷达1经DeChirp处理后的目标回波信号之间的相位差为:
其中,εiq=Δτiq-Δτ1q为第q个目标等效散射中心的相对延时差,βi=αi-α1为雷达i与雷达1之间的初始相位差,其中α1是雷达1的初始相位。在实际观测中,每个雷达选取的参考目标可能并非对应同一个参考目标。因此,此时的相位误差就变为:
其中,δi是第i部雷达采用参考距离不一致导致的一次相位误差系数,且其是与目标等效散射中心q无关的相位项,因此其对于所有目标等效散射中心而言是公共项,是可以估计出来并补偿掉的;而εiq与目标等效散射中心q是耦合的,不同的散射中心对应的值不同,因此在处理中无法完全补偿,但由于在雷达相邻配置时,相位往往很小,对成像质量影响近似可以忽略,所以相位差最终简化为:
其中,βi为初始相位差,为线性相位差。
由于参考目标点处的目标回波信号的相位为雷达的初始相位,因此可以通过计算所述参考目标散射中心点处目标回波信号的频率响应的相位,得到雷达之间的初始相位差。然后以其中一个雷达接收的回波信号为参考对另外的回波信号进行相干补偿,具体过程如下介绍。
对式(3)进行全相位傅里叶变换(apFFT)处理,可得:
其中,k=0,...,Ni-1,ηiq=-NiKriTsiΔτiq,Tsi为雷达i的采样周期,Ni=Ti/Tsi为对应的采样点数。从式(8)可以看出,对于位于参考目标点处的目标,其回波信号经Dechirp处理后对应的频率为零,并且apFFT相位谱的相位值恒等于αi。因此可以通过计算apFFT相位谱的零频谱线的相位值得到雷达i的初始相位值,进而得到雷达i与雷达1之间的初始相位差βi。
采用apFFT的方法保证了相位估计的准确性,进一步提高了初始相位差补偿的精度。但实际应用中,参考目标点处可能不存在目标,因此在使用apFFT方法前需要对采样数据进行预处理,进而得到线性相位差,具体过程如下介绍。
在两部雷达得到的一维距离像中,由于线性相位差的存在,同一散射中心会被映射到不同的距离单元。因此可以利用互相关方法得到同一散射中心对应在两部雷达一维距离像中的距离单元差,然后由所得距离单元差估计出线性相位差。假定N1=Ni,即雷达1与雷达i的采样点数相同,由式(3)、式(8)可得雷达i(i=2,3,…,L)与雷达1之间经过去调频处理后的目标回波信号的关系:
Sri[k]=Sr1[k]exp{j(βi-2πkδi)} (9)
其中,k=0,1,2,…,Ni-1,计算去调频处理后的目标回波信号之间的互相关系数序列:
其中,n=0,1,…,Ni-1,符号(·)*为共轭算子。则线性相位误差δi可以由下式估计:
δi=-nimax·2π/Ni
其中,nimax为ri[n]序列中最大值对应的索引号。
步骤S103:对所述初始相位差和线性相位差进行相干补偿,获得相干补偿后的多频带雷达目标回波信号。
得到初始相位差βi和线性相位差δi后,以雷达1为参考,对其他雷达回波信号进行相干补偿,得到相干补偿后的多频带雷达目标回波信号表示向量s1,s2,…,sL。
在步骤S2中,所述多频带雷达信号融合模型具体如下:
s=A·y
其中,相干补偿后的多频带雷达目标回波信号L为雷达数目,上标T表示转置操作;第i部雷达经过相干补偿后的采样数据 Ni为采样点数,且回波总点数 y=([y0,y1,...,yK-1]T)K×1表示融合后的多频带雷达目标回波信号的谱响应序列,K为谱响应序列的长度。表示多频带雷达信号的导向矩阵,为第i部雷达对应目标回波信号的导向矩阵,为对应的导向矢量, 为采样时刻。
在步骤S3中,根据步骤2得到的多频带雷达信号融合模型,将基于迭代自适应(IAA)的多频带融合问题转化为如下的加权最小二乘WLS问题,即通过最小化下式,得到目标谱响应序列中的各个元素:
其中,干扰协方差矩阵Q(ωk)=RM-pka(ωk)aH(ωk),pk=|yk|2表示目标回波信号在频率ωk处的功率, IAA协方差矩阵RM=APKAH,PK=diag{p0,p1,p2,…,pK-1},符号diag{·}表示对角矩阵,符号H表示共轭转置。最小化式(10)可得:
针对成像场景的稀疏特性,可采用l1范数作为迭代终止条件,因此基于IAA的多频带融合过程为:
步骤S301:初始化IAA协方差矩阵:IM×M为行、列数均为M的单位阵,置当前迭代次数i=0。
步骤S302:计算第i次迭代时,y中yk的估计值(k=0,1,2,…,K-1),进而得到y的估计值
步骤S303:判断若i≥2且满足且则终止迭代,输出多频带雷达融合信号式中的l1范数定义为否则,判断若i大于或等于最大迭代次数,则终止迭代,输出多频带雷达融合信号否则,令i递增1,更新IAA矩阵其中∑为噪声协方差矩阵,若噪声协方差矩阵∑无法确定,则转至步骤S302进行下次迭代。
下面通过仿真的方法对本发明进行验证。仿真时采用两部雷达验证多频带雷达信号融合算法,假设雷达1工作载频为9GHz,发射线性调频信号带宽为500MHz,雷达2工作载频为10GHz,发射线性调频信号带宽为500MHz。目标为4个点散射中心,相对雷达观测中心的位置分别为-1.3m、0m、1.2m、1.4m。考虑到两部雷达对应的距离分辨率均为0.3m,因此常规的单雷达成像算法并不能将目标3和目标4有效的分开。
首先假设两部雷达的回波信号相位误差为零,图2(a)给出了回波信号信噪比为20dB时,基于IAA的多频带信号融合的目标一维距离像。从图中可以看出,多频带融合方法能够有效的提高雷达成像的分辨率。图2(b)给出了回波信噪比为0dB时,基于IAA的多频带信号融合的目标一维距离像。
然后,设置两部雷达回波信号的初始相位误差为0.3π,测距误差为0.5m,图3给出了信噪比分别为20dB,0dB时频带融合的目标一维距离像。由图2(a)和图3(a)可知,在无噪声情况下,apFFT方法能够精确的估计各回波的相位误差。由图2(b)和图3(b)可知,apFFT方法具有较高的鲁棒性。图4给出了采样数据信噪比为20dB时,采用apFFT方法估计的初始相位精度。图5给出了两部雷达初始相位差为0.2π时,apFFT方法估计的初始相位差随信噪比的变化情况。
以上所述的具体实施例,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (10)
1.一种基于迭代自适应方法的多频带融合方法,其包括:
步骤S1:接收多频带雷达目标回波信号,利用全相位快速傅里叶变换方法估计所述多频带雷达目标回波信号间的初始相位差及线性相位差,进而对所述多频带雷达目标回波信号间的初始相位差和线性相位差进行相干补偿;
步骤S2:根据相干补偿后的多频带雷达目标回波信号,建立多频带雷达信号融合模型;
步骤S3:根据所述多频带雷达信号融合模型,应用迭代自适应方法获得融合后的多频带雷达融合信号。
2.如权利要求1所述的方法,其特征在于,步骤1具体包括以下步骤:
步骤S101:接收多频带雷达目标回波信号,并对所述多频带雷达目标回波信号进行去调频处理;
步骤S102:对所述去调频处理后的多频带雷达目标回波信号进行全相位傅里叶变换,估计获得多频带雷达目标回波信号间的初始相位差及线性相位差;
步骤S103:对所述初始相位差和线性相位差进行相干补偿,获得相干补偿后的多频带雷达目标回波信号。
3.如权利要求2所述的方法,其特征在于,步骤S101中,采用DeChirp处理对多频带雷达目标回波信号进行去调频,去调频后的多频带雷达目标回波信号如下表示:
其中,fci、Ti和Kri分别为雷达i的发射信号载频、发射信号脉冲宽度和调频率,i=1,2,…,L,L为雷达数目;rect(·)为矩形窗函数;σq为目标的等效散射中心q的后向散射系数,Q为目标的等效散射中心数目;τiq为目标的等效散射中心q到雷达i的时间延迟,τis为雷达i接收到的目标回波信号采用DeChirp处理的参考延时,Δτiq=τiq-τis为第q个目标等效散射中心的相对延时。
4.如权利要求3所述的方法,其特征在于,步骤S102中对去调频后的多频带雷达目标回波信号进行全相位傅里叶变换得:
其中,k=0,…,Ni-1,ηiq=-NiKriTsiΔτiq,Tsi为雷达i的采样周期,Ni=Ti/Tsi为对应的采样点数,αi为雷达i的初始相位;
通过零频谱线的相位值得到多频带雷达目标回波信号的初始相位值,进而得到所述多频带雷达目标回波信号间的初始相位差;
所述线性相位差如下计算:
δi=-nimax·2π/Ni
其中,nimax为ri[n]序列中最大值对应的索引号。
5.如权利要求2所述的方法,其特征在于,步骤S103中,在得到初始相位差和线性相位差后,以其中一部雷达目标回波信号为参考,对其它雷达目标回波参考信号进行相位补偿,进而得到相干补偿后的多频带雷达目标回波信号。
6.如权利要求1所述的方法,其特征在于,步骤S2中所建立的多频带雷达信号融合模型如下表示:
s=A·y
其中,相干补偿后的多频带雷达目标回波信号L为雷达数目,上标T表示转置操作;第i部雷达经过相干补偿后的目标回波信号的采样数据Ni为采样点数,且多频带雷达目标回波信号的回波总点数y=([y0,y1,…,yK-1]T)K×1表示融合后的多频带雷达目标回波信号的谱响应序列,K为谱响应序列的长度;
表示多频带雷达目标回波信号的导向矩阵,为第i部雷达对应的目标回波信号的导向矩阵, 为对应的导向矢量, 为采样时刻。
7.如权利要求6所述的方法,其特征在于,在步骤S3中,所述多频带雷达目标回波信号的融合通过将迭代自适应问题转化为最小二乘问题来实现,所述最小二乘问题的表示模型如下:
其中,干扰协方差矩阵Q(ωk)=RM-pka(ωk)aH(ωk),pk=|yk|2表示目标回波信号在频率ωk处的功率,迭代自适应协方差矩阵RM=APKAH,PK=diag{p0,p1,p2,…,pK-1},K为谱响应序列的长度,符号diag{·}表示对角矩阵,符号H表示共轭转置。
8.如权利要求7所述的方法,其特征在于,通过求解上述最小二乘问题的表示模型,得到融合后的多频带雷达目标回波信号的谱响应序列表示模型:
9.如权利要求6所述的方法,其特征在于,步骤S3具体包括如下步骤:
步骤S301:置当前迭代次数i=0,初始化迭代自适应协方差矩阵IM×M为行、列数均为M的单位阵;
步骤S302:进入第i次迭代,计算融合后的多频带雷达目标回波信号的谱响应序列y中各元素yk的估计值:
其中,
步骤S303:判断若i≥2且满足且则终止迭代,并输出融合后的多频带雷达目标回波融合信号式中的l1范数定义为否则,判断若i大于或等于最大迭代次数,则终止迭代,并输出融合后的多频带雷达目标回波融合信号否则,令i递增1,更新迭代自适应矩阵其中∑为噪声协方差矩阵,并转至步骤S302进行下次迭代。
10.一种基于迭代自适应方法的多频带融合装置,其包括:
相干补偿装置,其接收多频带雷达目标回波信号,利用全相位快速傅里叶变换方法估计所述多频带雷达目标回波信号间的初始相位差及线性相位差,进而对所述多频带雷达目标回波信号间的初始相位差和线性相位差进行相干补偿;
融合模型建立装置,其根据相干补偿后的多频带雷达目标回波信号,建立多频带雷达信号融合模型;
回波信号融合装置,其根据所述多频带雷达信号融合模型,应用迭代自适应方法获得融合后的多频带雷达融合信号。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310038953.2A CN103630878B (zh) | 2013-01-31 | 2013-01-31 | 一种基于迭代自适应方法的频带融合方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310038953.2A CN103630878B (zh) | 2013-01-31 | 2013-01-31 | 一种基于迭代自适应方法的频带融合方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103630878A CN103630878A (zh) | 2014-03-12 |
CN103630878B true CN103630878B (zh) | 2015-08-12 |
Family
ID=50212129
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310038953.2A Active CN103630878B (zh) | 2013-01-31 | 2013-01-31 | 一种基于迭代自适应方法的频带融合方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103630878B (zh) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104182951A (zh) * | 2014-08-15 | 2014-12-03 | 张建伟 | 一种基于双摄像头全景的多频段图像缝合方法 |
CN104331857A (zh) * | 2014-10-28 | 2015-02-04 | 南京理工大学 | 光强传输方程相位恢复中的相位差异迭代补偿方法 |
CN104459635B (zh) * | 2014-12-08 | 2017-03-01 | 西安科技大学 | 基于迭代收缩加权融合的自适应空滤滤波方法 |
CN105785123B (zh) * | 2016-03-22 | 2018-04-06 | 电子科技大学 | 一种基于apFFT相位差的雷达信号频率计算方法 |
CN107272000B (zh) * | 2017-07-20 | 2020-08-07 | 中国科学院电子学研究所 | 滑动散射中心方位向相位误差的计算方法 |
CN109031289B (zh) * | 2018-07-13 | 2020-06-09 | 清华大学 | 认知捷变频雷达波形设计方法和装置 |
CN110794401A (zh) * | 2019-11-20 | 2020-02-14 | 长沙理工大学 | 多频带雷达融合成像方法、装置、设备及可读存储介质 |
CN112596051B (zh) * | 2020-12-11 | 2023-08-15 | 中国科学院空天信息创新研究院 | 基于改进松弛算法的多带融合算法 |
CN113376634A (zh) * | 2021-05-25 | 2021-09-10 | 西安电子科技大学 | 基于gniaa的舰船目标isar成像方法 |
CN114942415A (zh) * | 2022-05-17 | 2022-08-26 | 哈尔滨工程大学 | 一种基于自编码器的多频段子带信号融合方法 |
CN115856883B (zh) * | 2023-02-14 | 2023-05-16 | 北京优诺信创科技有限公司 | 一种基于互补随机波形的双基雷达协同成像系统 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2009008451A (ja) * | 2007-06-26 | 2009-01-15 | Toyota Motor Corp | 広帯域レーダ装置 |
CN101614810A (zh) * | 2008-06-25 | 2009-12-30 | 电子科技大学 | 一种线阵三维成像合成孔径雷达分辨率融合方法 |
-
2013
- 2013-01-31 CN CN201310038953.2A patent/CN103630878B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2009008451A (ja) * | 2007-06-26 | 2009-01-15 | Toyota Motor Corp | 広帯域レーダ装置 |
CN101614810A (zh) * | 2008-06-25 | 2009-12-30 | 电子科技大学 | 一种线阵三维成像合成孔径雷达分辨率融合方法 |
Non-Patent Citations (5)
Title |
---|
《多波段雷达融合一维超分辨率成像技术研究》;付耀文 等;《自然科学进展》;20061030;第16卷(第10期);第1310-1316页 * |
基于稀疏贝叶斯学习的多频带雷达信号融合;叶钒等;《电波科学学报》;20101030;第25卷(第05期);第990-994页 * |
梁福来.《多频段合成成像技术研究》.《中国优秀硕士学位论文全文数据库(信息科技辑)》.2010,(第05期), * |
王海峰.《多频段合成雷达成像技术》.《中国优秀硕士学位论文全文数据库(信息科技辑)》.2012,(第07期), * |
稀疏子带的多频段雷达信号融合超分辨距离成像;王成等;《电子学报》;20060625;第34卷(第06期);第985-990页 * |
Also Published As
Publication number | Publication date |
---|---|
CN103630878A (zh) | 2014-03-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103630878B (zh) | 一种基于迭代自适应方法的频带融合方法及装置 | |
CN111736131B (zh) | 一种剔除一比特信号谐波虚假目标的方法及相关组件 | |
CN103760558B (zh) | 一种太赫兹雷达isar成像方法 | |
CN112882005B (zh) | 一种基于贝叶斯学习的otfs雷达目标参数估计方法 | |
CN107797099B (zh) | 一种多通道数字接收机实时内定标处理方法及装置 | |
CN103901429A (zh) | 基于稀疏孔径的机动目标逆合成孔径雷达成像方法 | |
CN109633538B (zh) | 非均匀采样系统的最大似然时差估计方法 | |
US11982698B2 (en) | Joint denoising and delay estimation for the extraction of pulse-width of signals in RF interference | |
CN106199600B (zh) | 基于多普勒估计的方位多通道合成孔径雷达成像方法 | |
CN109471095A (zh) | 基于快速迭代插值的fmcw雷达距离估计算法 | |
CN103616687A (zh) | 分段线性估计的多项式拟合isar包络对齐方法 | |
US10061014B2 (en) | Radar signal processing method and apparatus for compensating for in-phase/quadrature-phase (I/Q) imbalance | |
CN104698431A (zh) | 模糊分量空间角的估计及多通道sar方位解模糊方法 | |
CN111935038B (zh) | 基于分数阶傅里叶变换的线性调频干扰消除方法 | |
CN105024951A (zh) | 一种功率时延谱pdp估计方法及装置 | |
CN106842198A (zh) | 一种基于凸优化理论的子带融合成像算法 | |
CN111781575B (zh) | 一比特信号单快拍波达方向估计方法及相关组件 | |
CN106646531B (zh) | 一种多星约束稳健空频抗干扰处理方法和装置 | |
CN111007457A (zh) | 一种基于块稀疏贝叶斯模型的辐射源直接定位方法 | |
CN103454615A (zh) | 稳健盲信号参数估计方法 | |
US7706426B2 (en) | Method for estimating one or more parameters of a ultra wideband signal and a receiver system for receiving ultra wideband signals | |
KR102076096B1 (ko) | 레이더 신호 처리 장치 및 방법 | |
CN113325382B (zh) | 基于全局最优化相位调制的ddma mimo雷达速度解模糊方法 | |
CN106534028B (zh) | 基于压缩感知的ofdm信道跟踪方法 | |
CN113359095B (zh) | 一种相干被动mimo雷达克拉美罗界的计算方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | 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 |