一种具有谱约束的连续相位调制波形的设计方法技术

本发明专利技术提出一种连续相位调制波形的设计方法，其频谱具有频带内可任意设置阻带和有效抑制频带外能量的特性。具有该特性的波形能够有效利用电磁频谱并增加功率利用率，在雷达和通信系统中有重要应用。该方法使用离散形式近似表示连续相位波形，并对相位序列进行优化，再通过插值产生连续相位调制波形。通过该方法设计出的优化波形具有良好的带内阻带、带外抑制和恒模特性，简称为恒模谱约束特性。

【技术实现步骤摘要】

本专利技术属于雷达波形设计领域。特别涉及连续相位调制波形设计。
技术介绍
在现代电磁波频谱日益拥挤的情况下，需要减小不同电磁设备之间的干扰、扩大系统带宽，因此雷达系统对发射波形的频谱利用有效性的需求日益增加。频谱利用有效性要求发射波形频谱的带内具有指定阻带且带外具有良好抑制特性，将频谱满足上述条件的称为满足谱约束的波形，简称为谱约束波形。另外，雷达系统的能量有效性通过发射恒模信号来实现，这样使发射机始终工作在最大功率。由于编码信号由于具有更多的设计自由度，提升了性能和灵活性，目前对基带离散编码信号做了广泛研究，已有多种算法设计满足带内阻带要求的编码序列。然而，瞬时相位变化导致具有sin(x)/x包络的频谱旁瓣的扩展，不能够严格抑制频带外的能量。因此在编码信号的基础上，减慢其瞬时相位变化，得到连续相位调制波形可以满足上述恒模谱约束要求。对连续相位波形已有一些研究，在文献“Highly bandlimited radar signals,proceedings of 2002IEEE Radar Conference,pp.220-226,2002.”和文献“A spectrally clean transmitting system for solid-state phased-array radars,proceedings of 2004IEEE Radar Conference,pp.140-144.2004”提出高斯sinc函数作为核函数来对离散编码信号进行卷积，设计带限信号，但设计出的信号不能满足恒模特性。在文献“CPM-based radar waveforms for efficiently bandlimiting a transmitted spectrum，proceedings of 2009 IEEE Radar Conference,pp.1-6,2009”中，提出对相位差卷积核函数后积分的方法得连续相位信号。但该方法执行较为复杂，且未考虑带内任意设置阻带约束。
技术实现思路
本专利技术所要解决的技术问题是，提供了一种连续相位波形设计方法，其频谱具有频带内可设置阻带和有效抑制频带外能量的特性。本专利技术为解决上述技术问题所采用的技术方案是，提出用长度大于波形时宽带宽积的序列来近似表示谱约束波形。在此基础上，提出对该序列的设计方法，然后通过插值得到在一定范围频谱内具有严格频带限制的连续相位调制信号；包括以下步骤：步骤1：通过虚拟采样，获取有N个采样点的相位调制信号s(n)；步骤2：对s(n)进行L倍过采样，获得由有NL个采样点的调制信号x(n)；步骤3：采用公式：对x(n)进行频谱优化；其中R为一个NL×NL维的矩阵，矩阵华中第m行n列上的元素为Rmn：(fk1,fk2)表示一个阻带，Ns表示阻带的个数，ωk表示第k(k＝1,…,Ns)个阻带对应的权值，，fk2表示第k个阻带的终止频率，fk1表示第k个阻带的起始频率；步骤4：采用循环迭代法对相位进行优化；步骤4.1：对步骤3获得的经过频谱优化的调制信号x(n)的相位进行初始化；步骤4.2：计算辅助变量α，α＝Bx(n)，其中B＝U(I-Σ)U*，将R进行特征值分解，即R＝UΣU*，其中Σ为将R进行分解后的的特征值构成的矩阵，U为将R进行分解后的特征向量；步骤4.3：采用公式x(n)new＝exp(j arg(α))计算出新的调制信号x(n)new；步骤4.4：如果||x(n)new-x(n)||2<ε，ε是根据实际情况设定的阈值，||·||2表示向量的2范数，则停止迭代，最优解x(n)opt＝x(n)new；否则令x(n)＝x(n)new，重复步骤4.2-4.4，直到满足条件为止；步骤5：对步骤4获得x(n)进行差值，获取最终的连续调制信号x(n)。本专利技术可以有效设计带内具有指定阻带带外具有良好抑制的发射波形。从图1和图2的对比可以看出，当过采样率L增加时，得到的优化波形的带内阻带就越低，带外抑制效果越好，即谱约束特性就越好。附图说明图1为优化得到的波形频谱图；图2为得到的优化波形插值后得到的波形频谱图。具体实施方式本专利技术主要采用计算机仿真的方法进行验证，所有步骤、结论都在MATLAB-R2010a上验证确认设波形的时宽带宽积为400，先用过采样的波形进行优化，设过采样的倍数为L，分别令L＝1,1.5,2,3。用循环迭代法进行优化，迭代次数为200次。具体实施步骤如下：步骤1初始化参数；通常情况下，编码信号的码片数量与时域信号的时宽带宽积相等，记为N。则长度为N的相位调制信号s，其各码片表示如下：其中φ(n)(0≤φ(n)<2π)是信号的第n个子脉冲的相位，其中n＝1,…,N。φ表示φ(n)组成的向量，即φ＝[φ(1),…,φ(N)]T。不失一般性，采用归一化的频率，将信号带宽记为1，则s(n)所在频带记为[-1/2,1/2]。步骤2过采样；为了在时域上给定的时间间隔内提供更平滑的相位变化，需要增加采样率。假设将采样率增大到原来的L倍，即将每个码片分为L个子码片。令过采样后的波形相位表示为对应的波形记为合理选择L的值与方法的性能和效率密切相关。当L选择较小时，不能够有效表示连续的相位；当L选择较大时，增加波形设计方法的计算量。我们将在仿真部分通过不同L的取值观察仿真结果给出合适的L值。步骤3构造目标函数；过采样后信号频谱范围为[-L/2,L/2]，需要考虑[-L/2,L/2]频率范围内的理想频谱特性。带外的频率需要全部抑制，带外的频率可记为ΩO＝(-L/2,-1/2)∪(1/2,L/2)带内指定的阻带频率为其中(fk1,fk2)表示一个阻带，Ns表示阻带的个数。因此需要抑制的频谱范围写为Ω＝ΩO∪ΩI。定义NL×NL维矩阵R，其中R矩阵的第m行n列上的元素为一般而言，R是一个秩亏缺的矩阵，即rank(R)<NL。将R进行特征值分解，即R＝UΣU*，其中Σ是R的特征值构成的矩阵。定义矩阵B，B＝U(I-Σ)U*。为使阻带Ω内的干扰能量最小，设计代价函数：s.t.|x(n)|＝1,n＝1,2,…,NL步骤4循环迭代法计算最优相位；4.1、l＝0，随机初始化恒模向量，对过采样相位进行初始化，如对的NL个相位独立随机分布。4.2、令l＝l+1，计算辅助变量α，α＝Bx；4.3、计算使代价函数最小的xnew，xnew＝exp(jarg(α))；4.4、计算||xnew-x||2，如果||xnew-x||2<ε(ε是一个给定的值，||·||2表示向量的2范数)，则停止迭代，最优解xopt＝xnew；否则令x＝xnew，重复步骤4.2-4.4，直到满足条件为止。步骤5插值；由于考虑到计算量的问题，L取值有限，因此上一步骤优化得到的相位不能保证可以近似表示连续相位信号，因此我们对得到的过采样相位α进行M倍插值。插值的方法有样条曲线 法、线性法等。本专利技术的优点是能够获得具有带内阻带和带外能量抑制的连续相位调制波形，增加功率利用率，减小不同电磁设备间的干扰。该方法简便，易于实现。本专利技术可应用于雷达连续相位调制波形设计等领域。上述仿真实验用线性调频信号作为初始信号得到的优化波形的幅度谱为图1；对得到的相位本文档来自技高网...

