基于gs算法的mimo雷达正交相位编码信号设计方法
基于gs算法的mimo雷达正交相位编码信号设计方法
【技术领域】
[0001] 本发明属于雷达信号设计技术领域，特别涉及一种基于GS算法的MM0雷达正交 相位编码信号设计方法，即基于盖师贝格 -撒克斯通（Gerchberg-Saxton，GS)算法的多输 入多输出（Multiple Input Multiple Output, MIM0)雷达正交相位编码信号设计方法，适 合用于多输入多输出（MM0)雷达正交相位编码信号的设计。
【背景技术】
[0002] 多输入多输出（Multiple Input Multiple Output，MIM0)雷达是一种新兴的有源 探测技术，是近年来雷达领域提出的一种全新的雷达体制，它能使雷达系统通过独特的时 间-能量管理技术实现多个独立宽波束同时照射，同时采用多信号发射、多信号接收，并且 多信号之间可以是时域、空域或极化域分离的。该技术具有处理维数极高、收发孔径利用更 充分、角分辨率更高的优点，目前已成为雷达技术领域的一个研宄热点。
[0003] 根据雷达天线的间距大小，多输入多输出（MM0)雷达可以分为分布式和集中式 两类。对于分布式多输入多输出（MIM0)雷达来说，由于各个雷达天线对目标有不同的观测 视角以及目标回波的独立性，使其具有空间分集的能力。在统计意义下，这类多输入多输出 (MM0)雷达可以通过大间距布阵形式获得空间分集增益，提高对起伏目标的检测性能；集 中式多输入多输出（MM0)雷达的天线间距较小，相对于相控阵雷达，其自由度明显增加， 可以利用波形分集进行处理，具有更高的角度分辨率，提高了雷达的抗截获性能和杂波背 景中探测低速、弱目标的能力，同时也可根据先验信息，灵活地设计发射方向图，提高雷达 资源利用率。
[0004] 由于多输入多输出（MM0)雷达发射信号是正交的，所以采用相位编码信号 设计多输入多输出（MM0)雷达信号波形，即雷达发射信号之间满足正交性；关于雷达 正交相位编码信号的设计方法，已有的方法有文献[Deng H.Polyphase code design for orthogonal netted radar systems[J]. IEEE Transactions on Signal Processi ng，2004,52(ll):3126-3135.]中应用的混合模拟退火算法，文献[Liu B，He Z S，Zeng J K, et al. . Polyphase orthogonal code design for MIM0 radar systems[C]. Proceedings of the International Conference on Radar, Shanghai, 2006:1-4.]中应用的混合遗传算 法和文献[Hu Liang-bing, Liu Hong-wei, and Wu Shun-jun. Orthogonal waveform design for MIM0 radar via constrained nonlinear programming[J]. Systems Engineering and Electronics, 2011，33(1) :64-68.]中应用的序列二次规划方法。这些文献采用算法设计均 能够得到性能较好的正交信号集，但是由于这些算法本身的限制，使得所得码集的非周期 自相关峰值旁辦（Autocorrelation Sidelobe Peak, ASP)、互相关峰值（Crosscorrelation Peak, CP)和积分旁瓣电平（Integrated Sidelobe Level, ISL)仍然较高。此外，由于这些 算法的复杂度较大，使得所得相位编码信号集大小受到一定限制，并且计算的存储量较大， 求解信号集所需时间较长，不利于工程中的实际应用。
【发明内容】
[0005] 本发明的目的在于克服上述已有方法的不足，提出一种基于GS算法的MM0雷达 正交相位编码信号设计方法，即基于盖师贝格-撒克斯通（Gerchberg-Saxton，GS)算法的 多输入多输出（Multiple Input Multiple 0utput，MIM0)雷达正交相位编码信号设计方 法，该方法对集中式多输入多输出（MIM0)雷达和分布式多输入多输出（MIM0)雷达均有效。
[0006] 本发明的实现思路是：先对多输入多输出（Multiple Input Multiple Output，MM0)雷达天线的发射信号构造代价函数，再把该代价函数转化为适合求解的极小 化函数，最后利用改进的盖师贝格-撒克斯通（GS)算法求解该极小化函数，得到要求解的 发射信号矩阵；具体来说，首先对多输入多输出（MIM0)雷达的发射信号集构造正交相位编 码信号的代价函数，其次利用傅立叶变化将该代价函数转化为适合求解的极小化函数，最 后利用改进的盖师贝格-撒克斯通（GS)算法求解其极小化函数，得到多输入多输出（MIM0) 雷达天线发射相位编码信号矩阵。
[0007] 为实现上述技术目的，本发明采用如下技术方案予以实现。
[0008] -种基于GS算法的MM0雷达正交相位编码信号设计方法，其特征在于，包括以下 步骤：
[0009] 步骤1，构造MM0雷达M个天线发射正交相位编码信号的代价函数E，
[0011] 其中，&表示MM0雷达M个天线发射正交相位编码信号的二阶中心矩，N表示MM0 雷达每个天线发射正交相位编码信号的码长，I为单位矩阵，R n表示MIM0雷达M个天线发 射正交相位编码信号的协方差矩阵，I I I I表示矩阵范数；
[0012] 步骤2,将MM0雷达M个天线发射正交相位编码信号的代价函数E先转化为原始 极
[0014]小化函数s. t. | sm(n)| =1，m G {1，2,…，M}，n G {1，2,…，N}，并将该原始极小 化函数转化为I|ak||2=1，k G {1，2,…，2N}
[0015]适合求解的极小化函数 |Ids' -Ql I2，
Q = [a1； a2, . . . , a2N]T,
[0016] 其中，D表示2NX2N维离散傅里叶变换矩阵，S'表示MMO雷达M个天线发射正 交相位编码信号矩阵的扩展矩阵，I I I I表示矩阵范数，Q表示辅助变量矩阵，ak表示辅助变 量，且ak是满足| |z(? k) | |2= 1的所有解的集合，z(? k)表示MM0雷达M个天线发射正 交相位编码信号的第n个码元7"的傅里叶变换，s m(n)表示MM0雷达第m个天线发射正交 相位编码信号的第n个码元的具体码元值；
[0017] 步骤3,利用GS算法求解步骤2中适合求解的极小化函数| |DS' -Q| |2,得到MIM0 雷达M个天线发射正交相位编码信号矩阵的扩展矩阵S '，进而得到MIM0雷达M个天线发 射正交相位编码信号矩阵S。
[0018] 优选的，步骤3中，所述GS算法求解步骤2中适合求解的极小化函数 |DS' -q||2,其具体子步骤如下：
[0019] 3.1求解频谱面输出函数Q，（i):
[0020] 利用[0, 2it]之间均匀分布的相位初始化NXM维的MM0雷达M个天线发射正交 相位编码信号矩阵S(i)，并将该正交相位编码信号矩阵S(i)扩展为
得到频 谱面输出函数Q' (i) = DS' (i)，其中，D表示2N X 2N维离散傅里叶变换矩阵，i表示迭代次 数；
[0021] 3. 2确定频谱面的幅度约束：
[0022] 通过限定离散傅里叶变换矩阵D的每一个行向量对应的2N个通道系数向量的模 值平方和为1，使得满足求解需要的约束条件成立，即频谱面施加的幅度约束为：
[0024] 其中，〇n)表示MM0雷达第m个天线的第k个通道的幅度归一化系数向量， m G {1，2,…，M}，M表示MM0雷达天线数，N表示MM0雷达每个天线发射的正交相位编码 信号码长，N也表示MIM0雷达每个天线的通道系数个数，i表示迭代次数。使用改进的约束 条件更适合求解本发明中的极小化函数；
[0025] 3. 3求解物面输入函数S"⑴：
[0026] 对频谱面输出函数Q' (i)的每一列进行逆离散傅里叶变换，得到物面输入函数 S"w，即S"《=DhQ' w，其中，D表示2NX2N维离散傅里叶变换（DFT)矩阵，{}H表示共 轭转置，i表示迭代次数，Q' (i)表示频谱面输出函数；
[0027] 3. 4求解物面的幅度约束：
[0028] 物面的幅度约束即为所求的MIM0雷达M个天线发射正交相位编码信号矩阵S，代 价函数E转化式中的约束条件|s m(n) | = l，sm(n)表示MIM0雷达第m个天线发射正交相位 编码信号的第n个码元的具体码元值，m G {1，2,…，M}，n G {1，2,…，N};
[0029] 约束条件|sjn)| =1表明MM0雷达每个天线发射正交相位编码信号为恒模信 号，将MM0雷达M个天线发射正交相位编码信号矩阵S (i)中的每个元素幅度限定为1，此处 定义A为NXM维全1矩阵，且截取矩阵
，则得到更新的MM0雷达M个天线发 射正交相位编码信号矩阵的扩展矩阵S'(i+1):
[0030] S，（1+1) =J?exp(jarg(S"（1)))
[0031] 其中，?表示Hadamard积，D表示2NX 2N维离散傅里叶变换（DFT)矩阵， arg(S" (i))表示取S" (i)的角度值，S"（i)表示物面输入函数，i表示迭代次数；
[0032] 3. 5重复子步骤3. 1~3. 4,对MM0雷达M个天线发射正交相位编码信号矩阵S(i) 进行迭代计算，直到两次迭代之间的误差小于预先设定值e，g卩| |s⑴-s(i+1)| I < e，停止 迭代，该第i次迭代所得的MM0雷达M个天线发射正交相位编码信号矩阵S(i)就是所求的 MM0雷达M个天线发射正交相位编码信号矩阵S。
[0033] 本发明具有以下有益效果：
[0034] (1)本发明设计的多输入多输出（MM0)雷达正交相位编码信号所得码集的平均 自相关峰值旁瓣、平均互相关峰值和积分旁瓣电平较低，而且该所得码集的积分旁瓣电平 接近于积分旁瓣电平的下界。
[0035] (2)本发明的主要步骤是基于快速傅里叶变换（FFT)和逆快速傅里叶变换（IFFT) 完成，效率较高，且优化过程耗时较短，比现有算法（比如混合遗传算法、混合模拟退火算 法）的效率高两个数量级（1〇 2)，而且计算存储量较小，可以设计较大的雷达信号波形集，如 雷达天线个数M最大可以设计到200,每个雷达天线发射的相位编码信号码长N最大可以设 计到1〇 4,有利于工程实现。
[0036] (3)本发明可以用于设计单个序列，且所设计单个序列的峰值旁瓣电平和积分旁 瓣电平比m序列、p4序列更低。
【附图说明】
[0037] 下面结合附图和【具体实施方式】对本发明作进一步详细说明。
[0038] 图1是本发明的基于盖师贝格-撒克斯通（GS)算法的多输入多输出（MIM0)雷达 正交相位编码信号设计方法流程图；
[0039] 图2是本发明中改进的盖师贝格-撒克斯通（GS)算法流程图；
[0040] 图3是使用本发明所得发射正交相位编码信号波形集的自相关函数图；
[0041] 图4是使用本发明所得发射正交相位编码信号波形集的互相关函数图；
[0042] 图5是使用本发明所得发射正交相位编码信号波形集的平均自相关峰值旁瓣 (ASP)和平均互相关峰值（CP)随多输入多输出（MMO)雷达天线数目M的变化关系示意图；
[0043] 图6是使用本发明所得发射正交相位编码信号波形集的平均自相关峰值旁瓣 (ASP)和平均互相关峰值（CP)随每个天线发射相位编码信号的码长N的变化关系示意图；
[0044] 图7是使用本发明所得单个序列分别与m序列、p4序列和随机相位序列的峰值旁 瓣电平的对比示意图；
[0045] 图8是使用本发明所得单个序列分别与m序列、p4序列和随机相位序列的积分旁 瓣电平的对比示意图；
【具体实施方式】
[0046] 参照图1，为本发明的基于GS算法的MMO雷达正交相位编码信号设计方法流程 图，该基于基于GS算法的MMO雷达正交相位编码信号设计方法包括以下步骤：
[0047] 步骤1，构造MMO雷达M个天线发射正交相位编码信号的代价函数E，
[0049] 其中，&表示MMO雷达M个天线发射正交编码信号的二阶中心矩，N表示MMO雷 达每个天线发射正交相位编码信号的码长，I为单位矩阵，R n表示MIMO雷达M个天线发射 相位编码信号的协方差矩阵，I I I I表示矩阵范数；
[0050] 步骤1的具体构造过程为：
[0051] 假定一个MIMO雷达具有M个天线，M为自然数，每个天线发射码长为N的正交相 位编码信号，MMO雷达M个天线发射正交相位编码信号矩阵S可以表示为：
[0052] S = [s" s2,…，sm，…，sM]NXM= [y " y2,…，yn，…，yN]T <1>
[0053] 其中，sm= [Sm(l)，sj2)，…，sm(N)]T表示MMO雷达第m个天线发射正交相位编 码信号，y n= [s ?，s2(n), ~sm(n),…，sM(n)]T表示MIMO雷达M个天线发射正交相位编 码信号的第n个码元，
表示MM0雷达第m个雷达天线发射正交相位编码信号 的第n个码元的具体码元值，(K(n)表示MIM0雷达的每个天线发射正交相位编码信号的相 位，且0彡伞m(n)彡2 31。
[0054] 进而，MM0雷达M个天线发射正交相位编码信号的协方差矩阵Rn表示为：
[0056]其中，『表示共轭转置，n G {1，2, ...，N}，m G {1，2, ...，M}，k G {1，2,…，2N}，yk =[Sl (k)，s2 (k)，…，sm (k)，…，sM(k) ]T表示M个雷达天线发射正交相位编码信号的第k个 码元。
[0057] 然后，得到MM0雷达M个天线发射正交相位编码信号的代价函数E，其表示式如下 所示：
[0059] 其中，&表示MM0雷达M个天线发射正交相位编码信号的二阶中心矩，N表示MM0 雷达每个天线发射正交相位编码信号的码长，I为单位矩阵，R n表示MIM0雷达M个天线发 射正交相位编码信号的协方差矩阵，I I I I表示矩阵范数，此处的矩阵范数为F-范数。如无 特别说明，本专利中矩阵范数为F-范数，向量范数为2-范数。
[0060] 步骤2,将MM0雷达M个天线发射正交相位编码信号的代价函数E先转化为原始 极
[0062]小化函数 s. t. | sm(n) | = 1，m G {1，2,…，M}，n G {1，2,…，N}，并将该原始极小 化函数转化为I |ak| |2= 1，k G {1，2,…，2N}
[0063]适合求解的极小化函数 | |DS' -Ql I2，
，Q = [ai, a2, . . .，a2N]T，
[0064] 其中，D表示2NX2N维离散傅里叶变换矩阵，S'表示MMO雷达M个天线发射正 交相位编码信号矩阵的扩展矩阵，I I I I表示矩阵范数，Q表示辅助变量矩阵，ak表示辅助变 量，且ak是满足| |z(? k) | |2= 1的所有解的集合，z(? k)表示MM0雷达M个天线发射正 交发射正交相位编码信号的第n个码元7"的傅里叶变换，s m(n)表示MM0雷达第m个天线 发射正交相位编码信号的第n个码元的具体码元值；
[0065] 步骤2的具体过程为：
[0066] 定义MM0雷达M个天线发射正交发射相位编码信号的第n个码元yn的傅里叶变 换为z (?)，且
，得到MM0雷达M个天线发射正交相位编码信号的第n个 码元的功率谱密度矩阵s(?):
[0068]由式〈3>和式〈4>得出如下两式等价：
[0070] 即MM0雷达M个天线发射正交相位编码信号的代价函数E = 0与功率谱密度矩 阵S(u) =1是等价的，其中I为单位矩阵。
[0071] 因此，MM0雷达M个天线发射正交相位编码信号的代价函数E可以转化为等价代 价函数e 1:
[0075] 忽略与变量无关的常数，式〈7>可以转化为原始的极小化函数：
[0077] s.t. |sm(n) | = l,m G {1, 2, - , M}, n G {1, 2, - ,N}. <8>
[0078] | |ak| |2= 1，k G {1，2,…，2N}
[0079] 其中，= eM'(n)表示MMO雷达第m个天线发射正交相位编码信号的第n个码 元的具体码元值，z〇k)表示MIMO雷达M个天线发射正交相位编码信号的第n个码元y n 的傅里叶变换，S(?k)表示功率谱密度矩阵，ak表示辅助变量，为满足| |z(? k) | |2= 1的 所有解的集合，下标k表示第几个辅助变量，| | | |表示F-范数，M表示MIMO雷达发射天线 数目，N表示每个天线发射正交相位编码信号的码长。
[0080] 为表示方便，定义2NX2N维离散傅里叶变换（DFT)矩阵D的第（k，n)个元素为 ，其中，令
，并将NXM维的MMO雷达发射正交相位编码信号 矩阵S扩展为'
，辅助变量矩阵Q = [ai，a2，…，ak，…，a2N]T，则式〈8>中的原始 的极小化函数可转化为适合求解的极小化函数，BP :
[0081] |DSr -Q| |2 <9>
[0082] 步骤3,利用GS算法 求解步骤2中适合求解的极小化函数I IDS' -Ql I2,得到MIMO 雷达M个天线发射正交相位编码信号矩阵的扩展矩阵S '，进而得到MIMO雷达M个天线发 射正交相位编码信号矩阵S。
[0083]本发明利用gs算法求解步骤2中适合求解的极小化函数|Ids'-Q||2。由于原始 盖师贝格-撒克斯通GS算法的极小化函数是物面和频谱面的一维变换，而本发明中的极小 化函数是多维的，故先将GS算法扩展至多维，并将其物面和频谱面的约束条件进行改进， 以满足本发明优化问题的求解条件。
[0084] 参照图2,式〈9>的求解过程可以分为以下循环迭代的四个子步骤：
[0085] 3.1求解频谱面输出函数Q'（i):
[0086] 利用[0, 2 JT ]之间均匀分布的相位初始化NXM维的MM0雷达M个天线发射正交 相位编码信号矩阵s(i)，并将该正交相位编码信号矩阵s(i)扩展为
，得到频 谱面输出函数Q' (i) = DS' (i)，其中，D表示2N X 2N维离散傅里叶变换矩阵，i表示迭代次 数；
[0087] 具体地，
[0088] 原始盖师贝格-撒克斯通（GS)算法利用交替进行傅里叶变换和傅里叶逆变换，对 物面和频谱面的一维函数进行求解；而本发明中的所求的MM0雷达M个天线发射正交相位 编码信号矩阵S为NXM维，如果要使用利用盖师贝格-撒克斯通（GS)算法求解该极小化 函数，必须将其物面和频谱面的约束条件进行改进，以满足本发明的求解条件。
[0089] 利用[0, 2 JT ]之间均匀分布的相位初始化NXM维的MM0雷达M个天线发射正交 相位编码信号矩阵S(i)(以后每次迭代保留步骤3. 1MM0雷达M个天线发射正交相位编码 信号矩阵S(i))，然后将其扩展为
；区别于原始GS算法进行的一维傅里叶变 换，本发明对MM0雷达M个天线发射正交相位编码信号矩阵的扩展矩阵S' (i)的每一列进 行离散傅立叶变换〇FT)(相当于M维），得到频谱面的输出函数Q' (i)=DS' (i);其中，D 表示2NX 2N维离散傅里叶变换矩阵，i表示迭代次数。
[0090] 3. 2确定频谱面的幅度约束：
[0091] 通过限定离散傅里叶变换矩阵D的每一个行向量对应的2N个通道系数向量的模 值平方和为1，使得满足求解需要的约束条件成立，即频谱面施加的幅度约束为：
[0093] 其中，〇〇表示MM0雷达第m个天线的第k个通道的幅度归一化系数向量， m G {1，2,…，M}，M表示MM0雷达天线个数，N表示MM0雷达每个天线发射正交相位编码 信号码长，N也表示MIM0雷达每个天线的通道系数个数，i表示迭代次数。使用改进的约束 条件更适合求解本发明中的极小化函数；
[0094] 具体地，
[0095] 令频谱面的输出函数
.其中，频谱面的输出函数 Q' (1)的行向量为M列，kG {1，2，...，2叫，]/[表示|/[頂0雷达天线个数，；[表示迭代次 数。
[0096]对M列的行向量>进行幅度归一化得：
[0098] 其中，巧1表示行向量 中第k个通道的幅度归一化系数向量，k G {1，2,… ，2N}，i表示迭代次数。
[0099] 频谱面的输出函数Q' (i)的列向量的2N个元素与离散傅里叶变换（DFT)矩阵D 的2N个通道系数向量一一对应；2N X 2N维离散傅里叶变换矩阵D表示行数为2N、列数为2N 的矩阵，且该离散傅里叶变换矩阵D的每一个行向量对应2N个通道系数向量。
[0100] 此时，MM0雷达每个天线发射码长为2N的正交相位编码信号，同时，MM0雷达每 个天线也对应2N个通道系数向量，通过限定离散傅里叶变换（DFT)矩阵D的每一个行向量 对应的2N个通道系数向量的模值平方和为1，使得满足求解需要的约束条件成立，即频谱 面施加的幅度约束为：
[0102] 其中，〇0表示MM0雷达第m个天线的第k个通道的幅度归一化系数向量， m G {1，2,…，M}，M表示MM0雷达天线数，N表示MM0雷达每个天线发射正交相位编码信 号码长，N也表示MIM0雷达每个天线的通道系数个数，i表示迭代次数。使用改进的约束条 件更适合求解本发明中的极小化函数。
[0103] 3. 3求解物面输入函数S"⑴：
[0104] 对频谱面输出函数Q' (i)的每一列进行逆离散傅里叶变换，得到物面输入函数 S"w，即S"《=DHQ' w，其中，D表示2NX2N维离散傅里叶变换（DFT)矩阵，{}H表示共 轭转置，i表示迭代次数，Q' (i)表示频谱面输出函数；
[0105] 3. 4求解物面的幅度约束：
[0106] 物面的幅度约束即为所求的MM0雷达M个天线发射正交相位编码信号矩阵S，代 价函数E转化式中的约束条件|s m(n) | = l，sm(n)表示MIM0雷达第m个天线发射正交相位 编码信号的第n个码元的具体码元值，m G {1，2,…，M}，n G {1，2,…，N};
[0107]约束条件| sjn) | = 1表明MM0雷达每个天线发射正交相位编码信号为恒模信 号，将MM0雷达M个天线发射正交相位编码信号矩阵S(i)中的每个元素幅度限定为1，此处 定义A为NXM维全1矩阵，且截取矩阵.
，则得到更新的MM0雷达M个天线发 射正交相位编码信号矩阵的扩展矩阵S'(i+1):
[0108] S'（1+1) = J ? exp(jarg(S"（1)))
[0109] 其中，?表示Hadamard积，D表示2NX 2N维离散傅里叶变换（DFT)矩阵， arg(S" (i))表示取S" (i)的角度值，S"（i)表示物面输入函数，i表示迭代次数；
[0110] 具体地，
[0111]由于MM0雷达M个天线发射正交相位编码信号矩阵S(i)是NXM维矩阵，而MM0 雷达M个天线发射正交相位编码信号矩阵的扩展矩阵S' (i+1)是2NXM维矩阵，因此需要对 更新的MMO雷达M个天线发射正交相位编码信号矩阵的扩展矩阵S'(i+1)的维数进行调整， 得到更新的MMO雷达M个天线发射正交相位编码信号矩阵S (i+1)。
[0112] 并且，将MM0雷达M个天线发射正交相位编码信号矩阵S(i)中的每个元素幅度限 定为1，此处定义A为NXM维全1矩阵，且截取矩阵
[0113] 3. 5重复子步骤3. 1~3. 4,对MM0雷达M个天线发射正交相位编码信号矩阵S(i) 进行迭代计算，直到两次迭代之间的误差小于预先设定值e，g卩| |S⑴-S(i+1)| | < e，停止 迭代，该第i次迭代所得的MM0雷达M个天线发射正交相位编码信号矩阵S(i)就是所求的 MMO雷达M个天线发射正交相位编码信号矩阵S。
[0114] 本发明的效果通过以下仿真对比试验进一步说明：
[0115] (1)仿真条件：
[0116] 本发明仿真采用收发共置的集中式MM0雷达系统，假定目标雷达散射截面 (Radar Cross Section，RCS)在积累周期内保持恒定，并忽略MM0雷达M个天线发射编码 信号内的多普勒频移和距离走动问题，仿真中如无特别说明，采用MM0雷达天线数为3,每 个天线发射正交相位编码信号的码长为40,即M = 3, N = 40。
[0117] (2)仿真内容：
[0118] MM0雷达天线个数为M，每个天线发射正交相位编码信号的码长为N的正交相位 编码信号的积分旁瓣电平（Integrated Sidelobe Level, ISL)下界为：
[0119] ISL ^ N2M(M-1) = Bisl〇
[0120] 其中，Bia表示平均积分旁瓣电平（ISL)的下界值。
[0121] 表1给出了本发明所得MM0雷达M个天线发射正交相位编码信号的波形和随机 相位序列信号波形的平均积分旁瓣电平（ISL)与平均积分旁瓣电平（ISL)下界的比较。
[0122] 由表1可以看出，对于各种不同大小的信号波形集，本发明所得多输入多输出 (MIM0)雷达M个天线发射正交相位编码信号波形的平均积分旁瓣电平（ISL)均接近于平均 积分旁瓣电平（ISL)的下界，而随机相位序列信号波形的平均积分旁瓣电平（ISL)则远远 高于平均积分旁瓣电平（ISL)的下界。
[0123] 表1本发明所得波形和随机相位序列的ISL与ISL下界的比较
[0125] 图3为MM0雷达天线数为3,每个天线发射正交相位编码信号的码长为40时，采 用本发明方案而产生的波形集自相关函数示意图，横坐标为时间/采样序号，纵坐标为归 一化幅度；图4为MM0雷达天线数为3,每个天线发射正交相位编码信号的码长为40时， 采用本发明方案而产生的波形集互相关函数示意图，横坐标为时间/采样序号，纵坐标为 归一化幅度。
[01 26] 由图3和图4可以看出，该波形集的平均自相关旁瓣峰值为-16. 6dB，平均互相关 峰值为-14. 3dB ;可见本发明所得波形集的平均自相关旁瓣峰值和平均互相关峰值较低， 相关性能较好。
[0127] 在算法执行效率方面，由于本发明算法采用的离散傅里叶变换（DFT)和逆离散傅 里叶变换（IDFT)，是基于快速傅里叶变换（FFT)和逆快速傅里叶变换（IFFT)进行处理的， 故本发明算法的效率较高，优化耗时较短。而且，当波形集较大时本发明算法仍然适用，如 MM0雷达天线个数M最大可设计到200,每个雷达天线发射正交相位编码信号的码长N最 大可设计到1〇 4。
[0128] 此外，已有文献中提出的混合模拟退火算法、混合遗传算法和序列二次规划算法 的复杂度都较高，优化耗时较长，计算存储量较大，仅适用于波形集较小的情况。
[0129] 表2给出了本发明算法和上述几种算法的执行时间对比，其所得时间均为波形集 大小为M = 3, N = 40时执行10次算法所需的时间平均值。为了使所述几种算法执行时间 对比具有说服力，本文仅列出时间的数量级（计算机硬件条件为：P entium(R)双核处理器， 主频分别是2. 50GHz和2. 29GHz，内存1. 98GB，采用Matlab语言编写程序）。
[0130] 由表2可知，已有算法中执行时间最少的是混合遗传算法，而设计同样大小的波 形集，本发明所需时间是混合遗传算法所需时间的1/100 ;由此可见，本发明算法在执行效 率上较已有算法具有很强的优势，这在工程应用中极为有利。
[0131] 表2本发明算法与文献[Deng Liu Hu]算法的执行时间对比
[0133] 图5给出了每个雷达天线发射正交相位编码信号的码长N固定为128时，波形集 自相关峰值旁辦（Autocorrelation Sidelobe Peak, ASP)、互相关峰值（Crosscorrelation Peak，CP)随MIM0雷达天线个数M的变化关系示意图；图6为MIM0雷达天线个数M固定为 4时，波形集平均自相关峰值旁瓣（ASP)和平均互相关峰值（CP)随每个雷达天线发射正交 相位编码信号的码长N变化的关系示意图。从图5可以看出，当码长固定时，波形集自相关 峰值旁瓣（ASP)和互相关峰值（CP)随MM0雷达天线个数M的增加而增大，并且当MM0雷 达天线个数M增大到一定程度时，波形集自相关峰值旁瓣（ASP)和互相关峰值（CP)的增大 变得不明显；由图6可知，当MM0雷达天线个数M固定时，波形集自相关峰值旁瓣（ASP)和 互相关峰值（CP)与log 2N近似成负线性关系，M表示MM0雷达天线数，N表示每个天线发 射正交相位编码信号的码长。
[0134] 当MM0雷达天线个数M = 1时，使用本发明算法获得的正交相位编码信号矩阵为 单个序列，因为m序列、p4序列的相关性能较好，而且已经广泛应用于工程领域中，而随机 序列是未经算法优化涉及的序列，因此，使用m序列、p4序列与和随机序列与本发明算法所 得单个序列做对比可以说明本算法的优化性能。图7和图8分别给出了本发明算法所得单 个序列与m序列波形、p4序列波形和随机相位序列波形的峰值旁瓣电平和积分旁瓣电平的 对比。
[0135] 由图7和图8得，本发明算法所得单个序列分别比m序列、p4序列和随机相位序 列的自相关峰值旁瓣电平和积分旁瓣电平都要低，随机相位序列波形的自相关峰值旁瓣电 平和积分旁瓣电平最高，P4序列波形比m序列波形的要低。而且，随着MM0雷达天线个数 的增加，使用本发明算法所得发射正交相位编码信号波形、m序列波形、p4序列波形和随机 相位序列波形的自相关峰值旁瓣电平均逐渐减小，且积分旁瓣电平也均逐渐增大。
【主权项】
1. 一种基于GS算法的MMO雷达正交相位编码信号设计方法，其特征在于，包括以下步 骤： 步骤1，构造 MMO雷达M个天线发射正交相位编码信号的代价函数E，其中，Rtl表示MMO雷达M个天线发射正交相位编码信号的二阶中心矩，N表示MMO雷 达每个天线发射正交相位编码信号的码长，I表示单位矩阵，Rn表示MIMO雷达M个天线发 射正交相位编码信号的协方差矩阵，I I I I表示矩阵范数； 步骤2,将MMO雷达M个天线发射正交相位编码信号的代价函数E先转化为原始极小 化函数并将该原始极小化函数转化为适 合求解的极小化函数其中，D表示2NX2N维离散傅里叶变换矩阵，表示MMO雷达M个天线发射正交相 位编码信号矩阵的扩展矩阵，I I I I表示矩阵范数，Q表示辅助变量矩阵，ak表示辅助变量， 且ak是满足I |ζ(ω k) I I2= 1的所有解的集合，ζ(ω k)表示MMO雷达M个天线发射相位编 码信号的第η个码元yn的傅里叶变换，s m (η)表示MMO雷达第m个天线发射正交相位编码 信号的第η个码元的具体码元值； 步骤3,利用GS算法求解步骤2中适合求解的极小化函数I IDS' -Ql |2,得到MIMO雷 达M个天线发射正交相位编码信号矩阵的扩展矩阵S '，进而得到MIMO雷达M个天线发射 正交相位编码信号矩阵S。2. 根据权利要求1所述的基于GS算法的MMO雷达正交相位编码信号设计方法，其特 征在于，步骤3中，所述GS算法求解步骤2中适合求解的极小化函数I IDS' -Ql 12,其具体 子步骤如下： 3. 1求解频谱面输出函数V (i): 利用[〇, 2 π ]之间均匀分布的相位初始化NXM维的MMO雷达M个天线发射正交相位 编码信号矩阵S(i)，并将该正交相位编码信号矩阵S(i)扩展为'得到频谱面 输出函数V (i)=DS< ω，其中，D表示2NX2N维离散傅里叶变换矩阵，i表示迭代次数； 3. 2确定频谱面的幅度约束： 通过限定离散傅里叶变换矩阵D的每一个行向量对应的2N个通道系数向量的模值平 方和为1，使得满足求解需要的约束条件成立，即频谱面施加的幅度约束为：其中，〇n)表示MMO雷达第m个天线的第k个通道的幅度归一化系数向量， m e {1，2,…，M}，M表示MMO雷达天线数，N表示MMO雷达每个天线发射正交相位编码信 号码长，N也表示MIMO雷达每个天线的通道系数个数，i表示迭代次数。使用改进的约束条 件更适合求解本发明中的极小化函数； 3.3求解物面输入函数5"(1): 对频谱面输出函数V (i)的每一列进行逆离散傅里叶变换，得到物面输入函数S"（i)， 即S"ω =D1V ω，其中，D表示2NX 2N维离散傅里叶变换（DFT)矩阵，『表示共轭转置， i表示迭代次数，(^ (1)表示频谱面输出函数； 3. 4求解物面的幅度约束： 物面的幅度约束即为所求的MIMO雷达M个天线发射正交相位编码信号矩阵S，代价函 数E转化式中的约束条件I sm (η) I = 1，sm (η)表示MIMO雷达第m个天线发射正交相位编码 信号的第η个码元的具体码元值，m e {1，2,…，M}，n e {1，2,…，N}; 约束条件I Sm(n) I = 1表明MMO雷达每个天线发射正交相位编码信号为恒模信号，将 MMO雷达M个天线发射正交相位编码信号矩阵S(i)中的每个元素幅度限定为1，此处定义 A为NXM维全1矩阵，且截取矩阵項U得到更新的MMO雷达M个天线发射正交 相位编码信号矩阵的扩展矩阵S' (i+1): S'（1+1)= J Θ exp(jarg(S" ω)) 其中，Θ表示Hadamard积，D表示2NX2N维离散傅里叶变换（DFT)矩阵，arg(S"(i))表 示取S" ω的角度值，S" ω表示物面输入函数，i表示迭代次数； 3. 5重复子步骤3. 1~3. 4,对MMO雷达M个天线发射正交相位编码信号矩阵S(i)进行 迭代计算，直到两次迭代之间的误差小于预先设定值ε，即I |S⑴_S(i+1)| I < ε，停止迭代， 该第i次迭代所得的MMO雷达M个天线发射正交相位编码信号矩阵S(i)就是所求的MMO 雷达M个天线发射正交相位编码信号矩阵S。
【专利摘要】本发明公开了一种基于GS算法的MIMO雷达正交相位编码信号设计，其主要实现思路是：先对多输入多输出(MIMO)雷达天线的发射信号构造代价函数，再把该代价函数转化为适合求解的极小化函数，最后利用改进的盖师贝格-撒克斯通(GS)算法求解该极小化函数，得到要求解的发射信号矩阵；具体来说，首先对多输入多输出(MIMO)雷达的发射信号集构造正交相位编码信号的代价函数，其次利用傅立叶变化将该代价函数转化为适合求解的极小化函数，最后利用改进的盖师贝格-撒克斯通(GS)算法求解其极小化函数，得到多输入多输出(MIMO)雷达天线发射正交相位编码信号矩阵。
【IPC分类】G01S7/02
【公开号】CN104898093
【申请号】CN201510288013
【发明人】苏涛, 杜晓林
【申请人】西安电子科技大学
【公开日】2015年9月9日
【申请日】2015年5月29日
【技术领域】
[0001] 本发明属于雷达信号设计技术领域，特别涉及一种基于GS算法的MM0雷达正交 相位编码信号设计方法，即基于盖师贝格 -撒克斯通（Gerchberg-Saxton，GS)算法的多输 入多输出（Multiple Input Multiple Output, MIM0)雷达正交相位编码信号设计方法，适 合用于多输入多输出（MM0)雷达正交相位编码信号的设计。
【背景技术】
[0002] 多输入多输出（Multiple Input Multiple Output，MIM0)雷达是一种新兴的有源 探测技术，是近年来雷达领域提出的一种全新的雷达体制，它能使雷达系统通过独特的时 间-能量管理技术实现多个独立宽波束同时照射，同时采用多信号发射、多信号接收，并且 多信号之间可以是时域、空域或极化域分离的。该技术具有处理维数极高、收发孔径利用更 充分、角分辨率更高的优点，目前已成为雷达技术领域的一个研宄热点。
[0003] 根据雷达天线的间距大小，多输入多输出（MM0)雷达可以分为分布式和集中式 两类。对于分布式多输入多输出（MIM0)雷达来说，由于各个雷达天线对目标有不同的观测 视角以及目标回波的独立性，使其具有空间分集的能力。在统计意义下，这类多输入多输出 (MM0)雷达可以通过大间距布阵形式获得空间分集增益，提高对起伏目标的检测性能；集 中式多输入多输出（MM0)雷达的天线间距较小，相对于相控阵雷达，其自由度明显增加， 可以利用波形分集进行处理，具有更高的角度分辨率，提高了雷达的抗截获性能和杂波背 景中探测低速、弱目标的能力，同时也可根据先验信息，灵活地设计发射方向图，提高雷达 资源利用率。
[0004] 由于多输入多输出（MM0)雷达发射信号是正交的，所以采用相位编码信号 设计多输入多输出（MM0)雷达信号波形，即雷达发射信号之间满足正交性；关于雷达 正交相位编码信号的设计方法，已有的方法有文献[Deng H.Polyphase code design for orthogonal netted radar systems[J]. IEEE Transactions on Signal Processi ng，2004,52(ll):3126-3135.]中应用的混合模拟退火算法，文献[Liu B，He Z S，Zeng J K, et al. . Polyphase orthogonal code design for MIM0 radar systems[C]. Proceedings of the International Conference on Radar, Shanghai, 2006:1-4.]中应用的混合遗传算 法和文献[Hu Liang-bing, Liu Hong-wei, and Wu Shun-jun. Orthogonal waveform design for MIM0 radar via constrained nonlinear programming[J]. Systems Engineering and Electronics, 2011，33(1) :64-68.]中应用的序列二次规划方法。这些文献采用算法设计均 能够得到性能较好的正交信号集，但是由于这些算法本身的限制，使得所得码集的非周期 自相关峰值旁辦（Autocorrelation Sidelobe Peak, ASP)、互相关峰值（Crosscorrelation Peak, CP)和积分旁瓣电平（Integrated Sidelobe Level, ISL)仍然较高。此外，由于这些 算法的复杂度较大，使得所得相位编码信号集大小受到一定限制，并且计算的存储量较大， 求解信号集所需时间较长，不利于工程中的实际应用。
【发明内容】
[0005] 本发明的目的在于克服上述已有方法的不足，提出一种基于GS算法的MM0雷达 正交相位编码信号设计方法，即基于盖师贝格-撒克斯通（Gerchberg-Saxton，GS)算法的 多输入多输出（Multiple Input Multiple 0utput，MIM0)雷达正交相位编码信号设计方 法，该方法对集中式多输入多输出（MIM0)雷达和分布式多输入多输出（MIM0)雷达均有效。
[0006] 本发明的实现思路是：先对多输入多输出（Multiple Input Multiple Output，MM0)雷达天线的发射信号构造代价函数，再把该代价函数转化为适合求解的极小 化函数，最后利用改进的盖师贝格-撒克斯通（GS)算法求解该极小化函数，得到要求解的 发射信号矩阵；具体来说，首先对多输入多输出（MIM0)雷达的发射信号集构造正交相位编 码信号的代价函数，其次利用傅立叶变化将该代价函数转化为适合求解的极小化函数，最 后利用改进的盖师贝格-撒克斯通（GS)算法求解其极小化函数，得到多输入多输出（MIM0) 雷达天线发射相位编码信号矩阵。
[0007] 为实现上述技术目的，本发明采用如下技术方案予以实现。
[0008] -种基于GS算法的MM0雷达正交相位编码信号设计方法，其特征在于，包括以下 步骤：
[0009] 步骤1，构造MM0雷达M个天线发射正交相位编码信号的代价函数E，
[0011] 其中，&表示MM0雷达M个天线发射正交相位编码信号的二阶中心矩，N表示MM0 雷达每个天线发射正交相位编码信号的码长，I为单位矩阵，R n表示MIM0雷达M个天线发 射正交相位编码信号的协方差矩阵，I I I I表示矩阵范数；
[0012] 步骤2,将MM0雷达M个天线发射正交相位编码信号的代价函数E先转化为原始 极
[0014]小化函数s. t. | sm(n)| =1，m G {1，2,…，M}，n G {1，2,…，N}，并将该原始极小 化函数转化为I|ak||2=1，k G {1，2,…，2N}
[0015]适合求解的极小化函数 |Ids' -Ql I2，
Q = [a1； a2, . . . , a2N]T,
[0016] 其中，D表示2NX2N维离散傅里叶变换矩阵，S'表示MMO雷达M个天线发射正 交相位编码信号矩阵的扩展矩阵，I I I I表示矩阵范数，Q表示辅助变量矩阵，ak表示辅助变 量，且ak是满足| |z(? k) | |2= 1的所有解的集合，z(? k)表示MM0雷达M个天线发射正 交相位编码信号的第n个码元7"的傅里叶变换，s m(n)表示MM0雷达第m个天线发射正交 相位编码信号的第n个码元的具体码元值；
[0017] 步骤3,利用GS算法求解步骤2中适合求解的极小化函数| |DS' -Q| |2,得到MIM0 雷达M个天线发射正交相位编码信号矩阵的扩展矩阵S '，进而得到MIM0雷达M个天线发 射正交相位编码信号矩阵S。
[0018] 优选的，步骤3中，所述GS算法求解步骤2中适合求解的极小化函数 |DS' -q||2,其具体子步骤如下：
[0019] 3.1求解频谱面输出函数Q，（i):
[0020] 利用[0, 2it]之间均匀分布的相位初始化NXM维的MM0雷达M个天线发射正交 相位编码信号矩阵S(i)，并将该正交相位编码信号矩阵S(i)扩展为
得到频 谱面输出函数Q' (i) = DS' (i)，其中，D表示2N X 2N维离散傅里叶变换矩阵，i表示迭代次 数；
[0021] 3. 2确定频谱面的幅度约束：
[0022] 通过限定离散傅里叶变换矩阵D的每一个行向量对应的2N个通道系数向量的模 值平方和为1，使得满足求解需要的约束条件成立，即频谱面施加的幅度约束为：
[0024] 其中，〇n)表示MM0雷达第m个天线的第k个通道的幅度归一化系数向量， m G {1，2,…，M}，M表示MM0雷达天线数，N表示MM0雷达每个天线发射的正交相位编码 信号码长，N也表示MIM0雷达每个天线的通道系数个数，i表示迭代次数。使用改进的约束 条件更适合求解本发明中的极小化函数；
[0025] 3. 3求解物面输入函数S"⑴：
[0026] 对频谱面输出函数Q' (i)的每一列进行逆离散傅里叶变换，得到物面输入函数 S"w，即S"《=DhQ' w，其中，D表示2NX2N维离散傅里叶变换（DFT)矩阵，{}H表示共 轭转置，i表示迭代次数，Q' (i)表示频谱面输出函数；
[0027] 3. 4求解物面的幅度约束：
[0028] 物面的幅度约束即为所求的MIM0雷达M个天线发射正交相位编码信号矩阵S，代 价函数E转化式中的约束条件|s m(n) | = l，sm(n)表示MIM0雷达第m个天线发射正交相位 编码信号的第n个码元的具体码元值，m G {1，2,…，M}，n G {1，2,…，N};
[0029] 约束条件|sjn)| =1表明MM0雷达每个天线发射正交相位编码信号为恒模信 号，将MM0雷达M个天线发射正交相位编码信号矩阵S (i)中的每个元素幅度限定为1，此处 定义A为NXM维全1矩阵，且截取矩阵
，则得到更新的MM0雷达M个天线发 射正交相位编码信号矩阵的扩展矩阵S'(i+1):
[0030] S，（1+1) =J?exp(jarg(S"（1)))
[0031] 其中，?表示Hadamard积，D表示2NX 2N维离散傅里叶变换（DFT)矩阵， arg(S" (i))表示取S" (i)的角度值，S"（i)表示物面输入函数，i表示迭代次数；
[0032] 3. 5重复子步骤3. 1~3. 4,对MM0雷达M个天线发射正交相位编码信号矩阵S(i) 进行迭代计算，直到两次迭代之间的误差小于预先设定值e，g卩| |s⑴-s(i+1)| I < e，停止 迭代，该第i次迭代所得的MM0雷达M个天线发射正交相位编码信号矩阵S(i)就是所求的 MM0雷达M个天线发射正交相位编码信号矩阵S。
[0033] 本发明具有以下有益效果：
[0034] (1)本发明设计的多输入多输出（MM0)雷达正交相位编码信号所得码集的平均 自相关峰值旁瓣、平均互相关峰值和积分旁瓣电平较低，而且该所得码集的积分旁瓣电平 接近于积分旁瓣电平的下界。
[0035] (2)本发明的主要步骤是基于快速傅里叶变换（FFT)和逆快速傅里叶变换（IFFT) 完成，效率较高，且优化过程耗时较短，比现有算法（比如混合遗传算法、混合模拟退火算 法）的效率高两个数量级（1〇 2)，而且计算存储量较小，可以设计较大的雷达信号波形集，如 雷达天线个数M最大可以设计到200,每个雷达天线发射的相位编码信号码长N最大可以设 计到1〇 4,有利于工程实现。
[0036] (3)本发明可以用于设计单个序列，且所设计单个序列的峰值旁瓣电平和积分旁 瓣电平比m序列、p4序列更低。
【附图说明】
[0037] 下面结合附图和【具体实施方式】对本发明作进一步详细说明。
[0038] 图1是本发明的基于盖师贝格-撒克斯通（GS)算法的多输入多输出（MIM0)雷达 正交相位编码信号设计方法流程图；
[0039] 图2是本发明中改进的盖师贝格-撒克斯通（GS)算法流程图；
[0040] 图3是使用本发明所得发射正交相位编码信号波形集的自相关函数图；
[0041] 图4是使用本发明所得发射正交相位编码信号波形集的互相关函数图；
[0042] 图5是使用本发明所得发射正交相位编码信号波形集的平均自相关峰值旁瓣 (ASP)和平均互相关峰值（CP)随多输入多输出（MMO)雷达天线数目M的变化关系示意图；
[0043] 图6是使用本发明所得发射正交相位编码信号波形集的平均自相关峰值旁瓣 (ASP)和平均互相关峰值（CP)随每个天线发射相位编码信号的码长N的变化关系示意图；
[0044] 图7是使用本发明所得单个序列分别与m序列、p4序列和随机相位序列的峰值旁 瓣电平的对比示意图；
[0045] 图8是使用本发明所得单个序列分别与m序列、p4序列和随机相位序列的积分旁 瓣电平的对比示意图；
【具体实施方式】
[0046] 参照图1，为本发明的基于GS算法的MMO雷达正交相位编码信号设计方法流程 图，该基于基于GS算法的MMO雷达正交相位编码信号设计方法包括以下步骤：
[0047] 步骤1，构造MMO雷达M个天线发射正交相位编码信号的代价函数E，
[0049] 其中，&表示MMO雷达M个天线发射正交编码信号的二阶中心矩，N表示MMO雷 达每个天线发射正交相位编码信号的码长，I为单位矩阵，R n表示MIMO雷达M个天线发射 相位编码信号的协方差矩阵，I I I I表示矩阵范数；
[0050] 步骤1的具体构造过程为：
[0051] 假定一个MIMO雷达具有M个天线，M为自然数，每个天线发射码长为N的正交相 位编码信号，MMO雷达M个天线发射正交相位编码信号矩阵S可以表示为：
[0052] S = [s" s2,…，sm，…，sM]NXM= [y " y2,…，yn，…，yN]T <1>
[0053] 其中，sm= [Sm(l)，sj2)，…，sm(N)]T表示MMO雷达第m个天线发射正交相位编 码信号，y n= [s ?，s2(n), ~sm(n),…，sM(n)]T表示MIMO雷达M个天线发射正交相位编 码信号的第n个码元，
表示MM0雷达第m个雷达天线发射正交相位编码信号 的第n个码元的具体码元值，(K(n)表示MIM0雷达的每个天线发射正交相位编码信号的相 位，且0彡伞m(n)彡2 31。
[0054] 进而，MM0雷达M个天线发射正交相位编码信号的协方差矩阵Rn表示为：
[0056]其中，『表示共轭转置，n G {1，2, ...，N}，m G {1，2, ...，M}，k G {1，2,…，2N}，yk =[Sl (k)，s2 (k)，…，sm (k)，…，sM(k) ]T表示M个雷达天线发射正交相位编码信号的第k个 码元。
[0057] 然后，得到MM0雷达M个天线发射正交相位编码信号的代价函数E，其表示式如下 所示：
[0059] 其中，&表示MM0雷达M个天线发射正交相位编码信号的二阶中心矩，N表示MM0 雷达每个天线发射正交相位编码信号的码长，I为单位矩阵，R n表示MIM0雷达M个天线发 射正交相位编码信号的协方差矩阵，I I I I表示矩阵范数，此处的矩阵范数为F-范数。如无 特别说明，本专利中矩阵范数为F-范数，向量范数为2-范数。
[0060] 步骤2,将MM0雷达M个天线发射正交相位编码信号的代价函数E先转化为原始 极
[0062]小化函数 s. t. | sm(n) | = 1，m G {1，2,…，M}，n G {1，2,…，N}，并将该原始极小 化函数转化为I |ak| |2= 1，k G {1，2,…，2N}
[0063]适合求解的极小化函数 | |DS' -Ql I2，
，Q = [ai, a2, . . .，a2N]T，
[0064] 其中，D表示2NX2N维离散傅里叶变换矩阵，S'表示MMO雷达M个天线发射正 交相位编码信号矩阵的扩展矩阵，I I I I表示矩阵范数，Q表示辅助变量矩阵，ak表示辅助变 量，且ak是满足| |z(? k) | |2= 1的所有解的集合，z(? k)表示MM0雷达M个天线发射正 交发射正交相位编码信号的第n个码元7"的傅里叶变换，s m(n)表示MM0雷达第m个天线 发射正交相位编码信号的第n个码元的具体码元值；
[0065] 步骤2的具体过程为：
[0066] 定义MM0雷达M个天线发射正交发射相位编码信号的第n个码元yn的傅里叶变 换为z (?)，且
，得到MM0雷达M个天线发射正交相位编码信号的第n个 码元的功率谱密度矩阵s(?):
[0068]由式〈3>和式〈4>得出如下两式等价：
[0070] 即MM0雷达M个天线发射正交相位编码信号的代价函数E = 0与功率谱密度矩 阵S(u) =1是等价的，其中I为单位矩阵。
[0071] 因此，MM0雷达M个天线发射正交相位编码信号的代价函数E可以转化为等价代 价函数e 1:
[0075] 忽略与变量无关的常数，式〈7>可以转化为原始的极小化函数：
[0077] s.t. |sm(n) | = l,m G {1, 2, - , M}, n G {1, 2, - ,N}. <8>
[0078] | |ak| |2= 1，k G {1，2,…，2N}
[0079] 其中，= eM'(n)表示MMO雷达第m个天线发射正交相位编码信号的第n个码 元的具体码元值，z〇k)表示MIMO雷达M个天线发射正交相位编码信号的第n个码元y n 的傅里叶变换，S(?k)表示功率谱密度矩阵，ak表示辅助变量，为满足| |z(? k) | |2= 1的 所有解的集合，下标k表示第几个辅助变量，| | | |表示F-范数，M表示MIMO雷达发射天线 数目，N表示每个天线发射正交相位编码信号的码长。
[0080] 为表示方便，定义2NX2N维离散傅里叶变换（DFT)矩阵D的第（k，n)个元素为 ，其中，令
，并将NXM维的MMO雷达发射正交相位编码信号 矩阵S扩展为'
，辅助变量矩阵Q = [ai，a2，…，ak，…，a2N]T，则式〈8>中的原始 的极小化函数可转化为适合求解的极小化函数，BP :
[0081] |DSr -Q| |2 <9>
[0082] 步骤3,利用GS算法 求解步骤2中适合求解的极小化函数I IDS' -Ql I2,得到MIMO 雷达M个天线发射正交相位编码信号矩阵的扩展矩阵S '，进而得到MIMO雷达M个天线发 射正交相位编码信号矩阵S。
[0083]本发明利用gs算法求解步骤2中适合求解的极小化函数|Ids'-Q||2。由于原始 盖师贝格-撒克斯通GS算法的极小化函数是物面和频谱面的一维变换，而本发明中的极小 化函数是多维的，故先将GS算法扩展至多维，并将其物面和频谱面的约束条件进行改进， 以满足本发明优化问题的求解条件。
[0084] 参照图2,式〈9>的求解过程可以分为以下循环迭代的四个子步骤：
[0085] 3.1求解频谱面输出函数Q'（i):
[0086] 利用[0, 2 JT ]之间均匀分布的相位初始化NXM维的MM0雷达M个天线发射正交 相位编码信号矩阵s(i)，并将该正交相位编码信号矩阵s(i)扩展为
，得到频 谱面输出函数Q' (i) = DS' (i)，其中，D表示2N X 2N维离散傅里叶变换矩阵，i表示迭代次 数；
[0087] 具体地，
[0088] 原始盖师贝格-撒克斯通（GS)算法利用交替进行傅里叶变换和傅里叶逆变换，对 物面和频谱面的一维函数进行求解；而本发明中的所求的MM0雷达M个天线发射正交相位 编码信号矩阵S为NXM维，如果要使用利用盖师贝格-撒克斯通（GS)算法求解该极小化 函数，必须将其物面和频谱面的约束条件进行改进，以满足本发明的求解条件。
[0089] 利用[0, 2 JT ]之间均匀分布的相位初始化NXM维的MM0雷达M个天线发射正交 相位编码信号矩阵S(i)(以后每次迭代保留步骤3. 1MM0雷达M个天线发射正交相位编码 信号矩阵S(i))，然后将其扩展为
；区别于原始GS算法进行的一维傅里叶变 换，本发明对MM0雷达M个天线发射正交相位编码信号矩阵的扩展矩阵S' (i)的每一列进 行离散傅立叶变换〇FT)(相当于M维），得到频谱面的输出函数Q' (i)=DS' (i);其中，D 表示2NX 2N维离散傅里叶变换矩阵，i表示迭代次数。
[0090] 3. 2确定频谱面的幅度约束：
[0091] 通过限定离散傅里叶变换矩阵D的每一个行向量对应的2N个通道系数向量的模 值平方和为1，使得满足求解需要的约束条件成立，即频谱面施加的幅度约束为：
[0093] 其中，〇〇表示MM0雷达第m个天线的第k个通道的幅度归一化系数向量， m G {1，2,…，M}，M表示MM0雷达天线个数，N表示MM0雷达每个天线发射正交相位编码 信号码长，N也表示MIM0雷达每个天线的通道系数个数，i表示迭代次数。使用改进的约束 条件更适合求解本发明中的极小化函数；
[0094] 具体地，
[0095] 令频谱面的输出函数
.其中，频谱面的输出函数 Q' (1)的行向量为M列，kG {1，2，...，2叫，]/[表示|/[頂0雷达天线个数，；[表示迭代次 数。
[0096]对M列的行向量>进行幅度归一化得：
[0098] 其中，巧1表示行向量 中第k个通道的幅度归一化系数向量，k G {1，2,… ，2N}，i表示迭代次数。
[0099] 频谱面的输出函数Q' (i)的列向量的2N个元素与离散傅里叶变换（DFT)矩阵D 的2N个通道系数向量一一对应；2N X 2N维离散傅里叶变换矩阵D表示行数为2N、列数为2N 的矩阵，且该离散傅里叶变换矩阵D的每一个行向量对应2N个通道系数向量。
[0100] 此时，MM0雷达每个天线发射码长为2N的正交相位编码信号，同时，MM0雷达每 个天线也对应2N个通道系数向量，通过限定离散傅里叶变换（DFT)矩阵D的每一个行向量 对应的2N个通道系数向量的模值平方和为1，使得满足求解需要的约束条件成立，即频谱 面施加的幅度约束为：
[0102] 其中，〇0表示MM0雷达第m个天线的第k个通道的幅度归一化系数向量， m G {1，2,…，M}，M表示MM0雷达天线数，N表示MM0雷达每个天线发射正交相位编码信 号码长，N也表示MIM0雷达每个天线的通道系数个数，i表示迭代次数。使用改进的约束条 件更适合求解本发明中的极小化函数。
[0103] 3. 3求解物面输入函数S"⑴：
[0104] 对频谱面输出函数Q' (i)的每一列进行逆离散傅里叶变换，得到物面输入函数 S"w，即S"《=DHQ' w，其中，D表示2NX2N维离散傅里叶变换（DFT)矩阵，{}H表示共 轭转置，i表示迭代次数，Q' (i)表示频谱面输出函数；
[0105] 3. 4求解物面的幅度约束：
[0106] 物面的幅度约束即为所求的MM0雷达M个天线发射正交相位编码信号矩阵S，代 价函数E转化式中的约束条件|s m(n) | = l，sm(n)表示MIM0雷达第m个天线发射正交相位 编码信号的第n个码元的具体码元值，m G {1，2,…，M}，n G {1，2,…，N};
[0107]约束条件| sjn) | = 1表明MM0雷达每个天线发射正交相位编码信号为恒模信 号，将MM0雷达M个天线发射正交相位编码信号矩阵S(i)中的每个元素幅度限定为1，此处 定义A为NXM维全1矩阵，且截取矩阵.
，则得到更新的MM0雷达M个天线发 射正交相位编码信号矩阵的扩展矩阵S'(i+1):
[0108] S'（1+1) = J ? exp(jarg(S"（1)))
[0109] 其中，?表示Hadamard积，D表示2NX 2N维离散傅里叶变换（DFT)矩阵， arg(S" (i))表示取S" (i)的角度值，S"（i)表示物面输入函数，i表示迭代次数；
[0110] 具体地，
[0111]由于MM0雷达M个天线发射正交相位编码信号矩阵S(i)是NXM维矩阵，而MM0 雷达M个天线发射正交相位编码信号矩阵的扩展矩阵S' (i+1)是2NXM维矩阵，因此需要对 更新的MMO雷达M个天线发射正交相位编码信号矩阵的扩展矩阵S'(i+1)的维数进行调整， 得到更新的MMO雷达M个天线发射正交相位编码信号矩阵S (i+1)。
[0112] 并且，将MM0雷达M个天线发射正交相位编码信号矩阵S(i)中的每个元素幅度限 定为1，此处定义A为NXM维全1矩阵，且截取矩阵
[0113] 3. 5重复子步骤3. 1~3. 4,对MM0雷达M个天线发射正交相位编码信号矩阵S(i) 进行迭代计算，直到两次迭代之间的误差小于预先设定值e，g卩| |S⑴-S(i+1)| | < e，停止 迭代，该第i次迭代所得的MM0雷达M个天线发射正交相位编码信号矩阵S(i)就是所求的 MMO雷达M个天线发射正交相位编码信号矩阵S。
[0114] 本发明的效果通过以下仿真对比试验进一步说明：
[0115] (1)仿真条件：
[0116] 本发明仿真采用收发共置的集中式MM0雷达系统，假定目标雷达散射截面 (Radar Cross Section，RCS)在积累周期内保持恒定，并忽略MM0雷达M个天线发射编码 信号内的多普勒频移和距离走动问题，仿真中如无特别说明，采用MM0雷达天线数为3,每 个天线发射正交相位编码信号的码长为40,即M = 3, N = 40。
[0117] (2)仿真内容：
[0118] MM0雷达天线个数为M，每个天线发射正交相位编码信号的码长为N的正交相位 编码信号的积分旁瓣电平（Integrated Sidelobe Level, ISL)下界为：
[0119] ISL ^ N2M(M-1) = Bisl〇
[0120] 其中，Bia表示平均积分旁瓣电平（ISL)的下界值。
[0121] 表1给出了本发明所得MM0雷达M个天线发射正交相位编码信号的波形和随机 相位序列信号波形的平均积分旁瓣电平（ISL)与平均积分旁瓣电平（ISL)下界的比较。
[0122] 由表1可以看出，对于各种不同大小的信号波形集，本发明所得多输入多输出 (MIM0)雷达M个天线发射正交相位编码信号波形的平均积分旁瓣电平（ISL)均接近于平均 积分旁瓣电平（ISL)的下界，而随机相位序列信号波形的平均积分旁瓣电平（ISL)则远远 高于平均积分旁瓣电平（ISL)的下界。
[0123] 表1本发明所得波形和随机相位序列的ISL与ISL下界的比较
[0125] 图3为MM0雷达天线数为3,每个天线发射正交相位编码信号的码长为40时，采 用本发明方案而产生的波形集自相关函数示意图，横坐标为时间/采样序号，纵坐标为归 一化幅度；图4为MM0雷达天线数为3,每个天线发射正交相位编码信号的码长为40时， 采用本发明方案而产生的波形集互相关函数示意图，横坐标为时间/采样序号，纵坐标为 归一化幅度。
[01 26] 由图3和图4可以看出，该波形集的平均自相关旁瓣峰值为-16. 6dB，平均互相关 峰值为-14. 3dB ;可见本发明所得波形集的平均自相关旁瓣峰值和平均互相关峰值较低， 相关性能较好。
[0127] 在算法执行效率方面，由于本发明算法采用的离散傅里叶变换（DFT)和逆离散傅 里叶变换（IDFT)，是基于快速傅里叶变换（FFT)和逆快速傅里叶变换（IFFT)进行处理的， 故本发明算法的效率较高，优化耗时较短。而且，当波形集较大时本发明算法仍然适用，如 MM0雷达天线个数M最大可设计到200,每个雷达天线发射正交相位编码信号的码长N最 大可设计到1〇 4。
[0128] 此外，已有文献中提出的混合模拟退火算法、混合遗传算法和序列二次规划算法 的复杂度都较高，优化耗时较长，计算存储量较大，仅适用于波形集较小的情况。
[0129] 表2给出了本发明算法和上述几种算法的执行时间对比，其所得时间均为波形集 大小为M = 3, N = 40时执行10次算法所需的时间平均值。为了使所述几种算法执行时间 对比具有说服力，本文仅列出时间的数量级（计算机硬件条件为：P entium(R)双核处理器， 主频分别是2. 50GHz和2. 29GHz，内存1. 98GB，采用Matlab语言编写程序）。
[0130] 由表2可知，已有算法中执行时间最少的是混合遗传算法，而设计同样大小的波 形集，本发明所需时间是混合遗传算法所需时间的1/100 ;由此可见，本发明算法在执行效 率上较已有算法具有很强的优势，这在工程应用中极为有利。
[0131] 表2本发明算法与文献[Deng Liu Hu]算法的执行时间对比
[0133] 图5给出了每个雷达天线发射正交相位编码信号的码长N固定为128时，波形集 自相关峰值旁辦（Autocorrelation Sidelobe Peak, ASP)、互相关峰值（Crosscorrelation Peak，CP)随MIM0雷达天线个数M的变化关系示意图；图6为MIM0雷达天线个数M固定为 4时，波形集平均自相关峰值旁瓣（ASP)和平均互相关峰值（CP)随每个雷达天线发射正交 相位编码信号的码长N变化的关系示意图。从图5可以看出，当码长固定时，波形集自相关 峰值旁瓣（ASP)和互相关峰值（CP)随MM0雷达天线个数M的增加而增大，并且当MM0雷 达天线个数M增大到一定程度时，波形集自相关峰值旁瓣（ASP)和互相关峰值（CP)的增大 变得不明显；由图6可知，当MM0雷达天线个数M固定时，波形集自相关峰值旁瓣（ASP)和 互相关峰值（CP)与log 2N近似成负线性关系，M表示MM0雷达天线数，N表示每个天线发 射正交相位编码信号的码长。
[0134] 当MM0雷达天线个数M = 1时，使用本发明算法获得的正交相位编码信号矩阵为 单个序列，因为m序列、p4序列的相关性能较好，而且已经广泛应用于工程领域中，而随机 序列是未经算法优化涉及的序列，因此，使用m序列、p4序列与和随机序列与本发明算法所 得单个序列做对比可以说明本算法的优化性能。图7和图8分别给出了本发明算法所得单 个序列与m序列波形、p4序列波形和随机相位序列波形的峰值旁瓣电平和积分旁瓣电平的 对比。
[0135] 由图7和图8得，本发明算法所得单个序列分别比m序列、p4序列和随机相位序 列的自相关峰值旁瓣电平和积分旁瓣电平都要低，随机相位序列波形的自相关峰值旁瓣电 平和积分旁瓣电平最高，P4序列波形比m序列波形的要低。而且，随着MM0雷达天线个数 的增加，使用本发明算法所得发射正交相位编码信号波形、m序列波形、p4序列波形和随机 相位序列波形的自相关峰值旁瓣电平均逐渐减小，且积分旁瓣电平也均逐渐增大。
【主权项】
1. 一种基于GS算法的MMO雷达正交相位编码信号设计方法，其特征在于，包括以下步 骤： 步骤1，构造 MMO雷达M个天线发射正交相位编码信号的代价函数E，其中，Rtl表示MMO雷达M个天线发射正交相位编码信号的二阶中心矩，N表示MMO雷 达每个天线发射正交相位编码信号的码长，I表示单位矩阵，Rn表示MIMO雷达M个天线发 射正交相位编码信号的协方差矩阵，I I I I表示矩阵范数； 步骤2,将MMO雷达M个天线发射正交相位编码信号的代价函数E先转化为原始极小 化函数并将该原始极小化函数转化为适 合求解的极小化函数其中，D表示2NX2N维离散傅里叶变换矩阵，表示MMO雷达M个天线发射正交相 位编码信号矩阵的扩展矩阵，I I I I表示矩阵范数，Q表示辅助变量矩阵，ak表示辅助变量， 且ak是满足I |ζ(ω k) I I2= 1的所有解的集合，ζ(ω k)表示MMO雷达M个天线发射相位编 码信号的第η个码元yn的傅里叶变换，s m (η)表示MMO雷达第m个天线发射正交相位编码 信号的第η个码元的具体码元值； 步骤3,利用GS算法求解步骤2中适合求解的极小化函数I IDS' -Ql |2,得到MIMO雷 达M个天线发射正交相位编码信号矩阵的扩展矩阵S '，进而得到MIMO雷达M个天线发射 正交相位编码信号矩阵S。2. 根据权利要求1所述的基于GS算法的MMO雷达正交相位编码信号设计方法，其特 征在于，步骤3中，所述GS算法求解步骤2中适合求解的极小化函数I IDS' -Ql 12,其具体 子步骤如下： 3. 1求解频谱面输出函数V (i): 利用[〇, 2 π ]之间均匀分布的相位初始化NXM维的MMO雷达M个天线发射正交相位 编码信号矩阵S(i)，并将该正交相位编码信号矩阵S(i)扩展为'得到频谱面 输出函数V (i)=DS< ω，其中，D表示2NX2N维离散傅里叶变换矩阵，i表示迭代次数； 3. 2确定频谱面的幅度约束： 通过限定离散傅里叶变换矩阵D的每一个行向量对应的2N个通道系数向量的模值平 方和为1，使得满足求解需要的约束条件成立，即频谱面施加的幅度约束为：其中，〇n)表示MMO雷达第m个天线的第k个通道的幅度归一化系数向量， m e {1，2,…，M}，M表示MMO雷达天线数，N表示MMO雷达每个天线发射正交相位编码信 号码长，N也表示MIMO雷达每个天线的通道系数个数，i表示迭代次数。使用改进的约束条 件更适合求解本发明中的极小化函数； 3.3求解物面输入函数5"(1): 对频谱面输出函数V (i)的每一列进行逆离散傅里叶变换，得到物面输入函数S"（i)， 即S"ω =D1V ω，其中，D表示2NX 2N维离散傅里叶变换（DFT)矩阵，『表示共轭转置， i表示迭代次数，(^ (1)表示频谱面输出函数； 3. 4求解物面的幅度约束： 物面的幅度约束即为所求的MIMO雷达M个天线发射正交相位编码信号矩阵S，代价函 数E转化式中的约束条件I sm (η) I = 1，sm (η)表示MIMO雷达第m个天线发射正交相位编码 信号的第η个码元的具体码元值，m e {1，2,…，M}，n e {1，2,…，N}; 约束条件I Sm(n) I = 1表明MMO雷达每个天线发射正交相位编码信号为恒模信号，将 MMO雷达M个天线发射正交相位编码信号矩阵S(i)中的每个元素幅度限定为1，此处定义 A为NXM维全1矩阵，且截取矩阵項U得到更新的MMO雷达M个天线发射正交 相位编码信号矩阵的扩展矩阵S' (i+1): S'（1+1)= J Θ exp(jarg(S" ω)) 其中，Θ表示Hadamard积，D表示2NX2N维离散傅里叶变换（DFT)矩阵，arg(S"(i))表 示取S" ω的角度值，S" ω表示物面输入函数，i表示迭代次数； 3. 5重复子步骤3. 1~3. 4,对MMO雷达M个天线发射正交相位编码信号矩阵S(i)进行 迭代计算，直到两次迭代之间的误差小于预先设定值ε，即I |S⑴_S(i+1)| I < ε，停止迭代， 该第i次迭代所得的MMO雷达M个天线发射正交相位编码信号矩阵S(i)就是所求的MMO 雷达M个天线发射正交相位编码信号矩阵S。
【专利摘要】本发明公开了一种基于GS算法的MIMO雷达正交相位编码信号设计，其主要实现思路是：先对多输入多输出(MIMO)雷达天线的发射信号构造代价函数，再把该代价函数转化为适合求解的极小化函数，最后利用改进的盖师贝格-撒克斯通(GS)算法求解该极小化函数，得到要求解的发射信号矩阵；具体来说，首先对多输入多输出(MIMO)雷达的发射信号集构造正交相位编码信号的代价函数，其次利用傅立叶变化将该代价函数转化为适合求解的极小化函数，最后利用改进的盖师贝格-撒克斯通(GS)算法求解其极小化函数，得到多输入多输出(MIMO)雷达天线发射正交相位编码信号矩阵。
【IPC分类】G01S7/02
【公开号】CN104898093
【申请号】CN201510288013
【发明人】苏涛, 杜晓林
【申请人】西安电子科技大学
【公开日】2015年9月9日
【申请日】2015年5月29日
最新回复(0)