基于雷达大脉宽信号的目标测速测距方法
【专利摘要】本发明公开了一种基于雷达大脉宽信号的目标测速测距方法，其主要思路为：雷达发射端发射N个大脉宽的脉冲信号s(t)，根据N个大脉宽的脉冲信号s(t)，计算雷达接收端接收的回波信号y(t)，所述雷达接收端接收的回波信号y(t)包含M个目标的信息，每一个目标的信息包含目标的斜距和目标的径向速度；根据雷达接收端接收的回波信号y(t)，计算回波信号y(t)的相关矩阵R；根据回波信号y(t)的相关矩阵R，计算得到M个目标各自的多普勒频率和M个目标各自的目标径向速度；根据所述M个目标各自的目标径向速度，对雷达接收端接收的回波信号y(t)进行多普勒频率补偿，计算多普勒频率补偿后M个目标的回波信号，然后对多普勒频率补偿后M个目标的回波信号进行脉冲压缩处理，计算M个目标各自的斜距。
【专利说明】
基于雷达大脉宽信号的目标测速测距方法
技术领域
[0001 ]本发明属于雷达发射信号中测量目标距离和目标速度的技术领域，特别涉及一种 基于雷达大脉宽信号的目标测速测距方法，通过雷达大脉宽信号实现目标速度和目标距离 的准确测量，适用于高分辨雷达观测远距离高速的空中飞行目标的距离及速度的准确测 量。
【背景技术】
[0002] 随着空警戒雷达以及各种飞行器技术的发展，使得雷达的作用距离、分辨能力和 测量精度等性能指标必须得到相应提尚；空基雷达或者天基雷达对空中尚速目标和低空远 距离飞行目标的探测分别也变得越来越重要，而且准确检测出目标速度和目标距离对后期 的目标跟踪也至关重要。
[0003] 传统的脉冲雷达为了提高距离分辨率以及减少脉冲遮挡问题，多发射窄脉宽信 号，且窄脉宽信号的脉宽为微秒级别，发射窄脉宽信号能量低;为了检测空中远距离目标时 需要长时间的积累检测，导致高速运动目标在长时间的积累检测时产生了不同程度的距离 走动，影响相参积累的效果；并且，根据脉冲重复频率（PRF)的高低，现代雷达有低重频 (LPRF)、中重频(MPRF)和高重频(HPRF)三种工作模式，传统的雷达测量多普勒时利用多个 脉冲做快速傅里叶变换（FFT)处理，且传统的雷达能够测量的无模糊多普勒的范围为
，fr为脉冲重复频率。当雷达工作在低重频、目标运动速度较快、且雷达多普勒 频率大于雷达的脉冲重复频率时目标会出现速度模糊;而雷达测量的最大不模糊的距离计
，(：为光速；当雷达工作在高重频时，由于脉冲的周期较短，容易造成目 标距离模糊;传统的脉冲雷达对于检测空中远距离高速运行目标时会同时存在测速与测距 模糊问题，因此传统测量目标速度与目标距离的方法多采用测量出目标的模糊距离和目标 的模糊速度后再采用多重频或者多载频的方法分别解除目标速度和目标距离中的模糊，从 而求得目标的真实速度和目标的真实距离；而传统的方法测量距离和速度的方法比较复 杂，并且误差大精度不高。
【发明内容】
[0004] 针对以上现有技术适用条件的不足，以及大脉宽信号对检测空中高速远距离目标 的优势，本发明的目的在于提出一种基于雷达大脉宽信号的目标测速测距方法，该种基于 雷达大脉宽信号的目标测速测距方法能够实现空中远距离高速飞行的目标距离和目标速 度的准确测量。
[0005] 本发明的主要思路是:考虑到目标的多普勒频率在脉冲内可以体现，所以选择雷 达的发射脉冲为大脉宽的脉冲信号，并获取雷达接收端接收的回波信号，然后对雷达接收 端接收的回波信号采用迭代自适应算法(IAA)估计出雷达接收端接收的回波信号的相关矩 阵R，再结合多信号分类(Multiple Signal Class if i cat ion，MUSIC)算法，对相关矩阵R进 行特征值分解，进而计算得到目标的多普勒频率，从而实现目标多普勒频率的高分辨处理， 得到目标的径向速度;最后利用目标的径向速度对雷达接收端接收的回波信号进行多普勒 频率补偿后再做脉冲压缩处理，获得目标斜距，进而达到检测出目标径向速度和目标斜距 的目的。
[0006] 为达到上述技术目的，本发明采用如下技术方案予以实现。
[0007] 一种基于雷达大脉宽信号的目标测速测距方法，包括以下步骤：
[0008] 步骤1，雷达发射端发射N个大脉宽的脉冲信号s(t)，脉冲信号s(t)的脉冲重复周 期为IV，脉冲信号s(t)的脉冲重复频率为
脉冲信号s(t)的脉冲宽度为τ，脉冲宽度τ 为毫米级时对应的脉冲为大脉宽，脉冲信号s(t)波长为λ〇:
为脉冲信号s(t)的中 心频率，C为光速，N为自然数；
[0009] 步骤2,根据N个大脉宽的脉冲信号s(t)，计算得到雷达接收端接收的回波信号y (t)，所述雷达接收端接收的回波信号y(t)包含Μ个目标的信息，每一个目标的信息包含目 标的斜距和目标的径向速度；
[0010]步骤3,根据雷达接收端接收的回波信号y(t)，估计得到回波信号y(t)的相关矩阵 R；
[0011] 步骤4,根据回波信号y(t)的相关矩阵R，依次计算得到Μ个目标各自的多普勒频率 和Μ个目标各自的径向速度；
[0012] 步骤5,根据所述Μ个目标各自的多普勒频率，对雷达接收端接收的回波信号y(t) 进行多普勒频率补偿，计算得到多普勒频率补偿后Μ个目标的回波信号，然后对多普勒频率 补偿后Μ个目标的回波信号进行脉冲压缩处理，计算得到Μ个目标各自的斜距。
[0013] 本发明有益效果：
[0014] 第一，本发明中的雷达采用脉冲发射方式发射大脉宽毫秒级的脉冲信号，该脉冲 信号能量大，能够观测远距离目标，并且大脉宽信号能够保证远距离空中的目标距离不存 在模糊，且多普勒频率在单个脉冲内就能够体现出来，使得直接采用单个脉冲就能够实现 目标测距和目标测速;另外大脉宽信号能量大，探测目标距离远；当主要探测空中远距离目 标时，忽略近距离中脉冲遮挡问题；
[0015] 第二，传统的对单脉冲的回波信号做快速傅里叶变换(FFT)处理的方法进行目标 测速时多普勒分辨率很低，即多普勒分辨率为脉宽的倒数;而本发明采用迭代自适应算法 (ΙΑΑ)结合多信号分类(Mul tip le Signal Class if i cat ion，MUSIC)算法的高分辨处理法， 能够使得单脉冲回波信号的多普勒分辨率约为快速傅里叶变换(FFT)方法分辨率的5~10 倍，提高目标测速的精度；
[0016] 第三，传统的目标测速方法中多普勒频率的不模糊区间为
为脉冲 重复频率;但本发明中测速方法，多普勒频率的不模糊区间大为
不8为采样频 率，Fs> >fr，采用本发明方法不用考虑目标距离和目标速度的模糊问题。
【附图说明】
[0017] 下面结合附图和【具体实施方式】对本发明作进一步详细说明。
[0018] 图1为本发明的一种基于雷达大脉宽信号的目标测速测距方法流程图；
[0019] 图2为本发明中雷达发射的大脉宽脉冲信号的示意图，其中，Tr为脉冲信号的脉冲 重复周期，τ为脉冲信号的脉冲宽度，N为雷达发射的大脉宽脉冲信号包含的脉冲个数；
[0020] 图3为脉冲信号的多普勒频率分别为4000Hz和0Hz时进行脉冲压缩处理后的对比 示意图；
[0021] 图4为雷达发射端发射的脉冲信号进行快速傅里叶变换(FFT)后的频谱示意图； [0022]图5(a)为两个目标的多普勒频率分别为2000Hz和4000Hz时，采用快速傅里叶变换 (FFT)方法测得的目标的多普勒频率示意图，
[0023]图5(b)为两个目标的多普勒频率分别为3000Hz和4000Hz时，采用快速傅里叶变换 (FFT)方法测得的目标的多普勒频率示意图；
[0024]图6(a)为设置两个目标的多普勒频率分别为[3000Hz 4000Hz]时采用IAA算法结 合MUSIC算法分别估计两个目标的多普勒频率示意图，
[0025]图6(b)为两个目标的多普勒频率分别为[3500Hz 4000Hz]时采用IAA算法结合 MUSIC算法分别估计两个目标的多普勒频率示意图，
[0026]图6(c)为两个目标的多普勒频率分别为[3800Hz 4000Hz]时采用IAA算法结合 MUSIC算法分别估计两个目标的多普勒频率示意图，
[0027]图6(d)为两个目标的多普勒频率分别为[3850Hz 4000Hz]时采用IAA算法结合 MUSIC算法分别估计两个目标的多普勒频率示意图，
[0028]图6(e)为两个目标的多普勒频率分别为[3900Hz 4000Hz]时采用IAA算法结合 MUSIC算法分别估计两个目标的多普勒频率示意图；
[0029] 图7(a)为目标存在多普勒频率，且雷达接收端接收的回波信号没有经过速度补 偿，直接经过脉冲压缩处理得到的结果示意图，
[0030] 图7(b)为1000Hz的目标多普勒频率补偿掉后再经过脉压处理后得到的结果示意 图，
[0031] 图7(c)为4000Hz的目标多普勒频率补偿掉后再经过脉压处理后得到的结果示意 图。
【具体实施方式】
[0032] 参照图1，为本发明的一种基于雷达大脉宽信号的目标测速测距方法流程图，该种 基于大脉宽信号优化设计的测速测距技术，包括以下步骤：
[0033] 步骤1，雷达发射端采用脉冲发射方式发射N个大脉宽的脉冲信号s(t)，脉冲信号s (t)的脉冲重复周期为IV，脉冲信号s(t)的脉冲重复频率为
，脉冲信号s(t)的脉冲宽 度为τ，脉冲宽度τ为毫米级时对应的脉冲为大脉宽，脉冲信号s(t)波长为λ〇，
，:fo为脉 冲信号s(t)的中心频率，C为光速，脉冲信号s(t)的初始相位为θ〇，Ν为自然数。
[0034]具体地，参照图2,为本发明中雷达发射的大脉宽脉冲信号的示意图，其中，Tr为脉 冲信号的脉冲重复周期，τ为脉冲信号的脉冲宽度，N为雷达发射的大脉宽脉冲信号包含的 脉冲个数;所述Ν个大脉宽的脉冲信号s (t )，其具体表达式为：
[0036] 其中，1￡{1，2，一，《4为雷达发射端发射大脉宽的个数，1^为脉冲信号8(〇的脉 冲重复周期，τ为脉冲信号s(t)的脉冲宽度，fo为脉冲信号s(t)的中心频率，θ〇为脉冲信号 s (t)的初始相位，rect( ·)为矩形函数，t为时间变量，且满足，_N为 自然数。
[0037] 步骤2,根据Ν个大脉宽的脉冲信号s(t)，计算得到雷达接收端接收的回波信号y (t)，所述雷达接收端接收的回波信号y(t)包含Μ个目标的信息，每一个目标的信息包含目 标的斜距和目标的径向速度。
[0038] 具体地，根据Ν个大脉宽的脉冲信号s(t)，计算得到雷达接收端接收的回波信号y (t)，其表达式为：
[0040] 其中，雷达接收端接收的回波信号y(t)中包含Μ个目标的信息，每一个目标的信息 包含目标的斜距和目标的径向速度，第i个目标的斜距为心，第i个目标的径向速度为Vn，S (t)为雷达发射端发射的N个大脉宽的脉冲信号，t为时间变量，1^{1，2，一，}，1^为第1个 目标的回波衰减，A tl为第i个目标的回波对应的延迟时间，
，(：为光速，fdl为第i个 目标的多普勒频率：
，λ〇为脉冲信号s (t)的波长。
[0041] 步骤3,对雷达接收端接收的回波信号y(t)，采用迭代自适应（IAA)算法计算得到 回波信号y(t)的相关矩阵R。
[0042] 具体地，假定雷达接收端接收的回波信号y(t)的采样频率为Fs，则回波信号y(t) 测量多普勒频率的不模糊区间为
，在不模糊区间
内等间隔取K个 频率点，K个频率点对应的多普勒频率向量为fd，fd = {fdl，fd2，……fdd ;考虑雷达接收端接 收的回波信号y(t)内存在噪声干扰，所以得到回波信号y(t)的离散形式y(n)，其表达式为：
[0043] y(n)=A(fd)Xs(n)+e(n)
[0044] s(n) = [si(n)，…，sk(n)，…，SK(n)]
[0045] A(fd) = [a(fdi)，…，a(fdk)，…，a(fdK)]
[0047]其中^^丨^^-丄丨^^一为冗个频率点对应的多普率频率构成的阵列流行矩 阵，L为回波信号y(t)脉冲内的采样点数，e(n)为K个频率点对应的回波信号中包含的噪声， 一般为高斯白噪声；s(n)为K个频率点对应的回波信号波形组成的矩阵，Sk(n)为第η个采样 点处第k个频率点对应的回波信号波形，a(fdk)为第k个频率点的对应多普勒频率组成的阵 列流行向量，fdk为第k个频率点的对应多普勒频率，匕为雷达接收端接收的回波信号y (t)的 采样频率，k= {1，2，…，K}，K为回波信号y (t)测量多普勒频率的不模糊区间 内等间隔取的频率点个数。
[0048] 对雷达接收端接收的回波信号y(t)，采用迭代自适应（IAA)算法计算得到回波信 号y (t)的相关矩阵R，其具体子步骤为：
[0049] 3.1初始化:构造一个KXK对角矩阵P，KXK对角矩阵P对角线上的元素为等间隔取 的Κ个频率点各自对应的回波功率，即Κ个频率点各自对应的回波功率，与Κ X Κ对角矩阵Ρ对 角线上的元素分别相同;k表示第k个频率点，k={l，2,……，1(}，1(\1(对角矩阵？的初始矩阵
为Pi，初始矩阵Pi对角线上包含K个元素值，其中第k个元素值为为， a(fdk)为第k个频率点的对应多普勒频率组成的阵列流行向量，y(n)为回波信号y(t)的离散 形式；回波信号y(t)的相关矩阵初始值为心，
表示定义;设定迭代总次 数为D，m表示第m次迭代，m初值为1。
[0050] 3.2计算得到第m次迭代后第k个频率点对应的信号
a(fdk)为第k个频率点的对应多普勒频率组成的阵列流行向量，fdk为第k个频率点的对应多 普勒频率，y(n)为回波信号y(t)的离散形式，上标Η为共辄转置，上标-1为求逆操作。
[0051 ] 3.3根据第m次迭代后第k个频率点对应的信号波形&⑷，计算得到第m次迭代后 第k个频率点对应的回波功率值
[0052] 3.4令k加1，返回子步骤3.2,直到得到第m次迭代后第K个频率点对应的回波功率 值?恤，并根据此时得到的第m次迭代后第1个频率点对应的回波功率值P ml到第m次迭代后第 K个频率点对应的回波功率值PmK，得到第m次迭代后的KXK对角矩阵?"，然后令m加1，返回子 步骤3.2,直到得到第D次迭代后的KXK对角矩阵P D，并将根据第D次迭代后得到KXK对角矩 阵Pd，计算得到第D次迭代后得到的回波信号y(t)的相关矩阵R D，并将第D次迭代后得到的回 波信号y(t)的相关矩阵Rd，作为回波信号y(t)的相关矩阵R。
[0053] 步骤4,根据回波信号y(t)的相关矩阵R，利用MUSIC算法依次计算得到Μ个目标各 自的多普勒频率和Μ个目标各自的径向速度。
[0054] 具体地，根据回波信号y(t)的相关函数R，利用MUSIC算法依次计算得到Μ个目标各 自的多普勒频率和Μ个目标各自的目标径向速度，其具体子步骤为：
[0055] 4.1利用MUSIC算法对回波信号y(t)的相关函数R进行特征值分解，即R = U2UH，得 到F个特征值以及F个特征值对应的特征向量U，Σ为F个特征值对应的对角矩阵;对F个特征 值进行从大到小排序，得到：
[0056] λχ^λ2^ · · · · · ·λρ
[0057] Vl , V2 , , VM, VM+1, , VF
[0058] 其中，F个特征值中包含Μ个大的特征值和F-M个接近零的小特征值，Μ个大的特征 值为XgljAM，F-M个接近零的小特征值为1(? +1到\^^到#为大到小排序后?个特征值各自对应 的特征向量。
[0059] 4.2将F-M个接近零的小特征值分别对应的特征向量构成回波信号y (t)的噪声子 5^iClG，G = span{vM+i，VM+2，."，VF};
[0060] 4.3使用K个频率点对应的多普率频率构成的阵列流行矩阵A (fd)对回波信号y (t) 的噪声子空间G作投影，计算得到回波信号y(t)的功率谱S(fd) 回波信号y (t)的功率谱S( fd)的谱峰处对应的多普勒频率向量，即为Μ个大的特征值对应的 多普勒频率向量，进而依次得到Μ个目标各自的多普勒频率和Μ个目标各自的径向速度;其 中第i个大的特征值对应的第i个目标的多普勒频率f dl，并根据第i个目标的多普勒频率 fdl，计算得到第i个目标的径向速度vri，
，λ〇为脉冲信号s(t)波长? = {1，2，···，Μ}。
[0061] 步骤5,根据所述Μ个目标各自的多普勒频率，对雷达接收端接收的回波信号y(t) 进行多普勒频率补偿，计算得到多普勒频率补偿后Μ个目标的回波信号，然后对多普勒频率 补偿后Μ个目标的回波信号进行脉冲压缩处理，计算得到Μ个目标各自的斜距。
[0062] 具体地，设定Μ个目标各自的多普勒频率补偿因子，然后根据所述Μ个目标各自的 多普勒频率，对雷达接收端接收的回波信号y(t)进行多普勒频率补偿，计算得到多普勒频 率补偿后Μ个目标的回波信号;在设定的Μ个目标各自的多普勒频率补偿因子中，第i个目标 的多普勒频率补偿因子为并根据雷达接收端接收的回波信号y(t)，计算得到多普 勒频率补偿后第i个目标的回波信号_r(〇xe~>w，进而计算得到多普勒频率补偿后Μ个目 标的回波信号Y(t)，其表达式为：
[0064]然后对多普勒频率补偿后Μ个目标的回波信号Y(t)进行脉冲压缩处理，计算得到 脉冲压缩处理后Μ个目标的回波信号B(t)，其表达式为：
[0066]
为与Μ个目标各自的斜距有关的分量，C为光速，R为Μ个目标各自的斜 距集合，Φ为多普勒频率补偿后Μ个目标的回波信号Y(t)中多普勒频率没有被补偿完的部 分；由于Μ个目标中的每一个目标都存在多普勒频率，导致脉冲压缩处理后Μ个目标的回波 信号B(t)恶化，由此说明脉冲信号s(t)是多普勒敏感的，进而说明没有被补偿完的部分Φ 对脉冲压缩处理后Μ个目标的回波信号B(t)影响不大，可以忽略中的R为Μ个目标各 自的斜距，即得到Μ个目标各自的斜距。
[0067]下面通过仿真实验对本发明作进一步验证说明。
[0068](一)实验条件
[0069] 本发明中雷达发射的脉冲信号为大脉宽信号s(t)，即脉宽为毫秒级，噪声为高斯 白噪声，采用本发明的方法进行空中高速远距离目标进行测速测距时，脉冲信号s(t)的脉 冲重复周期T r = 3ms，脉冲信号s(t)的脉冲宽度为T = lms，雷达接收端接收的回波信号y(t) 的采样频率Fs = 0.2MHz，目标个数设定为2个，2个目标对应的多普勒频率fd设置分别为 [2000Hz 4000Hz],[3000Hz 4000Hz],[3500Hz 4000Hz],[3800Hz 4000Hz],[3850Hz 4000Hz]，[3900Hz 4000Hz];2个目标的距离R设置为[100km 200km]，2个目标对应的距离门
尤=3\108为光速，行1(>)为取小于*的最大整数。通过上述公式计算 目标的距离门，则2个目标对应的距离门个数分别为Δ n i和Δ n 2，
[0070] 采用MUSIC算法获得的2个目标的搜索矢量a(fd)中fd的搜索范围为-5000ΗΖ~ 5000Hz，2个目标对应的信噪比 SNR1 = 10dB，SNR2 = 10dB。
[0071] (二)实验内容
[0072]参照图3,为脉冲信号的多普勒频率分别为4000Hz和0Hz时进行脉冲压缩处理后的 对比示意图；参照图4,为雷达发射端发射的脉冲信号进行快速傅里叶变换(FFT)后的频谱 示意图。
[0073]图5(a)为两个目标的多普勒频率分别为2000Hz和4000Hz时，采用快速傅里叶变换 (FFT)方法测得的目标的多普勒频率示意图，图5(b)为两个目标的多普勒频率分别为 3000Hz和4000Hz时，采用快速傅里叶变换(FFT)方法测得的目标的多普勒频率示意图；图6 (a)为设置两个目标的多普勒频率分别为[3000Hz 4000Hz]时采用IAA算法结合MUSIC算法 分别估计两个目标的多普勒频率示意图，图6(b)为两个目标的多普勒频率分别为[3500Hz 4000Hz]时采用IAA算法结合MUSIC算法分别估计两个目标的多普勒频率示意图，图6(c)为 两个目标的多普勒频率分别为[3800Hz 4000Hz]时采用IAA算法结合MUSIC算法分别估计两 个目标的多普勒频率示意图，图6(d)为两个目标的多普勒频率分别为[3850Hz 4000Hz]时 采用IAA算法结合MUSIC算法分别估计两个目标的多普勒频率示意图，图6(e)为两个目标的 多普勒频率分别为[3900Hz 4000Hz]时采用IAA算法结合MUSIC算法分别估计两个目标的多 普勒频率示意图。
[0074]图7(a)为目标存在多普勒频率，且雷达接收端接收的回波信号没有经过速度补 偿，直接经过脉冲压缩处理得到的结果示意图，图7(b)为1000Hz的目标多普勒频率补偿掉 后再经过脉压处理后得到的结果示意图，图7(c)为4000Hz的目标多普勒频率补偿掉后再经 过脉压处理后得到的结果示意图；第一个目标的斜距Ri = 100km，第二个目标的斜距R2 = 200km，也就是分别对应第133个和266个距离门的信号速度补偿前的脉压结果和速度补偿 掉后的脉压结果。
[0075](三)结果分析
[0076]由图3可以看出，当脉冲信号的多普勒频率不为0时，脉压结果差，无法获取目标的 斜距，可以得出脉冲信号是对多普勒频率比较敏感的信号，故符合雷达发射端发射要选的 信号标准；图4为发射脉冲的频谱，脉冲信号的脉冲宽度τ = lms，故理论上主瓣宽度为
，也就是采用快速傅里叶变换(FFT)方法测量目标的多普勒分辨率为1000Hz，实 际为(833.8HZ+166.8) = 1000.6Hz。图5(a)是采用FFT的分析方法测量目标的多普勒频率， 两个目标的多普勒频率相差为2000Hz，所以可以分辨开来，但图5(b)是采用FFT分析的方法 测量目标的多普勒频率，因为两个目标的多普勒频率相差为1000Hz，已经无法分辨开，可知 采用传统的快速傅里叶变换（FFTFFT)的方法测量目标的多普勒频率时最大的分辨率为
，分辨率比较低。
[0077]故在相同的参数条件下，采用IAA算法结合MUSIC算法的高分辨处理方法测量目标 的多普勒频率，试验结果为图6(a)到图6(e);可以看出，采用IAA算法结合MUSIC算法的高分 辨处理方法的分辨率在150Hz左右，能够达到FFT方法的5~10倍，分辨率高，测量多普勒频 率的精度高，测量准确。
[0078] 对目标的多普勒频率为4000Hz，斜距R = 200km，也就是第266个距离门的信号，图7 (a)-7(c)实验结果为速度补偿前的脉压结果和速度补偿掉后的脉压结果对比，对应第133 个和266个距离门的信号速度补偿前的脉压结果和速度补偿掉后的脉压结果由图7(b)可以 看出脉压后目标的距离门约为第267个，故通过检测出目标的速度后通过先补偿掉多普勒 频率再进行脉压处理来测量目标的斜距。
[0079]综上所述，仿真实验验证了本发明的正确性，有效性和可靠性。
[0080]显然，本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精 神和范围；这样，倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围 之内，则本发明也意图包含这些改动和变型在内。
【主权项】
1. 一种基于雷达大脉宽信号的目标测速测距方法，其特征在于，包括W下步骤： 步骤1，雷达发射端发射N个大脉宽的脉冲信号s(t)，脉冲信号s(t)的脉冲重复周期为 Tr，脉冲信号s(t)的脉冲重复频率为乂 =^，脉冲信号s(t)的脉冲宽度为τ，脉冲宽度τ为毫 米级时对应的脉冲为大脉宽，脉冲信号s(t)波长为λ〇,為=^，C为光速，fo为脉冲信号s(t) 的中屯、频率，N为自然数； 步骤2,根据N个大脉宽的脉冲信号s(t)，计算得到雷达接收端接收的回波信号y(t)，所 述雷达接收端接收的回波信号y(t)包含Μ个目标的信息，每一个目标的信息包含目标的斜 距和目标的径向速度； 步骤3,根据雷达接收端接收的回波信号y(t)，计算得到回波信号y(t)的相关矩阵R; 步骤4,根据回波信号y(t)的相关矩阵R，依次计算得到Μ个目标各自的多普勒频率和Μ 个目标各自的径向速度； 步骤5，根据所述Μ个目标各自的多普勒频率，对雷达接收端接收的回波信号y(t)进行 多普勒频率补偿，计算得到多普勒频率补偿后Μ个目标的回波信号，然后对多普勒频率补偿 后Μ个目标的回波信号进行脉冲压缩处理，计算得到Μ个目标各自的斜距。2. 如权利要求1所述的一种基于雷达大脉宽信号的目标测速测距方法，其特征在于，在 步骤1中，所述Ν个大脉宽的脉冲信号s(t)，其具体表达式为：其中，le{l，2，…，N}，N为雷达发射端发射大脉宽的个数，Tr为脉冲信号s(t)的脉冲重 复周期，τ为脉冲信号s(t)的脉冲宽度，fo为脉冲信号s(t)的中屯、频率，θ〇为脉冲信号s(t)的 初始相位，rect( ·)为矩形函数，t为时间变量，N为自然数。3. 如权利要求1所述的一种基于雷达大脉宽信号的目标测速测距方法，其特征在于，在 步骤2中，所述雷达接收端接收的回波信号y(t)，其表达式为：其中，雷达接收端接收的回波信号y(t)中包含Μ个目标的信息，每一个目标的信息包含 目标的斜距和目标的径向速度，第i个目标的斜距为Ri，第i个目标的径向速度为vri，i e {1， 2,…，Μ}，ki为第i个目标的回波衰减，s(t)为雷达发射端发射的Ν个大脉宽的脉冲信号，t为 时间变量，Ati为第i个目标的回波对应的延迟时间：，(：为光速，fdi为第i个目标的 多普勒频率，I为脉冲信号S (t)的波长。4. 如权利要求1所述的一种基于雷达大脉宽信号的目标测速测距方法，其特征在于，在 步骤3中，所述回波信号y (t)的相关矩阵R，其计算过程为： 3.1初始化:构造一个KXK对角矩阵Ρ,ΚΧΚ对角矩阵P对角线上的元素为等间隔取的K 个频率点各自对应的回波功率，目化个频率点各自对应的回波功率，与Κ X Κ对角矩阵Ρ对角 线上的元素分别相同;k表示第k个频率点，k={l，2,……，Κ}，ΚΧΚ对角矩阵P的初始矩阵为 Pi，初始矩阵Pi对角线上包含Κ个元素值，其中第k个元素值为为，为第k个频率点的对应多普勒频率组成的阵列 流行向量，y (η)为回波信号y(t)的离散形式；回波信号y(t)的相关矩阵初始值为Ri，良示定义;设定迭代总次数为D，m表示第m次迭代，m初值为1; 3.2计算得到第m次迭代后第k个频率点对应的信号波形a(fdk)为第k个频率点的对应多普勒频率组成的阵列流行向量，fdk为第k个频率点的对应多 普勒频率，y(n)为回波信号y(t)的离散形式，上标Η为共辆转置，上标-1为求逆操作； 3.3根据第m次迭代后第k个频率点对应的信号波形，计算得到第m次迭代后第k个 频率点对应的回波功率值3.4令k加1，返回子步骤3.2，直到得到第m次迭代后第K个频率点对应的回波功率值Ρλ， 并根据此时得到的第m次迭代后第1个频率点对应的回波功率值Pmi到第m次迭代后第Κ个频 率点对应的回波功率值Ι\κ，得到第m次迭代后的KXK对角矩阵Pm，然后令m加1，返回子步骤 3.2，直到得到第D次迭代后的K X K对角矩阵时，并将根据第D次迭代后得至化X K对角矩阵时， 计算得到第D次迭代后得到的回波信号y(t)的相关矩阵化，并将第D次迭代后得到的回波信 号y(t)的相关矩阵化，作为回波信号y(t)的相关矩阵R。5.如权利要求1所述的一种基于雷达大脉宽信号的目标测速测距方法，其特征在于，在 步骤4中，所述Μ个目标各自的多普勒频率和Μ个目标各自的径向速度，其计算过程为： 4.1对回波信号y(t)的相关函数R进行特征值分解，即R=USUH，得到F个特征值W及F个 特征值对应的特征向量U，Σ为F个特征值对应的对角矩阵；对F个特征值进行从大到小排 序，得到： λ? ^ λ2 ^ · · · ^ Αμ+1 ^ · · · λρ VI, V2, , VM,VM+1 , ,VF 其中，F个特征值中包含Μ个大的特征值和F-M个接近零的小特征值，Μ个大的特征值为λι 到Am，F-M个接近零的小特征值为λΜ+1到λρ;ν?到VF为大到小排序后F个特征值各自对应的特 征向量； 4.2将F-M个接近零的小特征值分别对应的特征向量构成回波信号y(t)的噪声子空间 G，G = span{vM+i,VM+2,... ,vf}; 4.3使用Κ个频率点对应的多普率频率构成的阵列流行矩阵A(fd)对回波信号y(t)的噪 声子空间G作投影，计算得到回波信号y(t)的功率谱S(fd回波 信号y(t)的功率谱S(fd)的谱峰处对应的多普勒频率向量，即为Μ个大的特征值对应的多普 勒频率向量，进而依次得到Μ个目标各自的多普勒频率和Μ个目标各自的径向速度;其中第i 个大的特征值对应的第i个目标的多普勒频率fdi，并根据第i个目标的多普勒频率fdi，计算 得到第i个目标的径向速度vr、〇为脉冲信号s(t)波长i = {l，2，一，M}，K = F。6. 如权利要求1所述的一种基于雷达大脉宽信号的目标测速测距方法，其特征在于，在 步骤5中，所述多普勒频率补偿后的回波信号，其计算过程为:设定Μ个目标各自的多普勒频 率补偿因子，然后根据所述Μ个目标各自的多普勒频率，对雷达接收端接收的回波信号y(t) 进行多普勒频率补偿，计算得到多普勒频率补偿后Μ个目标的回波信号；在设定的Μ个目标 各自的多普勒频率补偿因子中，第i个目标的多普勒频率补偿因子为了/*'，并根据雷达接 收端接收的回波信号y(t)，计算得到多普勒频率补偿后第i个目标的回波信号 如)><€~^2^'"''，进而计算得到多普勒频率补偿后1个目标的回波信号￥4)，其表达式为：其中，ki为第i个目标的回波衰减，Ati为第i个目标的回波对应的延迟时间，s(t)为雷 达发射端发射的N个大脉宽的脉冲信号，t为时间变量，fdi为第i个目标的多普勒频率，W为 第111 '个目标的回波衰减，上标^表示共辆。7. 如权利要求6所述的一种基于雷达大脉宽信号的目标测速测距方法，其特征在于，所 述对多普勒频率补偿后Μ个目标的回波信号进行脉冲压缩处理，具体为： 对多普勒频率补偿后Μ个目标的回波信号Y(t)进行脉冲压缩处理，计算得到脉冲压缩 处理后Μ个目标的回波信号B(t)，其表达式为：其中为与Μ个目标各自的斜距有关的分量，C为光速，R为Μ个目标各自的斜距集 合，Φ为多普勒频率补偿后Μ个目标的回波信号Y(t)中多普勒频率没有被补偿完的部分;ki 为第i个目标的回波衰减，Ati为第i个目标的回波对应的延迟时间，s(t)为雷达发射端发 射的N个大脉宽的脉冲信号，t为时间变量，fdi为第i个目标的多普勒频率，W为第m'个目标 的回波衰减，上标^表示共辆。
【文档编号】G01S13/58GK106093927SQ201610369076
【公开日】2016年11月9日
【申请日】2016年5月30日
【发明人】廖桂生, 王英, 杨志伟, 黄鹏辉, 刘永军, 许华健, 徐雪菲, 郭学敏, 袁赛
【申请人】西安电子科技大学