一种多输入多输出雷达波形设计方法
一种多输入多输出雷达波形设计方法
【技术领域】
[0001] 本发明属于雷达通信技术领域，特别涉及MM0雷达满足时域频域任务需求的雷 达波形设计方法。
【背景技术】
[0002] 当今军事或民用需求使得雷达任务场景日趋复杂。为应对复杂变化的场景， 雷达自适应性显得尤为重要。其中，波形分集可视为自适应性的一种体现且成为目前 石开宄的热点（见文献：Cognitive radar: a way of the future, Haykin S. ；Signal Processing Magazine, IEEE, 2006, 23(1):30-40 ；Adaptive radar phase-coded waveform design, J. D. Zhang, X. H. Zhu, H. Q. Wang ；Electron. Lett. 2009, 45 (20) 1052 - 1053 ；Spatial diversity in radars-modeIs and detection performance, Eran, F. , Alexander, H. , &Ri ck，S.B. ; IEEE Transactions on Signal Processing, 2006, 54 (3) ,823 - 838)。无论是相 控阵雷达还是MM0雷达，波形分集既能提高目标识别能力又能提高分辨力，已然成为提高 雷达自适应性的重要切入点。以MIM0雷达为例，考虑到接收端滤波器可视为波形的相关函 数，结合波形设计与复杂任务场景分类：避免强散射体距离旁瓣遮蔽、防止不同阵元间波形 互扰、躲避特定频带电磁干扰，这三者即成为波形设计的基本要求。此外，为充分利用发射 机功率，恒模特性亦成为客观要求。
[0003] 避免强散射体遮蔽效应需要发射波形具备较低的自相关旁瓣；消除不同阵元相 互影响要求发射波形正交化；抑制特定频带电磁干扰要求发射波形频谱具备稀疏特性； 另外，恒模约束往往使得设计问题非凸。以上这些要求使得波形设计问题较为困难。针 对恒模波形提出的循环算法族（Cyclic Algorithm New，CAN)，其思路是基于相关特性 与功率谱之间的傅里叶对应关系，逐步迭代以优化所选目标函数（见文献：Waveform design for active sensing systems:a computational approach, He H, Li J, Stoica P. ;Cambridge University Press, 2012)。但该算法耗时过长、相关旁瓣电平幅度较差， 难以满足工程实践。但其给后续研宄提供了众多启示，以交替投影为机制的改进循环算 法框架，获得了比CAN算法族更明显的效果（见文献：零自相关区相位编码波形设计， 李风从，赵宜楠，乔晓林；电子学报，2013, 12:2499-2502;基于秩亏傅里叶变换的交替 投影编码波形设计，赵宜楠，李风从，王军，乔晓林；电子学报，2014, 06:1216-1219; Computational design of optimal waveforms for MIM0 radar via multi-dimensional iterative spectral approximation, Zhao Y N, Li F C, Zhang T ;Multidm syst sign process，11045-01400288-1，2014);但算法初始化导致非凸问题求解稳定性较差。除此 之外，利用遗传算法对MIM0雷达正交波形进行设计（见文献：Polyphase orthogonal code design for MIM0 radar systems, Liu B, He Z, Zeng J ；2006 CIE international conference on radar, 2006:1-4);但该算法在迭代优化的后期往往出现群体多样性较差， 导致局部最优解而陷入停滞，另外算法规模也成为限制其效率的因素。但针对实际MIM0雷 达场景，此类算法可能耗时过长且难以收敛。利用序列二次规划与加权迭代方法求解波形 设计问题，权值步长的搜索机制难以自适应；此类基于梯度矩阵的算法在多阵元情形计算 复杂度较高，且可能出现高阶矩阵奇异而难以收敛（见文献：一种空时联合优化的MMO雷 达波形设计方法，王旭，周生华，刘宏伟，保铮；西安电子科技大学学报，2014, 03:41-48)。
[0004] 因恒模导致的非凸问题往往收敛缓慢且难以得到全局最优解，而交替投 影算法对于凸范数求解亦可能存在局部停滞而难以收敛（见文献：〇n pro jection algorithms for solving convex feasibility problems,Bauschke H H,Borwein J M,SIAM review,1996, 38(3):367-426 ；An alternating projection that does not converge in norm, Hundal H S，Nonlinear Analysis:Theory,Methods&Applicatio ns，2004, 57 (1) : 35-61)，而松弛交替投影算法作为一种高效且耗时较少算法（见文献： Relaxed alternating projection methods, A. Cegielski, A. Suchocka, SIAM Journal on Optimization, vol. 19(3)，2008, pp. 1093-1106)，结合迭代谱近似思想可用来解决以恒模 波形设计为代表的非凸问题。
【发明内容】
[0005] 本发明的目的是提供一种具有更低的相关旁瓣和频谱抑制深度，且效率高、耗时 少、可有效避免局部最优，并具有较高的鲁棒性、具有良好时频抗干扰性能的多输入多输出 雷达波形设计方法。
[0006] 本发明为解决上述技术问题采取的技术方案是：
[0007] -种多输入多输出雷达波形设计方法，包括以下步骤：
[0008] 步骤1:根据雷达场景中强散射体与待测目标的相对位置，预估自相关旁瓣抑制 模糊区间，进而构造相应的目标函数：
[0009] 设M个发射阵元的MIM0雷达系统，其第m个阵元窄带恒模相位编码波形序列可表 不为
，N表征信号码长，itnG [0, 2 JT ] 为第nth时刻编码相位；因接收端滤波器可视为发射波形的相关函数，假定感兴趣的待测目 标在第pth距离单元；那么在第Ith单元的强散射体势必要对pth单元目标造成影响；对于第 p th距离单元，为最小化临近距离单元中强散射体产生的距离旁瓣对它的干扰，应使发射波 形1的自相关函数a (Xm)满足如下条件：
[0011]其中，a (xm) = [ a 0(xm)…an-JxJ0 a _N+1 (xm) ??? a < (xm) ]T表示自相关函数序 列；毛，表示理想波形；Qi= {± |z-p| |z G ZiUm-L-N+l}表示旁瓣干扰区间，&表示距离 向分布的强散射体位置集合；
[0012] 对于强散射体遮蔽干扰问题，发射波形Xm应使得目标函数J(Xm)最小化，即：
[0014]其中N ? | |表示欧式范数，视为距离度量函数；由Parseval等价性可知：
[0016]其中，
表示离散傅里叶变换矩阵，
表示自相 关函数的傅里叶变换，即功率谱；如果令?/" :和fm分别表征理想波形和设计波形的频谱， 即
，C为扩展矩阵，其通过扩充零元素使得矩阵运算时维度一致；那么易得：
I其中?表示按元素的Hadamard积，（? r表示取共轭操作；目标函数 写为：
[0018]步骤2 :分析MM0雷达波形正交性约束，构造满足正交性约束的目标函数：
[0019] 对于MM0发射波形矩阵
而言，波形之间正交性体现在互相关 函数上；第nth时刻的发射波形相关矩阵Rn(X)的非主对角线部分尽可能逼近零，主对角元 素为各阵元发射波形的自相关表示：
[0023] 其中diag(a"(l))表示nth时刻由理想波形矩阵.文的自相关形成的对角矩阵；U(n) 表示偏移矩阵，I表示单位矩阵；
[0024] 为解决不同阵元波形间互扰构造目标函数J2 (X)，同样考虑Parseval等价性，得到 如下：
[0026] 其中| | ? | |F表示矩阵的Frobenius范数；
[0027] 步骤3 :根据场景先验信息预估频域干扰模糊频带区间，进而构造相应目标函数：
[0028] 假定可能存在电磁干扰的频段区间集合为
，那么发射波形频谱应 在此区间内设为阻带，由傅里叶变换的性质可知，频谱为：
[0030]其中：NS表示相应频段数目，表示单位离散傅里叶变换矩阵，且有 为使得频谱离散采样足够密以表征区间信息，设定及= 2#;由上式可知，如果从离散傅里 叶变换矩阵& I的行空间中抽取特定行来构成矩阵P，且所抽取行序号与D在频谱W中频 段&相对应，用向量?斤表征&中稀疏凹口序列且满足
fn为频 谱W中频率的采样，从而使得：
[0032]只要# [幅度足够低就满足频谱稀疏特性，即目标函数J3(X)最小化，从而达到任 务要求：
[0034] 步骤4:构造恒模相位编码波形约束条件，即最大化发射功率要求恒模约束，也 即：
[0035] |Xn'J = 1，n = 1，…N. m = 1，? ?，M
[0036] 步骤5 :构造松弛交替投影算法框架：
[0037]假定集合A代表约束集合（即满足恒模约束），集合B代表目标集合（即满足子任 务要求），构造变量X，y到集合A,B的投影ProjA(x)和projB(y)，如下所示：
[0040] 其中d(i，j)= I Ii-jI I表示变量i，j的距离函数。
[0041] 那么交替投影机制可简述为：迭代映射T:A - B - A。
[0042] 然后，引入松弛交替投影框架：
[0043] Ts(.)，入x = ProjA(x+人? 8 (x) ? (ProjA(ProjB(x))-x))
[0044]xk+1= Ts(.)，Xxk
[0045] 上式中，A表示松弛算子；S (?)表示加速因子，有如下不同形式：
[0048] 此算法框架命名为RISAAP。
[0049] 步骤6、上述步骤1、2、3分别完成了三个子任务目标函数的构建，以适应雷达多任 务场景的时变性，某时刻场景可能涉及上述子任务目标中的一个或多个，根据步骤5所提 松弛交替投影算法框架求解波形设计，分别给出如下三种波形优化输出方式：
[0050] (1)表示满足抑制旁瓣遮蔽效应的算法步骤；
[0051] (2)表示获得波形正交性的算法步骤；
[0052] (3)表示联合抑制特定频带干扰及获得良好相关特性的算法步骤；
[0053] (1)对于抑制旁瓣遮蔽效应的算法，过程如下；
[0054] 以满足步骤1目标函数条件的集合为目标集合B ;以满足步骤4条件集合为约束 集合A ;(?) (1/2)表示按元素的平方根操作，返回值维度与其参数一致；
[0055] 那么获得较好时域自相关旁瓣抑制性能的恒模波形设计方法可归纳为：
[0056] a)初始化波形序列：x=[x(l). . . x(n) x(N)]TG A，其中
itnG [0, 2it]为随机相位；
[0057] b)迭代优化过程：
[0058] 由公式
[0060]获得f G B，
[0061] 再由
[0063]获得x G A，
[0064] 经公式
[0066]获得i并更新X = X;
[0067] c)相邻两次迭代产生J差值小于给定误差e或算法迭代次数大于给定总数K，则 算法停止并输出，否则继续b)，c);
[0068] (2)对于获得波形正交性的算法设计，过程如下：
[0069] 以满足步骤2目标函数条件的集合，即为目标集合B ;以满足步骤4条件的集合为 约束集合A ;那么获得较好时域相关旁瓣抑制性能的恒模波形设计方法可归纳为：
[0070] a)初始化波形矩阵：X= [Xl.? ? xM] G A，其中
，步[0, 2 n]为 随机相位；
[0071] b)迭代优化过程：
[0072] 由公式
[0073]
[0074] 获得 Fd
[0075] 再由
[0077] 获得 X G A，
[0078] 经公式
[0079] Ts(.)，入x = ProjA(x+人? 8 (x) ? (ProjA(ProjB(x))-x))
[0081] 获得X并更新X二X ;
[0082] c)若相邻两次迭代产生f差值小于给定误差e或算法迭代次数大于给定总数K， 则算法停止并输出，否则继续b)，c);
[0083] (3)对于联合抑制特定频带干扰及获得良好相关特性的算法设计，过程如下：
[0084] 以满足步骤1、2目标函数条件的集合，即为目标集合B1;以满足步骤3的集合，即 为目标集合B 2;以满足步骤4条件的集合为约束集合A ;那么获得较好时域相关旁瓣抑制性 能及频域特性的恒模波形设计方法可归纳为：
[0085] a)初始化波形矩阵：X = [Xl. . . xM] G A，其中
，步 n，me [0, 2 n ]为 随机相位；
[0086] b)迭代优化过程：
[0087] 由公式
[0089] 获得片e^，
[0090]由公式
[0092] 获得尽e /七，由
[0094]获得F e B，
[0095] 再由
[0097]获得 X e A，
[0098] 经公式
[0099] Ts(.)，入x = ProjA(x+人? 8 (x) ? (ProjA(ProjB(x))-x))
[0101] 获得1并更新x =
[0102] c)若相邻两次迭代产生差值小于给定误差e或算法迭代次数大于给定总数 K，则算法停止并输出，否则继续b)，c)。
[0103] 所述步骤1中时域干扰区间确定可根据先验信息，或通过事先发射传统波形以确 定大致区间。
[0104] 所述步骤2中多输入多输出雷达可包含多个阵元。
[0105] 所述步骤3中频域干扰区间的确定可根据先验信息、或参考通用电子手册确定大 致区间。
[0106] 所述步骤5中松弛算子A的取值范围是：〇〈 A〈2 ;Parteo权值y的取值范围是： 0〈丫〈1〇
[0107] 本发明的有益效果是：
[0108] 本发明采用松弛投影算法来设计多输入多输出雷达波形，相比采用现有的波形在 线设计方法，本发明具有更低的相关旁瓣和频谱抑制深度，且效率高、耗时少、可有效避免 局部最优，并具有较高的鲁棒性。应用本发明的松弛交替投影恒模波形编码设计，可使MIM0 雷达具有更好的检测性能。
[0109] 本发明所设计波形可有效避免强散射体遮蔽干扰、阵元波形间互扰及特定频带干 扰，更有利于目标的检测与跟踪，其优点如下：
[0110] (1)就其采用的恒模相位编码产生方式而言，本发明提出的MIM0雷达波形具有精 确的表达式及具体的步骤优化程序，避免了现有技术中涉及的基于梯度的高维非线性搜索 问题。
[0111] (2)就理想波形的构建而言，采用设计波形与其渐进拟合的方式，理想波形建立在 迭代设计的基础之上，避免了先期考虑的过多因素。
[0112] (3)就波形优化效果而言，本发明提出的MIM0雷达波形设计方法耗时少、抑制自 相关旁瓣及特定频带区间干扰远强于Li Jian等人所提出的CAN循环算法；更便于在线实 时设计。
[0113] 其他性能指标（归一化自相关幅值电平：ACL ;阻带峰值功率：PSP ;)与下述已有 技术相比，完全优于已有技术的相应指标。
[0114] 已有技术：以文献(waveform design for active sensing systems:a computational approach, He H, Li J, Stoica P. Cambridge University Press, 2012) 提出的CAN、WeCAN、SCAN算法族，包括Multi-CAN、WeCAN、SCAN ;及以文献（零自相关 区相位编码波形设计，李风从，赵宜楠，乔晓林；电子学报，2013, 12:2499-2502;基于 秩亏傅里叶变换的交替投影编码波形设计，赵宜楠，李风从，王军，乔晓林；电子学 报，2014,06:1216-1219 ;Computational design of optimal waveforms for MIM0 radar via multi-dimensional iterative spectral approximation, Zhao Y N, Li F C, Zhang T ； Multidm syst sign process, 11045-01400288-1，2014))提出的 ISAA、MDISAA 等。但这些 技术耗时过长，波形优化设计效果欠佳，不适合工程应用。
【附图说明】
[0115] 图1是本发明实施方式中，在特定的参数关系下，MM0雷达阵元信号在单个特定 区域自相关旁瓣抑制示意图。仿真参数：码长N = 150 ;算法迭代数K = 10000 ;迭代误差 阈值为e = 10、松弛因子A = 〇. 7。
[0116] 图2是本发明实施方式中，在特定的参数关系下，MM0雷达阵元信号在多个特定 区域自相关旁瓣抑制示意图。仿真参数：码长N = 150 ;算法迭代数K = 10000 ;迭代误差 阈值为e = 10、松弛因子A = 〇. 7。
[0117] 图3是本发明实施方式中，在特定的参数关系下，MM0雷达阵元间信号自相关旁 瓣抑制示意图。仿真参数：码长N = 64 ;阵元数M = 3 ;算法迭代数K = 10000 ;迭代误差阈 值e = 1(T4;松弛因子A = 〇? 7。
[0118] 图4是本发明实施方式中，在特定的参数关系下，MM0雷达阵元间信号互相关旁 瓣抑制示意图。仿真参数：码长N = 64 ;阵元数M = 3 ;算法迭代数为K = 10000 ;迭代误差 阈值为e = 10、松弛因子A = 〇. 7。
[0119] 图5是本发明实施方式中，在特定的参数关系下，MM0雷达信号在单个频带区域 抑制示意图。仿真参数：码长N = 100 ;阵元数M = 1 ;脉冲持续时间Tdming= 200 y s ;采样 频率fsample= 1000KHZ ;算法迭代数为K = 10000 ;迭代误差阈值为e = 10_4;松弛因子入 =0? 7 ;Pareto 权值 y = 0? 6〇
[0120] 图6是本发明实施方式中，在特定的参数关系下，MM0雷达信号在多个频带区域 抑制示意图。仿真参数：码长N = 100 ;阵元数M = 1 ;脉冲持续时间Tdming= 200 y s ;采样 频率fsample= 1000KHZ ;算法迭代数为K = 10000 ;迭代误差阈值为e = 10_4;松弛因子入 =0? 7 ;Pareto 权值 y = 0? 6〇
【具体实施方式】
[0121] 本实施方式所述的一种多输入多输出雷达波形设计方法包括以下步骤：
[0122] 步骤1 :根据雷达场景中强散射体与待测目标的相对位置，预估自相关旁瓣抑制 模糊区间，进而构造相应的目标函数：
[0123] 设M个发射阵元的MM0雷达系统，其第m个阵元窄带恒模相位编码波形序列可表 不为
，N表征信号码长，itnG[0, 2 31 ] 为第nth时刻编码相位；因接收端滤波器可视为发射波形的相关函数，假定感兴趣的待测目 标在第pth距离单 元；那么在第Ith单元的强散射体势必要对pth单元目标造成影响；对于第 pth距离单元，为最小化临近距离单元中强散射体产生的距离旁瓣对它的干扰，应使发射波 形1的自相关函数a(Xm)满足如下条件：
[0125]其中，a(xm) = [a0(xm)…an-JxJ0a_N+1 (xm) ???a< (xm) ]T表示自相关函数序 列；毛表示理想波形说={± |z_p||zGZJVN-l,-N+l}表示旁瓣干扰区间，表示距 离向分布的强散射体位置集合；
[0126] 对于强散射体遮蔽干扰问题，发射波形Xm应使得目标函数J(Xm)最小化，即：
[0128]其中N? | |表示欧式范数，视为距离度量函数；由Parseval等价性可知：
[0130]其中，CA>XA>表示离散傅里叶变换矩阵，
；表示自相 关函数的傅里叶变换，即功率谱；如果令大和分别表征理想波形和设计波形的频谱， 即
C为扩展矩阵，其通过扩充零元素使得矩阵运算时维度一致；那么易得：
，其中?表示按元素的Hadamard积，（?r表示取共轭操作；目标函数 写为：
[0132] 步骤2:分析MM0雷达波形正交性约束，构造满足正交性约束的目标函数：
[0133] 对于MM0发射波形矩阵：
而言，波形之间正交性体现在互相关 函数上；第nth时刻的发射波形相关矩阵Rn(X)的非主对角线部分尽可能逼近零，主对角元 素为各阵元发射波形的自相关表示：
[0137] 其中diag(a"(；e))表示nth时刻由理想波形矩阵x的自相关形成的对角矩阵；U(n) 表示偏移矩阵，I表示单位矩阵；
[0138]为解决不同阵元波形间互扰构造目标函数J2 (X)，同样考虑Parseval等价性，得到 如下：
[0140] 其中| | ? | |F表示矩阵的Frobenius范数；
[0141] 步骤3:根据场景先验信息预估频域干扰模糊频带区间，进而构造相应目标函数：
[0142] 假定可能存在电磁干扰的频段区间集合为
，那么发射波形频谱应 在此区间内设为阻带，由傅里叶变换的性质可知，频谱为：
[0144]其中：NS表示相应频段数目，巧表示单位离散傅里叶变换矩阵，且有
为使得频谱离散采样足够密以表征区间信息，设定及= 2W;由上式可知，如果从离散傅里 叶变换矩阵的行空间中抽取特定行来构成矩阵P，且所抽取行序号与D在频谱W中频 段相对应，用向量必表征&中稀疏凹口序列且满足
fn为频 谱W中频率的采样，从而使得：
[0146]只要# f幅度足够低就满足频谱稀疏特性，即目标函数J3(X)最小化，从而达到任 务要求：
[0148] 步骤4:构造恒模相位编码波形约束条件，即最大化发射功率要求恒模约束，也 即：
[0149] |Xn'J= 1，n= 1，…N.m= 1，??，M
[0150] 步骤5:构造松弛交替投影算法框架：
[0151] 假定集合A代表约束集合（即满足恒模约束），集合B代表目标集合（即满足子任 务要求），构造变量X，y到集合A,B的投影ProjA (x)和projB (y)，如下所示：
[0154] 其中d(i,j)= | | i-j| |表示变量i,j的距离函数。
[0155] 那么交替投影机制可简述为：迭代映射T:A-B-A。
[0156] 然后，引入松弛交替投影框架：
[0157] T5(.)；Xx=ProjA(x+A? 8 (x) ? (ProjA(ProjB(x))-x))
[0158]xk+1=Ts(.)，Xxk
[0159] 上式中，X表示松弛算子；S(?)表示加速因子，有如下不同形式：
[0162] 此算法框架命名为RISAAP。
[0163] 步骤6、上述步骤1、2、3分别完成了三个子任务目标函数的构建，以适应雷达多任 务场景的时变性，某时刻场景可能涉及上述子任务目标中的一个或多个，根据步骤5所提 松弛交替投影算法框架求解波形设计，分别给出如下三种波形优化输出方式：
[0164] (1)表示满足抑制旁瓣遮蔽效应的算法步骤；
[0165] (2)表示获得波形正交性的算法步骤；
[0166] (3)表示联合抑制特定频带干扰及获得良好相关特性的算法步骤；
[0167](1)对于抑制旁瓣遮蔽效应的算法，过程如下；
[0168] 以满足步骤1目标函数条件的集合为目标集合B;以满足步骤4条件集合为约束 集合A;(?)(1/2)表示按元素的平方根操作，返回值维度与其参数一致；
[0169] 那么获得较好时域自相关旁瓣抑制性能的恒模波形设计方法可归纳为：
[0170] a)初始化波形序列：x= [x(l). . .x(n)x(N)]TGA，其中x(?) =ew", itnG[0, 2it]为随机相位；
[0171]b)迭代优化过程：
[0172] 由公式
[0174]获得fGB，
[0175] 再由
[0177]获得xEA，
[0178] 经公式
[0179] Ts(.)，入x=ProjA(x+人? 8 (x) ? (ProjA(ProjB(x))_x))
[0181] 获得；并更新x= x;
[0182]c)相邻两次迭代产生差值小于给定误差e或算法迭代次数大于给定总数K，则 算法停止并输出，否则继续b)，c);
[0183] (2)对于获得波形正交性的算法设计，过程如下：
[0184] 以满足步骤2目标函数条件的集合，即为目标集合B;以满足步骤4条件的集合为 约束集合A;那么获得较好时域相关旁瓣抑制性能的恒模波形设计方法可归纳为：
[0185]a)初始化波形矩阵：X= [Xl...xM]GA，其中i
，步 n，me[0, 2n]为 随机相位；
[0186]b)迭代优化过程：
[0187] 由公式
[0189]获得FeB,
[0190] 再由
[0192]获得XGA，
[0193] 经公式
[0194] T5(.)；Ax=ProjA(x+A?8 (x)?(ProjA(ProjB(x))-x))
[0195] ^ -T&'pJ
[0196] 获得X并更新
[0197]c)若相邻两次迭代产生1差值小于给定误差e或算法迭代次数大于给定总数 K，则算法停止并输出，否则继续b)，c);
[0198] (3)对于联合抑制特定频带干扰及获得良好相关特性的算法设计，过程如下：
[0199] 以满足步骤1、2目标函数条件的集合，即为目标集合B1;以满足步骤3的集合，即 为目标集合B2;以满足步骤4条件的集合为约束集合A;那么获得较好时域相关旁瓣抑制性 能及频域特性的恒模波形设计方法可归纳为：
[0200] a)初始化波形矩阵：X = [Xl. ? ? xM] G A，其中
，步n，me[0, 2 n ]为 随机相位；
[0201] b)迭代优化过程：
[0202] 由公式
[0204] 获得巧e马，
[0205] 由公式
[0207]获得 e/i:，由
[0209]获得 FeB,
[0210] 再由
[0212] 获得XGA，
[0213] 经公式
[0214] Ts(.)，入x=ProjA(x+人? 8 (x) ? (ProjA(ProjB(x))_x))
[0216] 获得X并更新
[0217]c)若相邻两次迭代产生X差值小于给定误差e或算法迭代次数大于给定总数 K，则算法停止并输出，否则继续b)，c)。
[0218] 所述步骤1中时域干扰区间确定可根据先验信息，或通过事先发射传统波形以确 定大致区间。所述步骤2中多输入多输出雷达可包含多个阵元。所述步骤3中频域干扰区 间的确定可根据先验信息、或参考通用电子手册确定大致区间。所述步骤5中松弛算子入 的取值范围是：〇〈入〈2;Parteo权值y的取值范围是：0〈y〈1。
[0219] 如图1~6所示，结合步骤6中所列（1)、⑵、（3)方式，给出了基于本发明的技术 方案，从而可设计出任意阵元数及任意码长的具有良好时频抗干扰性能的离散恒模相位编 码信号的实施例：
[0220] 1.单阵元150个子码的抗单、多旁瓣干扰抑制波形设计：
[0221] 利用本技术方案中步骤6中（1)所述，假定强散射体可能位于单旁瓣干扰区 间％= [1:30]及多旁瓣干扰区间Q1= [20:30] U [40:60];初始化发射波形序列为：
，步ne [0,231]为随机相位且由随机数 产生；码长N = 150,算法迭代数K = 10000,迭代误差阈值e = 1(T4,松弛因子A = 〇. 7。 基于步骤6中（1)方式所列公式求解，直到相邻两次迭代产生分差值小于e或算法迭代 数大于K，则算法停止并输出。方案所提RISAAP算法能够在单旁瓣区间及多旁瓣区间获得 313dB和230dB的优化效果且耗时仅为0. 08s和0. 16s，其良好的抗干扰性能如表1、图1及 图2所示。
[0222] 2. 3阵元64码长阵元波形间正交性优化设计：
[0223] 利用本技术方案步骤6中（2)所述，初始化发射波形矩阵为：
其中
，[0,2 3!]为随机相位且由随机数产生；取码长N = 64,阵元数 M = 3 ;算法迭代数K = 10000 ;迭代误差阈值e = KT4;松弛因子A = 〇. 7。基于步骤6 中（2)方式所列公式求解，直到相邻两次迭代产生i差值小于e或算法迭代数大于K，则 算法停止并输出。作图仅以第一个阵元波形为例，RISAAP算法优化自相关幅度和互相关幅 度（第一、三阵元波形）仿真结果见图3、4,各算法随码长增加耗时比较见表二。
[0224 ] 3.单阵元100码长的联合抗旁瓣干扰及抑制特定频带电磁干扰的波形设计：
[0225] 利用本技术方案中步骤6中（3)所述，假设电磁干扰可能位于单频段区间 [200, 300] (KHz)和多频段区间[200, 300] U [500, 600] (KHz)，初始化发射波形矩阵为：
，步[0,231]为随机相位且由随机数产生；取 码长N = 100,阵元数M = 1 ;算法迭代数K = 10000 ;迭代误差阈值e = 1(T4;松弛因子入 =0. 7, Pareto权值y = 0. 6。基于步骤6中（3)方式所列公式求解，直到相邻两次迭代 产生差值小于e或算法迭代数大于K，则算法停止并输出。方案所提RISAAP算法能够 在单频段干扰区间及多频段干扰区间获得约14dB和13dB的优化效果且耗时仅为0. 06s和 0. 07s (如表3所示），其良好的抗干扰性能图5、6所示。
[0226] 表1不同算法抑制特定区间旁瓣比较
[0228] 表2不同序列长度算法耗时比较
[0229]
[0230] 表3不同频段干扰区间算法性能比较
[0232] 以上所述，仅为本发明较佳的【具体实施方式】，但本发明的保护范围并不局限于此， 任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，根据本发明的技术方案及其 发明构思加以等同替换或改变，都应涵盖在本发明的保护范围之内。
【主权项】
1. 一种多输入多输出雷达波形设计方法，其特征是包括以下步骤： 步骤1 :根据雷达场景中强散射体与待测目标的相对位置，预估自相关旁瓣抑制模糊 区间，进而构造相应的目标函数： 设M个发射阵元的MIMO雷达系统，其第m个阵元窄带恒模相位编码波形序列可表示为 Xm =[?)... Xm(W) xm(iV)]T e Cam，Xw(W) = ejv" , N 表征信号码长，Φ ne [〇, 2 π ]为第 nth时刻编码相位；因接收端滤波器可视为发射波形的相关函数，假定感兴趣的待测目标在 第Pth距离单元；那么在第Ith单元的强散射体势必要对P th单元目标造成影响；对于第P th距离单元，为最小化临近距离单元中强散射体产生的距离旁瓣对它的干扰，应使发射波形 Xm的自相关函数a (Xm)满足如下条件：其中，a (xm) = [ a Q(Xm)…a N_i (Xm) 〇 a _N+1 (xm)…α (xm) ]T表示自相关函数序列；· 表示理想波形；Qi={± |z-p| |z e Z1UiN-I,-Ν+1}表示旁瓣干扰区间，Z1表示距离向分布的 强散射体位置集合； 对于强散射体遮蔽干扰问题，发射波形^应使得目标函数J(Xm)最小化，即： min·/(λ:," H|a(xm)-Ct(^)II2 Xm 其中N · I I表示欧式范数，视为距离度量函数；由Parseval等价性可知：其中，表示离散傅里叶变换矩阵，# = 2#; FxrOt(Xm)eC2wxl表示自相关 函数的傅里叶变换，即功率谱；如果令/"和乙分别表征理想波形和设计波形的频谱，即 九= i^cT·^，C为扩展矩阵，其通过扩充零元素使得矩阵运算时维度一致；那么易得： Αα(毛,）= ./；"〇_/；二其中Θ表示按元素的Hadamard积，（·）*表示取共轭操作；目标函数 写为：步骤2 :分析MMO雷达波形正交性约束，构造满足正交性约束的目标函数： 对于MMO发射波形矩阵X=K..xM]e Cyxw而言，波形之间正交性体现在互相关函数 上；第nth时刻的发射波形相关矩阵Rn(X)的非主对角线部分尽可能逼近零，主对角元素为 各阵元发射波形的自相关表示：其中diag(a"(；e))表示nth时刻由理想波形矩阵X的自相关形成的对角矩阵；U(n)表示 偏移矩阵，I表示单位矩阵； 为解决不同阵元波形间互扰构造目标函数J2(X),同样考虑Parseval等价性，得到如 下：其中I I · I If表示矩阵的Frobenius范数； 步骤3 :根据场景先验信息预估频域干扰模糊频带区间，进而构造相应目标函数： 假定可能存在电磁干扰的频段区间集合为，那么发射波形频谱应在此 区间内设为阻带，由傅里叶变换的性质可知，频谱为： 其中：NS表示相应频段数目，G表示单位离散傅里叶变换矩阵，且有为使 得频谱离散采样足够密以表征区间信息，设定A = 27V;由上式可知，如果从离散傅里叶变 换矩阵的行空间中抽取特定行来构成矩阵P，且所抽取行序号与Ω在频谱W中频段& 相对应，用向量命表征Ω中稀疏凹口序列且满足也=l,/"e 亡，fn为频谱W 中频率的采样，从而使得： W = PC 'X 只要#幅度足够低就满足频谱稀疏特性，即目标函数J3(X)最小化，从而达到任务要 求： JAX) = m\nPC'X J X 步骤4 :构造恒模相位编码波形约束条件，即最大化发射功率要求恒模约束，也即： |Xn，J=l，n=l, ...N.m=l，· ·，M 步骤5 :构造松弛交替投影算法框架： 假定集合A代表约束集合（即满足恒模约束)，集合B代表目标集合（即满足子任务要 求)，构造变量X，y到集合A, B的投影Pro jA (X)和pro jB (y)，如下所示：其中d(i，j) = | |i-j| I表示变量i，j的距离函数。 那么交替投影机制可简述为：迭代映射T:A - B - A。 然后，引入松弛交替投影框架： Ts(.)，x=Pr〇jA(x+X · δ (X) · (ProjA(ProjB(x))-x)) Xk+l T δ ( . )，λ Xk 上式中，λ表示松弛算子；δ ( ·)表示加速因子，有如下不同形式：此算法框架命名为RISAAP。 步骤6、上述步骤1、2、3分别完成了三个子任务目标函数的构建，以适应雷达多任务场 景的时变性，某时刻场景可能涉及上述子任务目标中的一个或多个，根据步骤5所提松弛 交替投影算法框架求解波形设计，分别给出如下三种波形优化输出方式： (1) 表示满足抑制旁瓣遮蔽效应的算法步骤； (2) 表示获得波形正交性的算法步骤； (3) 表示联合抑制特定频带干扰及获得良好相关特性的算法步骤； (1)对于抑制旁瓣遮蔽效应的算法，过程如下； 以满足步骤1目标函数条件的集合为目标集合B ;以满足步骤4条件集合为约束集合 A ; ( · )(1/2)表示按元素的平方根操作，返回值维度与其参数一致； 那么获得较好时域自相关旁瓣抑制性能的恒模波形设计方法可归纳为： a) 初始化波形序列：x=[x(l)··· x(n) x(N)]TeA，其中 X(?) = ejv", Φηε[〇,2π] 为随机相位； b) 迭代优化过程：由公式 获得f G B， 再由获得X e A， 经公式获得4并更新.V = C)相邻两次迭代产生i差值小于给定误差ε或算法迭代次数大于给定总数K，则算法 停止并输出，否则继续b)，c); (2) 对于获得波形正交性的算法设计，过程如下： 以满足步骤2目标函数条件的集合，即为目标集合B ;以满足步骤4条件的集合为约束 集合A ;那么获得较好时域相关旁瓣抑制性能的恒模波形设计方法可归纳为： a) 初始化波形矩阵A=Iix1... xM] e A，其中XmO?) = , Φ[〇, 2 π ]为随机相 位； b) 迭代优化过程：获得X e A，经公式Ts(.)，x=Pr〇jA(x+X · δ (X) · (ProjA(ProjB(x))-x))λ - 由公式 获得FeB, 再由 获得X并更新X=义； C)若相邻两次迭代产生i差值小于给定误差ε或算法迭代次数大于给定总数K，则算 法停止并输出，否则继续b)，c); (3) 对于联合抑制特定频带干扰及获得良好相关特性的算法设计，过程如下： 以满足步骤1、2目标函数条件的集合，即为目标集合B1;以满足步骤3的集合，即为目 标集合B2;以满足步骤4条件的集合为约束集合A ;那么获得较好时域相关旁瓣抑制性能及 频域特性的恒模波形设计方法可归纳为： a)初始化波形矩阵：X=[Xl... XM] e A，其中 Xw(W) = ejv"' Φn，me [〇, 2 JT ]为随机相 位； b) 迭代优化过程：由公式 获得F1 ￡尽， 由公式 获得由 获得F e/i , 再由 获得X e A， 经公式获得1并更新 c) 若相邻两次迭代产生义差值小于给定误差ε或算法迭代次数大于给定总数K，则算 法停止并输出，否则继续b)，c)。2. 根据权利要求1所述的一种多输入多输出雷达波形设计方法，其特征是，所述步骤1 中时域干扰区间确定可根据先验信息，或通过事先发射传统波形以确定大致区间。3. 根据权利要求2所述的一种多输入多输出雷达波形设计方法，其特征是，所述步骤2 中多输入多输出雷达可包含多个阵元。4. 根据权利要求3所述的一种多输入多输出雷达波形设计方法，其特征是，所述步骤3 中频域干扰区间的确定可根据先验信息、或参考通用电子手册确定大致区间。5. 根据权利要求1、2、3或4所述的一种多输入多输出雷达波形设计方法，其特征是，所 述步骤5中松弛算子λ的取值范围是：〇〈λ〈2 ;Parteo权值γ的取值范围是：0〈γ〈1。
【专利摘要】一种多输入多输出雷达波形设计方法，属于雷达通信技术领域。本发明的目的是提供一种具有更低的相关旁瓣和频谱抑制深度，且效率高、耗时少、并具有较高的鲁棒性、具有良好时频抗干扰性能的设计方法。根据雷达场景中强散射体与待测目标的相对位置，预估自相关旁瓣抑制模糊区间，进而构造相应的目标函数；分析MIMO雷达波形正交性约束，构造满足正交性约束的目标函数；根据场景先验信息预估频域干扰模糊频带区间，进而构造相应目标函数；构造恒模相位编码波形约束条件；构造松弛交替投影算法框架；根据所提松弛交替投影算法框架求解波形设计，给出三种波形优化输出方式。采用松弛交替投影恒模波形编码设计，可使MIMO雷达具有更好的检测性能。
【IPC分类】G01S13/00
【公开号】CN104898113
【申请号】CN201510346063
【发明人】赵宜楠, 赵占锋, 冯翔, 周志权
【申请人】哈尔滨工业大学
【公开日】2015年9月9日
【申请日】2015年6月19日
【技术领域】
[0001] 本发明属于雷达通信技术领域，特别涉及MM0雷达满足时域频域任务需求的雷 达波形设计方法。
【背景技术】
[0002] 当今军事或民用需求使得雷达任务场景日趋复杂。为应对复杂变化的场景， 雷达自适应性显得尤为重要。其中，波形分集可视为自适应性的一种体现且成为目前 石开宄的热点（见文献：Cognitive radar: a way of the future, Haykin S. ；Signal Processing Magazine, IEEE, 2006, 23(1):30-40 ；Adaptive radar phase-coded waveform design, J. D. Zhang, X. H. Zhu, H. Q. Wang ；Electron. Lett. 2009, 45 (20) 1052 - 1053 ；Spatial diversity in radars-modeIs and detection performance, Eran, F. , Alexander, H. , &Ri ck，S.B. ; IEEE Transactions on Signal Processing, 2006, 54 (3) ,823 - 838)。无论是相 控阵雷达还是MM0雷达，波形分集既能提高目标识别能力又能提高分辨力，已然成为提高 雷达自适应性的重要切入点。以MIM0雷达为例，考虑到接收端滤波器可视为波形的相关函 数，结合波形设计与复杂任务场景分类：避免强散射体距离旁瓣遮蔽、防止不同阵元间波形 互扰、躲避特定频带电磁干扰，这三者即成为波形设计的基本要求。此外，为充分利用发射 机功率，恒模特性亦成为客观要求。
[0003] 避免强散射体遮蔽效应需要发射波形具备较低的自相关旁瓣；消除不同阵元相 互影响要求发射波形正交化；抑制特定频带电磁干扰要求发射波形频谱具备稀疏特性； 另外，恒模约束往往使得设计问题非凸。以上这些要求使得波形设计问题较为困难。针 对恒模波形提出的循环算法族（Cyclic Algorithm New，CAN)，其思路是基于相关特性 与功率谱之间的傅里叶对应关系，逐步迭代以优化所选目标函数（见文献：Waveform design for active sensing systems:a computational approach, He H, Li J, Stoica P. ;Cambridge University Press, 2012)。但该算法耗时过长、相关旁瓣电平幅度较差， 难以满足工程实践。但其给后续研宄提供了众多启示，以交替投影为机制的改进循环算 法框架，获得了比CAN算法族更明显的效果（见文献：零自相关区相位编码波形设计， 李风从，赵宜楠，乔晓林；电子学报，2013, 12:2499-2502;基于秩亏傅里叶变换的交替 投影编码波形设计，赵宜楠，李风从，王军，乔晓林；电子学报，2014, 06:1216-1219; Computational design of optimal waveforms for MIM0 radar via multi-dimensional iterative spectral approximation, Zhao Y N, Li F C, Zhang T ;Multidm syst sign process，11045-01400288-1，2014);但算法初始化导致非凸问题求解稳定性较差。除此 之外，利用遗传算法对MIM0雷达正交波形进行设计（见文献：Polyphase orthogonal code design for MIM0 radar systems, Liu B, He Z, Zeng J ；2006 CIE international conference on radar, 2006:1-4);但该算法在迭代优化的后期往往出现群体多样性较差， 导致局部最优解而陷入停滞，另外算法规模也成为限制其效率的因素。但针对实际MIM0雷 达场景，此类算法可能耗时过长且难以收敛。利用序列二次规划与加权迭代方法求解波形 设计问题，权值步长的搜索机制难以自适应；此类基于梯度矩阵的算法在多阵元情形计算 复杂度较高，且可能出现高阶矩阵奇异而难以收敛（见文献：一种空时联合优化的MMO雷 达波形设计方法，王旭，周生华，刘宏伟，保铮；西安电子科技大学学报，2014, 03:41-48)。
[0004] 因恒模导致的非凸问题往往收敛缓慢且难以得到全局最优解，而交替投 影算法对于凸范数求解亦可能存在局部停滞而难以收敛（见文献：〇n pro jection algorithms for solving convex feasibility problems,Bauschke H H,Borwein J M,SIAM review,1996, 38(3):367-426 ；An alternating projection that does not converge in norm, Hundal H S，Nonlinear Analysis:Theory,Methods&Applicatio ns，2004, 57 (1) : 35-61)，而松弛交替投影算法作为一种高效且耗时较少算法（见文献： Relaxed alternating projection methods, A. Cegielski, A. Suchocka, SIAM Journal on Optimization, vol. 19(3)，2008, pp. 1093-1106)，结合迭代谱近似思想可用来解决以恒模 波形设计为代表的非凸问题。
【发明内容】
[0005] 本发明的目的是提供一种具有更低的相关旁瓣和频谱抑制深度，且效率高、耗时 少、可有效避免局部最优，并具有较高的鲁棒性、具有良好时频抗干扰性能的多输入多输出 雷达波形设计方法。
[0006] 本发明为解决上述技术问题采取的技术方案是：
[0007] -种多输入多输出雷达波形设计方法，包括以下步骤：
[0008] 步骤1:根据雷达场景中强散射体与待测目标的相对位置，预估自相关旁瓣抑制 模糊区间，进而构造相应的目标函数：
[0009] 设M个发射阵元的MIM0雷达系统，其第m个阵元窄带恒模相位编码波形序列可表 不为
，N表征信号码长，itnG [0, 2 JT ] 为第nth时刻编码相位；因接收端滤波器可视为发射波形的相关函数，假定感兴趣的待测目 标在第pth距离单元；那么在第Ith单元的强散射体势必要对pth单元目标造成影响；对于第 p th距离单元，为最小化临近距离单元中强散射体产生的距离旁瓣对它的干扰，应使发射波 形1的自相关函数a (Xm)满足如下条件：
[0011]其中，a (xm) = [ a 0(xm)…an-JxJ0 a _N+1 (xm) ??? a < (xm) ]T表示自相关函数序 列；毛，表示理想波形；Qi= {± |z-p| |z G ZiUm-L-N+l}表示旁瓣干扰区间，&表示距离 向分布的强散射体位置集合；
[0012] 对于强散射体遮蔽干扰问题，发射波形Xm应使得目标函数J(Xm)最小化，即：
[0014]其中N ? | |表示欧式范数，视为距离度量函数；由Parseval等价性可知：
[0016]其中，
表示离散傅里叶变换矩阵，
表示自相 关函数的傅里叶变换，即功率谱；如果令?/" :和fm分别表征理想波形和设计波形的频谱， 即
，C为扩展矩阵，其通过扩充零元素使得矩阵运算时维度一致；那么易得：
I其中?表示按元素的Hadamard积，（? r表示取共轭操作；目标函数 写为：
[0018]步骤2 :分析MM0雷达波形正交性约束，构造满足正交性约束的目标函数：
[0019] 对于MM0发射波形矩阵
而言，波形之间正交性体现在互相关 函数上；第nth时刻的发射波形相关矩阵Rn(X)的非主对角线部分尽可能逼近零，主对角元 素为各阵元发射波形的自相关表示：
[0023] 其中diag(a"(l))表示nth时刻由理想波形矩阵.文的自相关形成的对角矩阵；U(n) 表示偏移矩阵，I表示单位矩阵；
[0024] 为解决不同阵元波形间互扰构造目标函数J2 (X)，同样考虑Parseval等价性，得到 如下：
[0026] 其中| | ? | |F表示矩阵的Frobenius范数；
[0027] 步骤3 :根据场景先验信息预估频域干扰模糊频带区间，进而构造相应目标函数：
[0028] 假定可能存在电磁干扰的频段区间集合为
，那么发射波形频谱应 在此区间内设为阻带，由傅里叶变换的性质可知，频谱为：
[0030]其中：NS表示相应频段数目，表示单位离散傅里叶变换矩阵，且有 为使得频谱离散采样足够密以表征区间信息，设定及= 2#;由上式可知，如果从离散傅里 叶变换矩阵& I的行空间中抽取特定行来构成矩阵P，且所抽取行序号与D在频谱W中频 段&相对应，用向量?斤表征&中稀疏凹口序列且满足
fn为频 谱W中频率的采样，从而使得：
[0032]只要# [幅度足够低就满足频谱稀疏特性，即目标函数J3(X)最小化，从而达到任 务要求：
[0034] 步骤4:构造恒模相位编码波形约束条件，即最大化发射功率要求恒模约束，也 即：
[0035] |Xn'J = 1，n = 1，…N. m = 1，? ?，M
[0036] 步骤5 :构造松弛交替投影算法框架：
[0037]假定集合A代表约束集合（即满足恒模约束），集合B代表目标集合（即满足子任 务要求），构造变量X，y到集合A,B的投影ProjA(x)和projB(y)，如下所示：
[0040] 其中d(i，j)= I Ii-jI I表示变量i，j的距离函数。
[0041] 那么交替投影机制可简述为：迭代映射T:A - B - A。
[0042] 然后，引入松弛交替投影框架：
[0043] Ts(.)，入x = ProjA(x+人? 8 (x) ? (ProjA(ProjB(x))-x))
[0044]xk+1= Ts(.)，Xxk
[0045] 上式中，A表示松弛算子；S (?)表示加速因子，有如下不同形式：
[0048] 此算法框架命名为RISAAP。
[0049] 步骤6、上述步骤1、2、3分别完成了三个子任务目标函数的构建，以适应雷达多任 务场景的时变性，某时刻场景可能涉及上述子任务目标中的一个或多个，根据步骤5所提 松弛交替投影算法框架求解波形设计，分别给出如下三种波形优化输出方式：
[0050] (1)表示满足抑制旁瓣遮蔽效应的算法步骤；
[0051] (2)表示获得波形正交性的算法步骤；
[0052] (3)表示联合抑制特定频带干扰及获得良好相关特性的算法步骤；
[0053] (1)对于抑制旁瓣遮蔽效应的算法，过程如下；
[0054] 以满足步骤1目标函数条件的集合为目标集合B ;以满足步骤4条件集合为约束 集合A ;(?) (1/2)表示按元素的平方根操作，返回值维度与其参数一致；
[0055] 那么获得较好时域自相关旁瓣抑制性能的恒模波形设计方法可归纳为：
[0056] a)初始化波形序列：x=[x(l). . . x(n) x(N)]TG A，其中
itnG [0, 2it]为随机相位；
[0057] b)迭代优化过程：
[0058] 由公式
[0060]获得f G B，
[0061] 再由
[0063]获得x G A，
[0064] 经公式
[0066]获得i并更新X = X;
[0067] c)相邻两次迭代产生J差值小于给定误差e或算法迭代次数大于给定总数K，则 算法停止并输出，否则继续b)，c);
[0068] (2)对于获得波形正交性的算法设计，过程如下：
[0069] 以满足步骤2目标函数条件的集合，即为目标集合B ;以满足步骤4条件的集合为 约束集合A ;那么获得较好时域相关旁瓣抑制性能的恒模波形设计方法可归纳为：
[0070] a)初始化波形矩阵：X= [Xl.? ? xM] G A，其中
，步[0, 2 n]为 随机相位；
[0071] b)迭代优化过程：
[0072] 由公式
[0073]
[0074] 获得 Fd
[0075] 再由
[0077] 获得 X G A，
[0078] 经公式
[0079] Ts(.)，入x = ProjA(x+人? 8 (x) ? (ProjA(ProjB(x))-x))
[0081] 获得X并更新X二X ;
[0082] c)若相邻两次迭代产生f差值小于给定误差e或算法迭代次数大于给定总数K， 则算法停止并输出，否则继续b)，c);
[0083] (3)对于联合抑制特定频带干扰及获得良好相关特性的算法设计，过程如下：
[0084] 以满足步骤1、2目标函数条件的集合，即为目标集合B1;以满足步骤3的集合，即 为目标集合B 2;以满足步骤4条件的集合为约束集合A ;那么获得较好时域相关旁瓣抑制性 能及频域特性的恒模波形设计方法可归纳为：
[0085] a)初始化波形矩阵：X = [Xl. . . xM] G A，其中
，步 n，me [0, 2 n ]为 随机相位；
[0086] b)迭代优化过程：
[0087] 由公式
[0089] 获得片e^，
[0090]由公式
[0092] 获得尽e /七，由
[0094]获得F e B，
[0095] 再由
[0097]获得 X e A，
[0098] 经公式
[0099] Ts(.)，入x = ProjA(x+人? 8 (x) ? (ProjA(ProjB(x))-x))
[0101] 获得1并更新x =
[0102] c)若相邻两次迭代产生差值小于给定误差e或算法迭代次数大于给定总数 K，则算法停止并输出，否则继续b)，c)。
[0103] 所述步骤1中时域干扰区间确定可根据先验信息，或通过事先发射传统波形以确 定大致区间。
[0104] 所述步骤2中多输入多输出雷达可包含多个阵元。
[0105] 所述步骤3中频域干扰区间的确定可根据先验信息、或参考通用电子手册确定大 致区间。
[0106] 所述步骤5中松弛算子A的取值范围是：〇〈 A〈2 ;Parteo权值y的取值范围是： 0〈丫〈1〇
[0107] 本发明的有益效果是：
[0108] 本发明采用松弛投影算法来设计多输入多输出雷达波形，相比采用现有的波形在 线设计方法，本发明具有更低的相关旁瓣和频谱抑制深度，且效率高、耗时少、可有效避免 局部最优，并具有较高的鲁棒性。应用本发明的松弛交替投影恒模波形编码设计，可使MIM0 雷达具有更好的检测性能。
[0109] 本发明所设计波形可有效避免强散射体遮蔽干扰、阵元波形间互扰及特定频带干 扰，更有利于目标的检测与跟踪，其优点如下：
[0110] (1)就其采用的恒模相位编码产生方式而言，本发明提出的MIM0雷达波形具有精 确的表达式及具体的步骤优化程序，避免了现有技术中涉及的基于梯度的高维非线性搜索 问题。
[0111] (2)就理想波形的构建而言，采用设计波形与其渐进拟合的方式，理想波形建立在 迭代设计的基础之上，避免了先期考虑的过多因素。
[0112] (3)就波形优化效果而言，本发明提出的MIM0雷达波形设计方法耗时少、抑制自 相关旁瓣及特定频带区间干扰远强于Li Jian等人所提出的CAN循环算法；更便于在线实 时设计。
[0113] 其他性能指标（归一化自相关幅值电平：ACL ;阻带峰值功率：PSP ;)与下述已有 技术相比，完全优于已有技术的相应指标。
[0114] 已有技术：以文献(waveform design for active sensing systems:a computational approach, He H, Li J, Stoica P. Cambridge University Press, 2012) 提出的CAN、WeCAN、SCAN算法族，包括Multi-CAN、WeCAN、SCAN ;及以文献（零自相关 区相位编码波形设计，李风从，赵宜楠，乔晓林；电子学报，2013, 12:2499-2502;基于 秩亏傅里叶变换的交替投影编码波形设计，赵宜楠，李风从，王军，乔晓林；电子学 报，2014,06:1216-1219 ;Computational design of optimal waveforms for MIM0 radar via multi-dimensional iterative spectral approximation, Zhao Y N, Li F C, Zhang T ； Multidm syst sign process, 11045-01400288-1，2014))提出的 ISAA、MDISAA 等。但这些 技术耗时过长，波形优化设计效果欠佳，不适合工程应用。
【附图说明】
[0115] 图1是本发明实施方式中，在特定的参数关系下，MM0雷达阵元信号在单个特定 区域自相关旁瓣抑制示意图。仿真参数：码长N = 150 ;算法迭代数K = 10000 ;迭代误差 阈值为e = 10、松弛因子A = 〇. 7。
[0116] 图2是本发明实施方式中，在特定的参数关系下，MM0雷达阵元信号在多个特定 区域自相关旁瓣抑制示意图。仿真参数：码长N = 150 ;算法迭代数K = 10000 ;迭代误差 阈值为e = 10、松弛因子A = 〇. 7。
[0117] 图3是本发明实施方式中，在特定的参数关系下，MM0雷达阵元间信号自相关旁 瓣抑制示意图。仿真参数：码长N = 64 ;阵元数M = 3 ;算法迭代数K = 10000 ;迭代误差阈 值e = 1(T4;松弛因子A = 〇? 7。
[0118] 图4是本发明实施方式中，在特定的参数关系下，MM0雷达阵元间信号互相关旁 瓣抑制示意图。仿真参数：码长N = 64 ;阵元数M = 3 ;算法迭代数为K = 10000 ;迭代误差 阈值为e = 10、松弛因子A = 〇. 7。
[0119] 图5是本发明实施方式中，在特定的参数关系下，MM0雷达信号在单个频带区域 抑制示意图。仿真参数：码长N = 100 ;阵元数M = 1 ;脉冲持续时间Tdming= 200 y s ;采样 频率fsample= 1000KHZ ;算法迭代数为K = 10000 ;迭代误差阈值为e = 10_4;松弛因子入 =0? 7 ;Pareto 权值 y = 0? 6〇
[0120] 图6是本发明实施方式中，在特定的参数关系下，MM0雷达信号在多个频带区域 抑制示意图。仿真参数：码长N = 100 ;阵元数M = 1 ;脉冲持续时间Tdming= 200 y s ;采样 频率fsample= 1000KHZ ;算法迭代数为K = 10000 ;迭代误差阈值为e = 10_4;松弛因子入 =0? 7 ;Pareto 权值 y = 0? 6〇
【具体实施方式】
[0121] 本实施方式所述的一种多输入多输出雷达波形设计方法包括以下步骤：
[0122] 步骤1 :根据雷达场景中强散射体与待测目标的相对位置，预估自相关旁瓣抑制 模糊区间，进而构造相应的目标函数：
[0123] 设M个发射阵元的MM0雷达系统，其第m个阵元窄带恒模相位编码波形序列可表 不为
，N表征信号码长，itnG[0, 2 31 ] 为第nth时刻编码相位；因接收端滤波器可视为发射波形的相关函数，假定感兴趣的待测目 标在第pth距离单 元；那么在第Ith单元的强散射体势必要对pth单元目标造成影响；对于第 pth距离单元，为最小化临近距离单元中强散射体产生的距离旁瓣对它的干扰，应使发射波 形1的自相关函数a(Xm)满足如下条件：
[0125]其中，a(xm) = [a0(xm)…an-JxJ0a_N+1 (xm) ???a< (xm) ]T表示自相关函数序 列；毛表示理想波形说={± |z_p||zGZJVN-l,-N+l}表示旁瓣干扰区间，表示距 离向分布的强散射体位置集合；
[0126] 对于强散射体遮蔽干扰问题，发射波形Xm应使得目标函数J(Xm)最小化，即：
[0128]其中N? | |表示欧式范数，视为距离度量函数；由Parseval等价性可知：
[0130]其中，CA>XA>表示离散傅里叶变换矩阵，
；表示自相 关函数的傅里叶变换，即功率谱；如果令大和分别表征理想波形和设计波形的频谱， 即
C为扩展矩阵，其通过扩充零元素使得矩阵运算时维度一致；那么易得：
，其中?表示按元素的Hadamard积，（?r表示取共轭操作；目标函数 写为：
[0132] 步骤2:分析MM0雷达波形正交性约束，构造满足正交性约束的目标函数：
[0133] 对于MM0发射波形矩阵：
而言，波形之间正交性体现在互相关 函数上；第nth时刻的发射波形相关矩阵Rn(X)的非主对角线部分尽可能逼近零，主对角元 素为各阵元发射波形的自相关表示：
[0137] 其中diag(a"(；e))表示nth时刻由理想波形矩阵x的自相关形成的对角矩阵；U(n) 表示偏移矩阵，I表示单位矩阵；
[0138]为解决不同阵元波形间互扰构造目标函数J2 (X)，同样考虑Parseval等价性，得到 如下：
[0140] 其中| | ? | |F表示矩阵的Frobenius范数；
[0141] 步骤3:根据场景先验信息预估频域干扰模糊频带区间，进而构造相应目标函数：
[0142] 假定可能存在电磁干扰的频段区间集合为
，那么发射波形频谱应 在此区间内设为阻带，由傅里叶变换的性质可知，频谱为：
[0144]其中：NS表示相应频段数目，巧表示单位离散傅里叶变换矩阵，且有
为使得频谱离散采样足够密以表征区间信息，设定及= 2W;由上式可知，如果从离散傅里 叶变换矩阵的行空间中抽取特定行来构成矩阵P，且所抽取行序号与D在频谱W中频 段相对应，用向量必表征&中稀疏凹口序列且满足
fn为频 谱W中频率的采样，从而使得：
[0146]只要# f幅度足够低就满足频谱稀疏特性，即目标函数J3(X)最小化，从而达到任 务要求：
[0148] 步骤4:构造恒模相位编码波形约束条件，即最大化发射功率要求恒模约束，也 即：
[0149] |Xn'J= 1，n= 1，…N.m= 1，??，M
[0150] 步骤5:构造松弛交替投影算法框架：
[0151] 假定集合A代表约束集合（即满足恒模约束），集合B代表目标集合（即满足子任 务要求），构造变量X，y到集合A,B的投影ProjA (x)和projB (y)，如下所示：
[0154] 其中d(i,j)= | | i-j| |表示变量i,j的距离函数。
[0155] 那么交替投影机制可简述为：迭代映射T:A-B-A。
[0156] 然后，引入松弛交替投影框架：
[0157] T5(.)；Xx=ProjA(x+A? 8 (x) ? (ProjA(ProjB(x))-x))
[0158]xk+1=Ts(.)，Xxk
[0159] 上式中，X表示松弛算子；S(?)表示加速因子，有如下不同形式：
[0162] 此算法框架命名为RISAAP。
[0163] 步骤6、上述步骤1、2、3分别完成了三个子任务目标函数的构建，以适应雷达多任 务场景的时变性，某时刻场景可能涉及上述子任务目标中的一个或多个，根据步骤5所提 松弛交替投影算法框架求解波形设计，分别给出如下三种波形优化输出方式：
[0164] (1)表示满足抑制旁瓣遮蔽效应的算法步骤；
[0165] (2)表示获得波形正交性的算法步骤；
[0166] (3)表示联合抑制特定频带干扰及获得良好相关特性的算法步骤；
[0167](1)对于抑制旁瓣遮蔽效应的算法，过程如下；
[0168] 以满足步骤1目标函数条件的集合为目标集合B;以满足步骤4条件集合为约束 集合A;(?)(1/2)表示按元素的平方根操作，返回值维度与其参数一致；
[0169] 那么获得较好时域自相关旁瓣抑制性能的恒模波形设计方法可归纳为：
[0170] a)初始化波形序列：x= [x(l). . .x(n)x(N)]TGA，其中x(?) =ew", itnG[0, 2it]为随机相位；
[0171]b)迭代优化过程：
[0172] 由公式
[0174]获得fGB，
[0175] 再由
[0177]获得xEA，
[0178] 经公式
[0179] Ts(.)，入x=ProjA(x+人? 8 (x) ? (ProjA(ProjB(x))_x))
[0181] 获得；并更新x= x;
[0182]c)相邻两次迭代产生差值小于给定误差e或算法迭代次数大于给定总数K，则 算法停止并输出，否则继续b)，c);
[0183] (2)对于获得波形正交性的算法设计，过程如下：
[0184] 以满足步骤2目标函数条件的集合，即为目标集合B;以满足步骤4条件的集合为 约束集合A;那么获得较好时域相关旁瓣抑制性能的恒模波形设计方法可归纳为：
[0185]a)初始化波形矩阵：X= [Xl...xM]GA，其中i
，步 n，me[0, 2n]为 随机相位；
[0186]b)迭代优化过程：
[0187] 由公式
[0189]获得FeB,
[0190] 再由
[0192]获得XGA，
[0193] 经公式
[0194] T5(.)；Ax=ProjA(x+A?8 (x)?(ProjA(ProjB(x))-x))
[0195] ^ -T&'pJ
[0196] 获得X并更新
[0197]c)若相邻两次迭代产生1差值小于给定误差e或算法迭代次数大于给定总数 K，则算法停止并输出，否则继续b)，c);
[0198] (3)对于联合抑制特定频带干扰及获得良好相关特性的算法设计，过程如下：
[0199] 以满足步骤1、2目标函数条件的集合，即为目标集合B1;以满足步骤3的集合，即 为目标集合B2;以满足步骤4条件的集合为约束集合A;那么获得较好时域相关旁瓣抑制性 能及频域特性的恒模波形设计方法可归纳为：
[0200] a)初始化波形矩阵：X = [Xl. ? ? xM] G A，其中
，步n，me[0, 2 n ]为 随机相位；
[0201] b)迭代优化过程：
[0202] 由公式
[0204] 获得巧e马，
[0205] 由公式
[0207]获得 e/i:，由
[0209]获得 FeB,
[0210] 再由
[0212] 获得XGA，
[0213] 经公式
[0214] Ts(.)，入x=ProjA(x+人? 8 (x) ? (ProjA(ProjB(x))_x))
[0216] 获得X并更新
[0217]c)若相邻两次迭代产生X差值小于给定误差e或算法迭代次数大于给定总数 K，则算法停止并输出，否则继续b)，c)。
[0218] 所述步骤1中时域干扰区间确定可根据先验信息，或通过事先发射传统波形以确 定大致区间。所述步骤2中多输入多输出雷达可包含多个阵元。所述步骤3中频域干扰区 间的确定可根据先验信息、或参考通用电子手册确定大致区间。所述步骤5中松弛算子入 的取值范围是：〇〈入〈2;Parteo权值y的取值范围是：0〈y〈1。
[0219] 如图1~6所示，结合步骤6中所列（1)、⑵、（3)方式，给出了基于本发明的技术 方案，从而可设计出任意阵元数及任意码长的具有良好时频抗干扰性能的离散恒模相位编 码信号的实施例：
[0220] 1.单阵元150个子码的抗单、多旁瓣干扰抑制波形设计：
[0221] 利用本技术方案中步骤6中（1)所述，假定强散射体可能位于单旁瓣干扰区 间％= [1:30]及多旁瓣干扰区间Q1= [20:30] U [40:60];初始化发射波形序列为：
，步ne [0,231]为随机相位且由随机数 产生；码长N = 150,算法迭代数K = 10000,迭代误差阈值e = 1(T4,松弛因子A = 〇. 7。 基于步骤6中（1)方式所列公式求解，直到相邻两次迭代产生分差值小于e或算法迭代 数大于K，则算法停止并输出。方案所提RISAAP算法能够在单旁瓣区间及多旁瓣区间获得 313dB和230dB的优化效果且耗时仅为0. 08s和0. 16s，其良好的抗干扰性能如表1、图1及 图2所示。
[0222] 2. 3阵元64码长阵元波形间正交性优化设计：
[0223] 利用本技术方案步骤6中（2)所述，初始化发射波形矩阵为：
其中
，[0,2 3!]为随机相位且由随机数产生；取码长N = 64,阵元数 M = 3 ;算法迭代数K = 10000 ;迭代误差阈值e = KT4;松弛因子A = 〇. 7。基于步骤6 中（2)方式所列公式求解，直到相邻两次迭代产生i差值小于e或算法迭代数大于K，则 算法停止并输出。作图仅以第一个阵元波形为例，RISAAP算法优化自相关幅度和互相关幅 度（第一、三阵元波形）仿真结果见图3、4,各算法随码长增加耗时比较见表二。
[0224 ] 3.单阵元100码长的联合抗旁瓣干扰及抑制特定频带电磁干扰的波形设计：
[0225] 利用本技术方案中步骤6中（3)所述，假设电磁干扰可能位于单频段区间 [200, 300] (KHz)和多频段区间[200, 300] U [500, 600] (KHz)，初始化发射波形矩阵为：
，步[0,231]为随机相位且由随机数产生；取 码长N = 100,阵元数M = 1 ;算法迭代数K = 10000 ;迭代误差阈值e = 1(T4;松弛因子入 =0. 7, Pareto权值y = 0. 6。基于步骤6中（3)方式所列公式求解，直到相邻两次迭代 产生差值小于e或算法迭代数大于K，则算法停止并输出。方案所提RISAAP算法能够 在单频段干扰区间及多频段干扰区间获得约14dB和13dB的优化效果且耗时仅为0. 06s和 0. 07s (如表3所示），其良好的抗干扰性能图5、6所示。
[0226] 表1不同算法抑制特定区间旁瓣比较
[0228] 表2不同序列长度算法耗时比较
[0229]
[0230] 表3不同频段干扰区间算法性能比较
[0232] 以上所述，仅为本发明较佳的【具体实施方式】，但本发明的保护范围并不局限于此， 任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，根据本发明的技术方案及其 发明构思加以等同替换或改变，都应涵盖在本发明的保护范围之内。
【主权项】
1. 一种多输入多输出雷达波形设计方法，其特征是包括以下步骤： 步骤1 :根据雷达场景中强散射体与待测目标的相对位置，预估自相关旁瓣抑制模糊 区间，进而构造相应的目标函数： 设M个发射阵元的MIMO雷达系统，其第m个阵元窄带恒模相位编码波形序列可表示为 Xm =[?)... Xm(W) xm(iV)]T e Cam，Xw(W) = ejv" , N 表征信号码长，Φ ne [〇, 2 π ]为第 nth时刻编码相位；因接收端滤波器可视为发射波形的相关函数，假定感兴趣的待测目标在 第Pth距离单元；那么在第Ith单元的强散射体势必要对P th单元目标造成影响；对于第P th距离单元，为最小化临近距离单元中强散射体产生的距离旁瓣对它的干扰，应使发射波形 Xm的自相关函数a (Xm)满足如下条件：其中，a (xm) = [ a Q(Xm)…a N_i (Xm) 〇 a _N+1 (xm)…α (xm) ]T表示自相关函数序列；· 表示理想波形；Qi={± |z-p| |z e Z1UiN-I,-Ν+1}表示旁瓣干扰区间，Z1表示距离向分布的 强散射体位置集合； 对于强散射体遮蔽干扰问题，发射波形^应使得目标函数J(Xm)最小化，即： min·/(λ:," H|a(xm)-Ct(^)II2 Xm 其中N · I I表示欧式范数，视为距离度量函数；由Parseval等价性可知：其中，表示离散傅里叶变换矩阵，# = 2#; FxrOt(Xm)eC2wxl表示自相关 函数的傅里叶变换，即功率谱；如果令/"和乙分别表征理想波形和设计波形的频谱，即 九= i^cT·^，C为扩展矩阵，其通过扩充零元素使得矩阵运算时维度一致；那么易得： Αα(毛,）= ./；"〇_/；二其中Θ表示按元素的Hadamard积，（·）*表示取共轭操作；目标函数 写为：步骤2 :分析MMO雷达波形正交性约束，构造满足正交性约束的目标函数： 对于MMO发射波形矩阵X=K..xM]e Cyxw而言，波形之间正交性体现在互相关函数 上；第nth时刻的发射波形相关矩阵Rn(X)的非主对角线部分尽可能逼近零，主对角元素为 各阵元发射波形的自相关表示：其中diag(a"(；e))表示nth时刻由理想波形矩阵X的自相关形成的对角矩阵；U(n)表示 偏移矩阵，I表示单位矩阵； 为解决不同阵元波形间互扰构造目标函数J2(X),同样考虑Parseval等价性，得到如 下：其中I I · I If表示矩阵的Frobenius范数； 步骤3 :根据场景先验信息预估频域干扰模糊频带区间，进而构造相应目标函数： 假定可能存在电磁干扰的频段区间集合为，那么发射波形频谱应在此 区间内设为阻带，由傅里叶变换的性质可知，频谱为： 其中：NS表示相应频段数目，G表示单位离散傅里叶变换矩阵，且有为使 得频谱离散采样足够密以表征区间信息，设定A = 27V;由上式可知，如果从离散傅里叶变 换矩阵的行空间中抽取特定行来构成矩阵P，且所抽取行序号与Ω在频谱W中频段& 相对应，用向量命表征Ω中稀疏凹口序列且满足也=l,/"e 亡，fn为频谱W 中频率的采样，从而使得： W = PC 'X 只要#幅度足够低就满足频谱稀疏特性，即目标函数J3(X)最小化，从而达到任务要 求： JAX) = m\nPC'X J X 步骤4 :构造恒模相位编码波形约束条件，即最大化发射功率要求恒模约束，也即： |Xn，J=l，n=l, ...N.m=l，· ·，M 步骤5 :构造松弛交替投影算法框架： 假定集合A代表约束集合（即满足恒模约束)，集合B代表目标集合（即满足子任务要 求)，构造变量X，y到集合A, B的投影Pro jA (X)和pro jB (y)，如下所示：其中d(i，j) = | |i-j| I表示变量i，j的距离函数。 那么交替投影机制可简述为：迭代映射T:A - B - A。 然后，引入松弛交替投影框架： Ts(.)，x=Pr〇jA(x+X · δ (X) · (ProjA(ProjB(x))-x)) Xk+l T δ ( . )，λ Xk 上式中，λ表示松弛算子；δ ( ·)表示加速因子，有如下不同形式：此算法框架命名为RISAAP。 步骤6、上述步骤1、2、3分别完成了三个子任务目标函数的构建，以适应雷达多任务场 景的时变性，某时刻场景可能涉及上述子任务目标中的一个或多个，根据步骤5所提松弛 交替投影算法框架求解波形设计，分别给出如下三种波形优化输出方式： (1) 表示满足抑制旁瓣遮蔽效应的算法步骤； (2) 表示获得波形正交性的算法步骤； (3) 表示联合抑制特定频带干扰及获得良好相关特性的算法步骤； (1)对于抑制旁瓣遮蔽效应的算法，过程如下； 以满足步骤1目标函数条件的集合为目标集合B ;以满足步骤4条件集合为约束集合 A ; ( · )(1/2)表示按元素的平方根操作，返回值维度与其参数一致； 那么获得较好时域自相关旁瓣抑制性能的恒模波形设计方法可归纳为： a) 初始化波形序列：x=[x(l)··· x(n) x(N)]TeA，其中 X(?) = ejv", Φηε[〇,2π] 为随机相位； b) 迭代优化过程：由公式 获得f G B， 再由获得X e A， 经公式获得4并更新.V = C)相邻两次迭代产生i差值小于给定误差ε或算法迭代次数大于给定总数K，则算法 停止并输出，否则继续b)，c); (2) 对于获得波形正交性的算法设计，过程如下： 以满足步骤2目标函数条件的集合，即为目标集合B ;以满足步骤4条件的集合为约束 集合A ;那么获得较好时域相关旁瓣抑制性能的恒模波形设计方法可归纳为： a) 初始化波形矩阵A=Iix1... xM] e A，其中XmO?) = , Φ[〇, 2 π ]为随机相 位； b) 迭代优化过程：获得X e A，经公式Ts(.)，x=Pr〇jA(x+X · δ (X) · (ProjA(ProjB(x))-x))λ - 由公式 获得FeB, 再由 获得X并更新X=义； C)若相邻两次迭代产生i差值小于给定误差ε或算法迭代次数大于给定总数K，则算 法停止并输出，否则继续b)，c); (3) 对于联合抑制特定频带干扰及获得良好相关特性的算法设计，过程如下： 以满足步骤1、2目标函数条件的集合，即为目标集合B1;以满足步骤3的集合，即为目 标集合B2;以满足步骤4条件的集合为约束集合A ;那么获得较好时域相关旁瓣抑制性能及 频域特性的恒模波形设计方法可归纳为： a)初始化波形矩阵：X=[Xl... XM] e A，其中 Xw(W) = ejv"' Φn，me [〇, 2 JT ]为随机相 位； b) 迭代优化过程：由公式 获得F1 ￡尽， 由公式 获得由 获得F e/i , 再由 获得X e A， 经公式获得1并更新 c) 若相邻两次迭代产生义差值小于给定误差ε或算法迭代次数大于给定总数K，则算 法停止并输出，否则继续b)，c)。2. 根据权利要求1所述的一种多输入多输出雷达波形设计方法，其特征是，所述步骤1 中时域干扰区间确定可根据先验信息，或通过事先发射传统波形以确定大致区间。3. 根据权利要求2所述的一种多输入多输出雷达波形设计方法，其特征是，所述步骤2 中多输入多输出雷达可包含多个阵元。4. 根据权利要求3所述的一种多输入多输出雷达波形设计方法，其特征是，所述步骤3 中频域干扰区间的确定可根据先验信息、或参考通用电子手册确定大致区间。5. 根据权利要求1、2、3或4所述的一种多输入多输出雷达波形设计方法，其特征是，所 述步骤5中松弛算子λ的取值范围是：〇〈λ〈2 ;Parteo权值γ的取值范围是：0〈γ〈1。
【专利摘要】一种多输入多输出雷达波形设计方法，属于雷达通信技术领域。本发明的目的是提供一种具有更低的相关旁瓣和频谱抑制深度，且效率高、耗时少、并具有较高的鲁棒性、具有良好时频抗干扰性能的设计方法。根据雷达场景中强散射体与待测目标的相对位置，预估自相关旁瓣抑制模糊区间，进而构造相应的目标函数；分析MIMO雷达波形正交性约束，构造满足正交性约束的目标函数；根据场景先验信息预估频域干扰模糊频带区间，进而构造相应目标函数；构造恒模相位编码波形约束条件；构造松弛交替投影算法框架；根据所提松弛交替投影算法框架求解波形设计，给出三种波形优化输出方式。采用松弛交替投影恒模波形编码设计，可使MIMO雷达具有更好的检测性能。
【IPC分类】G01S13/00
【公开号】CN104898113
【申请号】CN201510346063
【发明人】赵宜楠, 赵占锋, 冯翔, 周志权
【申请人】哈尔滨工业大学
【公开日】2015年9月9日
【申请日】2015年6月19日
最新回复(0)