CN105637824B - 从接收信号中恢复稀疏通信信号的方法 - Google Patents
从接收信号中恢复稀疏通信信号的方法 Download PDFInfo
- Publication number
- CN105637824B CN105637824B CN201380078991.9A CN201380078991A CN105637824B CN 105637824 B CN105637824 B CN 105637824B CN 201380078991 A CN201380078991 A CN 201380078991A CN 105637824 B CN105637824 B CN 105637824B
- Authority
- CN
- China
- Prior art keywords
- signal
- communication
- indicate
- collection
- coefficient
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L25/00—Baseband systems
- H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
- H04L25/03—Shaping networks in transmitter or receiver, e.g. adaptive shaping networks
- H04L25/03987—Equalisation for sparse channels
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B17/00—Monitoring; Testing
- H04B17/20—Monitoring; Testing of receivers
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Power Engineering (AREA)
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Mobile Radio Communication Systems (AREA)
- Complex Calculations (AREA)
Abstract
本发明提供一种从接收信号中恢复稀疏通信信号的方法(100),所述接收信号为所述稀疏通信信号的信道输出版本,所述信道包括用于形成信道矩阵的信道系数,所述方法包括从所述信道矩阵和所述接收信号,确定(101)支持集,所述支持集指示一组非零通信信号系数的第一指数;基于所述支持集、所述信道矩阵以及所述接收信号,确定(103)所述稀疏通信信号的估计;确定(105)所述支持集未指示的通信信号系数的第二指数;基于所述支持集、所述稀疏通信信号的估计、所述第二指数以及所述信道矩阵,确定(107)所述稀疏通信信号。
Description
技术领域
本发明涉及稀疏信号恢复领域,尤其涉及性能和质量方面的问题。
背景技术
许多不同的应用中都出现了稀疏信号的恢复问题。在通信系统中从接收信号中恢复稀疏通信信号已引起了人们的广泛关注。
一个n维通信信号或者向量,若其全部n个元素中仅包括少量几个非零元素,则称为稀疏。与总维数n相比,可以通过较少数量的测量来恢复稀疏通信信号。稀疏性可出现在不同的基,所以通信信号可例如在时域或者频域中表现为稀疏通信信号。
稀疏通信信号的恢复的计算量通常比最小二乘法中涉及到的矩阵求逆或者伪逆的计算量大。此外,由于系统欠定,矩阵通常是不可逆的。
普通的稀疏通信信号恢复方法存在高计算复杂度和/或低性能的问题。这些方法通常不适合在实时系统中使用。
E.J.Candes,和Candes在“线性程序译码”(IEEE信息传输理论,2005年12月,第51卷,第12期,4203-4215页)中提出并分析了一种线性程序方法。
T.Cai,L.Wang在“用于带有噪音的稀疏信号恢复的正交匹配追踪”(IEEE信息传输理论,第57卷,第7期,4680-4688页)中使用一种正交匹配追踪方法讨论了稀疏信号恢复中的不完整测量问题。
发明内容
本发明的目的是提供一种从接收信号中恢复稀疏通信信号的有效方案。
该目的通过独立权利要求的特征来实现。进一步的实施形式通过从属权利要求、说明书和附图显而易见。
本发明是基于以下发现:为了确定稀疏通信信号,可以对稀疏通信信号的估计进行细化。
按照本发明的第一方面,本发明提供一种从接收信号中恢复稀疏通信信号的方法,所述接收信号为所述稀疏通信信号的信道输出版本,所述信道包括用于形成信道矩阵的信道系数,所述方法包括从所述信道矩阵和所述接收信号,确定支持集,所述支持集指示一组非零通信信号系数的第一指数;基于所述支持集、所述信道矩阵以及所述接收信号,确定所述稀疏通信信号的估计;确定所述支持集未指示的通信信号系数的第二指数;基于所述支持集、所述稀疏通信信号的估计、所述第二指数以及所述信道矩阵,确定所述稀疏通信信号。因此,稀疏通信信号可以被有效的恢复。
稀疏通信信号可以表现为一个向量。稀疏通信信号可以包括通信信号系数。该通信信号的系数可以是实数,如1.5或2.3,或复数,例如1+j或5+3j。
接收信号可以表现为一个向量。接收信号可以包括接收信号系数。接收信号系数可以是实数,如1.3或2.7,或复数,例如2+3j或1-5j。
信道可以定义稀疏通信信号和接收信号之间的线性关系。该通道可包括附加的噪声。信道系数可以是实数,如0.3或1.5,或复数,例如3-2j或1+4j。
支持集可指示一组非零通信信号系数的指数。
根据第一方面的第一种实施方式中,该方法进一步包括基于所述通信信号系数的第二指数,更新所述支持集。因此,该支持集可以被细化以进行进一步处理。
支持集的更新可以涉及到为支持集增加通信信号系数的第二指数,或从支持集中减去通信信号系数的第二指数。
根据第一方面或第一方面的第一种实施方式,在第二种实施方式中,该方法还包括:确定所述支持集或所述更新的支持集的基数,所述基数用于指示所述支持集或所述更新的支持集的元素数量。因此,可以得到稀疏通信信号的稀疏度度量。
根据第一方面或第一方面的任一前述实施方式,在第三种实施方式中,所述稀疏通信信号包括通信信号系数,所述接收信号包括接收信号系数,其中,所述通信信号系数的数量大于所述接收信号系数的数量。因此,数量减少的接收信号系数可以用于执行该方法。
根据第一方面或第一方面的任一前述实施方式,在第四种实施方式中,所述确定所述支持集或者所述确定所述稀疏通信信号的估计使用正交匹配追踪法。因此,用于确定所述支持集和/或确定所述稀疏通信信号的估计的有效方法可以被应用。
正交匹配追踪方法可以适用于优化受线性等式约束和/或线性不等式约束的线性目标函数。正交匹配追踪方法例如由T.Cai,L.Wang在“用于带有噪音的稀疏信号恢复的正交匹配追踪”(IEEE信息传输理论,第57卷,第7期,4680-4688页)中被描述。
根据第一方面或第一方面的任一前述实施方式,在第五种实施方式中,所述确定通信信号系数的第二指数包括最小化对偶线性规划问题的残差,所述对偶线性规划问题通过以下公式定义:
-e≤HTz≤e,
yTz→min,
其中H表示所述信道矩阵,y表示所述接收信号,e表示一个全单元向量,z表示一个对偶变量。因此,可以以一种优化的方式确定通信信号系数的第二指数。
根据第一方面或第一方面的任一前述实施方式,在第六种实施方式中,所述确定通信信号系数的第二指数包括最小化对偶线性规划问题的残差,所述残差定义为:
或者
其中
若
若
其中S*表示根据向量的元素的符号被划分为两个子集的所述支持集,表示所述支持集S*的对应于的正元素的一个子集,表示所述支持集S*的对应于的负元素的一个子集,表示一个辅助向量,表示索引为i的辅助向量的系数,表示包括支持集S*所指示的列的信道矩阵,表示对偶线性规划问题的一个估计解,Hj表示包括索引j所指示的列的信道矩阵,或表示所述残差。因此,可以以一种优化的方式确定通信信号系数的第二指数。
根据第一方面或第一方面的任一前述实施方式,在第七种实施方式中,所述确定通信信号系数的第二指数包括最小化对偶线性规划问题的残差,所述确定通信信号系数的第二指数按照以下公式执行:
或
其中或者表示对偶线性规划问题的残差,j(1)表示所述通信信号系数的第二指数。因此,可以以一种优化的方式确定通信信号系数的第二指数。
根据第一方面或第一方面的任一前述实施方式,在第八种实施方式中,所述确定所述稀疏通信信号按照以下公式执行:
其中表示包括支持集S*所指示的列的信道矩阵,表示包括通信信号系数的第二指数j(1)所指示的列的信道矩阵,α表示一个预先确定的值,表示所述稀疏通信信号的估计,表示所述稀疏通信信号。因此,稀疏通信信号可以被有效地确定。
预先确定的值α可以是一个实数,如0.2,或复数,如1+j。
根据第一方面的第一种实施方式,在第九种实施方式中,所述更新所述支持集按照以下公式执行:
其中j(1)表示通信信号系数的第二指数,表示待从所述支持集中移除的通信信号系数的指数,i(1)表示待从所述支持集中移除的元素的数量,S*表示所述支持集,以及S(1)表示所述更新后的支持集。因此,所述支持集可以以一种优化的方式被细化。
根据第一方面或第一方面的任一前述实施方式,在第十种实施方式中,依次重复执行以下步骤:确定所述支持集,确定所述稀疏通信信号的估计,确定所述通信信号系数的第二指数,确定所述稀疏通信信号,更新所述支持集以及确定所述基数,直至达到停止标准。因此,该方法的性能可以被改善。
根据第一方面的第十种实施方式,在第十一种实施方式中,对应于具有最小基数的支持集的所述稀疏通信信号作为所述恢复后的稀疏通信信号被提供。因此,恢复的稀疏通信信号的稀疏度可以被优化。
根据第一方面的第十种实施方式或第一方面的第十一种实施方式,在第十二种实施方式中,所述停止标准为对偶线性规划问题的所有残差的值为正数或零,其中残差被定义为:
或者
其中
若
若
其中S*表示根据向量的元素的符号被划分为两个子集的所述支持集,表示所述支持集S*的对应于的正元素的一个子集,表示所述支持集S*的对应于的负元素的一个子集,表示一个辅助向量,表示索引为i的辅助向量的系数,表示包括支持集S*所指示的列的信道矩阵,表示对偶线性规划问题的一个估计解,Hj表示包括索引j所指示的列的信道矩阵,或表示残差。因此,该停止标准可以基于性能标准。
根据第一方面的第十种实施方式到第一方面的第十二种实施方式,在第十三种实施方式中,所述停止标准为到达一个预定的时间间隔长度。因此,该停止标准使得该方法可以应用于实时系统中。
预定的时间间隔长度可例如为1ms、5ms或10ms。
根据第二方面,本发明提供一种计算机程序,所述计算机程序在计算机上运行时用于执行第一方面或第一方面的任一前述实施方式所述的方法。
该计算机程序可以以一个机器可读代码的形式被提供。该计算机程序可以包括一系列用于计算机的处理器的指令。
该计算机可以包括处理器、存储器、输入接口和/或输出接口。计算机的处理器可用于执行计算机程序。
本发明可以通过硬件和/或软件来实现。
附图说明
下文中结合以下附图对本发明的进一步实施例进行描述,其中:
图1示出了一种从接收信号中恢复稀疏通信信号的方法的示意图;以及
图2示出了一种用于从接收信号中恢复稀疏通信信号的时隙分配的示意图。
具体实施方式
图1示出了一种从接收信号中恢复稀疏通信信号的方法100的示意图。
所述接收信号为所述稀疏通信信号的信道输出版本,所述信道包括用于形成信道矩阵的信道系数。
该方法100包括:101,从所述信道矩阵和所述接收信号,确定支持集,所述支持集指示一组非零通信信号系数的第一指数;103,基于所述支持集、所述信道矩阵以及所述接收信号,确定所述稀疏通信信号的估计;105,确定所述支持集未指示的通信信号系数的第二指数以及基于所述支持集、所述稀疏通信信号的估计、所述第二指数以及所述信道矩阵;以及107,确定所述稀疏通信信号。
稀疏通信信号可以表现为一个向量。稀疏通信信号可以包括通信信号系数。该通信信号的系数可以是实数,如1.5或2.3,或复数,例如1+j或5+3j。
接收信号可以表现为一个向量。接收信号可以包括接收信号系数。接收信号系数可以是实数,如1.3或2.7,或复数,例如2+3j或1-5j。
信道可以定义稀疏通信信号和接收信号之间的线性关系。该信道可包括附加的噪声。信道系数可以是实数,如0.3或1.5,或复数,例如3-2j或1+4j。
支持集可指示一组非零通信信号系数的指数。
图2示出了一种用于从接收信号中恢复稀疏通信信号的时隙分配的示意图。
使用OMP法及其变型能够解决从不完备测量恢复稀疏信号的问题。此外,也有线性规划方法LP。然而,没有将OMP和LP结合到一个数值化方案的尝试。
若n维向量0的全部个数n中仅包括少量几个非零元素,则称为稀疏。与总维数n相比,可以通过较少数量的测量来恢复稀疏可压缩向量。稀疏性可出现在不同的基,所以信号可在时域或者频域中表现为稀疏向量。
稀疏通信信号的重构的计算量通常比最小二乘法中涉及到的矩阵求逆或者伪逆的计算量大。此外,由于系统欠定,矩阵通常是不可逆的。
数值有效重构算法适合实时应用。线性系统
Hx=y, (1)
其中,若m<n,则H∈Rm×n,x∈Rn,y∈Rm欠定。更确切地说,若H的秩小于n(rank(H)<n),则该系统(1)欠定。这种系统有无穷多解。该系统的图案可以被示意性地表示为如下描述的矩阵向量图案:
用于减少测量数量的折衷是向量x的稀疏性。稀疏意味着零主导向量x的元素。在系统(1)的解中搜索最稀疏的解的方法是压缩感测(compressive sensing,CS)理论要考虑的关键问题。
(CS)理论要考虑的关键问题可以转化为重建对应于最稀疏向量的正确的支持集,即有效非零元素的集合:
Supp(x)={i:xi≠0}
可以看出当
是确定最稀疏解的问题的正确公式,其中||x||0=Supp(x),该问题的松弛形式为
也是充分的,其中是该向量的lp-标准,p>0。这个问题的非凸性使得它很难为数值求解。一个常用的方法用于确定最稀疏解,即具有最大数量的零元素的解,包括于线性规划问题(linear program,LP)的问题(P1)的解中:
其中,||x||1为向量的l1-标准,它是压缩感测(compressivesensing,CS)理论中使用的一种标准方法。
相比于问题(P0)的可能最优、但NP困难的解,其优势是包括诸如单纯形法的线性规划方法在内的所有已知凸优化方法的凸性和适应性,以及良好的性能。
一种实时实施更为有效的解决方案是正交匹配追踪(orthogonal matchingpursuit,OMP)算法。它的核心思想用公式表示如下。
步骤0:初始设定残差r(0)=y,初始设定支持集的估计
步骤k=1,2,…。找出解决问题的矩阵H中的列
并且更新S(k)=S(k-1)∪{j(k)},其中,为包括该些具有指数j∈S(k)的Hj的H的子矩阵。
重复执行该些步骤直到达到停止标准(stopping criterion,SC)。
这是一种所谓的贪婪算法,已纳入该支撑集中的指数不再从该集中取出。与其它方法相比,该OMP算法的优点是它的高计算复杂性。如果与通过(P1)得到的结果进行比较,该OMP算法在重建稀疏信号时显示出较低的性能结果。
所提到的问题可以通过将OMP方法和线性规划方法(P1)组合成一个数值算法来解决。
图2示出了该问题的解决方案的时隙分配的结构。假设分配给稀疏恢复问题的时隙为T∑,OMP迭代占用如图所示的小于T∑的TOMP。因此,通过引入OMP解决方案作为初始步骤,我们有时间T∑-TOMP进行改进。
在时隙T∑用尽后立即终止该过程,与常规的OMP相比,在TOMP和T∑之间获得了尽可能多的改进。
考虑到由l1最小化导致的LP问题,其中l1最小化可以等效地被表示为
其中,e为包括所有单元的向量,有合适的大小。此处x=u-v。该线性规划问题将被称为初始LP(primary LP,PLP)。与(PLP)一起考虑对偶线性规划问题(dual LP,DLP):
其中,z表示对偶变量。
假设完成OMP,结果为:
S*-通过OMP得到的该支持集的最后估计,l*≡card(S*)
集合S*可以进一步根据向量的元素的符号被划分为两个子集:
假设定义:
若
若
通常l*≤m,所以该矩阵不为正方形。
与基S*相对应的DLP的解的当前估计,可以按照欠定系统的最小范数解计算
其具有闭合解
这个向量对应于问题(DLP),其对于线性不等式组可能是可行的也可能是不可行的。其对于(DLP)可以是最优或非最优的。由于在第一阶段(OMP)的最后一个步骤中获得集合S*,保持最优性条件。所以,由于(SC),在OMP的最后一次迭代中获得的向量x*的可行性条件是“几乎满意”的。但(PLP)与(DLP)对偶,并且(PLP)的可行性条件是(DLP)最优性条件。
因此,仅检查向量的可行性条件。检查不等式的(DLP)的残差,以查看是否存在违反条件或类似具有不正确符号的残差的情况。
(DLP)的残差为
和对于所有j都可为非负。如果满足该条件,解和分别能够解决(PLP)和(DLP)的问题。令该些残差其中之一为负。不失一般性地,令
并且
则可以生成一个新的解
其通过采用如下l1改进性能生成:
i(1)为待从基(序列1,2,…,l*中的序列号)中移除的元素的数量。
增强的OMP算法可以总结如下:
1)k=0;
当(分配的时隙未用尽)
若n维向量x的全部个数n中仅包括少量几个非零元素,则为稀疏。
与总维数n相比,可以通过较少数量的测量来恢复稀疏可压缩向量。
稀疏性可出现在不同的基,所以信号可在时域或者频域中表现为稀疏向量。
稀疏信号的重构的计算量通常比仅测量矩阵求逆的计算量大。此外,由于系统欠定,矩阵通常是不可逆的。
理想的情况是,数值有效重构算法适合实时应用。
可以执行从欠定或不完备测量模型中的稀疏信号恢复。线性系统
Hx=y
其中,H∈Rm×n,y∈Rm,x∈Rn,m≤n欠定。具有图案
其具有无穷多解:Νull(H)+x0是与每一个单一解x0一起的一组解。为找出最稀疏的解,即,零元素数量最多的解,可以作为l1最小值被找出:
||x||1→min受限制于Hx=y
其为一个线性规划(linear program,LP)。它是压缩感测(compressive sensing,CS)理论中使用到的标准问题设置。
LP的计算量较大。正交匹配追踪(orthogonal matching pursuit,OMP)算法更易于实时实施。它的核心思想用公式表示如下。
步骤0:初始设定残差r0=y以及支持集的估计
步骤k=1,….:找出解决问题的矩阵H中的列
更新Sk=Sk-1∪{jk},rk=y-Hxk,其中 表示对应于该组指数的矩阵H的子矩阵。
重复执行该些步骤直到达到停止标准。
假设该测量模型通过公式y=Hx+ε,x∈Rn,y∈Rm描述,其中测量的数量m远小于二进制字x的长度。
假设向量x中不多于k个元素为非零。换言之n,
xi=0 forxi≠0 for i∈supp(x),card(supp(x))=k.
此外,矩阵H满足有限等距性(restricted isometry property,RIP)性能并且保持以下条件
OMP的停止标准可被选择为
该核心思想是结合OMP和LP方法的优势,将OMP作为加强的OMP(enhanced OMP,EOMP)的第一阶段,在(SC)满足之后,对改进单纯形法进行多次迭代,将Sk作为初始基,尝试按顺序修改其内容。
该结合算法是贪婪的:该基仅在第一OMP阶段扩展;在第二阶段中从该基移除的指数不再次插入。在执行第二阶段时,使用矩阵分解的低秩变型来修改伪逆矩阵。
考虑到由l1最小化导致的LP问题,其可以等效地被表示为
Hu-Hv=y,u≥0,v≥0,
eTu+eTv→min (PLP)
其中,e为包括所有单元的向量,有合适的大小。此处x=u-v。该LP将被称为初始LP(primary LP,PLP)。与(PLP)一起考虑对偶LP(dual LP,DLP):
-e≤HTz≤e,
yTz→min, (DLP)
其中,z表示对偶变量。
假设完成OMP,结果为:
S*-通过OMP得到的该支持集的最后估计,l*≡card(S*)
集合S*可以进一步根据向量的元素的符号被划分为两个子集:
假设定义:
若
若
与基S*相对应的DLP的解的当前估计,可以按照欠定系统的最小范数解计算
其具有闭合解
这个向量对应于它中的问题(DLP),其对于线性不等式组可能是可行的也可能是不可行的。其对于(DLP)可以是最优或非最优的。
由于在第一阶段(OMP)的最后一个步骤中获得集合S*,保持最优性条件。所以,由于(SC),可行性条件是“几乎满意”的。但(PLP)与(DLP)对偶,并且(PLP)的可行性条件是(DLP)的最优条件。因此,仅能检查可行性条件。
令为与基S*对应的矩阵H的列指数集。检查不等式的(DLP)的残差。
(DLP)的残差为
和都可为非负。如果是这样,解和分别能够解决(PLP)和(DLP)的问题。令该些残差其中之一为负。不失一般性地,令可以生成一个新的解
其通过采用如下l1改进性能生成:
i(1)为待从基中移除的元素的数量。
最终,该算法可以总结如下:
当(DLP约束违反并且迭代次数不超过限制)
在一个实施方式中,本发明涉及一种用于稀疏信号恢复的加强的正交匹配追踪法。
在一个实施方式中,本发明涉及稀疏信号检测,尤其涉及其性能质量方面的问题。
在一个实施方式中,本发明涉及结合OMP与一个或多个l1改进步骤的稀疏信号恢复方法,该方法包括:执行OMP,根据方程式(2)、(3)计算问题(DLP)的解,检查残差(4)、(5),做出改进(6)、(7),保持该记录解。
在一个实施方式中,对基数来说,该记录解是最好的。
在一个实施方式中,对l1标准值来说,该记录解是最好的。
在一个实施方式中,终止条件是(DLP)(5)的可行。
在一个实施方式中,终止条件是为该问题分配的时间的耗尽。
在一个实施方式中,与OMP和/或LP相比,本发明的性能更好。
在一个实施方式中,可在各步骤使用低秩矩阵因式分解以减少计算负荷。
在一个实施方式中,本发明涉及一种结合OMP和LP方法的优势的稀疏信号恢复方法。
Claims (12)
1.一种从接收信号中恢复稀疏通信信号的方法(100),其特征在于,所述接收信号为所述稀疏通信信号的信道输出版本,所述信道包括用于形成信道矩阵的信道系数,所述方法包括:
从所述信道矩阵和所述接收信号,确定(101)支持集,所述支持集指示一组非零通信信号系数的第一指数;
基于所述支持集、所述信道矩阵以及所述接收信号,确定(103)所述稀疏通信信号的估计;
确定(105)所述支持集未指示的通信信号系数的第二指数;
基于所述支持集、所述稀疏通信信号的估计、所述第二指数以及所述信道矩阵,确定(107)所述稀疏通信信号;
所述确定(105)通信信号系数的第二指数包括最小化对偶线性规划问题的残差,
所述对偶线性规划问题通过以下公式定义:
-e≤HTz≤e
z=arg min(yTz)
其中H表示所述信道矩阵,y表示所述接收信号,e表示一个全单元向量,z表示一个对偶变量;
所述残差定义为:
或者
其中
若
若
其中S*表示根据向量的元素的符号被划分为两个子集的所述支持集,表示所述支持集S*的对应于的正元素的一个子集,表示所述支持集S*的对应于的负元素的一个子集,表示一个辅助向量,表示索引为i的辅助向量的系数,表示包括支持集S*所指示的列的信道矩阵,表示对偶线性规划问题的一个估计解,Hj表示包括索引j所指示的列的信道矩阵,或表示所述残差。
2.根据权利要求1所述的方法(100),其特征在于,还包括:
基于所述通信信号系数的第二指数,更新所述支持集。
3.根据权利要求1或2所述的方法(100),其特征在于,还包括:
确定所述支持集或更新的支持集的基数,所述基数用于指示所述支持集或所述更新的支持集的元素数量。
4.根据权利要求1或2所述的方法(100),其特征在于,所述稀疏通信信号包括通信信号系数,所述接收信号包括接收信号系数,其中,所述通信信号系数的数量大于所述接收信号系数的数量。
5.根据权利要求1或2所述的方法(100),其特征在于,所述确定(101)所述支持集或者所述确定(103)所述稀疏通信信号的估计使用正交匹配追踪法。
6.根据权利要求1或2所述的方法(100),其特征在于,所述确定(105)通信信号系数的第二指数包括最小化对偶线性规划问题的残差,所述确定通信信号系数的第二指数按照以下公式执行:
或
其中或者表示对偶线性规划问题的残差,j(1)表示所述通信信号系数的第二指数。
7.根据权利要求1或2所述的方法(100),其特征在于,所述确定(107)所述稀疏通信信号按照以下公式执行:
其中表示包括支持集S*所指示的列的信道矩阵,表示包括通信信号系数的第二指数j(1)所指示的列的信道矩阵,α表示一个预先确定的值,表示所述稀疏通信信号的估计,表示所述稀疏通信信号。
8.根据权利要求2所述的方法(100),其特征在于,所述更新所述支持集按照以下公式执行:
其中j(1)表示通信信号系数的第二指数,表示待从所述支持集中移除的通信信号系数的指数,i(1)表示待从所述支持集中移除的元素的数量,S*表示所述支持集,以及S(1)表示所述更新后的支持集。
9.根据权利要求1或2所述的方法(100),其特征在于,依次重复执行以下步骤:所述确定(101)所述支持集,所述确定(103)所述稀疏通信信号的估计,所述确定(105)所述通信信号系数的第二指数,所述确定(107)所述稀疏通信信号,更新所述支持集以及确定基数,直至达到停止标准。
10.根据权利要求9所述的方法(100),其特征在于,对应于具有最小基数的支持集的所述稀疏通信信号作为所述恢复后的稀疏通信信号被提供。
11.根据权利要求9所述的方法(100),其特征在于,所述停止标准为对偶线性规划问题的所有残差的值为正数或零。
12.根据权利要求9所述的方法(100),其特征在于,所述停止标准为到达一个预定的时间间隔长度。
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/RU2013/000974 WO2015065230A1 (en) | 2013-11-01 | 2013-11-01 | Method for recovering a sparse communication signal from a receive signal |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105637824A CN105637824A (zh) | 2016-06-01 |
CN105637824B true CN105637824B (zh) | 2018-11-09 |
Family
ID=50842308
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201380078991.9A Active CN105637824B (zh) | 2013-11-01 | 2013-11-01 | 从接收信号中恢复稀疏通信信号的方法 |
Country Status (3)
Country | Link |
---|---|
US (1) | US9882750B2 (zh) |
CN (1) | CN105637824B (zh) |
WO (1) | WO2015065230A1 (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101995882B1 (ko) * | 2017-10-11 | 2019-07-04 | 서울대학교산학협력단 | 시분할 이중통신 시스템에서 상향링크 전송 방법 및 장치 |
CN114696946B (zh) * | 2020-12-28 | 2023-07-14 | 郑州大学 | 数据编码及解码方法、装置、电子设备及存储介质 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101795501A (zh) * | 2009-10-30 | 2010-08-04 | 中国科学院上海微系统与信息技术研究所 | 传感器网络物理信号协同压缩感知系统及感知方法 |
CN101895297A (zh) * | 2010-07-30 | 2010-11-24 | 哈尔滨工业大学 | 一种面向压缩感知的块稀疏信号重构方法 |
CN101908889A (zh) * | 2010-07-30 | 2010-12-08 | 哈尔滨工业大学 | 一种块稀疏度未知的稀疏信号压缩感知重构方法 |
CN102244624A (zh) * | 2011-07-12 | 2011-11-16 | 深圳数字电视国家工程实验室股份有限公司 | 基于正交匹配追踪的稀疏信道估计方法 |
CN102833020A (zh) * | 2012-09-10 | 2012-12-19 | 杭州电子科技大学 | 认知无线网络中基于自适应测量的贝叶斯压缩宽带频谱检测方法 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8855190B2 (en) * | 2009-10-21 | 2014-10-07 | Futurewei Technologies, Inc. | Communication system with compressive sensing |
US8648742B2 (en) * | 2011-02-25 | 2014-02-11 | California Institute Of Technology | Systems and methods for acquiring and decoding signals using compressed sensing |
US20120259590A1 (en) * | 2011-04-11 | 2012-10-11 | Jong Chul Ye | Method and apparatus for compressed sensing with joint sparsity |
KR101209908B1 (ko) * | 2011-08-04 | 2012-12-11 | 광주과학기술원 | 희소 신호 전송 방법 및 장치, 그리고 희소 신호 복구 방법 및 장치 |
US9178590B2 (en) * | 2011-12-27 | 2015-11-03 | Industrial Technology Research Institute | Channel information feedback method and wireless communication device using the same |
US8929390B2 (en) * | 2012-09-24 | 2015-01-06 | Alcatel Lucent | Methods and apparatuses for channel estimation in wireless networks |
-
2013
- 2013-11-01 WO PCT/RU2013/000974 patent/WO2015065230A1/en active Application Filing
- 2013-11-01 CN CN201380078991.9A patent/CN105637824B/zh active Active
-
2016
- 2016-04-29 US US15/142,731 patent/US9882750B2/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101795501A (zh) * | 2009-10-30 | 2010-08-04 | 中国科学院上海微系统与信息技术研究所 | 传感器网络物理信号协同压缩感知系统及感知方法 |
CN101895297A (zh) * | 2010-07-30 | 2010-11-24 | 哈尔滨工业大学 | 一种面向压缩感知的块稀疏信号重构方法 |
CN101908889A (zh) * | 2010-07-30 | 2010-12-08 | 哈尔滨工业大学 | 一种块稀疏度未知的稀疏信号压缩感知重构方法 |
CN102244624A (zh) * | 2011-07-12 | 2011-11-16 | 深圳数字电视国家工程实验室股份有限公司 | 基于正交匹配追踪的稀疏信道估计方法 |
CN102833020A (zh) * | 2012-09-10 | 2012-12-19 | 杭州电子科技大学 | 认知无线网络中基于自适应测量的贝叶斯压缩宽带频谱检测方法 |
Non-Patent Citations (2)
Title |
---|
Ordered Orthogonal Matching Pursuit;Deepak Baby, et al;《Communications (NCC), 2012 National Conference on》;20120403;第1-5页 * |
Subspace Pursuit Embeded in Orthogonal Matching Pursuit;Sooraj K. Ambat, ET AL;《TENCON 2012- 2012 IEEE Regin 10 Conference》;20121122;第1-5页 * |
Also Published As
Publication number | Publication date |
---|---|
US9882750B2 (en) | 2018-01-30 |
WO2015065230A1 (en) | 2015-05-07 |
US20160248611A1 (en) | 2016-08-25 |
CN105637824A (zh) | 2016-06-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Servadei et al. | A BREZIS-NIRENBERG RESULT FOR NON-LOCAL CRITICAL EQUATIONS IN LOW DIMENSION. | |
Needell et al. | CoSaMP: iterative signal recovery from incomplete and inaccurate samples | |
Li et al. | Iterative row sampling | |
CN107832406B (zh) | 海量日志数据的去重入库方法、装置、设备及存储介质 | |
Brustein et al. | Ratio of shear viscosity to entropy density in generalized theories of gravity | |
US9209833B1 (en) | Methods and apparatus for rational compression and decompression of numbers | |
WO2005119507A1 (ja) | 行列の高速高精度特異値分解方法、プログラムおよび装置 | |
CN106503268B (zh) | 数据对比方法、装置和系统 | |
Faucheaux et al. | Higher-order diagrammatic vibrational coupled-cluster theory | |
US8571255B2 (en) | Scalable media fingerprint extraction | |
CN105637824B (zh) | 从接收信号中恢复稀疏通信信号的方法 | |
Wang et al. | A Runge–Kutta type scheme for nonlinear stochastic partial differential equations with multiplicative trace class noise | |
Batenkov et al. | On the global-local dichotomy in sparsity modeling | |
Das et al. | Random convolutional coding for robust and straggler resilient distributed matrix computation | |
Huai et al. | Zerobn: Learning compact neural networks for latency-critical edge systems | |
CN109241498B (zh) | Xml文件处理方法、设备和存储介质 | |
Altschuler et al. | Shifted composition I: Harnack and reverse transport inequalities | |
CN110708160B (zh) | 基于sm2算法标量乘法编码的抗侧信道攻击方法及系统 | |
CN115982310B (zh) | 一种自带验证功能的链表生成方法及电子设备 | |
Xu | The performance of orthogonal multi-matching pursuit under RIP | |
Papapantoleon et al. | Stability results for martingale representations: the general case | |
US10474666B2 (en) | Metadata validation tool | |
CN105095513B (zh) | 一种破解Oracle数据库账户密码的方法和装置 | |
Wen et al. | Robust ensemble feature extraction for solving conditional nonlinear optimal perturbation | |
Eigel et al. | Local equilibration error estimators for guaranteed error control in adaptive stochastic higher-order Galerkin FEM |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |