本发明专利技术公开了一种稀疏信号的重构方法，主要解决从观测向量重构原始稀疏信号速率低的问题。该方法利用约束目标函数的可分解性，将优化约束目标函数问题分解为一系列小的约束目标函数来优化，提高重构速率。它包括获取原始数据部分、重构原始稀疏信号部分和重构准确度评价部分，其中：获取原始数据部分包括产生原始稀疏信号、观测矩阵和观测向量；重构原始稀疏信号部分主要包括设定无约束目标函数，推导约束目标函数、使用序列最小优化方法分解约束目标函数、计算重构信号和对重构信号进行除偏处理；重构准确度评价是比较均方误差的大小。本发明专利技术的重构方法可在保证重构准确率的前提下，提高重构速率，可应用于压缩感知等领域的稀疏信号重构问题。
【技术实现步骤摘要】
本专利技术属于信号处理
，特别是涉及稀疏信号重构，可用于压縮感知 或者相反问题上。
技术介绍
随着数字化信息时代的发展，信息来源的模拟化和信息处理工具的数字化越 来越占有主导地位，信号的采样是两者联系的桥梁。傅立叶变换和奈奎斯特采样 定理表明，欲无失真从离散信号中恢复信号，其采样速率必须是信号带宽的2倍。 然而随着人们对信息需求量的增加，携带信息的信号的带宽越来越宽，据此理论 为基础的信号处理框架，要求的采集速率和处理速度越来越高，对宽带信号处理 的困难也越来越大，给信号采集、存储、传输和处理等带来了巨大的压力。近几年来出现的一种新颖的理论——Compressed sensing,也称为Compressive sampling,得到了快速发展。该理论异于近代奈奎斯特采样定理，它指出利用随 机观测矩阵可以把一个稀疏或可压縮的高维信号投影到低维空间上，然后再利用 这些少量的投影通过解一个优化问题就可以以高概率重构原始稀疏信号，并证明 了这样的随机投影包含了原始稀疏信号的足够信息。于是一个可压缩信号能够用 比奈奎斯特采样要求的信号长度小得多的信号来准确表示。该理论利用其它变换 空间描述信号，建立新的信号描述和处理的理论框架，使得在保证信息不损失的 条件下，用远低于奈奎斯特采样定理要求的速率采样信号而又可以完全恢复信， 可极大地降低信号的采样频率、极大地减少数据存储和传输代价、极大地减少信 号处理时间和计算成本。它带领信号处理进入一个新的革命时代。在压縮感知理论中，由于观测数量《远小于信号长度iV，因此不得不面对求 解欠定方程组y-Axe7^,《〈iV的问题。乍一看，求解欠定方程组似乎是无望的， 但是，人们发现由于信号jce及w是稀疏的或可压縮的，也就是说，它在某个固定基V/下可以精确地写成少量的基向量的叠加。这个前提从根本上改变了问题，使 得求解压縮感知理论中欠定方程组问题可解，而观测矩阵具有RIP性质也为从《 个观测值中精确恢复信号提供了理论保证。目前，存在很多压縮感知中重构原始 稀疏信号的方法，例如GPSR、 state-of-the-art、 1ST、 /1—"和OMP等方法。IST最初应用于反巻积方面，对于解决压縮感知问题，观测矩阵的行数一般小于列数，ist并不有效。njs是基于共轭梯度下降算法提出的一种方法，每一次循环都需要寻找下降方向，浪费了时间。OMP方法虽然节约了大量时间，但是对 于压縮感知中不一致情况，它对噪声的鲁棒性太低。GPSR和state-of-the-art当循 环次数多时，浪费了时间。
技术实现思路
本专利技术的目的在于克服上述技术的不足，提出一种稀疏信号重构方法，以提高 重构原始稀疏信号的速率和准确度。为实现上述目的，本专利技术包括如下步骤(1) 利用计算机得到原始稀疏信号fei^、观测矩阵Aei "w和观测向量y"、(2) 从含有K个元素的观测向量y中，按如下步骤重构含有W个元素的原始 稀疏信号f :2a)设定重构原始稀疏信号f的无约束目标函数<formula>formula see original document page 6</formula>式中，xe及w为重构信号，s为惩罚因子，通过调节e的大小，控制重构原始稀疏信号f的准确度；2b)将无约束目标函数中重构信号^的^个元素用^-《,/ = 1，...,〃代替，得到约束目标函数，令》1 =p;，i = l，...,W，其中nei 2w，初始化向量fi的所有元素都为零；2c)求取约束目标函数的Lagrange对偶函数；2d)由约束目标函数和Lagrange对偶函数得到满足Karush-Kuhn-Tucker ，即KKT条件的最优解；2e)以KKT条件作为重构稀疏信号的终止条件判断，判断当前向量fi中所有元素是否满足KKT条件，若有违反KKT条件的元素，执行步骤2f)，否则跳至步 骤2h);2f)设定一个工作集woA，使得工作集woA中元素个数为miwA，并且选用 当前向量H中最大违反KKT条件的前"woA个元素作为工作集woA中元素；2g)计算工作集MwVt中元素的最优解，更新向量fi中与工作集woA中元素相 对应的元素；2h)计算重构信号的当前值、i,-"^,"l,…,iV;2i)对重构信号进行除偏处理，得到有除偏处理的重构信号X'; (3)计算重构原始稀疏信号f的准确度，用5-(l/A0lx-fl;评价重构信号的无除偏处理准确度，其中符号3表示重构信号的无除偏处理均方误差，用重构信号的有除偏处理均方误差m = (l/iV)||X'-fl;评价重构信号的有除偏处理重构准确度。本专利技术与其他技术相比具有以下优点1、 本专利技术在无约束目标函数的基础上，推导出一种约束目标函数，并且通过 调整参数e使重构信号的准确度和重构速率达到预想结果。2、 本专利技术考虑到约束目标函数的可分解性，将约束目标函数分解为一系列小的线性方程组来解决，节约了存储空间。3、 本专利技术将序列最小优化方法应用到重构原始稀疏信号的约束目标函数上， 并且使用KKT条件作为重构稀疏信号的终止条件判断，具有很好的收敛性和快速 性。4、 本专利技术在重构原始稀疏信号的过程中，通过人工设定当前工作集中元素的 个数，控制重构原始稀疏信号过程的循环次数，提高了重构速率。附图说明图1是本专利技术对稀疏信号重构的过程图2是用本专利技术方法和三种GPSR方法进行目标函数值随时间的仿真曲线图； 图3是用本专利技术方法和一种GPSR方法进行目标函数值和均方误差值随时间 的仿真曲线图4是用本专利技术方法、GPSR方法和两种OMP方法进行均方误差值和时间随 原始稀疏信号中包含非零元素个数的仿真曲线图。具体实施例方式参照图1，本专利技术包括原始数据获取过程、重构原始稀疏信号过程和准确度评价过程。具体实现如下步骤l:利用计算机得到原始稀疏信号f ei 4W、观测矩阵Ae W和观测 向量yeiT4。(la)利用计算机产生包含^个土1元素的原始稀疏信号feW，其他4096-0个元素为零；(lb)利用计算机产生满足均值为0、方差为1和独立同分布的1024 x 4096高 斯矩阵AA ，将矩阵AA按行进行正交化处理，得到观测矩阵A e W ;(lc)根据原始稀疏信号f和观测矩阵A，利用计算机按照等式y-Af + n计算观测向量y e及1024,其中n是方差为cr2 = 10—的高斯白噪声。步骤2:从含有1024个元素的观测向量y中精确重构含有4096个元素的原始稀疏信号f 。(2a)设定重构原始稀疏信号f的无约束目标函数； 1min5 ||y - Ax||2 + f ||xa)式中，xei ，为重构信号，s为惩罚因子，通过调节s的大小，控制重构原始稀疏信号f的准确度；(2b)将无约束目标函数中重构信号x的4096个元素用A - ， / = 1，4096代替，初始化向量fi的所有元素都为零，得到约束目标函数；mmW ^ 2 0,b)H -H -H Hc =—-yrA + e—，P =艾!' = 1，.."8192厂H 一H"l式中，H》(2c)求约束目标函数b)的Lagrange对偶函数;fd =^ H^ + c ji + ^y y-ra n",20，<y,々0，！' = 1，..本文档来自技高网...
【技术保护点】
一种稀疏信号的重构方法，包括如下步骤： （１）利用计算机得到原始稀疏信号ｆ∈Ｒ↑［Ｎ］、观测矩阵Ａ∈Ｒ↑［Ｋ×Ｎ］和观测向量ｙ∈Ｒ↑［Ｋ］； （２）从含有Ｋ个元素的观测向量ｙ中，按如下步骤重构含有Ｎ个元素的原始稀疏信号ｆ： ２ａ）设定重构原始稀疏信号ｆ的无约束目标函数： ｍｉｎ１／２‖ｙ－Ａｘ‖↓［２］↑［２］＋ε‖ｘ‖↓［１］ （１Ａ） 式中，ｘ∈Ｒ↑［Ｎ］为重构信号，ε为惩罚因子，通过调节ε的大小，控制重构原始稀疏信号ｆ的准确度； ２ ｂ）将无约束目标函数中重构信号ｘ的Ｎ个元素用β↓［ｉ］－β↓［ｉ］↑［＊］，ｉ＝１，．．．，Ｎ代替，得到约束目标函数，令μ＝［＊＊＊］，ｉ＝１，．．．，Ｎ，其中μ∈Ｒ↑［２Ｎ］，初始化向量μ的所有元素都为零； ２ｃ）求取约束目标函数的 Ｌａｇｒａｎｇｅ对偶函数； ２ｄ）由约束目标函数和Ｌａｇｒａｎｇｅ对偶函数得到满足Ｋａｒｕｓｈ－Ｋｕｈｎ－Ｔｕｃｋｅｒ，即ＫＫＴ条件的最优解； ２ｅ）以ＫＫＴ条件作为重构稀疏信号的终止条件判断，判断当前向量μ中所有元素是否满足Ｋ ＫＴ条件，若有违反ＫＫＴ条件的元素，执行步骤２ｆ），否则跳至步骤２ｈ）； ２ｆ）设定一个工作集ｗｏｒｋ，使得工作集ｗｏｒｋ中元素个数为ｎｗｏｒｋ，并且选用当前向量μ中最大违反ＫＫＴ条件的前ｎｗｏｒｋ个元素作为工作集ｗｏｒｋ中元素； ２ｇ）计算工作集ｗｏｒｋ中元素的最优解，更新向量μ中与工作集ｗｏｒｋ中元素相对应的元素； ２ｈ）计算重构信号的当前值ｘ↓［ｉ］＝ｕ↓［ｉ］－ｕ↓［ｉ＋Ｎ］，ｉ＝１，．．．，Ｎ； ２ｉ）对重构信号进行除偏处理，得到有除偏处理的 重构信号ｘ↑［＊］； （３）计算重构原始稀疏信号ｆ的准确度，用δ＝（１／Ｎ）‖ｘ－ｆ‖↓［２］↑［２］评价重构信号的无除偏处理准确度，其中符号δ表示重构信号的无除偏处理均方误差，用重构信号的有除偏处理均方误差ｍ＝（１／Ｎ）‖ｘ↑［＊］ －ｆ‖ ↓［２］↑［２］评价重构信号的有除偏处理重构准确度。...
【技术特征摘要】
【专利技术属性】
技术研发人员：张莉，陈桂荣，焦李成，周宏杰，
申请(专利权)人：西安电子科技大学，
类型：发明
国别省市：87[中国|西安]
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。