一种基于广义Jaccard...AMP压缩信号快速恢复算法_时天昊.pdf

第 卷 第期空军工程大学学报 年月 收稿日期：基金项目：国家自然科学基金（）作者简介：时天昊（），男，山东菏泽人，硕士生，研究方向为民航通讯频谱态势感知。：引用格式：时天昊，白银山，沈堤，等一种基于广义 系数的改进 压缩信号快速恢复算法空军工程大学学报，（）：，（）：一种基于广义 系数的改进 压缩信号快速恢复算法时天昊，白银山，沈堤，沈志远（南京航空航天大学民航学院，南京，；部队，新疆和田，；空军工程大学空管领航学院，西安，）摘要压缩感知算法作为一种信号处理方法，可以解决机场终端区实时频谱监测的问题。基于稀疏度自适应匹配追踪（）信号重构算法，引入了广义 系数、平均相关系数、变步长思想，提出了 算法。在算法的原子筛选部分引入广义 系数可以减少原子混淆导致的精度下降问题，平均相关系数的引入可以避免测量矩阵 系数的计算，降低了算法的复杂度，变步长思想中的大步长迭代，小步长靠近的步骤使得算法的效率及精度都大大提升。采用一维高斯随机稀疏信号作为测量信号进行仿真，可以有效的模拟机场终端区经过能量检测后的测量信号，经过仿真，在不同的测量数、稀疏度情况下的算法重构成功率的表现明显优于传统压缩感知算法，且与 算法相比，在重构误差和算法时间方面的表现均有显著提升。关键词压缩感知；信号重构；广义 系数；变步长 中图分类号 ；文献标志码文章编号 （），（，；，；，），（），；根据“十四五”民用航空发展规划，预计到 年，民用运输机场数量达到 个以上，保障起降架次 万。民航的快速发展导致机场终端区的频谱资源需求及安全隐患大大增加，对频谱的快速高效监测可以指导频谱的智能管理。在机场终端管制区认知无线电技术可以包括：利用信号能量检测实现信号状态感知及稀疏化处理，将电磁信号转换成二进制信号，并将二进制信号作为测量信号，引入 宽 带 压 缩 频 谱 感 知 模 型，利 用 压 缩 感 知（，）算法进行信号传递后的重构工作。压缩感知可以在较小的采样频率的情况下对频谱进行高精度监测，同时在广播式自动相关监视（，）的来波信号估计、干扰抑制、空管数据处理等方面也有着广泛的应用。压缩感知算法是一种高效的信号采集及重构算法，该算法在采样频率远低于奈奎斯特频率情况下依然有着优秀的重构精度，大大降低了信号的传输量，在图像重建、远距离通信等领域有广泛的应用前景。对于压缩感知的研究主要集中在测量矩阵构造、信号稀疏表示、重构算法个方面。重构算法的选择直接影响着信号恢复的精度及速度。重构算法主要有组合类算法、凸优化类算法和贪婪追踪类算法种。贪婪类算法凭借复杂度低，重构效果好，收敛速度快等特点得到了广泛应用。贪婪算法需要较多的先验信息，然后在每次迭代过程中，不断地将迭代目标像原始信号靠近，并通过设定一个收敛条件来结束算法，从而达到具有较高重构精度的信号。常用的贪婪类算法包括匹配追踪、正交匹配追踪、正则化正交匹配追踪、分段正交匹配追踪、压 缩 采 样 匹 配 追 踪、子 空 间 追 踪 算法。但是上述介绍的算法一般需要先验参数稀疏度来确定迭代的次数，在实际运用中不具备较高的应用性，为了解决这个问题，文献 提出了稀疏度自适应匹配追踪算法（，），该算法在稀疏度未知的情况下依旧具有较高的重构精度。算法在迭代过程中的原子筛选标准均采用的是内积匹配准则，有时会导致个相似原子信息的丢失，降低重构精度。算法 在原子筛选阶段采用了正则化回溯的思想，用正则化进行原子的二次筛选，大大提高了重构精度，但是由于增加了二次筛选，其算法运行时间也大大增加。文献 提出的 算法用 系数代替了 算法中的内积匹配准则进行原子筛选，系数可以更好地突出残差信号中的重要元素组成部分，更准确的选择符合信号重构的原子，大大提高了重构精度。同时，算法在运行前需要进行一个固定步长设置，过大的步长设置会导致算法的精度达无法达到理想效果，而过小的步长设置会降低算法的运行效率。文献 将变步长的思想引入到算法之中，作者利用 函数作为变步长的判断条件，利用大步长逼近，小步长渐进的方法提高了算法的重构精度及运行效率。不仅采用了正则化回溯的思想，还将抛物线函数和变步长思想结合，函数的变化率与算法迭代过程中步长的大小变化呈正相关，取得了良好的效果。在变步长思想中增加初始稀疏度估计有助于大小步长的设置，降低算法的迭代次数，根据定理：稀疏度小于等于测量数的时，信号可以完全准确地重建。将初始稀疏度大小设置为测量数的进行取根号处理，与未设置初始稀疏度情况相比重构效率大大提高。由于民航终端区频谱资源态势复杂，以上方法在解决频谱感知等问题时的精度及效率还未达到要求，为解决上述问题，本文提出了一种基于 算法改进的算法：（代表广义 系数，代表平均相关系数，代表变步长思想）。该算法在原子筛选过程中，利用广义 系数代替了 算法中的内积匹配准则，并引入平均相关系数对初始稀疏度进行预测来确定算法迭代的初始步长，同时引入了变步长的思想来优化算法的重构效率和重构精度。空军工程大学学报 年压缩感知 压缩感知的概述压缩感知的流程如图 所示，主要分为 步：信号稀疏表示、信号采样、信号重构。图压缩感知流程图常用的信号稀疏化处理包括小波变换、曲波变换、离散傅里叶变换 等。图中，是维稀疏表示阵，是长度为的一维信号，经过稀疏表示阵处理后，信号就具备了稀疏性，数学表达为：（）式中：为压缩信号；表示测量矩阵；为的感知矩阵。如果一个信号本来就具备稀疏性且为一维的，那就可以直接使用一个测量矩阵来对信号进行测量，数学表达为：（）式中：，为测量矩阵；是被测量矩阵测量后的测量信号。在信号的采样过程中，会丢弃原始信号大量的无用信息，因此测量信号会有一定程度上的能量缺失，和 证明了稀疏度为的测量信号要从个测量值准确的进行信号重构，就要保证测量矩阵满足有限等距特性（，）。信号重构过程即对式（）的逆向求解过程。在式中，因此的解会有无穷多个，考虑到具有稀疏性，可以把信号重构过程转换成范数问题进行计算，信号重构的数学模型可表示为：.（）式中：是一个非凸函数，代表输入信号中非零元素数量，对于公式的求解属于（，）难问题，无法直接求解，文献 将 范数最小化问题转换成了更加简单的 范数最小化问题，实现了从非凸向凸的转变，计算式如下：.（）算法 算法是贪婪类算法的一种，该算法不需要稀疏度等先验信息，更符合自然界信号的实际情况。算法步骤流程如下：表 算法步骤输入参数：测量信号；感知矩阵；步长；输出参数：恢复信号?初始化：估计支撑集，初始残差；支撑集长度，迭代次数；.计算 ，将排序后符合条件的个原子的索引加入索引集；.确定候选集：令；.用 最 小 二 乘 法 求 解 得 到 重 构 信 号?；.将重构信号?里个最大原子索引储存在集合中；.更新残差?；.停止 迭 代 条 件：如 果 ，转 步 骤；如 果 ，转步骤；若上述个条件均不满足，则，如果终止迭代进入步骤，否则转步骤；.更新步长，转步骤；.输出信号估计值?。算法 算法改进 基于广义 系数的原子筛选传统的 算法采用了内积匹配准则 来作为原子筛选标准，原子与迭代过程中的残差匹配度越高，其内积计算结果越大，将该原子加入到支撑集中就可以更快地逼近原始信号，可以达到较为优秀的重构精度，内积匹配准则表达式为：（，）（）式中：（，）；，（）。由于内积匹配准则更多地考虑了原子和信号残差之间角度的差距，缺乏对信号本身的关注，这会导致在原子筛选过程中会忽略个相似度较高的原子，造成信号丢失，降低重构质量。为了避免上述情况的出现，引用了广义 系数来代替内积匹配准则，计算式如下：（，）（）由式（）可以看出，广义 系数可以反应第期时天昊，等：一种基于广义 系数的改进 压缩信号快速恢复算法出个向量之间的相似度。分母上放大了个原子的差异部分，减去了两者的相似部分，使得原子不容易混淆，减少了重要信息的丢失，使得重构精度大大提升。初始稀疏度的预估计准确合适的初始稀疏度估计可以指导设置初始步长，有效地减少迭代次数，提高运行效率，本文采用的稀疏估计策略基于如下性质：当满足 条件的测量矩阵在稀疏度下有以下规定时：（）式中：，则为真。式（）中 系数计算复杂，不容易通过计算得出准确结果。与矩阵的正交性呈正相关，其值的大小与信号采样过程中保留的原始信号信息有关，与越靠近代表保存的原始信息更完整。根据上述特点，引入了平均相关系数来代替 系数，可以在保证稀疏度估计精确度的同时大大降低计算复杂度，如式（）所示：（）（）（）式中：（）是测量矩阵的平均相关系数。平均相关系数作为测量矩阵的评价指标，定义如下：（），（，）（，），（，）（）式中：矩阵，（，）为矩阵中第行第列的元素，是一个特定的阈值，它的大小由矩阵结构决定，设置阈值的大小为（）。由此，可以对初始稀疏度进行估计，初始化，代入式（）进行验证，如果满足条件，则，继续代入公式验证，否则就结束迭代并输出。迭代过程中的变步长处理在测量矩阵满足 条件的前提下，测量信号与原始稀疏信号满足以下关系：.（）在算法每次迭代过程中重构信号?都会逐渐逼近原始信号，?的数值会逐渐增大，由于作为一个定值，测量信号与重构信号之间的能量比逐渐降低，最后的能量比会落在公式所规定的区间。因此，在.时，重构信号与测量信号之间还有较大的差距，此时选择较大步长可以增加算法运行效率，大步长的变化规则采用（），为步长变化参数，数值大于，可以表示在大 步 长 迭 代 阶 段 保 证 较 高 的 运 行 效 率；在 .时，进行小步长选择阶段，小步长采用可以确保较高的重构精度。虽然考虑了大小步长变化的方法来平衡重构效率及重构精度之间的平衡，但是由于稀疏度的未知性，无法准确地确定算法迭代的次数，若在最后一次迭代计算时，步长仍然有着较大的值，可能会导致过估计 的现象，为了避免出现上述情况，结合回溯的思想，本文提出了一种步长的处理方法。若最后一次迭代过程中，步长，就不进行额外处理；若，则将迭代中的参数回溯至上次迭代的数据，然后以步长继续迭代，直到再一次满足迭代结束条件，经过回溯处理后，重构信号的误差会大大降低。算法步骤根据以上改进措施，本文在 算法的基础上提出了 算法，步骤流程如下：表 算法步骤输入参数：测量信号；测量矩阵；输出参数：恢复信号?初始化：稀疏度；估计支撑集，残差；计算测量矩阵的平均相关系数（）；判断：如果（）（），则，转步骤；.计算残差；.设置迭代初始步长，支撑集长度，迭代次数；.计算 ，（），将符合条件的个原子得索引加入索引集，令；.用 最 小 二 乘 法 求 解 得 到 重 构 信 号?；.将重构信号?里个最大原子索引储存在集合中；.更新残差?；.停止迭代条件：设定阀值，如果 且，则停止迭代；如果 且，则转步骤；若上述条件均不满足，转步骤；.支撑集回溯，步长，迭代次数，转到步骤；.如 果 ，转 步 骤 ，否 则，转 步骤；.如果?.，（），。否则，；，迭代次数，转步骤；.，迭代次数，转步骤。空军工程大学学报 年仿真结果及对比分析为了比较本算法和其他算法重构精度及重构效率上的优越性，本文选择一维高斯随机稀疏信号来模拟算法所需的测量信号，该信号可以有效的模拟机场终端区信号经 过 能量 检测 处 理 后 的 信 道 状态。测量矩阵选择为高斯随机矩阵进行信号采样，实验中所采用的对比算法为实验选用、，参数设置见表。表仿真参数表参数类型参数值测量信号长度 固定步长，稀疏度 ，重构成功阈值 测量数 ，本文所有实验运行的平台为 ，计算机配置为 ，。为避免额外因素干扰，程序运行期间保证后台其余程序全部关闭。不同观测值下的算法重构成功率设定采样信号的稀疏度为固定值，测量数的取值以 为起始，位步长向上增大直至 。设定当?时，即信号重构成功，重构成功率的计算如下所示：重构成功率信号重构成功次数实验总次数 每种算法实验重复 次，运行结果见图。图不同观测值信号重构成功率从图中的曲线我们可以看出，随着测量数的的增加，各个算法的信号重构成功率都有较大的提升。时，相比于 有着更为高的信号重构成功率，这是因为采用广义 系数保留了更多有用的原子，小步长阶段的算法迭代过程带来了更高的重构精度，当 时，除了 的其余算法信号重构成功率基本达到了 。不同稀疏度下的算法重构成功率本实验探究随着稀疏度的增加，信号重构成功率的变化。采样信号固定测量数 ，将稀疏度从 开始，步长为增加到。设定当?时，即信号重构成功，每种算法实验重复 次，实验结果见图。图不同稀疏度信号重构成功率从图可以看出，在测量值为定值的情况下，随着稀疏度的增加，所有的算法信号重构成功率均呈现不同程度的下降。算法下降速度远远慢于其它算法，在稀疏度 的时，依然保持着一定的信号重构成功率，而 算法的重构成功率已经接近。不同算法的重构误差比较为了验证 算法及 算法在信号重构成功情况下其重构精度的优劣，我们引入了重构误差这个指标来表示重构信号和原始信号的相似性，重构误差越大，代表信号重构效果越差，计算公式如下：?（）选择 ，共个稀疏度。算法的步长设置由 ，选取步长为和两种情况进行对比分析。对每个特定值条件下仿真 次，统计仿真结果计算信号的平均重构误差，见表。表重构误差对比表方法稀疏度重构误差 由表格可得，随着稀疏的的增加，、和 的重构误差都会逐渐增大，相比于，采用了更小的步长设置，在算法迭代过程中具有更高的精度，因此其重构误 差 小 于 的 重 构 误 差。的重构误差略优于 算法，这是因为 算法利用广义 系数筛选原子、小步长逼近等措施，使得该算法具有较高的精度。第期时天昊，等：一种基于广义 系数的改进 压缩信号快速恢复算法 不同信号采样率情况下的重构精度分析信号采样率的大小是测量数与信号长度的比值，代表了采集过程中的所需的信号采集数量。为了证明算法在不同信号采样率情况下的重构性能，引入均方误差（）作为评价指标，其计算公式如下：（?）（）在压缩感知中，测量信号是稀疏信号完全恢复的充分条件。因此在本次仿真实验中，设置测量信号稀疏度。采用 算法作为对比算法，其步长设置分为和两种情况，探讨采样率为，情况下的改进算法与对比算法的重构精度，结果见图。图不同采样率下的算法均方误差比较图从图可以看出，随着采样率的增加，种算法的均方误差数值在逐渐减小，当采样率大于 时，所有的算法均可以完成高精度信号重构，但在低采样率的信号重构过程中，算法相比于 算法，具有更高的误差控制能力，表明了算法在低采样率情况下也具有较高的算法稳定性。算法重构时间的比较为了验证 算法的重构效率的优越性，使用 软件进行算法用时统计，与步长设置为、的 、算法进行比较。采样信号固定测量数 ，设置稀疏度从开始，步长为 增加到。进行循环 次处理计算平均运行时间，仿真结果见图。图平均运行时间对比图由图可以看出，随着稀疏的 的增加，种算法的平均运行时间都在增加，这是由于信号重构的内容增多，其耗费的时间也在逐渐增加。的平均运行时间大于其余两种算法，是由于其步长设置较小导致的效率偏低。算法在效率方 面 略 优 于 算 法，这 是 因 为 算法对初始稀疏度进行了预测，且采用了变步长的思想，在大步长阶段的快速迭代，大大降低了算法所需时间。结语本文基于的 算法，通过融合广义 系数、平均相关系数、变步长思想，提出了 算 法。相 比 于，等传统压缩感知算法，在使用二进制信号模拟 信 号 能 量 检 测 后 的 机 场 电 磁 信 号 时，算法在同稀疏度不同观测值、不同稀疏度同观测值下的情况下，算法重构的成功率表现更为优秀，且与其原算法 相比，在重构误差、低采样率信号稳定性、重构时间等方面均有了明显的提高。该改进算法未来可在机场终端区频谱分析、波达方向估计等领域广泛应用。参考文献 ，（）：毕号旗，向新，李娜，等 系统中基于压缩感知的慢衰落航空稀疏信道估计空军工程大学学报（自然科学版），（）：，（）：，：，：，（）：，空军工程大学学报 年 ：，（）：，（）：，（）：，：，（）：，（）：，：，：，（）：，：，：，：，（）：，：杜秀丽，胡兴，顾斌斌，等基于变步长的正则回溯 压缩感 知 重 构 算 法 计 算 机 应 用 研 究，（）：，（）：，：黎浩民基于深度学习的图像压缩感知重建算法研究 广州：广东工业大学，苏扬基于稀疏正则化的电磁勘探数据反演研究 长春：吉林大学，曹静杰，杨志权，杨歧焱一种基于压缩感知的地震数据重建方法及其在城市活断层地震勘探中的应用 地球物理学进展，（）：任杰基于深度学习的压缩感知技术成都：电子科技大学，（）：张连娜，张慧萍，李荣鹏，等基于二次筛选的回溯广义正交匹配追踪算法的稀疏信号重构计算机与现代化，（）：，（）：张永立，黄 芳，范 志勇压缩感知重构算法研 究科技视界，（）：，：程鑫淼基于 的终端管制区频谱分析研究南京：南京航空航天大学，（），（）：（编辑：徐楠楠）第期时天昊，等：一种基于广义 系数的改进 压缩信号快速恢复算法