[发明专利]一种非连续谱的连续调相信号设计方法

说明书

技术领域

本发明涉及一种信号的设计方法。

背景技术

雷达信号设计在雷达系统中起着至关重要的作用，其回波中包含着大量的目标信息，包括空间位置，运动情况等等。雷达工作频率所处的频段往往存在严重的窄带干扰，且干扰具有时变特性。采用稀疏频率信号可避开干扰频段，将发射能量分配到各个离散的“寂静”频带来获得大信干噪比。但频谱的稀疏性会引起距离旁瓣的抬高，容易出现强目标旁瓣淹没弱目标主瓣等不利情况，从而影响目标的检测。由于信干噪比与距离旁瓣以互相冲突的形式耦合，信号设计过程中必定要综合考虑两个指标，在设计过程中可以把信号设计问题转化为一个多目标优化问题。

之前有研究人员针对该问题，通过划分频谱通阻带与人为加权，同时优化功率谱与自相关。但其不足之处在于设计波形时，人为对通阻带进行加权，没有考虑每个带内的干扰起伏及各个频点的细节，不能与多变的干扰环境充分匹配，不利于信号信干比的进一步提高。

发明内容

本发明为了解决现有的波形设计方法中，并没有对每个通阻带的各个频点的细节进行考虑，导致信号设计时不能与多变的干扰环境充分匹配的问题和设计过程中信干比与旁瓣指标之间的比重权衡设计不够灵活的问题。

一种非连续谱的连续调相信号设计方法，包括以下步骤：

步骤1、设采样后的单脉冲内离散信号形式为

x＝[x₁；x₂；...；x_N](1)

其中，N为序列长度；

由于变模信号会使得发射机放大器难以一直工作在饱和状态，将对发射信号的威力范围及信噪比产生影响，因此一般采用恒模信号，在单位发射能量的假设条件下，信号形式由相位决定；信号的第n个采样值x_n表示为

其中，n＝1,2,...,N；Ψ_n为第n个采样值的相位；j表示复数单位；

步骤2、针对信号设计总目标函数

P_Total＝λx^HKx+(1-λ)(2N-x^Hb-b^Hx)(3)

其中，λ表示权重；K为干扰功率谱的相关矩阵；N＝x^HA₁₁x，A₁₁是一个对角线均为1，其他元素趋近于0的N×N矩阵；b＝[b₁；b₂；...；b_N]，b₁～b_N是向量b中的元素；

步骤3、利用梯度迭代法求解波形最优解：

构造迭代格式

其中，第p次迭代的第n个采样点的相位为初相在0到2π随机取值；μ为迭代步长；

有

其中，P_SIR为信干比(SIR)目标函数，P_ACF为距离旁瓣目标函数；Im[·]表示对括号中的值取虚部；Ψ_m表示第m个采样点的相位，每个采样点的相位梯度变化和其他采样点是有相关性的；b(n)是中间变量，表示向量b中的第n个元素；K(n,m)表示干扰功率谱的相关矩阵K中第n行第m列的值；

设置权重λ、迭代步长μ、期望值和收敛容限；按照公式(4)进行迭代，

当总目标函数达到预先的期望值或者前后两次迭代的总目标函数值的差值小于收敛容限时，终止迭代；得到当前权重下的最优信号，即最优相位组合，以及其对应的信干比与旁瓣指标值，记录P_SIR与P_ACF；

步骤4：改变权重λ，重复步骤1-3，得到另一权重下的最优相位组合以及信干比与旁瓣指标；从1到0遍历权重λ，得到一系列信干比与旁瓣的值，它们是在不同权重下一一对应的；以其中一个指标为横轴，另一指标为纵轴，得到的曲线即为Pareto最优曲线，该曲线上反应了不同目标函数间的代价变化关系；

将Pareto最优曲线作为多目标优化中辅助选取目标函数间最优权衡的工具设计非连续谱的连续调相信号。

本发明具有以下有益效果：