脉冲体制1-比特雷达非线性目标重构问题降维方法.pdf

《脉冲体制1-比特雷达非线性目标重构问题降维方法.pdf》由会员分享，可在线阅读，更多相关《脉冲体制1-比特雷达非线性目标重构问题降维方法.pdf（13页完成版）》请在专利查询网上搜索。

1、(19)中华人民共和国国家知识产权局 (12)发明专利申请 (10)申请公布号 (43)申请公布日 (21)申请号 201910284211.5 (22)申请日 2019.04.10 (71)申请人 南京航空航天大学 地址 211106 江苏省南京市江宁区将军大 道29号 (72)发明人 晋本周李建峰李泽昊 (74)专利代理机构 南京苏高专利商标事务所 (普通合伙) 32204 代理人 王安琪 (51)Int.Cl. G01S 7/41(2006.01) (54)发明名称 一种脉冲体制1-比特雷达非线性目标重构 问题降维方法 (57)摘要 本发明公开了一种脉冲体制1-比特雷达非 线性目标重构问。

2、题降维方法， 包括如下步骤： 接 收端， 在模拟域对每个脉冲的接收回波分别进行 移相， 实现对目标多普勒的搬移； 对单帧采样后 得到的1-比特数据立方体进行线性处理， 分别在 三个域完成对目标信号能量的积累， 得到线性处 理后的数据立方体； 对线性处理后的数据立方 体， 保留目标所在多普勒区间内的数据， 完成非 匹配谐波抑制； 对非匹配谐波抑制后的数据立方 体进行三维恒虚警预检测， 得到预检测目标点迹 信息； 基于预检测点迹所在的多普勒、 空域和距 离单元编号， 计算降维观测矩阵， 代替原始观测 矩阵， 获得1-比特雷达降维观测模型。 本发明解 决由于信号观测模型维度过高常规非线性重构 算法无。

3、法求解的问题。 权利要求书2页 说明书7页 附图3页 CN 110007284 A 2019.07.12 CN 110007284 A 1.一种脉冲体制1-比特雷达非线性目标重构问题降维方法， 其特征在于， 包括如下步 骤： (1)接收端， 在模拟域对每个脉冲的接收回波分别进行移相， 实现对目标多普勒的搬 移， 改变目标的多普勒区间； (2)对单帧回波采样后得到的1-比特数据立方体进行线性处理， 在慢时间域、 空域和快 时间域三个域完成对目标信号能量的积累， 提升目标信噪比SNR， 得到线性处理后的数据立 方体； (3)对线性处理后的数据立方体， 保留目标所在多普勒区间内的数据， 剔除其它数据。

4、， 实现对非匹配谐波的抑制； (4)对非匹配谐波抑制后的数据立方体进行三维恒虚警CFAR预检测， 得到预检测目标 点迹信息， 包括所在的多普勒单元、 空域单元和距离单元编号； (5)基于预检测点迹所在的多普勒、 空域和距离单元编号， 计算降维观测矩阵， 用降维 观测矩阵代替原始观测矩阵， 获得1-比特雷达降维观测模型。 2.如权利要求1所述的脉冲体制1-比特雷达非线性目标重构问题降维方法， 其特征在 于， 步骤(1)中， 对每个脉冲的接收回波进行移相， 移相值为2 (k-1)fd， max， 其中， 0kK-1， k为当前帧脉冲序号， K为当前帧脉冲数， fd， max为移相前目标多普勒的上界。

5、， 移相后目标多普 勒区间由(-fd， max， fd， max)变为(0， 2fd， max)。 3.如权利要求1所述的脉冲体制1-比特雷达非线性目标重构问题降维方法， 其特征在 于， 步骤(1)中， 对每个脉冲的接收回波进行移相， 移相值也可以为-2 (k-1)fd， max， 此时， 移 相后目标多普勒区间由(-fd， max， fd， max)变为(-2fd， max， 0)。 4.如权利要求1所述的脉冲体制1-比特雷达非线性目标重构问题降维方法， 其特征在 于， 步骤(2)中， 线性处理在空域、 快时间域和慢时间域分别对应同时多波束合成、 脉冲压缩 和快速傅里叶变换FFT。 5.如权。

