太赫兹ISAR成像系统及基于后向投影的图像重建方法.pdf

(19)中华人民共和国国家知识产权局(12)发明专利申请(10)申请公布号CN108983234A(43)申请公布日2018.12.11(21)申请号201810691237.7(22)申请日2018.06.28(71)申请人山西大学地址030006山西省太原市小店区坞城路92号(72)发明人梁美彦张存林(74)专利代理机构太原市科瑞达专利代理有限公司14101代理人卢茂春(51)Int.Cl.G01S13/90(2006.01)G01S7/41(2006.01)权利要求书3页说明书6页附图4页(54)发明名称太赫兹ISAR成像系统及基于后向投影的图像重建方法(57)摘要太赫兹ISAR成像系统及基于后向投影的图像重建方法，所述步进太赫兹ISAR系统，包括太赫兹频率步进逆合成孔径雷达成像系统、二维高精度电控旋转台、步进电机控制器；步进电机控制器控制二维高精度电控旋转台旋转；所述太赫兹频率步进雷达成像系统包括跳频模块、射频前端、中频模块、喇叭天线模块、信号采集与处理模块；所述二维高精度电控旋转台包括超精密轴系、精密涡轮、蜗杆传动、光栅度盘、液晶屏以及标配计算机通讯接口；所述步进电机控制器包括控制单元、电机驱动器、稳压电路。CN108983234ACN108983234A权利要求书1/3页1.一种太赫兹ISAR成像系统，其特征是包括太赫兹频率步进雷达成像系统、二维高精度电控旋转台、步进电机控制器；在步进电机控制器的控制下，二维高精度电控旋转台旋转，实现二维扫描；所述太赫兹频率步进雷达成像系统包括频综模块、射频前端、中频模块、喇叭天线模块、信号采集处理模块、上位机控制单元、图像显示单元。2.基于后向投影的图像重建方法，其特征是利用权利要求1所述的太赫兹ISAR成像系统，操作下述步骤：(1)太赫兹雷达雷达系统工作时发射214GHz-226GHz频率步进信号s(t)，可以表示为其中，rect(·)表示矩形包络，τ为脉冲串中子脉冲的宽度，Tr是脉冲重复周期，f0为脉冲的起始频率，Δf为脉冲间频率步进量，每一帧信号中第n个子脉冲的频率fn＝f0+nΔf，其中，n＝0,1…N-1，N为一帧信号中子脉冲的个数。t表示快时间，表示虚数单位。在雷达发送信号的同时，二维高精度电控旋转台在步进电机的驱动下二维步进旋转扫描，二维高精度电控旋转台的步进转角用i来表示，M为ISAR成像时，一幅图像所需信号的帧数，雷达发射信号帧数与转台的步进转角i一一对应，因此，i＝1…M。设在二维高精度电控旋转台旋转到步进转角i时，二维高精度电控旋转台上散射点p到雷达的距离用矢量表示，则该散射点产生的雷达回波时延为：其中，c为光速。对于转台确定的转角i，散射点p产生的雷达回波时延为ti为常数。(2)转台在步进转角i时，位于O'处的天线接收到p点反射的回波信号为：其中，为回波的幅度因子，s(t-ti)为发射信号s(t)经历时延ti后的回波，t仍为快时间。计算二维高精度电控旋转台在步进转角i上，雷达接收到的回波信号为：转台扫描的同时，雷达将接收到的回波信号采样并I/Q解调，二维扫描完成后，得到一个N×M的回波矩阵，该矩阵为复矩阵，矩阵的实部为I路采样的回波幅值，矩阵的虚部为Q路采样的回波幅值。将采样的数据根据公式计算回波矩阵距离方向上，每帧信号子脉冲频点的相位；2CN108983234A权利要求书2/3页然后，计算帧内子脉冲采样点的相位差并将这些相位差求平均，得到帧内平均相位差将多次采样计算出的帧内平均相位差求和再平均得到的值，为理想相位差θ；求出每帧采样信号上各频点间的相位差对理想相位差的偏差作为相位补偿值，叠加在原始信号的相位上，得到相位补偿后的回波信号用表示；调用Matlab中的IFFT命令，将距离向回波数据变换到时域得到一个N×M的二维时域矩阵；在成像平面选取a×a的方形成像区域D，采用Matlab中的x＝[-a/2:Δε:a/2]；y＝[-a/2:Δε:a/2]'；来对二维成像区域D进行网格划分，其中，Δε为网格划分的精度；(3)二维高精度电控旋转台位于步进转角i时，成像区域D网格上的任意像素点q到雷达天线的距离为：计算二维高精度电控旋转台位于步进转角i时，网格上所有的点到雷达天线的距离；其中，rq是像素点q所处位置二维高精度电控旋转台的回转半径，R是雷达天线O'到二维高精度电控旋转台中心O的距离，这个距离为常量，θqi是二维高精度电控旋转台位于步进转角i时，rq与R的夹角；计算对于不同的二维高精度电控旋转台步进转角i上，雷达天线到成像区域D上所有网格的距离，再将步骤2中得到的N×M的时域矩阵反向投影到成像平面网格上，得到M个后向投影的子矩阵，M个投影矩阵相干叠加后，得到太赫兹二维ISAR图像位于处的信号强度是转台在步进转角i处，所接收到的回波信号沿距离向相位补偿并逆傅