基于频谱优化建模的双基地前视sar频域成像方法
【专利摘要】本发明公开了一种基于频谱优化建模的双基地前视SAR频域成像方法，包括以下步骤：S1、进行数据预处理，得到回波信号；S2、对回波信号进行二维傅立叶变换，建立系统的点目标响应最优二维频谱模型，并通过差分进化优化方法求解点目标响应精确二维频谱；S3、对点目标响应精确二维频谱进行粗匹配聚焦；S4、对二维频谱进行距离频率变换；S5、对距离频率变换后的二维频谱进行二维逆傅立叶变换，得到最终的复图像。本发明克服了现有成像算法无法处理双基地前视SAR距离向非线性空变的问题，通过差分进化优化方法得到精确的双基地前视SAR二维频谱，根据点目标响应精确的二维频谱在二维频域对回波信号进行成像处理，可实现双基地前视SAR原始数据的精确聚焦。
【专利说明】
基于频谱优化建模的双基地前视SAR频域成像方法
技术领域
[0001] 本发明属于雷达技术领域，特别涉及一种基于频谱优化建模的双基地前视SAR频 域成像方法。
【背景技术】
[0002] 双基地前视SAR成像是一种具有重要军事应用和民用价值的新概念、新体制SAR技 术，是指:发射站安装在飞行器之外的平台上（高地、飞机、卫星），飞行器只承载接收站，收 发波束指向接收站所在飞行器的正前方区域，利用收发站之间的协同，通过成像处理，实现 飞行器正前方的高分辨雷达成像。其应用背景广泛，有自主着陆、自主导航、物资和人员空 投、前视侦查、景象匹配末制导等。
[0003] 与单基地SAR类似，双基地SAR可以在方位频域进行成像处理，以加快其处理速度。 文南犬"B? Liu，T? Wang，Q.Wu，et al.Bistatic SAR Data Focusing Using an Omega-K Algorithm Based on Method of Series Reversion[J].Geoscience and Remote Sensing，IEEE Transactions on，2009,47(8) :2899_2912?"中提出了一种基于序列反转方 法计算点目标响应二维频谱的频域成像算法。文献"Junjie Wu，Zhongyu Li，Yulin Huang， J i anyu Yang ? Qing Huo Liu.An0mega-K Algorithm for Translational InvariantBi stati c SAR based on Generalized Loffeld'sBistatic Formula， Geoscience and Remote Sensing, IEEE Transactions on，2014,6699_6714,中提出了一 种基于广义Loffeld模型的频域成像算法。然而上述算法均假设双基地SAR距离向线性空 变。然而在双基地前视SAR中，接收机前视这一特殊几何结构导致距离向线性空变的假设并 不适用。因此，上述算法均无法直接应用于双基地前视SAR成像模式中。
【发明内容】
[0004] 本发明的目的在于克服现有成像算法无法处理双基地前视SAR距离向非线性空变 的问题，提供了一种通过差分进化优化方法得到精确的双基地前视SAR二维频谱，根据点目 标响应精确的二维频谱，在二维频域对回波信号进行成像处理，可实现双基地前视SAR原始 数据的精确聚焦的基于频谱优化建模的双基地前视SAR频域成像方法。
[0005] 本发明的目的是通过以下技术方案来实现的：基于频谱优化建模的双基地前视 SAR频域成像方法，包括以下步骤：
[0006] S1、进行数据预处理，得到回波信号；
[0007] S2、对回波信号进行二维傅立叶变换，建立系统的点目标响应最优二维频谱模型， 并通过差分进化优化方法求解点目标响应精确二维频谱；
[0008] S3、对步骤S2得到的点目标响应精确二维频谱进行粗匹配聚焦；
[0009] S4、对步骤S3进行粗匹配聚焦后的二维频谱进行距离频率变换；
[0010] S5、对距离频率变换后的二维频谱进行二维逆傅立叶变换，得到最终的复图像。
[0011] 具体地，步骤S1包括以下子步骤：
[0012] SI 1、初始化成像系统参数;记发射站的位置坐标为(xt，yT，hT)，接收站零时刻位置 坐标为(xr，yR，hR)，其中，xt、yT和hT分别为发射站的x轴、y轴和z轴坐标;xr、yR和hR分别为接 收站的x轴、y轴和z轴坐标;场景中心坐标为(0,0,0);波速中心位于场景坐标原点处时记为 零时刻，平台速度为V，场景中任一点目标的位置坐标为P(x，y);
[0013] S12、双基地前视SAR回波信号建模；
[0014]设发射信号为：
[0015] st(i) = so(x)exp( j23if〇x)
[0016] 其中，fQ为发射信号的载频;t为距离向时间；SQ(t)为发射信号St(i)的包络，是一 个调频斜率为K的线性调频信号:sob) = exp( jJiKi2);
[0017] 发射站和接收站的距离历史分别为：
[0020] rT，rR分别为发射站和接收站与点目标P(x，y)的最近斜距，
:;_0ST，0SR分别为零时刻发射站与接收站的斜视角;n表不方位向慢时间；
[0021] 双基距离和为：
[0022] Rbi(n;x，y) =RT(n;x，y)+RR(n;x，y)
[0023] 发射信号经过点目标P(x，y)散射后返回接收站的延迟时间为：
[0025]式中c为电磁波传播速度；
[0026]回波信号表示为：
[0027] sr(x,n；x,y) = exp( j3TK(x-Td)2)exp[ j23Tfo(x-Td)]〇
[0028] 具体地，步骤S2包括以下子步骤：
[0029] S21、对回波信号进行二维傅立叶变换，建立点目标响应精确二维频谱模型：
[0030] 对回波信号进行二维傅立叶变换，得到回波信号二维频域表达式为：
[0031] S2df(f，fn;x，y) = So(f)J'exp{-j 伞（n，fn;x，y)}dq
[0032] 其中，So(f)为sq(t)傅立叶变换;上式中：
[0034] 利用驻定相位定理，对巾（11，匕;1，7)进行泰勒级数展开及序列反转，求得系统的 点目标响应的精确二维频谱的表达式为：
[0035] Sc(f ,fn；x) = So(f )exp{-j 〇c(f ,fn；x)}
[0036] 其中，
[0038] 其中，ki(x)、k2(x)、k3(x)、k4(x)分别为双基距离和Rbi的一阶导、二阶导、三阶导、 四阶导在零多普勒时刻n = 〇时的取值，表达式如下：[0039] ki(x) =-vsin9sT-vsin9sR
[0043] S22、建立系统的点目标响应精确二维频谱的相位模型：
[0045]其中，f为距离向频率;fn为多普勒频率;Rbi (〇)为接收机在零时刻场景中心点到接 收站和发射站的双基距尚和;ai，a2，a3，bo，bi，b2，b3，b4，bs，b6，b7，bs，bg，bio为待求解的参数； [0046] S23、建立目标函数:基于系统的点目标响应最优二维频谱的相位模型与精确二维 频谱的相位模型，构造目标函数：
[0048]其中，Xl代表距离向地距坐标采样点，一共采样了 i个点；
[0049] 上式中，(了)(,(/，./:/;\.)-(1 )<,.(/，./:/3,)表示在地距坐标为11时，点目标响应精确二维 频谱与最优二维频谱之间的差值，卜％(/,/^)|| 2代表对该差值求二范数，
代表对i个地距坐标采样点上的点目标响应精确二维频谱 与最优二维频谱之间的差值的二范数进行求和；
[0050] S24、采用差分进化优化方法求解参数:差分进化方法的种群由若干个体组成，每 个个体代表优化问题的一个潜在解;差分进化方法的优化机制是根据不同个体之间的距离 和方向信息来生成新的候选个体，实现群体进化;差分进化方法采用变异、交叉、选择这三 个典型进化算子对种群进行更新；
[0051] 设种群为P种群规模为NP，解的维数为N，第G代个体i可以表示为实数解空间中的 一个向量二以卜/心,We.丨丨，2，"?，八;
[0052] 用差分进化优化方'法求解参数ai，a2，a3，bo，bi，b2，b3，b4，bs，b6，b7，bs，bg，bio，即转 化为解优化问题：
[0053] minF(pii，pi2, ? ? ?，piN)
[0054] S.t Pij min^Pij^Pij max, j = l,2, . . . ,N
[0055]其中，解的维数N为待求解参数的个数14 ;Pij为个体Pi的第j维元素 ，Pij min和Pij max 分别为PU的下界和上界；
[0056] 求解上述优化问题包括以下步骤：
[0057] S241、初始化种群，初始种群 ^ | 尸-S 巧,/ = U，:…,/VP: / = 1,2,...，iV)随 机产生：
[0058] /<=/，,,_ 十"臟/(〇,丨
[0059] 其中rand(0，l)表示在(0，1)区间均匀分布的随机数；
[0060] S242、判断进化是否终止，判断是否达到终止条件或进化代数达到最大；若是，则 进化终止，将此时的最佳个体作为解输出；否则继续操作；
[0061 ] S243、变异操作，从当代种群中随机选择3个个体作为父代个体进行变异，生成变 异个体〃f，变异公式如下所示：
[0062] < =Ph./+/7-(/^.,-/?；；.,)
[0063] 其中，<是变异个体<的第j维元素；/为第G代种群中3个互异的个体， 并且rl辛r2辛r3辛i;F是实数常量，称为"缩放因子"，其取值范围为FG [0,2];
[0064] S244、交叉操作，将多个父代个体按照二项式交叉组合，生成新的候选解#，称为 "试验个体"，二项式交叉公式如下所示：
[0066]其中vf为候选解vf的第j维元素;CR是[0，1]之间的常数，称之为"交叉概率"；Ir是 [1，N]之间的随机整数；
[0067] S245、选择操作，用贪婪选择策略来对种群进行更新，通过比较新生成的试验个体 与当代种群中对应个体的优劣，选择适应值更优的个体作为子代进入新种群，选择公式如 下所示：
[0069] 具体地，步骤S3的实现方法为:根据S2步骤，得到为实现粗匹配聚焦的参考函数的 表达式为：
[0070] 心财(/,/";:%) = cxpL./(\(/，.
[0071 ]其中，xo为场景中心点的x坐标，参考相位为：
[0073]将SAR回波信号进行二维傅里叶变换得到的信号S2df(f，fn;x，y^S RFM(f，fn;xo^ 乘，实现了粗匹配聚焦，移除二维频谱中的空不变项，得到的残余相位为：
[0075]其中 Ax = Xi_xo。
[0076]具体地，步骤S4的具体实现方法为：利用下列距离频率变换公式对步骤S3得到的 二维频谱进行距离频率变换：
[0078]进行了距离频率变换后的二维频谱的残余相位为：
[0080]本发明的有益效果是:利用差分进化优化方法求解双基地前视SAR点目标响应二 维频谱，首先建立系统的点目标响应最优二维频谱模型，然后建立目标函数，再通过差分进 化优化方法求解得到系统的点目标响应精确二维频谱，最后根据点目标响应精确二维频 谱，在二维频域对回波信号进行成像处理，解决了现有成像算法无法处理双基地前视SAR距 离向非线性空变的问题，从而实现双基地前视SAR在接收机正前方区域的成像;本发明对双 基地前视SAR点目标响应二维频谱建模和求解，在二维频域对回波信号进行成像处理，可实 现双基地前视SAR原始数据的精确聚焦，且聚焦精度高。
【附图说明】
[0081 ]图1为本发明的成像方法流程框图；
[0082]图2为本发明具体实施例采用的双基地前视SAR几何结构图，其中横坐标表示距离 维，纵坐标表示多普勒维；
[0083]图3为本发明具体实施例中采用的目标场景布置图；
[0084]图4为本发明具体实施例的点目标的成像结果图。
【具体实施方式】
[0085]本发明主要采用计算机仿真的方法进行验证，所有步骤、结论都在MATLAB2013上 验证正确。下面结合附图和具体实施例进一步说明本发明的技术方案。
[0086]如图1所示，基于频谱优化建模的双基地前视SAR频域成像方法，包括以下步骤：
[0087] S1、进行数据预处理，得到回波信号;具体包括以下子步骤：
[0088] S11、初始化成像系统参数;本实施例采用的双基地前视SAR集合结构如图2所示。 其参数如表一所示。本实施例采用的目标场景如图3所示，图中的黑色圆点为布置于地面上 的5个点目标，这5个点分布于x轴和y轴上，沿x方向（切航迹）间隔50m，沿y方向（沿航迹）间 隔50m，平台沿y轴运动。
[0089] 表一
[0091]波速中心位于场景坐标原点处时记为零时刻，平台速度为v发射站的位置坐标为 (-800，-500，1000)m，接收站零时刻位置坐标为(-600,-300，800)m，场景中心坐标为(0,0， 〇);波速中心位于场景坐标原点处时记为零时刻，平台速度为50m/s，场景中任一点目标的 位置坐标为P (x，y);用MATLAB仿真出回波数据。
[0092] S12、双基地前视SAR回波信号建模；
[0093]设发射信号为：
[0094] st(x) = so(x)exp( j23if〇x)
[0095] 其中，fo为发射信号的载频p为距离向时间；sob)为发射信号St(i)的包络，是一 个调频斜率为K的线性调频信号:sob) = exp( jJiKi2);
[0096]发射站和接收站的距离历史分别为：
[0099] rT，rR分别为发射站和接收站与点目标P(x，y)的最近斜距
;Qst，Qsr分别为零时刻发射站与接收站的斜视角;n表不方位向慢时间；
[0100] 双基距离和为：
[0101 ] Rbi(n；x,y) =Ri(n；x,y)+RR(n；x,y)
[0102]假设首发站能够很好的同步，发射信号经过点目标P(x，y)散射后返回接收站的延 迟时间为：
[0104] 式中c为电磁波传播速度；
[0105] 回波信号表示为：
[0106] Sr(x,n；x,y) = exp( j3TK(x-Td)2)exp[ j23Tfo(x-Td)] 〇
[0107] S2、对回波信号进行二维傅立叶变换，建立系统的点目标响应最优二维频谱模型， 并通过差分进化优化方法求解点目标响应精确二维频谱;具体包括以下子步骤：
[0108] S21、对回波信号进行二维傅立叶变换，建立点目标响应精确二维频谱模型：
[0109]对回波信号进行二维傅立叶变换，得到回波信号二维频域表达式为：
[0110] S2df(f，fn;x，y) = So(f)i'exp{-j 伞（n，fn;x，y)}dn
[0111] 其中，So(f)为sq(t)傅立叶变换;上式中：
[0113]利用驻定相位定理，对巾（11，匕;1，7)进行泰勒级数展开及序列反转，求得系统的 点目标响应的精确二维频谱的表达式为：
[0114] Sc(f ,fn；x) = So(f )exp{-j 〇c(f ,fn；x)}
[0115] 其中，
[0117] 其中，ki(x)、k2(x)、k3(x)、k4(x)分别为双基距离和Rbi的一阶导、二阶导、三阶导、 四阶导在零多普勒时刻n = 〇时的取值，表达式如下：
[0118] ki(x) =-vsin9sT-vsin9sR
[0122] S22、建立系统的点目标响应精确二维频谱的相位模型：
[0124] 其中，f为距离向频率;fn为多普勒频率;Rbi (0)为接收机在零时刻场景中心点到接 收站和发射站的双基距尚和;ai，a2，a3，bo，bi，b2，b3，b4，bs，b6，b7，bs，bg，bio为待求解的参数；
[0125] S23、建立目标函数:基于系统的点目标响应最优二维频谱的相位模型与精确二维 频谱的相位模型，构造目标函数：
[0127] 其中，Xl代表距离向地距坐标采样点，一共采样了i个点；
[0128] 上式中，私,(/，./:/4)-(1)(乂/，./:_,4)表示在地距坐标为1 1时，点目标响应精确二维 频谱与最优二维频谱之间的差值，||K/,/；;x,)(/，/,4)|2代表对该差值求二范数，
代表对i个地距坐标采样点上的点目标响应精确二维频谱 与最优二维频谱之间的差值的二范数进行求和；
[0129] S24、采用差分进化优化方法求解参数:差分进化方法的种群由若干个体组成，每 个个体代表优化问题的一个潜在解;差分进化方法的优化机制是根据不同个体之间的距离 和方向信息来生成新的候选个体，实现群体进化;差分进化方法采用变异、交叉、选择这三 个典型进化算子对种群进行更新；
[0130] 设种群为P种群规模为NP，解的维数为N，第G代个体i可以表示为实数解空间中的 一个向量/广），We '丨 12,.…八;
[0131 ] 用差分进化优化方'法求角军参数ai，a2，a3，bo，bi，b2，b3，b4，bs，b6，b7，bs，bg，bio，即车专 化为解优化问题：
[0132] min F(pn，pi2, ? ? ?，piN)
[01 33] S . t Pij min^sPij ^；Pij max , j = 1,2 , . . . , N
[0134]其中，解的维数N为待求解参数的个数14;Pij为个体Pi的第j维元素 ，Pij min和Pij max 分别为PU的下界和上界；
[0135] 求解上述优化问题包括以下步骤：
[0136] S241、初始化种群，初始种群^ I巧_ < ，/ = 1,2,…,1，2,…，iV)随机 产生：
[0137] 七騰/((U)?(八
[0138] 其中rand(0，l)表示在(0，1)区间均匀分布的随机数；
[0139] S242、判断进化是否终止，判断是否达到终止条件或进化代数达到最大；若是，则 进化终止，将此时的最佳个体作为解输出；否则继续操作；
[0140] S243、变异操作，从当代种群中随机选择3个个体作为父代个体进行变异，生成变 异个体，变异公式如下所示：
[01 41 ] t4=P^.J+F'(Pr2.J-Pri,j)
[0142] 其中，4是变异个体的第j维元素；为第G代种群中3个互异的个体， 并且rl辛r2辛r3辛i;F是实数常量，称为"缩放因子"，其取值范围为FG [0,2];
[0143] S244、交叉操作，将多个父代个体按照二项式交叉组合，生成新的候选解f，称为 "试验个体"，二项式交叉公式如下所示：
[0145] 其中 <为候选解vf的第j维元素;CR是[0，1]之间的常数，称之为"交叉概率"；Ir是 [1，N]之间的随机整数；
[0146] S245、选择操作，用贪婪选择策略来对种群进行更新，通过比较新生成的试验个体 与当代种群中对应个体的优劣，选择适应值更优的个体作为子代进入新种群，选择公式如 下所示：
[0148] S3、对步骤S2得到的点目标响应精确二维频谱进行粗匹配聚焦;其具体实现方法 为:根据S2步骤，得到为实现粗匹配聚焦的参考函数的表达式为：
[0149] SmM (/, fq;x0) = exp[70c(f, fit;x0)]
[0150] 其中，xo为场景中心点的x坐标，参考相位为：
[0152]将SAR回波信号进行二维傅里叶变换得到的信号S2df(f，fn;x，y^S RFM(f，fn;xo^ 乘，实现了粗匹配聚焦，移除二维频谱中的空不变项，得到的残余相位为：
[0154] 其中 Ax = Xi_xo。
[0155] S4、对步骤S3进行粗匹配聚焦后的二维频谱进行距离频率变换;其具体实现方法 为:利用下列距离频率变换公式对步骤S3得到的二维频谱进行距离频率变换：
[0157]在频率变换后，残余相位就与新的频域坐标成线性关系，进行了距离频率变换后 的二维频谱的残余相位为：
[0159] S5、对距离频率变换后的二维频谱进行二维逆傅立叶变换，得到最终的复图像。
[0160]本实施例采用MATLAB2013进行仿真，其点目标的成像结果如图4所示。由图4可以 看出，本发明能实现双基地前视合成孔径雷达成像，并得到良好聚焦。
[0161] 通过本发明的具体实施可以看出，本发明建立了距离非线性空变的双基地前视 SAR二维频谱模型，并利用差分进化优化方法求解精确二维频谱。最后，本方法在其二维频 域对回波信号进行成像处理，实现了双基地前视SAR原始回波数据的精确聚焦。
[0162] 本领域的普通技术人员将会意识到，这里所述的实施例是为了帮助读者理解本发 明的原理，应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。本领域的 普通技术人员可以根据本发明公开的这些技术启示做出各种不脱离本发明实质的其它各 种具体变形和组合，这些变形和组合仍然在本发明的保护范围内。
【主权项】
1. 基于频谱优化建模的双基地前视SAR频域成像方法，其特征在于，包括以下步骤： 51、 进行数据预处理，得到回波信号； 52、 对回波信号进行二维傅立叶变换，建立系统的点目标响应最优二维频谱模型，并通 过差分进化优化方法求解点目标响应精确二维频谱； 53、 对步骤S2得到的点目标响应精确二维频谱进行粗匹配聚焦； 54、 对步骤S3进行粗匹配聚焦后的二维频谱进行距离频率变换； 55、 对距离频率变换后的二维频谱进行二维逆傅立叶变换，得到最终的复图像。2. 根据权利要求1所述的基于频谱优化建模的双基地前视SAR频域成像方法，其特征在 于，所述步骤S1包括以下子步骤： SI 1、初始化成像系统参数;记发射站的位置坐标为(xt，yT，hT)，接收站零时刻位置坐标 为(XR，yR，hR)，其中，XT、yT和hT分另U为发射站的x轴、y轴私轴坐标;XR、yR和hR分另lj为接收站 的x轴、y轴和z轴坐标;场景中心坐标为(0，0，0);波速中心位于场景坐标原点处时记为零时 亥IJ，平台速度为V，场景中任一点目标的位置坐标为P(x，y); S12、双基地前视SAR回波信号建模； 设发射信号为： st ⑴=so ⑴ exp( j23if〇T) 其中，f〇为发射信号的载频；T为距离向时间；SQ(T)为发射信号St(T)的包络，是一个调 频斜率为K的线性调频信号:so⑴=exp (jJTKT2); 发射站和接收站的距离历史分别为：rT，rR分别为发射站和接收站与点目标P(x，y)的最近斜距;._0st，0sr分别为零时刻发射站与接收站的斜视角;n表不方位向慢时间； 双基距离和为： Rbi(n；x,y)=Ri(n；x,y)+RR(n；x,y) 发射信号经过点目标P(x，y)散射后返回接收站的延迟时间为：式中C为电磁波传播速度； 回波信号表不为： Sr(T，n;x，y) = exp( j3iK(T-Td)2)exp[ j23if〇(T-Td)] 〇3. 根据权利要求2所述的基于频谱优化建模的双基地前视SAR频域成像方法，其特征在 于，所述步骤S2包括以下子步骤： S21、对回波信号进行二维傅立叶变换，建立点目标响应二维频谱模型： 对回波信号进行二维傅立叶变换，得到回波信号二维频域表达式为： S2df(f ,fn;x,y)=So(f)/exp{-j <i> (q,fn；x,y)}dn 其中，So(f)为sob)傅立叶变换;上式中：利用驻定相位定理，对巾（n，fn;x，y)进行泰勒级数展开及序列反转，求得系统的点目标 响应最优二维频谱的表达式为： Sc(f ,fn;x) = So(f )exp{-j〇c(f ,fn；x)} 其中，其中，ki(x)、k2(x)、k3(x)、k4(x)分别为双基距离和Rbi的一阶导、二阶导、三阶导、四阶 导在零多普勒时刻n = 〇时的取值，表达式如下： ki(x) =-vsin9sT_vsin9sR522、 建立系统的点目标响应精确二维频谱的相位模型：其中，f为距离向频率;fn为多普勒频率;Rbi (0)为接收机在零时刻场景中心点到接收站 和发射站的双基距尚和;ai，a2，a3，bo，bi，b2，b3，b4，bs，b6，b7，bs，bg，bio为待求解的参数； 523、 建立目标函数:基于系统的点目标响应最优二维频谱的相位模型与精确二维频谱 的相位模型，构造目标函数：其中，Xl代表距离向地距坐标采样点，一共采样了 i个点； 上式中，屯,(/，./?<1^(/，/^6)表示在地距坐标为11时，点目标响应精确二维频谱 与最优二维频谱之间的差值，|%(/，./；4)(/,./；A)||2代表对该差值求二范数，;代表对i个地距坐标采样点上的点目标响应精确二维频谱 与最优二维频谱之间的差值的二范数进行求和； S24、采用差分进化优化方法求解参数:差分进化方法的种群由若干个体组成，每个个 体代表优化问题的一个潜在解;差分进化方法的优化机制是根据不同个体之间的距离和方 向信息来生成新的候选个体，实现群体进化;差分进化方法采用变异、交叉、选择这三个典 型进化算子对种群进行更新； 设种群为P种群规模为NP，解的维数为N，第G代个体i可以表示为实数解空间中的一个 向量，Vz'e{.l，2,; 用差分进化优化方法求解参数ai, a2, a3, bo, bi, b2, b3, b4, b5, b6, b7, bs, bg, bio，即转化为 解优化问题： min F(pu，pi2, ? ? ?，piN) S . tpijmin^Pij ^Pijmax , j - 1,2 , . . . , N 其中，解的维数N为待求解参数的个数14 ; Pij为个体Pi的第j维元素，Pijmin和Pijnmx分别为 Pij的下界和上界； 求解上述优化问题包括以下步骤： 3241、初始化种群，初始种群{只°|汽*43巧臟,，.=1，2,~八沪;./ = 1,2，..,，叫随机产 生：其中rand(0，l)表示在(0，1)区间均匀分布的随机数； 5242、 判断进化是否终止，判断是否达到终止条件或进化代数达到最大；若是，则进化 终止，将此时的最佳个体作为解输出；否则继续操作； 5243、 变异操作，从当代种群中随机选择3个个体作为父代个体进行变异，生成变异个 体，变异公式如下所示：其中是变异个体的第j维元素；为第G代种群中3个互异的个体，并且rl 辛r2辛r3#i;g是实数常量，称为"缩放因子"，其取值范围为ge[0,2]; 5244、 交叉操作，将多个父代个体按照二项式交叉组合，生成新的候选解vf，称为"试验 个体"，二项式交叉公式如下所示：其中vf为候选解vf的第j维元素;CR是[0，1]之间的常数，称之为"交叉概率"；Ir是[1， N]之间的随机整数； 5245、 选择操作，用贪婪选择策略来对种群进行更新，通过比较新生成的试验个体与当 代种群中对应个体的优劣，选择适应值更优的个体作为子代进入新种群，选择公式如下所 示：4. 根据权利要求3所述的基于频谱优化建模的双基地前视SAR频域成像方法，所述步骤 S3具体实现方法为:根据S2步骤，得到为实现粗匹配聚焦的参考函数的表达式为：其中，xo为场景中心点的x坐标，参考相位为：将SAR回波信号进行二维傅里叶变换得到的信号52<?^仏3,7)与51^^仏3())相乘， 移除二维频谱中的空不变项，得到的残余相位为：其中 A x = Xi_xo。5. 根据权利要求4所述的基于频谱优化建模的双基地前视SAR频域成像方法，其特征在 于，所述步骤S4的具体实现方法为:利用下列距离频率变换公式对步骤S3得到的二维频谱 进行距离频率变换：进行了距离频率变换后的二维频谱的残余相位为：
【文档编号】G01S13/90GK106054190SQ201610579683
【公开日】2016年10月26日
【申请日】2016年7月21日 公开号201610579683.X, CN 106054190 A, CN 106054190A, CN 201610579683, CN-A-106054190, CN106054190 A, CN106054190A, CN201610579683, CN201610579683.X
【发明人】武俊杰, 蒲巍, 孙稚超, 任建宇, 黄钰林, 杨建宇, 杨海光
【申请人】电子科技大学