6、利要求1所述的脉冲体制1-比特雷达非线性目标重构问题降维方法， 其特征在 于， 步骤(4)中， 对于线性处理后的数据立方体进行3维CFAR预检测， 检测器选择单元平均法 CA-CFAR、 排序式恒虚警OS-CFAR方法， 检测门限由虚警率确定， 而虚警率则根据信号观测模 型维度需要降低的量级来决定， 通常设置为10-210-6。 6.如权利要求1所述的脉冲体制1-比特雷达非线性目标重构问题降维方法， 其特征在 于， 步骤(5)中， 基于预检测到的点迹信息， 对单帧的接收数据观测模型进行近似； 基于预检 测点迹所在的空域、 多普勒和距离单元编号， 计算降维观测矩阵， 用降维观测矩阵代替原始 观测。

7、矩阵， 进而获得降维观测模型， 具体包括如下步骤： (501)将单帧接收到的数据立方体用向量的形式表示； 对于第k个脉冲， 所对应的接收 数据可用矩阵表示， 其中N表示快时间域采样点数， L表示天线单元数， 将矩阵Yk按 列排列成一个向量1帧内， K个脉冲对应K个观测向量， 可构成矩阵Rr0， .， rK-1， 将矩阵R按列排列成向量r， 即为观测向量； (502)根据CFAR检测器预检测结果构造降维观测矩阵， 假设非匹配谐波抑制后空域波 束数目为Csp， 慢时间域多普勒单元数为Cd， 快时间域距离单元数为Cr， CFAR预检测器检测到 的点迹数目为Ipd， 第ipd个点迹所对应的空域频率、 。

8、多普勒频率和距离单元分别记为 权利要求书 1/2 页 2 CN 110007284 A 2 和由此可以得到降维观测矩阵为： Aa(0)， .， a(ipd)， .， a(Ipd-1)， (1) 其中， 符号表示Kronecker积，和分别为： 其中， Np为1个脉冲宽度持续时间内的采样点数，为 其中， Ts为采样时间间隔， 为调频斜率， Nq为1个脉冲宽度持续时间内的采样点数； (503)建立降维观测模型； 降维前， 观测模型可以表示为： rcsignAx+w， (7) 其中， A表示降维前的观测矩阵， w为噪声向量， 为目标向量， 其 元素x(mdMspMr+mspMr+mr)表示多普勒频率。

