一种压缩感知的伪随机等效采样信号重构方法技术

本发明专利技术公开了一种压缩感知的伪随机等效采样信号重构方法，属于信息处理技术领域，利用伪随机等效采样方法对周期信号进行采样；将采样信号在某个变换域内进行稀疏表示或近似稀疏表示；构造观测矩阵；为了重构原始信号需求解式min||α||0s.t.y＝Φx。本发明专利技术利用伪随机等效采样方法对信号进行采样，最终采用OMP算法恢复原始信号，解决了ADC实时采样速率不足问题，降低了控制复杂度，且能够高精度恢复原信号；利用压缩感知理论中的OMP算法恢复原始信号，在不增加任何硬件条件的情况可通过较少采样点实现信号重构，在采样点较少时可精确重构原始信号，同时较基于压缩感知的随机等效采样重构方法具有更高的重构成功率。

A Pseudo-Random Equivalent Sampling Signal Reconstruction Method for Compressed Sensing

【技术实现步骤摘要】
一种压缩感知的伪随机等效采样信号重构方法
本专利技术属于信息处理
，具体涉及一种压缩感知的伪随机等效采样信号重构方法。
技术介绍
在传统的实时采样技术中，模数转换器(AnalogDigitalConverter,ADC)以信号最高频率两倍以上的转换速率进行采样，对于时域上极窄的宽带信号要求ADC具有很高的采样速率，因此极大地增加了硬件设计实施难度。然而等效采样方法通过以远低于奈奎斯特采样率对信号不同周期多次采样得到的采样点进行重构信号波形，从而能够达到较高的等效采样率，因此已成为学术界研究的热点。传统的等效采样方法分为顺序等效采样和随机等效采样。顺序等效采样方法按时间顺序依次获得采样点，每一次触发时增加一个小的增量Δt，直到采样点填满时间窗时按时间顺序恢复重构原始信号，然而该方法仅能在触发点后采样，且当增量Δt较小时难以实现。吴兵等人在论文“一种超宽带等效采样接收机的设计与实现”(雷达科学与技术，2017(4).15(4):443-448)中提出了一种改进的顺序等效采样方法，即在一个触发周期内采集多个数据，并通过在FPGA上内置延迟线控制波形产生系统，从而消除数据对齐问题并能达到较高的等效采样速率。随机等效采样方法通过以触发点为参考基准对信号进行采样，从而能够在若干个周期内采集足够的数据，但该方法难以精确测量出触发点与采样时钟的时间间隔Δt，且采样点分布的不均匀性导致大量的采样点仍无法填充所有的时间窗，因此该方法在采集数据不足时无法有效重构出原始信号。基于并行结构的随机等效采样方法是通过多路ADC对信号进行时间交替采样达到更高的等效采样速率，然而存在各路ADC增益和偏置不一致以及时间非均匀采样问题。针对上述两种方法存在的问题，车俐等人在论文“伪随机等效采样法在超宽带接收机中的应用”(现代雷达，2014,36(7):62-64))中提出了一种伪随机等效采样方法，该方法利用采样周期数和采样点数间的互质关系，使得各采样点均匀复现于同一周期内，因此该方法具有较高的采样重构性能。谢跃雷等人在论文“一种基于FPGA的超宽带雷达数字接收机”(现代雷达，2014，36(1):62-65.))中设计了一种并行的雷达数字接收机，通过FPGA控制驱动四路ADC实现伪随机等效采样，等效采样速率最高可达10Gs/s。实际应用过程中，为了精确重构出原始信号通常需要增加采样时间以获得大量的采样数据，从而导致上述方法采样时间过长，因此限制了其在实际工程中的应用。压缩感知(CompressedSensing,CS)理论是一种全新的采样恢复理论，将传统信号处理中的采样和压缩过程相结合。当信号满足稀疏性时，可以以远低于奈奎斯特采样速率对信号进行采样，并能通过最终的采样值精确重构原始信号，从而能够解决采样数据不足时无法精确重构原始信号的问题。Zhao等人在论文“Compressedsensingenhancedrandomequivalentsampling”(IEEEInternationalConferenceonAcoustics.IEEE,2012.3445-3448)中提出了一种基于压缩感知的随机等效采样重构方法，该方法构造了随机等效采样观测矩阵，通过转换求解l0范数最优化问题重构出原始信号，但由于随机等效采样点分布不均匀性导致了信号重构成功率低。因此，针对信号采样重构方法存在的问题，研究一种重构精度高、实时性强的采样恢复信号的方法是十分有必要的。
技术实现思路
专利技术目的：本专利技术的目的在于提供一种压缩感知的伪随机等效采样信号重构方法，该专利技术通过构造伪随机等效采样观测矩阵，并选择离散傅里叶变换基建立稀疏重构模型，然后利用压缩感知中的正交匹配追踪(OrthogonalMatchingPursuit,OMP)算法求解该模型，从而重构出原始信号。通过验证可知，本专利技术能够有效重构出原始信号且具有更好的重构性能。技术方案：为实现上述目的，本专利技术采用如下技术方案：一种压缩感知的伪随机等效采样信号重构方法，括如下步骤：步骤1)利用伪随机等效采样方法对周期信号进行采样获得采样信号；步骤2)将步骤1)中采样信号在离散傅里叶变换域内进行稀疏表示；步骤3)构造伪随机等效采样观测矩阵，建立稀疏重构模型；步骤4)利用OMP算法转换求解该稀疏重构模型，重构谱域后做逆傅里叶变换得到时域信号。进一步地，所述的步骤1)中，设被采样信号周期为T，根据等效采样处理的信号周期性以第一个采样点的时刻为起始时刻，求得之后的N-1个采样点相对于起始时刻为基准的排列时间点大小，N为传统伪随机等效采样重构信号一个周期所需要的采样点数，从而各采样点按时间顺序由小到大排列。进一步地，步骤2)中，所述的变换域选取快速傅里叶变换域或离散小波变换域。进一步地，所述的步骤2)中，设原始信号x∈RN在变换基Ψ的变换系数为α，其中Ψ＝[ψ1,ψ2,...,ψN]T，R为实数集，[·]T表示转置运算，ψ1,ψ2,...,ψN为变换基Ψ的列向量，则原始信号x可稀疏表示为其中α＝[α1,α2,...,αN]，i为正整数且i∈[1,N]，若α中的非零系数个数为K时，其中K＜＜N，则称信号x在变换基Ψ内是稀疏的，稀疏度为K。进一步地，所述的步骤3)中，所述的原始信号x与观测值y关系表示为y＝Φx，其中Φ为观测矩阵。进一步地，设A＝ΦΨ∈RM×N为压缩感知等效观测矩阵，其中M＜＜N，M为重构具有长度为N的信号波形所需的采样值序列长度，则观测值y∈RM，y的维度远小于x的维度，式y＝Φx＝ΦΨα没有唯一解，原始信号x在表示基Ψ内具有稀疏性，通过l0范数最优化问题由观测值y重构原始信号x，建立稀疏重构模型：min||α||0s.t.y＝Φx，式中，||·||0为向量的l0范数，表示向量α中非零元素个数。进一步地，所述的步骤4)中，重构原始信号需求解式min||α||0s.t.y＝Φx，采用凸优化松弛算法或贪婪算法求解l0范数。而求解l0范数是一个NP-hard(Non-deterministicPolynomialHard)问题，需将此问题进行转换处理。目前所提出的算法有凸优化松弛算法(ConvexOptimizationRelaxtion)和贪婪算法(GreedyAlgorithm)，凸优化松弛算法通过用l1范数替换l0范数转换为求解凸优化问题，然而贪婪算法比凸优化松弛算法工程实现更简单、计算复杂度低，因此本专利技术采用OMP算法进行信号重构，首先在投影前对搜寻到的原子进行正交化过程，然后在原子集上投影获得残差并利用最小二乘法进行更新，最终重构谱域后做逆傅里叶变换得到时域信号。有益效果：传统的实时采样技术中，模数转换器(AnalogDigitalConverter,ADC)以信号最高频率两倍以上的转换速率进行采样，对于时域上极窄的信号要求ADC具有很高的采样速率；本专利技术的一种压缩感知的伪随机等效采样信号重构方法，利用伪随机等效采样方法对信号进行采样，最终采用OMP算法恢复原始信号，解决了ADC实时采样速率不足问题，降低了控制复杂度，且能够高精度恢复原信号；现有的信号重构多是利用采样点重排列方法，本专利技术利用压缩感知理论中的OMP算法恢复原始信号，在不增加任何硬件条件的情况可通过较少采样点实现信号重构，在采本文档来自技高网...