【技术保护点】
一种具有谱约束的连续相位调制波形的设计方法，包括以下步骤：步骤1：通过虚拟采样，获取有N个采样点的相位调制信号s(n)；步骤2：对s(n)进行L倍过采样，获得由有NL个采样点的调制信号x(n)；步骤3：采用公式：对x(n)进行频谱优化；其中R为一个NL×NL维的矩阵，矩阵华中第m行n列上的元素为Rmn：Rmn=Σk=1Nsωk∫fk1fk2ej2πf(m-n)df=Σk=1Nsωkej2πfk2(m-n)-ej2πfk1(m-n)j2π(m-n),m≠nfk2-fk1,m=n]]>(fk1,fk2)表示一个阻带，Ns表示阻带的个数，ωk表示第k(k＝1,…,Ns)个阻带对应的权值，，fk2表示第k个阻带的终止频率，fk1表示第k个阻带的起始频率；步骤4：采用循环迭代法对相位进行优化；步骤4.1：对步骤3获得的经过频谱优化的调制信号x(n)的相位进行初始化；步骤4.2：计算辅助变量α，α＝Bx(n)，其中B＝U(I‑Σ)U*，将R进行特征值分解，即R＝UΣU*，其中Σ为将R进行分解后的的特征值构成的矩阵，U为将R进行分解后的特征向量；步骤4.3：采用公式x(n)new＝exp(jarg(α))计算出新的调制信号x(n)new；步骤4.4：如果||x(n)new‑x(n)||2<ε，ε是根据实际情况设定的阈值，||·||2表示向量的2范数，则停止迭代，最优解x(n)opt＝x(n)new；否则令x(n)＝x(n)new，重复步骤4.2‑4.4，直到满足条件为止；步骤5：对步骤4获得x(n)进行差值，获取最终的连续调制信号x(n)。...

【技术特征摘要】
1.一种具有谱约束的连续相位调制波形的设计方法，包括以下步骤：步骤1：通过虚拟采样，获取有N个采样点的相位调制信号s(n)；步骤2：对s(n)进行L倍过采样，获得由有NL个采样点的调制信号x(n)；步骤3：采用公式：对x(n)进行频谱优化；其中R为一个NL×NL维的矩阵，矩阵华中第m行n列上的元素为Rmn： R m n = Σ k = 1 N s ω k ∫ f k 1 f k 2 e j 2 π f ( m - n ) d f = Σ k = 1 N s ω k e j 2 πf k 2 ( m - n ) - e ...

【专利技术属性】
技术研发人员：崔国龙，杨亚，陈健，李溯琪，黎亮，姬红敏，易伟，孔令讲，杨晓波，
申请(专利权)人：电子科技大学，
类型：发明
国别省市：四川;51