9、、 空域频率和距离单元分别为和的 目标复幅度， Md， Msp和Mr分别表示慢时间域、 空域和快时间域所划分的格点数； MdCd， Msp Csp， MrCr， 为了讨论方便， 假设MdCd， MspCsp， MrCrN， 且目标位于格点上， csign() sign(Re()+jsign(Im()， Re()和Im()分别表示取变量的实部和虚部， sign ()表示取变量的符号， 基于降维观测矩阵， 可得如下降维模型： rcsign(Ax+w). (8) 其中，x(ipd)表示目标真实的复幅度， 其多 普勒频率、 空域频率和距离单元分别为和 权利要求书 2/2 页 3 CN 11000728。

10、4 A 3 一种脉冲体制1-比特雷达非线性目标重构问题降维方法 技术领域 0001 本发明涉及雷达探测技术领域， 尤其是一种脉冲体制1-比特雷达非线性目标重构 问题降维方法。 背景技术 0002 随着超高精度、 超低速度等新的雷达探测需求的出现， 宽带化、 全数字化成为重要 的发展趋势， 进而对模数转换器(ADC)采样率和数据传输、 存储、 处理能力提出了更高要求。 这一要求与器件水平不足、 数据传输与存储等资源受限之间的矛盾成为常规雷达系统面临 的挑战。 例如， 谷歌 “Soli” 项目， 将毫米波雷达搭载于手表等设备， 利用雷达捕捉人的手势， 对机器进行隔空操作。 为了实现超高距离精度， 。

11、所需探测带宽可达数GHz甚至更大， 同时要 求小体积和低功耗。 但是， 当前应用较多的主流模数转换器(ADC)采样率通常小于3GHz， 更 高速率的ADC短时间内尚难以广泛应用。 另一方面， 采样率1GHz时， ADC芯片功耗大。 比如， ADI公司AD9625-2000芯片， 采样率为2GHz、 12比特量化， 其功耗达3.48W。 0003 1-比特雷达针对上述挑战提供了新的解决途径。 与常规采用高精度量化ADC的雷 达不同， 1-比特雷达基于1-比特ADC对接收信号进行采样和量化， 1-比特ADC采样得到的信 号也称之为1-比特感知信号。 与常规雷达相比， 在二者均以奈奎斯特采样率进行采。

12、样的条 件下， 1-比特雷达可以大幅降低阵面数据产生速率， 进而降低数据传输、 存储能力需求。 另 一方面， 1-比特ADC仅需一个比较器即可实现， 在易于实现超高速率采样的同时保持低功 耗。 0004 然而， 由于1-比特ADC是一个高度非线性器件， 基于1-比特感知信号的目标检测面 临挑战。 1-比特量化使得各散射中心回波产生复杂耦合关系， 量化后的信号除包含各散射 中心自身分量外， 还将产生自身高阶谐波分量和各散射中心交叉谐波分量。 高阶谐波即包 含与脉压参考函数相匹配的谐波， 也包含非匹配谐波。 传统线性信号处理方法， 如匹配滤波 等， 难以完全抑制复杂谐波分量， 需要采用非线性处理方。

13、法。 但是在雷达中， 由于雷达观测 模型中观测矩阵维度巨大， 非线性处理方法难以求解。 以全数字化阵列雷达为例， 假设天线 单元数为100， 单帧脉冲个数为25， 快时间域采样点数为400， 目标检测时需要联合处理的采 样点数为10025400106。 进一步， 假设待重构目标向量维度与采样点数相同， 那么， 对 应的信号模型中观测矩阵维度为106106。 如果观测矩阵每个复元素用4个字节表示， 仅观 测矩阵的所需存储空间就高达3T Byte以上， 现有硬件实时存储、 读取和处理能力难以满足 需求。 因此， 高维度观测模型的降维是非线性处理方法在1-比特雷达中应用的关键。 发明内容 0005 。

14、本发明所要解决的技术问题在于， 提供一种脉冲体制1-比特雷达非线性目标重构 问题降维方法， 能够大幅降低信号观测型的维度， 解决由于观测模型维度过高常规非线性 重构算法无法求解的问题。 0006 为解决上述技术问题， 本发明提供一种脉冲体制1-比特雷达非线性目标重构问题 说明书 1/7 页 4 CN 110007284 A 4 降维方法， 包括如下步骤： 0007 (1)接收端， 在模拟域对每个脉冲的接收回波分别进行移相， 实现对目标多普勒的 搬移， 改变目标的多普勒区间； 0008 (2)对单帧回波采样后得到的1-比特数据立方体进行线性处理， 在慢时间域、 空域 和快时间域三个域完成对目标信。

15、号能量的积累， 提升目标信噪比SNR， 得到线性处理后的数 据立方体； 0009 (3)对线性处理后的数据立方体， 保留目标所在多普勒区间内的数据， 剔除其它数 据， 实现对非匹配谐波的抑制； 0010 (4)对非匹配谐波抑制后的数据立方体进行三维恒虚警CFAR预检测， 得到预检测 目标点迹信息， 包括所在的多普勒单元、 空域单元和距离单元编号； 0011 (5)基于预检测点迹所在的多普勒、 空域和距离单元编号， 计算降维观测矩阵， 用 降维观测矩阵代替原始观测矩阵， 获得1-比特雷达降维观测模型。 0012 优选的， 步骤(1)中， 对每个脉冲的接收回波进行移相， 移相值为2 (k-1)fd。

16、， max， 其 中， 0kK-1， k为当前帧脉冲序号， K为当前帧脉冲数， fd， max为移相前目标多普勒的上界， 移相后目标多普勒区间由(-fd， max， fd， max)变为(0， 2fd， max)。 0013 优选的， 步骤(1)中， 对每个脉冲的接收回波进行移相， 移相值也可以为-2 (k-1) fd， max， 此时， 移相后目标多普勒区间由(-fd， max， fd， max)变为(-2fd， max， 0)。 0014 优选的， 步骤(2)中， 线性处理在空域、 快时间域和慢时间域分别对应同时多波束 合成、 脉冲压缩和快速傅里叶变换FFT。 0015 优选的， 步骤(4。

17、)中， 对于线性处理后的数据立方体进行3维CFAR预检测， 检测器选 择单元平均法CA-CFAR、 排序式恒虚警OS-CFAR方法， 检测门限由虚警率确定， 而虚警率则根 据信号观测模型维度需要降低的量级来决定， 通常设置为10-210-6。 0016 优选的， 步骤(5)中， 基于预检测到的点迹信息， 对单帧的接收数据观测模型进行 近似； 基于预检测点迹所在的空域、 多普勒和距离单元编号， 计算降维观测矩阵， 用降维观 测矩阵代替原始观测矩阵， 进而获得降维观测模型， 具体包括如下步骤： 0017 (501)将单帧接收到的数据立方体用向量的形式表示； 对于第k个脉冲， 所对应的 接收数据可用。

18、矩阵表示， 其中N表示快时间域采样点数， L表示天线单元数， 将矩 阵Yk按列排列成一个向量1帧内， K个脉冲对应K个观测向量， 可构成矩阵R r0， .， rK-1， 将矩阵R按列排列成向量r， 即为观测向量； 0018 (502)根据CFAR检测器预检测结果构造降维观测矩阵， 假设非匹配谐波抑制后空 域波束数目为Csp， 慢时间域多普勒单元数为Cd， 快时间域距离单元数为Cr， CFAR预检测器检 测到的点迹数目为Ipd， 第ipd个点迹所对应的空域频率、 多普勒频率和距离单元分别记为 和由此可以得到降维观测矩阵为： 0019 Aa(0)， .， a(ipd)， .， a(Ipd-1)， 。

19、(1) 0020 其中， 0021 0022符号表示Kronecker积，和分别为： 说明书 2/7 页 5 CN 110007284 A 5 0023 0024 0025 0026其中， Np为1个脉冲宽度持续时间内的采样点数，为 0027 0028 其中， Ts为采样时间间隔， 为调频斜率， Nq为1个脉冲宽度持续时间内的采样点数； 0029 (503)建立降维观测模型； 降维前， 观测模型可以表示为： 0030 rcsignAx+w， (7) 0031 其中， A表示降维前的观测矩阵， w为噪声向量， 为目标向量， 其 元素x(mdMspMr+mspMr+mr)表示多普勒频率、 空域频率。

20、和距离单元分别为和的 目标复幅度， Md， Msp和Mr分别表示慢时间域、 空域和快时间域所划分的格点数； MdCd， Msp Csp， MrCr， 为了讨论方便， 假设MdCd， MspCsp， MrCrN， 且目标位于格点上， csign() sign(Re()+jsign(Im()， Re()和Im()分别表示取变量的实部和虚部， sign ()表示取变量的符号。 基于降维观测矩阵， 可得如下降维模型： 0032 rcsign(Ax+w). (8) 0033其中，x(ipd)表示目标真实的复幅 度， 其多普勒频率、 空域频率和距离单元分别为和 0034 本发明的有益效果为： 本发明通过分。

21、别在模拟域对目标进行多普勒频移、 线性信 号处理、 非匹配谐波抑制和恒虚警预检测， 提出一种脉冲体制1-比特雷达信号观测模型降 维方法， 根据实际对降维程度的需求， 可选择相应的CFAR预检测门限， 进而调整预检测器获 得的预检测目标点迹数目； 通常， 可将CFAR预检测器的虚警率可设置为10-210-6， 对应信 号观测模型中观测矩阵列的维度可降低102106量级， 进而解决由于观测矩阵维度过高常 规非线性重构算法无法求解的问题。 附图说明 0035 图1为本发明的1-比特雷达接收数据排列为数据立方体的示意图。 0036 图2为本发明的降维处理过程示意图。 0037 图3(a)为常规线性处理。

22、后目标所在波束的距离多普勒谱示意图。 0038 图3(b)为本发明移相后进行线性处理所获得的目标所在波束距离多普勒谱示意 图。 0039 图4为本发明基于降维模型进行非线性目标重构后的仿真结果示意图。 0040 图5为本发明基于非线性目标重构结果恢复的距离多普勒谱示意图。 说明书 3/7 页 6 CN 110007284 A 6 具体实施方式 0041 一种脉冲体制1-比特雷达非线性目标重构问题降维方法其特征在于， 包括以下步 骤： 0042 (1)接收端， 在模拟域对每个脉冲的接收回波分别进行移相， 实现对目标多普勒的 搬移， 改变目标的多普勒区间； 0043 (2)对单帧回波采样后得到的1。

23、-比特数据立方体进行线性处理， 在慢时间域、 空域 和快时间域三个域完成对目标信号能量的积累， 提升目标信噪比(SNR)， 得到线性处理后的 数据立方体； 0044 (3)对线性处理后的数据立方体， 保留目标所在多普勒区间内的数据， 剔除其它数 据， 实现对非匹配谐波的抑制； 0045 (4)对非匹配谐波抑制后的数据立方体进行三维恒虚警(CFAR)预检测， 得到预检 测目标点迹信息， 包括所在的多普勒单元、 空域单元和距离单元编号； 0046 (5)基于预检测点迹所在的多普勒、 空域和距离单元编号， 计算降维观测矩阵， 用 降维观测矩阵代替原始观测矩阵， 获得1-比特雷达降维观测模型。 004。

24、7 图1是本发明脉冲体制1-比特雷达单帧接收数据排列为数据立方体的示意图。 每 个天线单元的接收信号首先在射频端进行移相， 改变目标的多普勒频率。 经射频前端后， 采 用1-比特ADC分别对同相和正交两路信号进行采样和1-比特量化。 0048 对于第k个脉冲回波， 移相值为2 (k-1)fd， max， 其中fd， max2vmax/ 为目标多普勒频 率的上界， vmax为目标最大速度， 为载频波长。 移相后， 目标的多普勒区间由(-fd， max， fd， max) 变为(0， 2fd， max)。 为了将真实目标多普勒区间与非匹配谐波对应的多普勒区间分离， 要求脉 冲重复频率PRF满足PR。

25、F8fd， max。 当PRF8fd， max时， 非匹配谐波所在的多普勒区间为- 4fd， max， 02fd， max， 4fd， max。 0049 以一均匀线阵为例进行讨论， 假设天线单元数为L。 单帧内的采样点， 即单帧观测 数据， 可以排列成一个数据立方体。 三个维度分别为空域、 慢时间和快时间域。 不考虑距离 多普勒耦合和距离走动， 数据立方体的任一元素r(k， l， n)可以表示为： 0050 0051 其中， k、 l、 n、 p分别表示脉冲、 天线单元、 距离单元和目标序号， P为目标个数。 p和 p分别为第p个目标复幅度和延时， p2(R0， p-kTIvp)/c2R0，。

26、 p/c， c为光速， R0， p为当前帧目 标起始距离。 fd， p和fsp， p分别表示第p个目标的多普勒频率和空域频率。 fd， p2vp/ ， vp为目 标径向速度， 为载频波长。 fsp， pdsinp/ ， d为天线单元间距， p为目标方位角。 为调频 斜率， Ts为采样间隔， Tp为脉冲宽度， TI为脉冲间隔， w(k， l， n)为噪声。 csign()sign(Re ()+jsign(Im()， Re()和Im()分别表示取变量的实部和虚部， sign()表示取 变量的符号。 0052 在空域， 假设空域频率区间0， 1被分为Msp个格点， 且目标位于格点上。 对于第msp 。

27、说明书 4/7 页 7 CN 110007284 A 7 个格点， 其对应的空域频率为0mspMsp-1。 通常， 格点数Msp满足且MspL。 构建一个矩阵其第msp列为： 0053 0054 对于特定的脉冲k和距离单元n， 空域观测向量可以表示为： 0055 rsp(k， n)csignAspxsp(k， n)+wsp(k， n)， (11) 0056其中， 0057 类似地， 慢时间域观测向量可以表示为： 0058 rd(l， n)csignAdxd(l， n)+wd(l， n)， (12) 0059其中，Md为慢时间域的 格点数， 即多普勒区间(0， 2fd， max)被分为Md个格点。

28、， PRF为脉冲重复频率。 同样， 假设目标位 于格点上， 对于第md个格点， 其目标对应的多普勒频率为0mdMd-1。 通常， 格点数Md满足(PRF/Md)fd， 即MdK， 其中fd1/(KTI)表示多普勒分辨率。 观测 矩阵的第md列为： 0060 0061 在快时间域， 距离区间cTp/2， cTI/2被划分成Mr个格点， 通常MrN， 为了讨论方 便， 假设MrN。 对于第mr个格点， 其对应的距离为(mrcTs+cTp)/2。 快时间域观测向量为： 0062 rr(k， l)csignArxr(k， l)+wr(k， l)， (14) 0063观测矩阵的 第mr列为： 0064 。

29、0065其中， Np为1个脉冲宽度持续时间内的采样点数，为 0066 0067 联立式(11)、 (12)和(14)， 可将单帧接收数据立方体用如下形式表示： 0068 rcsignAx+w， (17) 0069其中， 其元素x(mdMspMr+mspMr+mr)表示多普勒频率、 空域频率和距离单元分别为和mr 的目标复幅度。通常， 目标在三维空间中是稀疏的， 向量 x的非零元素数目满足|x|0PMdMspMr。 0070 式(17)中矩阵A和向量x的维度非常高， 常规非线性重构算法无法求解。 本发明提 供了一种降低矩阵A和向量x维度的方法。 0071 图2是本发明降维处理过程， 包括如下步骤。

30、： 0072 S501、 每个天线单元的接收信号在射频端进行移相， 改变目标的多普勒。 对于第k 说明书 5/7 页 8 CN 110007284 A 8 个脉冲回波， 移相值为2 (k-1)fd， max， 其中fd， max2vmax/ 为目标多普勒频率的上界， vmax为 目标最大速度， 为载频波长。 移相后， 目标的多普勒区间由(-fd， max， fd， max)变为(0， 2fd， max)。 为了将真实目标多普勒区间与谐波对应的多普勒区间分离， 要求脉冲重复频率PRF满足PRF 8fd， max。 当PRF8fd， max时， 非匹配谐波所在的多普勒区间为-4fd， max， 0。

31、2fd， max， 4fd， max。 0073 对单帧回波采样后得到的1-比特数据立方体进行线性处理， 在慢时间域、 空域和 快时间域三个域完成对目标信号能量的积累， 提升目标信噪比(SNR)， 得到线性处理后的数 据立方体。 0074 图3(a)为常规线性处理后得到的目标所在波束距离多普勒谱。 图3(b)为移相后进 行线性处理， 所获得的目标所在波束距离多普勒谱。 仿真中， 加入2个目标， 归一化多普勒分 别为-0.06和0.1， 且均位于主瓣， ADC端两个目标信噪比(SNR)均为-3dB。 线性处理后， 取主 瓣所对应的波束数据， 两种处理方法所获得的距离多普勒谱分别如图3(a)和图3。

32、(b)所示。 从仿真结果可以看出， 通过移相， 可以将非匹配谐波多普勒区间与真实目标、 匹配谐波分 离。 0075 S502、 对线性处理后的数据立方体， 保留多普勒区间(0， 2fd， max)内的数据， 剔除其 它数据， 实现对非匹配谐波的抑制； 抑制后， 假设数据立方体的多普勒、 空域和距离单元数 目分别为记为Cd， Csp和Cr。 0076 S503、 非匹配谐波抑制后的数据立方体， 进行三维CFAR预检测。 参考单元、 保护单 元数目可根据常规CFAR方法选择， 检测门限基于给定的虚警率PFA确定。 PFA则根据降维需求 确定， 例如， 假设维度需要降低103量级， 则令PFA10-。

33、3， 基于该虚警率选择检测门限。 通常， 可根据实际需求， 令PFA10-210-6。 对于每个过检测门限的点迹， 记录其所对应的慢时间 域、 空域和快时间域单元编号。 0077 S504、 基于预检测点迹所在的多普勒、 空域和距离单元编号， 计算降维观测矩阵， 用降维观测矩阵代替原始观测矩阵， 获得1-比特雷达降维观测模型。 包括如下过程： 0078 首先， 将单帧接收到的数据立方体用向量的形式表示。 对于第k个脉冲， 所对应的 接收数据可用矩阵表示， 其中N表示快时间采样点数， L表示天线单元数， 将矩阵Yk 按列排列成一个向量1帧内， K个脉冲对应K个观测向量可构成矩阵Rr0， .， r。

34、K-1， 将矩阵R按列排列成向量r， 即为观测向量； 0079 接着， 构造降维观测矩阵。 非匹配谐波抑制后， 每一个单元进行CFAR预检测， 检测 到的点迹数目为Ipd， 第ipd个点迹所对应的空域频率、 多普勒频率和距离单元分别记为 和由此可以得到降维观测矩阵为： 0080 Aa(0)， .， a(ipd)， .， a(Ipd-1)， (18) 0081 其中， 0082 0083符号 表示Kronecker积，和分别为： 0084 说明书 6/7 页 9 CN 110007284 A 9 0085 0086 0087其中， Np为1个脉冲宽度持续时间内的采样点数，为 0088 0089 。

35、其中， Ts为采样时间间隔， 为调频斜率。 0090 最后， 建立降维观测模型： 0091 rcsign(Ax+w). (23) 0092其中，为降维后的待重构目标向量。 0093 图4是基于公式(23)所述降维观测模型， 采用近似消息传递算法(GAMP)算法重构 结果。 仿真中， 假设回波信号包含两个目标， 信噪比(SNR)相同， 且均为-3dB。 空域频率相同， 均为0， 且均位于主瓣。 归一化多普勒频率分别为-0.06和0.08。 天线单元数L100， 脉冲数K 50， 距离单元数N350， 脉内采样点数Np100， CFAR预检测器检测门限12dB。 0094 从图4可以看出， 该仿真。

36、中， CFAR后过检测门限点迹数目为22， 对应降维后观测矩 阵A的列数为22。 而降维前观测矩阵A的列数为100503501.75106， 因此， 观测模型 维度得到了大幅度降低。 从重构结果看， 真实目标信号得到保留， 谐波导致的虚假目标的强 度得到有效抑制。 0095 图5进一步对基于降维模型的重构效果进行了分析。 基于降维模型进行重构， 可以 得到重构向量进一步地， 可以基于恢复观测信号， 即将 重新排列为如图1 所述数据立方体， 然后进行常规的线性处理， 取真实目标所在空域单元数据， 可得到基于重 构结果恢复的距离多普勒谱， 如图5所示。 与图3(a)相比， 真实目标得到了恢复， 同时， 谐波 分量得到了有效抑制。 0096 本实施例公开的脉冲体制1-比特雷达非线性目标重构问题降维方法， 有效降低了 信号观测模型的维度， 使得将非线性重构算法应用于脉冲体制1-比特雷达目标重构问题成 为可能。 说明书 7/7 页 10 CN 110007284 A 10 图1 图2 说明书附图 1/3 页 11 CN 110007284 A 11 图3(a) 图3(b) 说明书附图 2/3 页 12 CN 110007284 A 12 图4 图5 说明书附图 3/3 页 13 CN 110007284 A 13 。