【技术保护点】
1.一种压缩感知的伪随机等效采样信号重构方法，其特征在于，包括如下步骤：步骤1)利用伪随机等效采样方法对周期信号进行采样获得采样信号；步骤2)将步骤1)中采样信号在离散傅里叶变换域内进行稀疏表示；步骤3)构造伪随机等效采样观测矩阵，建立稀疏重构模型；步骤4)利用OMP算法转换求解该稀疏重构模型，重构谱域后做逆傅里叶变换得到时域信号。

【技术特征摘要】
1.一种压缩感知的伪随机等效采样信号重构方法，其特征在于，包括如下步骤：步骤1)利用伪随机等效采样方法对周期信号进行采样获得采样信号；步骤2)将步骤1)中采样信号在离散傅里叶变换域内进行稀疏表示；步骤3)构造伪随机等效采样观测矩阵，建立稀疏重构模型；步骤4)利用OMP算法转换求解该稀疏重构模型，重构谱域后做逆傅里叶变换得到时域信号。2.根据权利要求1所述的一种压缩感知的伪随机等效采样信号重构方法，其特征在于：所述的步骤1)中，设被采样信号周期为T，根据等效采样处理的信号周期性以第一个采样点的时刻为起始时刻，求得之后的N-1个采样点相对于起始时刻为基准的排列时间点大小，N为传统伪随机等效采样重构信号一个周期所需要的采样点数，从而各采样点按时间顺序由小到大排列。3.根据权利要求1所述的一种压缩感知的伪随机等效采样信号重构方法，其特征在于：步骤2)中，所述的变换域选取快速傅里叶变换域或离散小波变换域。4.根据权利要求3所述的一种压缩感知的伪随机等效采样信号重构方法，其特征在于：所述的步骤2)中，设原始信号x∈RN在变换基Ψ的变换系数为α，其中Ψ＝[ψ1,ψ2,...,ψN]T，R为实数集，[·]T表示转置运算，ψ1,...

【专利技术属性】
技术研发人员：李家强，董石明，陈焱博，陈金立，
申请(专利权)人：南京信息工程大学，
类型：发明
国别省市：江苏,32