基于多发多收阵列的太赫兹近场稀疏成像方法.pdf

(19)中华人民共和国国家知识产权局(12)发明专利申请(10)申请公布号CN106772370A(43)申请公布日2017.05.31(21)申请号201611180361.4(22)申请日2016.12.19(71)申请人上海理工大学地址200093上海市杨浦区军工路516号(72)发明人丁丽伍斯璇丁茜叶阳阳王喜旺朱亦鸣(74)专利代理机构上海申汇专利代理有限公司31001代理人吴宝根(51)Int.Cl.G01S13/89(2006.01)权利要求书3页说明书9页附图3页(54)发明名称基于多发多收阵列的太赫兹近场稀疏成像方法(57)摘要本发明涉及一种基于多发多收阵列的太赫兹近场稀疏成像方法，经过成像场景设置、获取太赫兹回波、回波离散化、通过Akaike信息量准则估计出散射点数目、最后通过稀疏重构算法估计目标点坐标，实现多发多收阵列的太赫兹近场稀疏成像。所提方法保留了目标信息在路径延迟中的非线性，利用观测矩阵列原子序号与目标场景之间的一一映射关系，基于稀疏重构技术的贪婪思想实现目标重构。其中，为了改善经典贪婪算法对目标稀疏度（即散射点个数）的依赖性，所提算法在每次迭代选择原子序列时，利用Akaike信息量准则，基于余量对目标的稀疏度进行自适应估计，从而进一步提高了目标重构算法的普适性和实用性。CN106772370ACN106772370A权利要求书1/3页1.一种基于多发多收阵列的太赫兹近场稀疏成像方法，其特征在于，具体包括如下步骤：1)成像场景设置：太赫兹多发多收阵列的M个发射机和N个接收机分布在二维平面的同一基线上，同时发射机、接收机和目标平面三者都位于同一平面XY，以成像场景中心为直角坐标原点O建立坐标系XOY，则第m个发射机和第n个接收机所在位置分别表示为和设目标散射点的直角坐标为r＝(x,y)，散射系数记为σ(r)，发射机和接收机天线组成的天线阵列基线中心到目标平面中心的距离为R0；2)获取太赫兹回波，对回波进行去载波、滤波处理；3)回波离散化：将获取的太赫兹回波根据成像场景离散化，构建观测矩阵，构建的观测矩阵中保留含有目标信息在延迟路径中的非线性，完整保留路径延迟项中的目标信息；4)通过Akaike信息量准则估计出散射点数目；5)通过贪婪算法估计目标点坐标，同时在每次迭代选择原子序列时，利用Akaike信息量准则，基于余量对目标的稀疏度进行自适应估计，实现多发多收阵列的太赫兹近场稀疏成像。2.根据权利要求1所述基于多发多收阵列的太赫兹近场稀疏成像方法，其特征在于，所述步骤2)具体步骤包括如下：第m个发射机发射的信号：其中um(t)是第m个发射机发射波形的复包络，其中m＝1,...,M，其带宽为Bm，fm为第m个发射机的发射载频，发射频率分集信号，此时，任意两个相邻发射载频之间的差值足够大，使发射信号在频谱上占有分离不相交的频带，实现频域正交；反射回波被第n个接收机接收的回波为：其中n＝1,...,N，σ(x,y)是目标散射系数，τm,n(x,y)是第m个发射机到位于(x,y)的目标散射点再到第n个接收机的路径延迟；该回波经去载波处理得到：对回波利用发射信号频率正交性实现通道分离输出的第(n,m)个通道信号为：并通过处理器得到：其中表示卷积操作，上标*表示取共轭操作；2CN106772370A权利要求书2/3页对做傅里叶变换，其频域形式为：其中F{·}表示傅里叶变换；Um(f)＝F{um(t)}，F{um(t)}是复包络um(t)的频谱，设置低通匹配滤波器Hlp,m(f)为：则经过频域滤波输出为：3.根据权利要求2所述基于多发多收阵列的太赫兹近场稀疏成像方法，其特征在于，所述步骤3)具体步骤包括如下：a:将获取的太赫兹回波根据成像场景离散化：将目标平面根据分辨率划分网格，对路径延迟项τm,n(x,y)进行展开：其中(xl,yl)表示第l个强散射点的直角坐标，假设目标成像区间内共有L个待测的散射点，第m个发射天线到第l个散射点的距离为第n个接收天线到第l个散射的的距离为c为光速；将获取的太赫兹回波根据目标平面划分的网格离散化处理得到：其中表示第(n,m)个通道在第k个频率采样点fk处的回波，表示第l个强散射点与第m个发射机和第n个接收机的距离，σl表示第l个强散射点的散射系数；b:构建观测矩阵：收集第(n,m)个通道的K个采样点构成向量为：3CN106772370A权利要求书3/3页其中上标T为转置操作，则收集MN个通道的回波可组成矩阵：其中为MNK×L维的观测矩阵，σ为目标场景对应的L×1维向量，实际目标强散射点的个数为L0，满足L0＜＜L。4.根据权利要求3所述基于多发多收阵列的太赫兹近场稀疏成像方法，其特征在于，所述步骤4)具体步骤包括如下：得到接收回波的协方差矩阵R＝Sr×(Sr)H，其中，上标H表示共