一种太赫兹频段的旋翼叶片逆合成孔径雷达成像算法制造技术

本发明专利技术属于太赫兹雷达成像技术，提供一种太赫兹频段旋翼叶片的逆合成孔径成像方法。旋翼叶片转速很快，不再满足距离‑多普勒成像算法在相干积累时间内转动的角度较小的条件。本发明专利技术的技术方案是将逆合成孔径成像矩阵与距离像矩阵，通过系数矩阵联系起来，建立线性方程组；采用短时傅里叶方法计算叶片的旋转周期,并基于Gauss‑Seidel迭代求解方法解线性方程，得到逆合成孔径成像矩阵。本发明专利技术提出了旋翼叶片新的成像方法，为旋翼无人机的成像识别提供了重要方法。

【技术实现步骤摘要】
一种太赫兹频段的旋翼叶片逆合成孔径雷达成像算法
本专利技术属于雷达成像技术，更进一步涉及太赫兹雷达成像技术中的一种基于太赫兹频段的旋翼叶片逆合成孔径雷达成像算法。
技术介绍
在对刚性目标进行成像时，经过运动补偿后，可以采用距离-多普勒算法得到较为清晰的二维像。如果目标不是刚性的，有些部件在旋转运动或机械振动，目标回波的频率会受到微动的调制，导致目标多普勒展宽，出现微多普勒效应。此时，如果仍然用距离-多普勒算法进行成像，得到的将是受到污染的像。对含有旋转部件的目标进行成像，必须要充分认识旋转部件的微多普勒效应，将旋转部件产生的回波提取出来，采取与目标主体不同的成像方法，实现对目标高质量成像。张群等人在其论文“Separationofmicro-DopplerSignalUsingExtendedHoughTransform”中根据旋转部件的微多普勒谱与目标主体谱的不同表现形式，将旋转部件的微多普勒谱与目标主体谱进行分离，获得了具有大旋翼一类目标的清晰像，同时也获得了目标旋转部件的一些运动及结构信息，但该方法无法提取运动幅度小于距离单元的旋转散射点的微动特征，也就无法消除其对目标主体成像的影响。何元桥等人在其论文“具有旋转部件雷达目标ISAR成像新方法”中根据具有旋转部件的运动目标回波信号的特点，提出了以单频信号和正弦调频信号为原子的双信号集匹配分解法，用来提取微多普勒特征，进而可以分离目标主体回波分量和旋转部件回波分量，从而获得聚焦清晰的目标主体像，但该方法把旋转部件回波抑制了，目标像中无法体现旋转部件。西安电子科技大学王敏提出的专利申请“基于调频步进频的太赫兹逆合成孔径雷达成像方法”(申请号：201310416221.2，公开号：CN103454637A)中公开了一种基于调频步进频的太赫兹逆合成孔径雷达成像方法。该方法基于调频步进频信号，采用将脉冲压缩后的信号进行频移合成大宽带技术、以及加窗修正方法，采用距离-多普勒成像算法，获得高分辨率二维逆合成孔径雷达成像图。该方法的不足之处是不能对含有旋转部件的目标进行成像，其成像应用范围受到距离-多普勒方法的限制。白雪茹和保铮等人提出的专利申请“空中微动旋转目标的二维ISAR成像方法”(申请号：201110257606.X，公开号：CN102426360A)中公开了一种空中微动旋转目标的二维ISAR成像方法。该方法采用低调频率匹配滤波的方法将旋转部件回波从总回波中分离出来，对目标主体回波用距离-多普勒算法进行成像，对旋转部件回波用逆Radon变换进行成像，得到的图像聚焦良好、位置估计精度高，克服了自适应chirplet分解成像方法计算量大、耗时多的缺点，也克服了EHT算法得到的图像旁瓣高、估计位置精度低的缺点。该方法具有较高效率的关键是选择合适的能量门限，若能量门限较大，则提取不到微多普勒特征，若能量门限较小，则主体成像模糊。
技术实现思路
本专利技术方法基于旋翼叶片已经经过运动补偿，等效为转台转动模型，然后根据代数重构方法建立距离向与图像矩阵之间的线性方程组，求解方程组就得到图像矩阵。本专利技术的基本思路是：将图像矩阵拉直成向量，构建一个系数矩阵，两者相乘等于距离像从而建立一个线性方程组，然后基于Kaczmarz迭代格式进行求解线性方程组，得到图像矩阵。本专利技术的具体步骤如下：S1、雷达发射线性调频脉冲信号，并接收回波；S2、接收的回波采用解线性调频的脉冲压缩方法；S3、采用短时傅里叶变换估计叶片旋转速度ω；S4、构建系数矩阵；S5、建立线性方程；S6、基于Gauss-Seidel迭代求解方法求解线性方程。优选地，所述步骤S1中，雷达发射信号的中心频率在110GHz～325GHz之间，调频斜率为Kr，脉冲宽度为Tp；其中，当-1/2≤t≤1/2时，rect(t)＝1。优选地，若脉冲信号的重复周期为T，则发射时刻包含慢时间和快时间；其中，τ＝mT,(m＝0,1,…,M-1)为慢时间，M为采样点数；将以发射时刻为起点的时间用t(0＜t＜Tp)表示，其为快时间；快时间用来计算电波传播时间，慢时间用来记录发射脉冲的时刻；快时间和慢时间与全时间的关系为：设时刻τ目标上任意散射中心m到雷达的距离为Rm(τ)，则该散射中心的回波信号可表示为：优选地，所述步骤S2中，选择目标中心点的回波信号作为参考信号，与回波做混频，在快时间域内做傅里叶变换，得到距离压缩后的信号为：上式即为高分辨率距离像(HRRP)序列，对于每一个慢时间的HRRP，它的实包络为sinc函数，且对应的极大值点位于：表明散射点对应距离像上的点随慢时间在各个距离单元上移动。在方位向上，散射点的多普勒相位表现以一个随时间变化的函数：优选地，所述步骤S3中，变换公式为：STFT(t,f)＝∫s(t')w(t'-t)exp(-j2πft')dt'其中，s(t')是信号，w(t)是窗函数。优选地，所述步骤S4中，图像矩阵P的旋转中心位于P的中心，按行把P拉直为一个向量：P＝(p1,p2,…,pJ),J＝N×N在成像时间内发射M个脉冲，构建M行N×N列的系数矩阵W(τ)：其中，I为距离向的插值倍数；m,n＝1,2,…,N；狄拉克函数δ满足：优选地，所述步骤S5中，在距离慢时间域上建立线性方程组：W(τ)·P＝abs(Y(τ))式中，Y(τ)表示τ时刻的距离像；且慢时间τ取值离散，分别为t1,t2,…,tN；W(τ)和Y(τ)则为Wk和Yk。优选地，所述步骤S6中，线性方程组表达式为：wi1p1+wi2p2+…+wiJpJ＝yi式中，i为1,2,…,M；在W(τ)中wii不为0，那么Gauss-Seidel迭代的迭代公式为：其中，k＝0,1,2,…,J。本专利技术的有益效果是在太赫兹频段下，相同方位分辨率下太赫兹频段的积累时间比微波频段更小，而叶片转速较大，超过了距离-多普勒成像算法所需要的角度限制，本专利技术提出了一种太赫兹频段的旋翼叶片逆合成孔径成像算法，克服了在相干积累时间内叶片转过的角度过大，导致无法成像的问题。附图说明图1为本专利技术的流程图；图2为发射信号波形；图3为发射信号频率随时间变化的关系图；图4为图像矩阵的拉直变换示意图。具体实施方式下面结合附图1，对本专利技术具体实施方式作进一步的详细描述：步骤一，如附图2、附图3所示，雷达发射线性调频脉冲，载频fc＝110GHz，脉冲持续时间Tp＝1.0×10-6s，带宽B＝0.5GHz，调频率K＝B/Tp，那么发射线性调频脉冲可以表示为：设时刻τ叶片上任意一点m到雷达的距离为Rm(τ)，则该散射中心的回波信号可表示为：每个脉冲内采样Nsimple＝512个点，采样时刻分别为：t1,t2,…,tNsimple。步骤二，对每一个脉冲回波采用解线性调频的脉冲压缩方法，选择目标中心点的回波信号作为参考信号，与回波做混频，得到的差频信号为：上式中有三个相位项，第一个相位项是快时间t的线性函数，表明混频后的信号是快时间t的单频脉冲，可以认为第一个相位项是包络。由于ΔRm(τ)在单个脉冲期间可认为是固定的，第二、第三个相位项在快时间域处理期间为常数。在慢时间域，ΔRm(τ)随慢时间τ变化，第二个相位项使回波在慢时间域发生多普勒频率偏移，这正是实现横向多普勒分辨率所需要的，而第三个相位项是解线性调频本文档来自技高网...

【技术保护点】
1.一种太赫兹频段的旋翼叶片逆合成孔径成像方法，其特征在于，包括以下步骤：S1、雷达发射线性调频脉冲信号，并接收回波；S2、接收的回波采用解线性调频的脉冲压缩方法；S3、采用短时傅里叶变换估计叶片旋转速度ω；S4、构建系数矩阵；S5、建立线性方程；S6、基于Gauss‑Seidel迭代求解方法求解线性方程。

【技术特征摘要】
1.一种太赫兹频段的旋翼叶片逆合成孔径成像方法，其特征在于，包括以下步骤：S1、雷达发射线性调频脉冲信号，并接收回波；S2、接收的回波采用解线性调频的脉冲压缩方法；S3、采用短时傅里叶变换估计叶片旋转速度ω；S4、构建系数矩阵；S5、建立线性方程；S6、基于Gauss-Seidel迭代求解方法求解线性方程。2.如权利要求1所述的一种太赫兹频段的旋翼叶片逆合成孔径成像方法，其特征在于，所述步骤S1中，雷达发射信号的中心频率在110GHz～325GHz之间，调频斜率为Kr，脉冲宽度为Tp；其中，当-1/2≤t≤1/2时，rect(t)＝1。3.如权利要求2所述的一种太赫兹频段的旋翼叶片逆合成孔径成像方法，其特征在于，若脉冲信号的重复周期为T，则发射时刻包含慢时间和快时间；其中，τ＝mT,(m＝0,1,…,M-1)为慢时间，M为采样点数；将以发射时刻为起点的时间用t(0＜t＜Tp)表示，其为快时间；快时间用来计算电波传播时间，慢时间用来记录发射脉冲的时刻；快时间和慢时间与全时间的关系为：设时刻τ目标上任意散射中心m到雷达的距离为Rm(τ)，则该散射中心的回波信号可表示为：4.如权利要求3所述的一种太赫兹频段的旋翼叶片逆合成孔径成像方法，其特征在于，所述步骤S2中，选择目标中心点的回波信号作为参考信号，与回波做混频，在快时间域内做傅里叶变换，得到距离压缩后的HRRP序列信号为：对于每一个慢时间的HRRP，它的实包络为sinc函数，...

【专利技术属性】
技术研发人员：冯明，贾洁姝，高鹏程，张元，宋诚，花鹏程，
申请(专利权)人：上海无线电设备研究所，
类型：发明
国别省市：上海,31