本发明专利技术属于信号检测与估计(signal detection and estimation)领域，尤其涉及一种在被检测信号具有稀疏性的前提下，对字典和稀疏信号进行联合估计的超分辨率算法。本发明专利技术提供一种只需少量的信号采样即可在有噪声情况下恢复与估计信号的方法，并提供高精度的保证。本发明专利技术通过使用压缩感知技术，只使用少量的信号便可对字典和稀疏信号进行精确的估计。通过联合估计，可以消除由于网格不匹配产生的误差，从而达到高精度估计。
【技术实现步骤摘要】
本专利技术属于信号检测与估计（signaldetectionandestimation)领域，尤其涉 及一种在被检测信号具有稀疏性的前提下，对字典和稀疏信号进行联合估计的超分辨率算 法。
技术介绍
传统的信号处理中，要求对信号的采样率至少为信号最高频率的两倍，也即奈奎 斯特采样率。在信号频率较低的情况下，对采样频率的这个要求是很容易达到的。然而， 现在很多应用的实现中，信号的频率很高，这种情况下，要取得奈奎斯特采样率是十分困难 的。压缩感知是一种可以突破奈奎斯特采样率的技术，其以远低于奈奎斯特采样率的采样 率对信号进行采样，并精确的恢复信号。为了达到这个目的，压缩感知技术要求原信号具 有稀疏性，不过实际应用中绝大多数信号都具有稀疏性，所以近年来压缩感知技术得到了 很广泛的应用。同样，在很多实际应用场景中，对于能耗有着严格的要求，比如传感器网络， 一般要求尽量以低功耗来处理信号，以延长传感器网络的使用寿命。以上均说明了，利用低 米样率来进行?目号检测和估计的必要性，同时，在低米样率的情况下，希望获得尚精度的估 计。目前已经有很多基于压缩感知的算法着重于解决对信号的估计问题。例如，交匹 配追踪算法、基追踪算法、稀疏贝叶斯算法等。所有这些算法，在求解过程中，都预定义了一 个离散化的字典，并假设估计值在格点上。但大多数时候，估计值是偏离格点的，这就造成 了不可避免的误差，使得高精度估计十分困难。本专利技术中，对信号和字典进行联合估计，在 整个可能的区间对信号估计值进行搜索，最终得到精确的估计值。
技术实现思路
本专利技术的目的在于提供一种只需少量的信号采样即可在有噪声情况下恢复与估 计信号的方法，并提供高精度的保证。本专利技术通过使用压缩感知技术，只使用少量的信号便 可对字典和稀疏信号进行精确的估计。通过联合估计，可以消除由于网格不匹配产生的误 差，从而达到高精度估计。 为了方便描述，首先对术语进行定义： 感知矩阵：用以对信号进行线性采样，起到降维的作用，把η维信号映射到m维空 间，通常m<<η。 稀疏性：信号可以用一组基或者一个字典中若干元素线性表示。当这种表示是精 确的，就称这个信号是稀疏的。大多数的高维信号所包含的信息远远低于它的维度所能包 含的，稀疏信号模型为这种高维信号提供了数学上的解释。 稀疏表示：信号可以用一组基中若干元素线性表示，就称这组基为稀疏基。稀疏的 信号在稀疏基下的表不即为信号的稀疏表不。 替代目标函数：原目标函数的一个替代函数，对这个替代函数进行优化可以达到 对原目标函数优化目的。 ，包括如下步骤： S1、构造具有随机采样性质的感知矩阵ψ，对信号进行采样得到yn，其中，η为自然 数，η= 1，2,…，Ν;S2、设置初始参数e、#'t=0,设置初始解其中，t为采样时 间，，为初始格点，e为迭代终止条件，是经验值，εw为初始的噪声方差和为初始的 拉格朗日系数；S3、迭代，具体为：S31、根据构建替代目标函娄，其 中，，Θ为格点，y为已知的初始信号， Α(θ)表示只取括号内Θ向量元素对应的列；S32、固定格点Θ，对S31所述替代目标函数求导，求得最优解表达式，将所述f(Θ)代入S31所述替代目标函数，替代 \ 人''J 目标函数转化为；):/，利用梯度下降法搜索 一个新的估计值分叫，所述分叫满足不等式/ X#1pS33、更新信号的稀疏表示，使;S34、如I则更新ε(t+1)=ε(t)/10 ; S35、如果ε (t)<〇·〇〇〇1，则更新S36、如果迭代过程满足终止条件<e，则停止迭代，否则进行下一次迭 代；S4、最终返回信号的稀疏表示 4灰复的信号#可以通过稀疏表示i和对应的频率 备求得。 进一步地，S2 所述e= 10 6。本专利技术的有益效果是： 利用本专利技术方法，在高带宽或者对能耗要求很严格的情况下，只需要少量的信号 采样样本，即可完成信号的参数估计和恢复，同时可以达到超分辨率的精度。相比于传统的 超分辨率信号估计算法，本专利技术具有更强的实际操作性。【附图说明】 图1为使用本专利技术方法应用于谱估计问题的流程图。 图2为测量次数分别与RSNR、成功率的关系，其中，（a)为测量次数与RSNR的关 系，（b)测量次数为与成功率的关系。 图3为稀疏度分别与RSNR、成功率的关系，其中，（a)为测量次数与RSNR的关系， (b)测量次数为与成功率的关系。 图4为PSNR分别与RSNR、成功率的关系，其中，（a)为测量次数与RSNR的关系， (b)测量次数为与成功率的关系。 图5为频点间距分别与RSNR、成功率的关系，其中，（a)为测量次数与RSNR的关 系，（b)测量次数为与成功率的关系。【具体实施方式】 下面结合实施例和附图，详细说明本专利技术的技术方案。 本专利技术实施可用于信号的谱估计问题，为描述方便，下述实施例将建立以下系统 模型来进行说明。 设谱线估计问题中，测量信号3ν··>ν]是K个复正弦信号的叠加，即,1 = 1，2…Τ，其中，I表示独立同分布的零均值高斯噪声。 仿真中，需要预设一个初始格点Θw，统一设置为Θw= (2JI/NK0 1…Ν-1] τ，并取N= 64,终止条件e= 1〇6。其中频率{c〇J是在[0,2 30均匀产生的，复幅度值 {ak}是在复单位圆上产生的。测量值y通过从yT中随机选取Μ个值得到。并假设T= 64,Κ= 3,Μ= 20,即需要在一个对长度64的信号进行采样得到一个长度为20的采样信 号，并由这个采样信号精确估计出3个频率{cok}和原信号^。噪声强度用信噪峰值比 (peak-signal-to-noise ration，PSNR)来表不，即PSNR = 101og10(l/ σ 2)。 基于上述所构建模型及定义，本专利技术提供了一种基于压缩感知的超分辨率算法进 行参数估计和信号恢复。具体为： S1、构造具有随机采样性质的感知矩阵Ψ，对信号进行采样得到yn，其中，η为自然 数，η= 1，2,…，Ν; S2、设置初始参数e、ε 〇、，>、，>、t= 0,设置初始解^，其中，t为采样时 间，少U|为初始格点，e为终止条件，是经验值； S3、迭代，具体为： S31、根据p)、构建替代目标函数，其 中，Θ为格点；S32、固定格点Θ，对S31所述替代目标函数求导，求得最优解表达式「，将所述f(Θ)代入S31所述替代目标函数，替代 目标函数转化为，利用梯度下降法搜索 一个新的估计值妒+1)，所述沪+"满足不等式;S33、更新信号的稀疏表示，使;S34、如』则更新ε(t+1)=ε(t)/l〇 ;S35、如果ε(t)< 〇· 〇〇〇1，则更fS36、如果迭代过程满足终止条件|#+1)-妒;| <e，则停止迭代，否则进行下一次迭 代；S4、最终返回信号的稀疏表示：，恢复的信号j>可以通过稀疏表示I和对应的频率 &求得。 经过上述操作，就完成了对信号的频率估计和恢复。 下面将利用了离散化格点压缩感知的相关算法同本专利技术方法的算法性能对比分 析，以进一步验证本专利技术的性能。 采用两种衡量指标来度量算法的性能。一个是用来衡量信号的恢复性能，叫做 重构信噪比（reconstructionsignal-to-noiser本文档来自技高网...
【技术保护点】
基于压缩感知的超分辨率信号估计方法，其特征在于，包括如下步骤：S1、构造具有随机采样性质的感知矩阵Ψ，对信号进行采样得到yn，其中，n为自然数，n＝1,2,…,N；S2、设置初始参数∈、ε(0)、t＝0，设置初始解其中，t为采样时间，为初始格点，∈为迭代终止条件，是经验值，ε(0)为初始的噪声方差和为初始的拉格朗日系数；S3、迭代，具体为：S31、根据构建替代目标函数P(x,θ|x^(t))=xHD(t)x+λ^(t)||y-A(θ)x||22,]]>其中，D(t)=Δdiag{1|x^1(t)|2+ϵ(t),1|x^2(t)|2+ϵ(t),...,1|x^N(t)|2+ϵ(t)},]]>θ为格点，y为已知的初始信号，A(θ)表示只取括号内θ向量元素对应的列；S32、固定格点θ，对S31所述替代目标函数求导，求得最优解表达式x*(θ)=(AH(θ)A(θ)+1λ(t)D(t))-1AH(θ)y,]]>将所述x*(θ)代入S31所述替代目标函数，替代目标函数转化为f(θ)=Δ-yHA(θ)(AH(θ)A(θ)+1λ(t)D(t))-1AH(θ)y,]]>利用梯度下降法搜索一个新的估计值所述满足不等式S33、更新信号的稀疏表示，使S34、如果则更新ε(t+1)＝ε(t)/10；S35、如果ε(t)＜0.0001，则更新S36、如果迭代过程满足终止条件则停止迭代，否则进行下一次迭代；S4、最终返回信号的稀疏表示ω^←{θ^i(t+1)||x^i(t+1)|>10-3,i=1,2,...N},]]>恢复的信号可以通过稀疏表示和对应的频率求得。...
【技术特征摘要】
【专利技术属性】
技术研发人员：方俊，崔星星，李靖，
申请(专利权)人：电子科技大学，
类型：发明
国别省市：四川;51
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。