低自相关和低峰值平均功率比的序列是通信或信号处理应用中所需要的。本文研究了在频谱约束下具有低自相关和低峰均比的序列的合成。频谱约束限制了每个子载波上的最大允许功率,以避免特定保留频带上的干扰。为了获得性能更好的频域受限序列,本文提出采用模拟退火算法来对序列进行数值优化。基本原理是为序列搜索建立适当的目标函数,调试出适当的退火和停止规则。数值实验结果表明,该算法在满足频谱约束的条件下,具有更好的抑制PAPR和非周期自相关的能力,且该算法收敛性良好。 Sequences with low Autocorrelation (AC) and low Peak to Average Power Ratio (PAPR) are desired in many communications or signal processing applications. This work investigates the synthesis of sequence with the desired low AC and low PAPR under spectral constraints. The spectral constraints limit the maximum allowable power on each subcarrier to avoid interference on particular reserved bands. In order to obtain a Spectrally-Constrained Sequence (SCS) with better performance, Simulated Annealing (SA) algorithm is proposed to optimize the sequence numerically. The basic principle is to establish the appropriate objective function for sequence search and debug the appropriate annealing and stopping rules. Numerical results indicate that this algorithm features better suppression capabilities for both PAPR and aperiodic auto-correlations under the condition of spectrum constraint, and its convergence is good.
低自相关和低峰值平均功率比的序列是通信或信号处理应用中所需要的。本文研究了在频谱约束下具有低自相关和低峰均比的序列的合成。频谱约束限制了每个子载波上的最大允许功率,以避免特定保留频带上的干扰。为了获得性能更好的频域受限序列,本文提出采用模拟退火算法来对序列进行数值优化。基本原理是为序列搜索建立适当的目标函数,调试出适当的退火和停止规则。数值实验结果表明,该算法在满足频谱约束的条件下,具有更好的抑制PAPR和非周期自相关的能力,且该算法收敛性良好。
频域受限序列,峰值平均功率比,非周期自相关函数,模拟退火算法
Lushuang Han, Xudong Li
School of Science, Xihua University, Chengdu Sichuan
Received: Feb. 16th, 2021; accepted: Mar. 11th, 2021; published: Mar. 18th, 2021
Sequences with low Autocorrelation (AC) and low Peak to Average Power Ratio (PAPR) are desired in many communications or signal processing applications. This work investigates the synthesis of sequence with the desired low AC and low PAPR under spectral constraints. The spectral constraints limit the maximum allowable power on each subcarrier to avoid interference on particular reserved bands. In order to obtain a Spectrally-Constrained Sequence (SCS) with better performance, Simulated Annealing (SA) algorithm is proposed to optimize the sequence numerically. The basic principle is to establish the appropriate objective function for sequence search and debug the appropriate annealing and stopping rules. Numerical results indicate that this algorithm features better suppression capabilities for both PAPR and aperiodic auto-correlations under the condition of spectrum constraint, and its convergence is good.
Keywords:Spectrally-Constrained Sequence (SCS), Peak-to-Average Power Ratio (PAPR), Aperiodic Auto-Correlation Function (AACF), Simulated Annealing (SA)
Copyright © 2021 by author(s) and Hans Publishers Inc.
This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).
http://creativecommons.org/licenses/by/4.0/
几十年来,具有良好相关性的序列设计一直是经典的研究课题 [
通常,在传统的序列设计中假定了整个频谱带的可用性,这意味着将序列能量分配给所有连续的载波。对于当今日益拥挤和碎片化的频谱,现代通信和雷达感测中已引入了许多依赖有限频谱资源的应用。认知无线电/雷达被认为是解决可用频谱稀缺的一种有前途的范例,该频谱能够感知和搜索周围环境中的频谱机会,并将其提供给次要用户进行传输。频谱机会被称为频谱限制,它是被许可用户未使用的频谱波段。在某些频谱约束中具有低频谱功率的序列称为频谱受限序列(SCS)。传统的序列通常不适用于频谱受限的系统,如 [
作为 [
该算法是一种全局优化算法,其原理与序列最优相位搜索有一定的内在联系。如果我们建立起SA与序列最优相位搜索之间的联系,就可以绘制SA算法中相似的目标函数(或代价函数)。一旦建立了合适的目标函数,就可以搜索目标函数的最小值,即目标函数的最优解。
基本原理模拟退火之所以如此命名,是因为它类似于固体的物理退火过程,在这一过程中,结晶固体被加热,然后缓慢冷却,直到它达到最规则的可能晶格构型(即,其最小晶格能态),从而没有晶体缺陷。如果冷却进度足够慢,最终的配置会产生具有如此优异结构完整性的固体。模拟退火建立了这种类型的热力学行为和对离散优化问题的全局最小值的搜索之间的联系。此外,它提供了利用这种联系的算法手段。
在应用于离散优化问题的模拟退火算法的每次迭代中,目标函数为两个解(当前解和新选择的解)生成比较值。改进的解决方案总是被接受,而一小部分非改进的(劣质的)解决方案也被接受,希望在寻找全局最优解时避开局部最优解。接受非改进解的概率取决于温度参数,该参数通常不随算法的每次迭代而增加。这样的一个过程叫做Metropolis过程。
简单描述一下这个过程: ( s , f ) 表示一个优化问题。i和j是这个问题的两个解, f ( i ) 和 f ( j ) 是成本因子。从i到j的接收标准由接收概率设计:
P = { 1 f ( j ) ≤ f ( i ) exp ( f ( i ) − f ( j ) T ) f ( j ) > f ( i ) (1)
其中T是控制参数,在控制参数逐渐减小的条件下,SA算法是Metropolis过程的迭代。每个还原过程将产生L个可能的解,以形成长度为L的马尔科夫链,因此每个控制参数(温度)都有L个Metropolis过程。随着控制参数的减小,搜索范围逐渐缩小,算法将停止直到满足收敛条件。
模拟退火的关键算法特征是它提供了一种通过允许爬山移动(即恶化目标函数值的移动)来逃避局部最优的方法。随着温度参数降低到零,爬山运动发生得越来越少,与模拟算法行为的非均匀马尔可夫链相关联的解分布收敛到一种形式,其中所有概率都集中在全局最优解集上。
假设整个频带分为N个频率槽。用 M = [ m 0 , m 1 , ⋯ , m N − 1 ] T 表示“频率槽标记向量”,给出所有N个频率槽的状态。具体而言,如果第k个频率槽有效,即可供使用,则 m k ( 0 ≤ k ≤ N − 1 ) 的值设置为1;否则, m k = 0 ,表示该频率槽被其他应用占用或保留。用 Ω 表示所有空频槽的位置,即 Ω = { k | m k = 0 } 。在本文中, Ω 也被称为“频谱空穴约束”。如果 Ω 是非空的,即 | Ω | > 0 ,就称频谱空穴约束是非平凡的。
在本文中,用矩阵 F N = [ f i , j ] i , j = 0 N − 1 表示N阶傅里叶变换矩阵,即:
f i , j = 1 N ω N − i j , 0 ≤ i , j ≤ N − 1 (2)
其中 ω N = exp ( − 12 π / N ) 。值得注意的是 F N 是酉矩阵,即 F N F N H = I N ,因此N阶离散傅里叶逆变换矩阵是 F N H
设
Ω = { i 0 , i 1 , ⋯ , i k − 1 } ⊂ { 0 , 1 , 2 , ⋯ , N − 1 } (3)
其中 0 ≤ i 0 < i 1 < ⋯ < i k − 1 ≤ N − 1 ,为了便于表达,令
Ω ¯ = { 0 , 1 , 2 , ⋯ , N − 1 } \ Ω = { j 0 , j 1 , ⋯ , j N − k − 1 } (4)
其中 0 ≤ j 0 < j 1 < ⋯ < j N − k − 1 ≤ N − 1 。又设 B = [ B 0 , B 1 , ⋯ , B μ , ⋯ , B N − 1 ] T 是长度为N的频域序列,如果 μ ∈ Ω ,则 B μ = 0 。基于B,我们有相应的时域序列b如下:
b = [ b 0 , b 1 , ⋯ , b μ , ⋯ , b N − 1 ] T = F N H B (5)
如果b的空频率槽特征为 Ω ,则b被称为SCS。在这个意义上, Ω 被称为b的频谱空穴约束。
时域序列 { b n , n = 0 , 2 , ⋯ , N − 1 } 的非周期( C b ( τ ) )自相关函数(Aperiodic Auto-correlation Function, AACF)定义如下:
C b ( τ ) = ∑ n = τ + 1 N b n b n − τ ∗ = r ( − τ ) ∗ , 0 ≤ τ ≤ ( N − 1 ) (6)
设 δ b 为时域序列b的峰值相关旁瓣(Peak Sidelobe Level, PSL),其定义如下:
δ b = max { | C b ( τ ) | | 1 ≤ τ ≤ N − 1 } (7)
时域序列b的峰均比(PAPR)定义如下:
PAPR ( b ) = 10 log 10 ( max 0 ≤ n ≤ N − 1 | b n | 2 1 N ∑ n = 0 N − 1 | b n | 2 ) (8)
在本篇文章中,我们通过引入权重因子 λ ∈ [ 0 , 1 ] 来控制时域序列b的PSL和PAPR的相对权重,因此我们优化问题的目标函数就可以表述为:
E = λ δ b + ( 1 − λ ) PAPR ( b ) (9)
式(1)中表示的Metropolis规则为接收规则,接收概率计算函数为:
P = exp ( − Δ E T k ) (10)
退火规则为:
T k + 1 = α ⋅ T k (11)
在以上设置中, Δ E 为目标函数的变化量; T k 是控制参数; α 为常数,在本次设计中, α 取0.99;k为Metropolis过程中的循环指数,为模拟退火的第k个阶段。
本文是在 [
在温度 T k ( k > 0 ) 下,序列不断“突变”并按(10)的概率接收,直到目标函数分布达到平衡状态。然后根据(11)将温度降至 T k + 1 ,重复“突变”过程,直到在更新后的温度下达到新的平衡状态。在本次设计中,我们选择将每个温度 T k 下的最小目标函数值对应的序列作为下一次温度 T k + 1 下的初始序列,这样能从一定程度上加快算法的收敛速度。
如果在连续三次温度降低中都没有“突变”的相位被接受或 T < ε (其中 ε 为最小停止温度),则退火过程停止。选择的 ε 值是0.001。
算法具体流程如下所示:
步骤1:加载出 [
步骤2:设置初始温度 T 0 、马尔科夫链长度L (固定)、外循环终止迭代条件 ε = 0.001 、降温速率 α ;
步骤3:初始化随机序列:初始序列X为一条单模频域序列,将 [
步骤4:序列“突变”:
a) 确定“突变”位置:在序列X中随机选择两个“突变”位置;
b) 确定“突变”相位值:在 [ 0 , 2 π ] 之间随机选择两个数去替代a中已选位置上的相位值,由此可得到新的单模频域序列 X ′ ,也得到了满足频域受限的新频域序列 B ′ = | B | ⋅ X ′ ,其中 | B | 是 [
步骤5:计算“突变”前后目标函数的变化量 Δ E ,以Metropolis规则决定是否接收 X ′ ;
步骤6:重复步骤4~5;当每个温度下的迭代次数达到阈值,则根据降温速率进行降温,并选取每个温度下使得目标函数E最小的单模频域序列作为下一个温度下的初始序列;
步骤7:重复步骤4~6;达到外循环终止迭代条件时,循环结束。
应用我们提出的算法并使用参考文献 [
Ω = { 14 , 15 , ⋯ , 19 } ∪ { 40 , 41 , ⋯ , 47 } (12)
图1(a)和图1(b)分别显示了 [
图1. 由本文算法和 [
为了进一步优化 [
本工作受到教育部春晖计划项目(No.Z2017065)的支持。从论文选题到完稿,刘子龙教授给我提出了宝贵的意见。
韩露霜,李旭东. 基于模拟退火算法的频域受限序列搜索算法Spectrally-Constrained Sequence Search Based on Simulated Annealing Algorithm[J]. 计算机科学与应用, 2021, 11(03): 543-548. https://doi.org/10.12677/CSA.2021.113055