CN107271977B - 基于移动激励源fdtd算法的高精度sar回波仿真方法 - Google Patents

基于移动激励源fdtd算法的高精度sar回波仿真方法 Download PDF

Info

Publication number
CN107271977B
CN107271977B CN201710608975.6A CN201710608975A CN107271977B CN 107271977 B CN107271977 B CN 107271977B CN 201710608975 A CN201710608975 A CN 201710608975A CN 107271977 B CN107271977 B CN 107271977B
Authority
CN
China
Prior art keywords
excitation source
target
point
far
field
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201710608975.6A
Other languages
English (en)
Other versions
CN107271977A (zh
Inventor
邹斌
陆健
张腊梅
史庆武
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Harbin Institute of Technology
Original Assignee
Harbin Institute of Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Harbin Institute of Technology filed Critical Harbin Institute of Technology
Priority to CN201710608975.6A priority Critical patent/CN107271977B/zh
Publication of CN107271977A publication Critical patent/CN107271977A/zh
Application granted granted Critical
Publication of CN107271977B publication Critical patent/CN107271977B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/40—Means for monitoring or calibrating
    • G01S7/4052—Means for monitoring or calibrating by simulation of echoes

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

基于移动激励源FDTD算法的高精度SAR回波仿真方法，属于回波信号的仿真技术领域。本发明是为了解决现有SAR回波仿真系统精度低，不能考虑电磁波与目标作用过程的变化的问题。它包括以下步骤：设计SAR系统的相关参数；计算FDTD仿真的相关参数；对选定采样位置，利用时域有限差分法FDTD对目标进行电磁计算，得到目标的远场电场数据；设置远场观察点，计算目标上所有可视点元胞与远场观察点之间的距离，并把可视点分为不同的等距离带；对目标的远场电场数据进行处理，得到移动激励源的飞行轨迹上所有采样位置点的回波数据，对采样位置点的回波数据进行处理获得目标的SAR原始回波数据。本发明用于SAR的回波仿真。

Description

基于移动激励源FDTD算法的高精度SAR回波仿真方法
技术领域
本发明涉及基于移动激励源FDTD算法的高精度SAR回波仿真方法，属于回波信号的仿真技术领域。
背景技术
合成孔径雷达SAR是一种安装在运动载体上的高分辨率微波成像雷达，能实现全天时、全天候、高分辨率、宽幅对地观测，因此在军事、海洋、农业和林业等诸多领域具有广阔的应用前景。回波信号建模仿真在合成孔径雷达的研制过程中起着十分重要的作用，它一方面能够弥补由于技术和经济等原因造成的SAR系统特性验证性实验较少的不足，另一方面还可以以用于检测合成孔径雷达系统性能，用来验证和评价各种SAR成像处理算法的性能，用于典型目标SAR数据获取以及识别检测等。
传统的SAR仿真系统分辨率较低，在回波仿真过程中忽略了雷达散射截面积在观测时间内随空间角度及信号频率的变化，无法反映目标在不同观测条件下的后向散射特性差异。但是对于高分辨率SAR数据，需要考虑到目标在不同角度来波照射下后向散射特性的差异，尤其对于棱角分明的建筑物和军事目标。同时传统的SAR回波仿真方法，不考虑电磁波与目标作用的具体过程变化，不能在根本上解决散射回波无法完全反映目标的结构材质等特征的问题，很难得到目标精确的细节特性。因此随着星载SAR技术的发展和分辨率的提高，特别是在高分辨率成像模式下，传统SAR回波仿真方法已经难以满足SAR系统的设计和回波成像需求。
发明内容
本发明目的是为了解决现有SAR回波仿真系统精度低，不能考虑电磁波与目标作用过程的变化的问题，提供了一种基于移动激励源FDTD算法的高精度SAR回波仿真方法。
本发明所述基于移动激励源FDTD算法的高精度SAR回波仿真方法，它包括以下步骤：
步骤一：设计SAR系统的相关参数；
步骤二：根据步骤一中的相关参数计算FDTD仿真的相关参数；
步骤三：对移动激励源飞行轨迹上的所有选定采样位置，利用时域有限差分法FDTD对目标进行电磁计算，得到目标的远场电场数据；
步骤四：在移动激励源处设置远场观察点，计算目标上所有可视点元胞与远场观察点之间的距离，并把可视点分为不同的等距离带；
步骤五：根据可视点所处的距离带，对目标的远场电场数据进行处理，得到移动激励源的飞行轨迹上所有采样位置点的回波数据，对采样位置点的回波数据进行处理获得目标的SAR原始回波数据。
本发明的优点：本发明方法利用时域有限差分法(FDTD)精确计算了目标在移动激励源下的电磁场分布状况，并由设置在目标可视点表面的惠更斯面外推出目标局部细节的远场电场值，之后根据目标与激励源的相对位置，将目标点划分为连续的等距离带，对等距离带内点及带与带之间的远场值做不同的处理，最终得到目标的SAR回波数据。
本发明充分利用了FDTD计算的精确性来计算电磁波与目标的相互作用过程，精确地获取目标表面细节散射信息，并且通过设置移动激励源的方式模拟真实卫星工作方式，真实的反映了目标在不同状态下的后向散射特性，得到了精确的SAR回波数据，并能满足各种实际应用需求。
附图说明
图1是本发明所述基于移动激励源FDTD算法的高精度SAR回波仿真方法的流程图；
图2是SAR系统空间几何模型示意图；
图3是目标可视点元胞等距离带划分示意图；
图4是SAR回波数据块叠加处理示意图；
图5是FDTD算法仿真流程图；
图6是FDTD算法计算区域各边界划分示意图；
图7是二面角目标仿真得到的原始回波示意图；
图8是根据二面角回波计算的RCS变化曲线图；
图9是直升飞机模型示意图；
图10是直升飞机的SAR回波成像结果示意图。
具体实施方式
具体实施方式一：下面结合图1至图10说明本实施方式，本实施方式所述基于移动激励源FDTD算法的高精度SAR回波仿真方法，它包括以下步骤：
步骤一：设计SAR系统的相关参数；
步骤二：根据步骤一中的相关参数计算FDTD仿真的相关参数；
步骤三：对移动激励源飞行轨迹上的所有选定采样位置，利用时域有限差分法FDTD对目标进行电磁计算，得到目标的远场电场数据；
步骤四：在移动激励源处设置远场观察点，计算目标上所有可视点元胞与远场观察点之间的距离，并把可视点分为不同的等距离带；
步骤五：根据可视点所处的距离带，对目标的远场电场数据进行处理，得到移动激励源的飞行轨迹上所有采样位置点的回波数据，对采样位置点的回波数据进行处理获得目标的SAR原始回波数据。
步骤一中SAR系统的相关参数为：雷达距离向分辨率ρr、雷达方位向分辨率ρa、天线波束下视角θ、天线波束方位角
Figure BDA0001359060690000031
雷达运动速度V、移动激励源与目标垂直距离R0、移动激励源信号的中心频率f0、调频率μ和脉冲重复频率fPRF。
步骤二中FDTD仿真的相关参数为：移动激励源原始信号表达式Ein(t)、移动激励源的位置、远场观察点位置和能接收到移动激励源照射的可视点位置Pvisual的坐标；同时设计相应的仿真总时间、仿真时间步长Δt和剖分网格大小。
步骤三中得到目标的远场电场数据的具体方法为：针对移动激励源的某个采样位置，计算移动激励源原始信号表达式；利用FDTD方法对目标进行仿真，得到所有可视点的远场电场值数据Efarfield，远场电场值数据Efarfield为M×N的矩阵，M行表示远场电场值在时间上以aΔt为间隔采样M次，其中a为采样时间步数，N列表示远场观察点处的目标可视点Pvisual是N个。这些数据包含了目标被照射区域的细节信息。
步骤五中目标的SAR原始回波数据的获得方法为：根据每个采样位置点与目标中心的斜距R(t)，计算回波延迟时间；对所有采样位置点的回波数据做时间延迟R(t)/2c后，并行排列为一个二维矩阵，即得到目标的SAR原始回波数据
Figure BDA0001359060690000032
其中c为电磁波在真空中传播速度。
计算FDTD仿真的相关参数的具体方法为：
根据雷达距离向分辨率ρr要求计算出移动激励源原始信号的带宽B及脉冲宽度τ：
Figure BDA0001359060690000033
移动激励源原始信号Ein(t)的表达式为：
Figure BDA0001359060690000034
式中A为信号幅度，rect(·)为矩形函数，τ是矩形函数宽度；
根据雷达方位向分辨率ρa的要求计算出SAR的合成孔径长度Lmax：
Figure BDA0001359060690000041
在移动激励源运行轨迹上以fPRF的频率进行采样，得到的采样点作为SAR回波数据获取点，结合移动激励源与目标的几何模型计算出所有采样点的位置；
根据SAR系统空间几何模型，各采样点天线波束方位角
Figure BDA0001359060690000042
变化范围为：
Figure BDA0001359060690000043
公式(4)中的三项分别表示天线波束方位角的初始角度、相邻两点角度间隔和终止角度；
球坐标系中，每个采样点处移动激励源点的球坐标为
Figure BDA0001359060690000044
转换成对应的直角坐标系中的位置坐标为(xin,yin,zin)，具体表达式如下：
Figure BDA0001359060690000045
把目标剖分成立方体元胞形式，根据远场观察点与目标之间的几何关系计算出所有可视点元胞Pvisual(n)的坐标[xvisual(n),yvisual(n),zvisual(n)]；
根据数值稳定性条件，移动激励源信号的中心波长λ0＝c/f0，FDTD计算中，元胞的大小限制为：
Figure BDA0001359060690000046
其中Δx、Δy、Δz分别为立方体元胞的三个边长，δ为设定边长，根据FDTD计算量与元胞体积成反比，设置
Figure BDA0001359060690000047
FDTD计算时间步数为
Figure BDA0001359060690000048
其中仿真时间步长Δt根据Courant稳定性条件计算得
Figure BDA0001359060690000049
计算得到时间步数Num为：
Figure BDA0001359060690000051
步骤三中得到目标的远场电场数据的具体方法为：
A.读取剖分后的目标结构及材质信息，建立空间坐标系使目标位于第一象限，此空间坐标系与SAR系统空间几何模型坐标系一致，根据时域有限差分算法，求得空间每一个元胞上的离散电场值
Figure BDA0001359060690000052
与磁场值
Figure BDA0001359060690000053
Figure BDA0001359060690000054
Figure BDA0001359060690000055
式和中
Figure BDA0001359060690000056
表示位于
Figure BDA0001359060690000057
处第n个时间步长z方向的电场值，其中(i,j,k)表示该位置与坐标原点的距离为(iΔx,jΔy,kΔz)，式中迭代系数CA(m)、CB(m)、CP(m)、CQ(m)是随着计算的场值坐标变化而变化的与场值大小无关的量，m表示的是电磁场节点的坐标，如式(10)所示：
Figure BDA0001359060690000061
Figure BDA0001359060690000062
式(10)中ε(m)为节点m处材质的介电常数，σ(m)为节点m处材质的电导率；
B.如图2所示，移动激励源与目标距离最近点记为B处，此时移动激励源表达式EBin(t)＝Ein(t)，移动激励源传播方向为射线CO，指向目标中心，当移动激励源位于运动轨迹上某一个采样点C处，根据移动激励源的初始位置和移动之后的相对位置可以求得移动激励源位于C点处相对于B点发射信号的时间延迟tC，如式(11)所示：
tC＝f(tB,RB,Δt)\*MERGEFORMAT(11)
其中f(·)是个根据位置求得的已知函数，tB为移动激励源在B点发射电磁波的时刻，RB是激励源在B点处于目标的距离，Δt是移动激励源从B移动到C处电磁波作用在目标上的时间；
进而求得移动激励源位于C点处激励源ECin(t)的表达式，如式(12)所示：
Figure BDA0001359060690000063
C.在FDTD计算区域的总场-散射场边界引入式的移动激励源，并在目标边界设置卷积完善匹配层CPML吸收电磁波，限制计算区域大小，然后根据式和代表的电磁场迭代公式在时间域迭代计算空间各节点处的离散电磁场值；迭代计算Num次；
D.在每个可视点元胞表面设置惠更斯面，根据等效原理计算惠更斯面上的表面电流，然后根据电磁流的辐射公式计算远场观察点处的场值，由此近远场外推方法计算的电场如式(13)所示：
Figure BDA0001359060690000071
Figure BDA0001359060690000072
Figure BDA0001359060690000073
Figure BDA0001359060690000074
其中Eθ,
Figure BDA0001359060690000075
为球坐标系下θ,
Figure BDA0001359060690000076
方向上的点场分量，
Figure BDA0001359060690000077
为计算过程中间矢量，
Figure BDA0001359060690000078
是远场观察点在目标坐标系下的位置矢量，
Figure BDA0001359060690000079
分别表示位置矢量为
Figure BDA00013590606900000710
的点在t时刻的电流密度和磁流密度，
Figure BDA00013590606900000711
为方向是
Figure BDA00013590606900000712
的电位电场矢量；
以aΔt为时间间隔对得到的远场电场值在时间上采样，得到M＝Numa个采样值，采样电场值如式(14)所示：
Figure BDA00013590606900000713
其中δ[t]为单位脉冲函数；可视点的远场电场值数据Efarfield包含了目标被照射区域的细节电磁散射信息。
把N个可视点电场值序列并行排列得到大小M×N的电场值矩阵Efarfield，该数据矩阵记录了目标局部电场变化情况及电磁散射信息。
步骤四计算所有可视点与远场观察点的距离，并把可视点分为不同的等距离带，具体包括以下几个步骤：
A.取远场观察点位置与移动激励源所在位置相同，计算可视点与远场观察点的距离：
Figure BDA00013590606900000714
B.根据Rvisual(n)值的大小，把可视点元胞划分为不同的等距离带，每个等距离带宽为(aΔt)·c，按照与远场观察点的距离从近到远对等距离带编号为1,2,3,...,K，同一个等距离带内的可视点元胞，其在移动激励源作用下的回波同时到达远场观察点，不同等距离带内可视点的回波间隔2k·(aΔt)抵达远场观察点。
步骤五获得目标的SAR原始回波数据的具体方法为：
A.对编号为k的等距离带内的Nk个可视点，他们的远场值数据是一个M×Nk的矩阵，等距离带内可视点回波同时抵达移动激励源，把这些可视点的远场电场值对应相加即矩阵每一行相加，得到大小为M×1的矩阵记作Mk，将Mk看作是等距离带k内目标的回波数据；
B.对编号为k+1的等距离带处理得到另一个M×1的矩阵Mk+1，相比于等距离带k，等距离带k+1上的可视点的回波抵达观察点的时间延迟了2×(aΔt)，把Mk+1的第m行数据与Mk的第m+2行数据相加，即在Mk+1第一行前插两个0构成大小为(M+2)×1新的矩阵M′k+1，然后再与Mk从第一行开始对应相加，Mk矩阵元素为空时补零，得到编号为k、k+1的等距离带内可视点回波数据；
C.同理对编号为k+l的等距离带处理，计算得到大小为M×1的矩阵Mk+l，在Mk+l第一行前插2l个0构成大小为(M+2l)×1的新矩阵M′k+l，然后再与Mk、M′k+1、...、M′k+l-1、M′k+l从第一行开始对应相加，矩阵元素为空时补零，得到编号为k、k+1、...、k+l-1、k+l的等距离带内可视点回波数据；
D.计算出所有等距离带的回波数据M1,M2,M3,...,MK，然后计算出他们插零之后的回波矩阵M1,M′2,M′3,...,M′K，对插零后的矩阵从第一行开始对应相加，矩阵元素为空时补零，最终得到在该采样点处的回波数据SEcho(t)。
本发明方法中，对于选定的所有采样点位置，分别进行计算，获得所有轨道采样点处回波数据
Figure BDA0001359060690000081
根据每个采样点处雷达与目标中心的斜距R(t)，确定回波延迟时间，对回波
Figure BDA0001359060690000082
做时间延迟R(t)/2c后，并行排列为一个二维矩阵如图4下半部分所示，即得到目标在移动激励源下的SAR原始回波数据
Figure BDA0001359060690000083
步骤一中，首先对SAR系统空间几何关系建模，如图2所示，
以目标为坐标原点建立空间坐标系，雷达沿着直线AC朝X轴正方向飞行，雷达与目标垂直距离OB为R0，在B点处雷达天线波束下视角为θ，雷达运动速度为V，激励源为线性调频信号，其中心频率是f0、调频率为μ，雷达信号脉冲重复频率为fPRF，为了得到高分辨率数据设置雷达距离向分辨率ρr、雷达方位向分辨率ρa。
步骤二中，由于合成孔径的长度远大于一般目标的尺寸，所以在SAR系统模型中，雷达的运动轨迹长度一般选择为Lmax，其运动轨迹为直线，方向平行于SAR系统空间几何模型坐标系的X轴。
步骤三中，设置雷达与目标距离最近点B处激励源表达式EBin(t)＝Ein(t)，经过tC时间雷达从B运行到C点，此时与目标的距离是RC为：
Figure BDA0001359060690000091
相对于目标来说，B、C两点发射的电磁波间隔Δt到达目标，如此可以得到tC的方程如下所示：
Figure BDA0001359060690000092
其中RB＝R0，Δt可以从FDTD仿真时间中获得，从而求得雷达位于C点处相对于B点的时间延迟tC。
按照图6所示划分FDTD计算区域的各个边界，FDTD计算流程如图5所示。
式(13)中，
Figure BDA0001359060690000093
步骤四中，如图3所示，根据Rvisual(n)值的大小，把可视点元胞划分为不同的等距离带。
步骤五中，根据可视点所处的距离带，对远场电场值数据处理，如图4所示，得到雷达位于某个采样点处的回波数据
Figure BDA0001359060690000094
具体实施例：为了说明本发明的有效性，进行典型目标仿真实验，仿真参数为：雷达距离向分辨ρr＝1m，方位向分辨率ρa＝1m，雷达与目标垂直距离R0＝600km，天线波束下视角θ＝33°，雷达运动速度V＝7000m/s，激励源信号中心频率f0＝1.3×109Hz，线性调频信号调频率μ＝1014s-2，脉冲重复频率fPRF＝1500Hz，由此计算得到信号带宽B＝1.79×108，脉冲宽度τ＝1.79μs，SAR的合成孔径长度Lmax＝6.9×104m，方位角
Figure BDA0001359060690000101
变化范围[-3.3°:(4.46×10-4)°:3.3°]，YEE元胞边长δ＝1.9×10-2m，时间步长Δt＝3.7×10-11s，迭代时间步数为Num＝4.8×104。
选取一个垂直于XY平面的二面角作为仿真目标，二面角每个面大小为0.5m*0.5m，材质为金属，两面之间夹角为90°，雷达垂直于二面角角平分面方向飞行，得到的回波信号幅值量化之后的结果如图7所示；在0°到360°变化的方位角上，每隔5°做一次SAR回波仿真并计算该方位角下的目标RCS，得到如图8所示的二面角RCS变化图，这与实际二面角的RCS变化曲线一致。为了进一步验证得到的SAR回波，通过本方法对图9所示的直升飞机模型做回波仿真，然后根据RD成像算法得到如图10所示的成像结果，可以清晰看出直升飞机的轮廓信息。
从仿真结果可得，本发明充分利用了FDTD计算的精确性，真实的反映了目标不同状态下后向散射特性，得到了精确的SAR回波数据，初步的成像结果表明本发明不仅可以满足SAR系统特性验证，还可以满足SAR高分辨率成像模式等不同实际应用对高精度仿真信号的需求。

Claims (9)

1.一种基于移动激励源FDTD算法的高精度SAR回波仿真方法，其特征在于，它包括以下步骤：
步骤一：设计SAR系统的相关参数；
步骤二：根据步骤一中的相关参数计算FDTD仿真的相关参数；
步骤三：对移动激励源飞行轨迹上的所有选定采样位置，利用时域有限差分法FDTD对目标进行电磁计算，得到目标的远场电场数据；
步骤四：在移动激励源处设置远场观察点，计算目标上所有可视点元胞与远场观察点之间的距离，并把可视点分为不同的等距离带；
步骤五：根据可视点所处的距离带，对目标的远场电场数据进行处理，得到移动激励源的飞行轨迹上所有采样位置点的回波数据，对采样位置点的回波数据进行处理获得目标的SAR原始回波数据。
2.根据权利要求1所述的基于移动激励源FDTD算法的高精度SAR回波仿真方法，其特征在于，步骤一中SAR系统的相关参数为：雷达距离向分辨率ρr、雷达方位向分辨率ρa、天线波束下视角θ、天线波束方位角
Figure FDA0002331683500000011
雷达运动速度V、移动激励源与目标垂直距离R0、移动激励源信号的中心频率f0、调频率μ和脉冲重复频率fPRF。
3.根据权利要求2所述的基于移动激励源FDTD算法的高精度SAR回波仿真方法，其特征在于，步骤二中FDTD仿真的相关参数为：移动激励源原始信号表达式Ein(t)、移动激励源的位置、远场观察点位置和能接收到移动激励源照射的可视点位置Pvisual的坐标；同时设计相应的仿真总时间、仿真时间步长Δt和剖分网格大小。
4.根据权利要求3所述的基于移动激励源FDTD算法的高精度SAR回波仿真方法，其特征在于，步骤三中得到目标的远场电场数据的具体方法为：针对移动激励源的某个采样位置，计算移动激励源原始信号表达式；利用FDTD方法对目标进行仿真，得到所有可视点的远场电场值数据Efarfield，远场电场值数据Efarfield为M×N的矩阵，M行表示远场电场值在时间上以aΔt为间隔采样M次，其中a为采样时间步数，N列表示远场观察点处的目标可视点Pvisual是N个。
5.根据权利要求4所述的基于移动激励源FDTD算法的高精度SAR回波仿真方法，其特征在于，
步骤五中目标的SAR原始回波数据的获得方法为：根据每个采样位置点与目标中心的斜距R(t)，计算回波延迟时间；对所有采样位置点的回波数据做时间延迟R(t)/2c后，并行排列为一个二维矩阵，即得到目标的SAR原始回波数据
Figure FDA0002331683500000021
其中c为电磁波在真空中传播速度。
6.根据权利要求5所述的基于移动激励源FDTD算法的高精度SAR回波仿真方法，其特征在于，
计算FDTD仿真的相关参数的具体方法为：
根据雷达距离向分辨率ρr要求计算出移动激励源原始信号的带宽B及脉冲宽度τ：
Figure FDA0002331683500000022
移动激励源原始信号Ein(t)的表达式为：
Figure FDA0002331683500000023
式中A为信号幅度，rect(·)为矩形函数；
根据雷达方位向分辨率ρa的要求计算出SAR的合成孔径长度Lmax：
Figure FDA0002331683500000024
在移动激励源运行轨迹上以fPRF的频率进行采样，得到的采样点作为SAR回波数据获取点，结合移动激励源与目标的几何模型计算出所有采样点的位置；
根据SAR系统空间几何模型，各采样点天线波束方位角
Figure FDA0002331683500000025
变化范围为：
Figure FDA0002331683500000026
公式(4)中的三项分别表示天线波束方位角的初始角度、相邻两点角度间隔和终止角度；
球坐标系中，每个采样点处移动激励源点的球坐标为
Figure FDA0002331683500000027
转换成对应的直角坐标系中的位置坐标为(xin,yin,zin)，具体表达式如下：
Figure FDA0002331683500000028
把目标剖分成立方体元胞形式，根据远场观察点与目标之间的几何关系计算出所有可视点元胞Pvisual(n)的坐标[xvisual(n),yvisual(n),zvisual(n)]；
根据数值稳定性条件，移动激励源信号的中心波长λ0＝c/f0，FDTD计算中，元胞的大小限制为：
Figure FDA0002331683500000031
其中Δx、Δy、Δz分别为立方体元胞的三个边长，δ为设定边长，根据FDTD计算量与元胞体积成反比，设置
Figure FDA0002331683500000032
FDTD计算时间步数为
Figure FDA0002331683500000033
其中仿真时间步长Δt根据Courant稳定性条件计算得
Figure FDA0002331683500000034
计算得到时间步数Num为：
Figure FDA0002331683500000035
7.根据权利要求6所述的基于移动激励源FDTD算法的高精度SAR回波仿真方法，其特征在于，步骤三中得到目标的远场电场数据的具体方法为：
A.读取剖分后的目标结构及材质信息，建立空间坐标系使目标位于第一象限，此空间坐标系与SAR系统空间几何模型坐标系一致，根据时域有限差分算法，求得空间每一个元胞上的离散电场值
Figure FDA0002331683500000036
与磁场值
Figure FDA0002331683500000037
Figure FDA0002331683500000038
Figure FDA0002331683500000041
式(8)和(9)中
Figure FDA0002331683500000042
表示位于
Figure FDA0002331683500000043
处第n个时间步长z方向的电场值，其中(i,j,k)表示该位置与坐标原点的距离为(iΔx,jΔy,kΔz)，式中迭代系数CA(m)、CB(m)、CP(m)、CQ(m)是随着计算的场值坐标变化而变化的与场值大小无关的量，m表示的是电磁场节点的坐标，如式(10)所示：
Figure FDA0002331683500000044
式(10)中ε(m)为节点m处材质的介电常数，σ(m)为节点m处材质的电导率；
B.移动激励源与目标距离最近点记为B处，此时移动激励源表达式EBin(t)＝Ein(t)，移动激励源传播方向指向目标中心，当移动激励源位于运动轨迹上某一个采样点C处，根据移动激励源的初始位置和移动之后的相对位置可以求得移动激励源位于C点处相对于B点发射信号的时间延迟tC，如式(11)所示：
tC＝f(tB,RB,Δt) (11)
其中f(·)是个根据位置求得的已知函数，tB为移动激励源在B点发射电磁波的时刻，RB是激励源在B点处于目标的距离，Δt是移动激励源从B移动到C处电磁波作用在目标上的时间；
进而求得移动激励源位于C点处激励源ECin(t)的表达式，如式(12)所示：
Figure FDA0002331683500000051
C.在FDTD计算区域的总场-散射场边界引入式(12)的移动激励源，并在目标边界设置卷积完善匹配层CPML吸收电磁波，限制计算区域大小，然后根据式(8)和(9)代表的电磁场迭代公式在时间域迭代计算空间各节点处的离散电磁场值；
D.在每个可视点元胞表面设置惠更斯面，根据等效原理计算惠更斯面上的表面电流，然后根据电磁流的辐射公式计算远场观察点处的场值，由此近远场外推方法计算的电场如式(13)所示：
Figure FDA0002331683500000052
其中Eθ,
Figure FDA0002331683500000053
为球坐标系下θ,
Figure FDA0002331683500000054
方向上的点场分量，
Figure FDA0002331683500000055
为计算过程中间矢量，
Figure FDA0002331683500000056
是远场观察点在目标坐标系下的位置矢量，
Figure FDA0002331683500000057
分别表示位置矢量为
Figure FDA0002331683500000058
的点在t时刻的电流密度和磁流密度，
Figure FDA0002331683500000059
为方向是
Figure FDA00023316835000000510
的电位电场矢量；
以aΔt为时间间隔对得到的远场电场值在时间上采样，得到M＝Num/a个采样值，采样电场值如式(14)所示：
Figure FDA00023316835000000511
其中δ[t]为单位脉冲函数；
把N个可视点电场值序列并行排列得到大小M×N的电场值矩阵Efarfield，该数据矩阵记录了目标局部电场变化情况。
8.根据权利要求7所述的基于移动激励源FDTD算法的高精度SAR回波仿真方法，其特征在于，所述步骤四具体包括以下几个步骤：
A.取远场观察点位置与移动激励源所在位置相同，计算可视点与远场观察点的距离：
Figure FDA0002331683500000061
B.根据Rvisual(n)值的大小，把可视点元胞划分为不同的等距离带，每个等距离带宽为(aΔt)·c，按照与远场观察点的距离从近到远对等距离带编号为1,2,3,...,K，同一个等距离带内的可视点元胞，其在移动激励源作用下的回波同时到达远场观察点，不同等距离带内可视点的回波间隔2k·(aΔt)抵达远场观察点。
9.根据权利要求8所述的基于移动激励源FDTD算法的高精度SAR回波仿真方法，其特征在于，步骤五获得目标的SAR原始回波数据的具体方法为：
A.对编号为k的等距离带内的Nk个可视点，他们的远场值数据是一个M×Nk的矩阵，等距离带内可视点回波同时抵达移动激励源，把这些可视点的远场电场值对应相加即矩阵每一行相加，得到大小为M×1的矩阵记作Mk，将Mk看作是等距离带k内目标的回波数据；
B.对编号为k+1的等距离带处理得到另一个M×1的矩阵Mk+1，相比于等距离带k，等距离带k+1上的可视点的回波抵达观察点的时间延迟了2×(aΔt)，把Mk+1的第m行数据与Mk的第m+2行数据相加，即在Mk+1第一行前插两个0构成大小为(M+2)×1新的矩阵M′k+1，然后再与Mk从第一行开始对应相加，Mk矩阵元素为空时补零，得到编号为k、k+1的等距离带内可视点回波；
C.同理对编号为k+l的等距离带处理，计算得到大小为M×1的矩阵Mk+l，在Mk+l第一行前插2l个0构成大小为(M+2l)×1的新矩阵M′k+l，然后再与Mk、M′k+1、...、M′k+l-1、M′k+l从第一行开始对应相加，矩阵元素为空时补零，得到编号为k、k+1、...、k+l-1、k+l的等距离带内可视点回波；
D.计算出所有等距离带的回波数据M1,M2,M3,...,MK，然后计算出他们插零之后的回波矩阵M1,M′2,M′3,...,M′K，对插零后的矩阵从第一行开始对应相加，矩阵元素为空时补零，最终得到在该采样点处的回波数据SEcho(t)。
CN201710608975.6A 2017-07-25 2017-07-25 基于移动激励源fdtd算法的高精度sar回波仿真方法 Active CN107271977B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710608975.6A CN107271977B (zh) 2017-07-25 2017-07-25 基于移动激励源fdtd算法的高精度sar回波仿真方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710608975.6A CN107271977B (zh) 2017-07-25 2017-07-25 基于移动激励源fdtd算法的高精度sar回波仿真方法

Publications (2)

Publication Number Publication Date
CN107271977A CN107271977A (zh) 2017-10-20
CN107271977B true CN107271977B (zh) 2020-04-24

Family

ID=60079091

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710608975.6A Active CN107271977B (zh) 2017-07-25 2017-07-25 基于移动激励源fdtd算法的高精度sar回波仿真方法

Country Status (1)

Country Link
CN (1) CN107271977B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109035186B (zh) * 2018-07-06 2021-03-26 电子科技大学中山学院 基于并行随机采样与传播的地形非特征区域合成方法、装置和电子设备
CN110426688A (zh) * 2019-07-02 2019-11-08 中国航空工业集团公司雷华电子技术研究所 一种基于地形背景目标的sar回波模拟方法
EP3764110B8 (en) 2019-07-11 2022-04-20 Rohde & Schwarz GmbH & Co. KG Method and system for near-field reconstruction in indirect far-field systems
CN110765687B (zh) * 2019-10-23 2022-06-21 哈尔滨工业大学 基于时域有限差分方法的涡旋波束源合成方法
CN113051793A (zh) * 2021-03-12 2021-06-29 西安电子科技大学 一种谐振辐射装置传输功率和比吸收率的仿真外推方法
CN114415181B (zh) * 2022-03-31 2022-07-12 南方海洋科学与工程广东省实验室(广州) 一种合成孔径雷达的原始回波生成方法和装置

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3824561B2 (ja) * 2002-07-05 2006-09-20 独立行政法人科学技術振興機構 電磁環境解析方法、電磁環境解析プログラム及びそのプログラムを記録した記録媒体
CN101576622A (zh) * 2009-06-12 2009-11-11 成都理工大学 一种超宽带电磁波的模拟方法
JP5865262B2 (ja) * 2010-01-18 2016-02-17 コーニンクレッカ フィリップス エヌ ヴェＫｏｎｉｎｋｌｉｊｋｅ Ｐｈｉｌｉｐｓ Ｎ．Ｖ． 電気特性断層画像化方法及びシステム
CN102207987B (zh) * 2011-05-31 2012-11-21 中国航天标准化研究所 基于OpenCL的GPU加速三维时域有限差分电磁场仿真的方法
CN102810126A (zh) * 2012-07-18 2012-12-05 上海交通大学 一种用于降低mimo-ue sar的方法及其仿真系统
CN103969627B (zh) * 2014-05-26 2016-07-06 苏州市数字城市工程研究中心有限公司 基于fdtd的探地雷达大规模三维正演模拟方法
CN104317984B (zh) * 2014-09-09 2017-06-27 中国舰船研究设计中心 基于分域建模的船舶电磁散射预测方法及系统
CN104992001B (zh) * 2015-06-19 2018-01-23 中国人民解放军信息工程大学 大规模mimo阵列天线远场辐射场的精确快速计算方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
FDTD analysis of EM field effected by the motion of source and vehicle;S.Sahrani, etal;《2014 International Conference on Electromagnetics in Advanced Applications (ICEAA)》;20140922;77-80 *
HIGH-RESOLUTION SAR SIGNAL SIMULATION USING PARALLEL FDTD METHOD;Xin pan, etal;《IGARSS 2014》;20141231;490-493 *
基于FDTD的高保真SAR回波信号仿真方法;张豪杰 等;《系统工程与电子技术》;20160131;第38卷(第1期);45-52 *

Also Published As

Publication number Publication date
CN107271977A (zh) 2017-10-20

Similar Documents

Publication Publication Date Title
CN107271977B (zh) 基于移动激励源fdtd算法的高精度sar回波仿真方法
CN106501802B (zh) 高分辨率多维协同昆虫迁飞雷达测量仪
CN107976660B (zh) 弹载多通道雷达超低空目标分析与多径回波建模方法
CN107037429B (zh) 基于门限梯度追踪算法的线阵sar三维成像方法
CN103713288A (zh) 基于迭代最小化稀疏贝叶斯重构线阵sar成像方法
He et al. Fast 3D-ISAR image simulation of targets at arbitrary aspect angles through nonuniform fast Fourier transform (NUFFT)
CN103018741A (zh) 一种基于后向投影的InSAR成像去平地一体化方法
Ren et al. 3D imaging algorithm for down-looking MIMO array SAR based on Bayesian compressive sensing
Zhou et al. Synthesis of conformal array antenna for hypersonic platform SAR using modified particle swarm optimisation
CN104101869B (zh) 一种极坐标下的地波雷达运动目标仿真建模方法
He et al. Fast ISAR image simulation of targets at arbitrary aspect angles using a novel SBR method
Dai et al. Scattering simulation and reconstruction of a 3-D complex target using downward-looking step-frequency radar
CN109884621B (zh) 雷达高度表回波相参积累方法
Liu et al. Efficient and precise frequency-modulated continuous wave synthetic aperture radar raw signal simulation approach for extended scenes
CN113447896B (zh) 一种基于动态遮挡判断的起伏地形sar回波仿真方法
Borden Some issues in inverse synthetic aperture radar image reconstruction
CN111859704A (zh) 一种分布式多视角下非刚体目标电磁散射建模方法
Liu et al. Bistatic FMCW SAR raw signal simulator for extended scenes
Wei et al. Sparse reconstruction for linear array SAR 3-D imaging based on Bayesian estimation
Fang et al. Three-dimensional near-field microwave imaging approach based on compressed sensing
Qingwu et al. High Precision SAR Echo Simulation Based on FDTD Algorithm with Mobile Excitation Source
CN115379560B (zh) 一种无线传感网络中仅有距离量测信息下的目标定位与跟踪方法
Vaupel et al. Validation of a 3-D near-field ISAR imaging technique with far-field RCS extraction by means of a hybrid GO-PO/PTD ray tracing algorithm
AU2021105419A4 (en) Highly-dynamic Radar Platform Echo Modeling Method Based on Space-time Decomposition
Peng et al. The simulation of SAR echo-wave signals for natural target

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant