本发明涉及一种基于压缩感知及∑-δ量化的信息平滑与重构方法,属于模拟信息转换以及压缩感知技术领域。
背景技术:
在对高频信号、大数据量、宽带扩频信号的采集时,传统的数据采集技术面临着巨大的挑战。信息技术的飞速发展导致数据量极具增加,经典的奈奎斯特采样定理大大限制了信息数据的采集与存储。
模拟到信息转换技术类似于模拟数字转换,将对模拟信号样点进行压缩并经跟踪算法恢复。压缩感知是缓解数据采集压力的有效方法之一,尤其当信号稀疏或可压缩的情况。基于压缩感知的模拟信息转换器则直接对模拟信号以远小于奈奎斯特采样频率进行采样,最后通过相关的重构算法对原信号进行精确重构。现有模拟信号的数字转换技术中很重要的过程是如何实现压缩量化,即用几个比特表达采集的信号,就引发一个问题,就是到底能压缩到多少?压缩感知是否能用于获取适当信号的压缩表达?由此带来信息压缩比率的不确定性以及压缩表达的不适当性,导致原始信息重构困难。
技术实现要素:
本发明的目的是针对现有模拟信息转换技术中原始信号准确恢复困难以及鲁棒性差的技术缺陷,提出了一种基于压缩感知及∑-δ量化的信息平滑与重构方法。
一种基于压缩感知及∑-δ量化的信息平滑与重构方法,包括如下步骤:
步骤一、压缩采样;
其中,压缩采样具体通过采样操作因子φ实现,具体为:
y=φx;
其中,x是采样点,该采样点x的维度是n;y是测量矩阵,该测量矩阵的维度是m;
n>m,且n至少是m的1.5倍;y中的元素是x的压缩采样点;采样操作因子φ的维度为m行n列;
采样点x经过与采样操作因子φ进行矩阵相乘后,x中的全部信息都会完整的得到保留并传递给测量矩阵y;
步骤二、对步骤一输出的测量矩阵进行去趋势化处理,输出去趋势化测量矩阵y1;
步骤三、求步骤二输出的去趋势化测量矩阵y1的标准差;
步骤四、以步骤三输出的标准差的k倍,以0为中心上下划分s+1个等分区间段,统计矩阵y1中的数据标准差位于s段中各段内的数量并计算概率;
其中,k的取值范围为2到20,s的取值范围为20到100;
步骤五、计算步骤四输出的s+1段中每一段的renyi值,并找出最大renyi值对应的去趋势化测量矩阵y1中的数据,具体包括如下子步骤:
步骤5.1初始化renyi值re为0;
步骤5.2基于公式(1)计算测量矩阵中对应到s段中各段内数据概率的renyi值;
其中,re为renyi值;p(i)对应s+1个区间段的概率,q为renyi熵指数;q的取值范围为0.3到3;
步骤六、将步骤五输出的去趋势化测量矩阵y1中的数据替换为前后x段内对应数据的均值,更新去趋势化测量矩阵y1;
其中,x的取值范围为1到3;
至此,从步骤二到步骤六,完成了压缩采样所得测量矩阵y的去干扰处理,输出结果yy;
步骤七、将步骤六输出的更新后去趋势化测量矩阵y1通过量化操作因子ω将步骤六得到的更新后的去趋势化测量矩阵y1量化为测量值q;
q=ω{yy};
其中,q∈c0;c0是一个有限集合,ω{yy}表示对测量矩阵yy基于量化操作因子进行量化,且量化后的测量值q用有限比特流c0=am表达;
其中,量化基于无记忆标量量化以及∑-δ调制技术:m为1时,即1比特的量化集合am={-1,1};
步骤八、编码,具体将步骤七量化输出的测量值q用特定比特的数字表达,输出编码结果;
特定比特的数字通过log2|c0|表示,具体通过编码图ε:c0-->c实现;其中,c是有限集合,称为码本;
编码通过离散的johnon-lindenstrauss嵌入实现;编码的目标是保证准确重构的基础上降低比特数,步骤八用如下公式(2)通过编码矩阵乘积表示:
qc=bq;(2)
其中,b为c0-->c转换的编码矩阵,b的维度为l行m列;qc为编码结果;
步骤九、重构,基于步骤八输出的编码结果恢复出步骤一中的原始采样信号x经去干扰处理后的近似值,具体基于映射δ:c→rn;
其中,r为重构比特率,表达为:r=log2|c|。
有益效果
一种基于压缩感知及∑-δ化的信息平滑与重构方法,与现有技术相比,具有如下有益效果:
1.所述的信息平滑与重构方法基于压缩采样技术,提供了对稀疏及压缩信号的近似优化编码;
2.所述的信息平滑与重构方法采用renyi熵可以有效去除原始信号中的小概率干扰;
3.所述的信息平滑与重构方法是kolmogorov意义上的可实现的编码与重构,能够实现在去除干扰基础上,高准确率的还原原始信号。
附图说明
图1是本发明一种基于压缩感知及∑-δ量化的信息平滑与重构方法的结构组成;
图2是本发明一种基于压缩感知及∑-δ量化的信息平滑与重构方法中的压缩采样结构示意图;
图3是本发明一种基于压缩感知及∑-δ量化的信息平滑与重构方法的流程图;
图4是本发明一种基于压缩感知及∑-δ量化的信息平滑与重构方法中的信息重构流程图;
图5是基于本发明一种基于压缩感知及∑-δ量化的信息平滑与重构方法对含噪声信号进行采集及重构的仿真结果图。
具体实施方式
下面结合附图和实施例对本发明一种基于压缩感知及∑-δ量化的信息平滑与重构方法做进一步说明和详细描述。
实施例1
本实施例叙述了本发明一种基于压缩感知及∑-δ量化的信息平滑与重构方法的系统结构组成及具体实施。
如图1所示,为一种基于压缩感知及∑-δ量化的信息平滑与重构方法的所依托系统的结构框图。根据该结构框图,制作了本实施例。
从图1可以看出,一种基于压缩感知及∑-δ量化的信息平滑与重构方法的所依托系统,简称模拟信息转换系统,包含采样模块、平滑模块、量化模块、编码模块和重构模块;
模拟信息转换系统中各模块的连接关系如下:
采样模块平滑模块相连,平滑模块与量化模块相连,量化模块与编码模块相连,编码模块和重构模块相连。
图2是本发明一种基于压缩感知及∑-δ量化的信息平滑与重构方法中的压缩采样结构示意图,具体为:针对可压缩信号x1,首先要进行稀疏变换,变换为信号x;再对x进行压缩变换,从而完成压缩采样;
具体实施时,假设被采集的信号是稀疏的,即具备稀疏度为kk的特性;那么被采集信号可以直接经压缩变化被采集,具体通过公式y=φx实现,其中,φ为采样操作因子;
其中,x是采样点,该被采集信号x的维度是n;y是测量矩阵,该测量矩阵的维度是m;
n>m,且n至少是m的1.5倍;y中的元素是x的压缩采样点;采样操作因子φ的维度为m行n列;具体实施时,n=1600;k=8;m=100;采样操作因子φ选取m行n列的高斯测量矩阵,矩阵中的每一个元素是均值为0方差为1的变量;
采样点x经过与采样操作因子φ进行矩阵相乘后,x中的全部信息都会完整的得到保留并传递给测量矩阵y。
图3是本发明一种基于压缩感知及∑-δ量化的信息平滑与重构方法的流程图;
从图3可以看出,本发明所述信息平滑与重构方法,包括如下步骤:
步骤1、基于稀疏变换及采集操作因子进行压缩采集,得到测量矩阵;
步骤2、对测量矩阵进行去趋势化处理;
其中,具体到本实施例,去趋势化处理采用去均值函数。
步骤3、求去趋势化处理后测量矩阵的标准差;
其中,标准差计算采用matlab软件中的std函数;
步骤4、求去趋势化处理后测量矩阵的标准差以标准差k倍、0为中心上下划分s+1个等分区间段,统计各段内的数量并计算概率;
具体到本实施例的实施,k取3,s取50;
步骤5、计算每一段的renyi值,并找出最大renyi值对应的去趋势化测量矩阵y1中的数据,具体为:初始化renyi值re为0,通过如下公式计算测量矩阵中对应到s段中各段内数据概率的renyi值;
其中,re为renyi值;p(i)对应s+1个区间段的概率,q为renyi熵指数;当q为2时,为信息熵;具体到本实施例,q值取0.7;
步骤6、将去趋势化测量矩阵y1中的数据替换为前后t段内对应数据的均值;
具体实施时,t取2;
至此,从步骤2到步骤6,完成了压缩采样所得测量矩阵y的去干扰处理,输出yy;
步骤7、通过量化操作因子ω将更新后的去趋势化测量矩阵y1量化为测量值q;
其中,量化操作因子,记为ω,量化通过公式q=ω{yy},实现,其中,“yy”为步骤6的输出;
其中,q∈c0;c0是一个有限集合,ω{yy}表示对测量矩阵yy基于量化操作因子进行量化,且量化后的测量值q用有限比特流c0=am表达;
其中,量化基于无记忆标量量化以及∑-δ调制技术:m为1时,即1比特的量化集合am={-1,1};
步骤8、测量值q用特定比特的数字表达,输出编码结果;
其中,特定比特的数字通过log2|c0|表示,具体通过编码图ε:c0-->c实现;其中,c是有限集合,称为码本;
编码通过离散的johnon-lindenstrauss嵌入实现;编码的目标是保证准确重构的基础上降低比特数,步骤三用如下公式通过编码矩阵乘积表示:
qc=bq;
其中,b为c0-->c转换的编码矩阵,b的维度为l行m列;qc为编码结果;
步骤9、重构,编码结果恢复出原始采样信号经去干扰处理后的近似值。
具体将步骤8输出的编码结果恢复出yy的近似值,具体基于映射δ:c→rn;
其中,r为重构比特率,表达为:r=log2|c|。
具体到本实施例,重构基于映射恢复原始采样信号,采用正交匹配追踪算法,其拥有较高的重构精度并且运算复杂度低,正交匹配追踪重构算法步骤如图4所示。步骤8输出编码结果,即序列y[m];重构矩阵t=采样操作因子φ×傅里叶正交反变换矩阵ψ;算法迭代次数gamma次,gamma应不小于信号的稀疏度k;残差值ri,其初值r0=y[m],残差值随着每一次迭代而更新;增量矩阵ωi,其初值为空矩阵;重构矩阵的列向量和残差的投影系数向量βi,其初值β0为零向量;最大投影系数对应的位置向量λi,以及由其构成的位置集合λi,其初值为空集;重构的频域系数向量s,含有n个元素,其对应位置的第λi个元素与βi的第i个元素相等,而其余的n-gamma个元素由零向量构成;最后输出重构信号。
图5是基于本发明一种基于压缩感知及∑-δ量化的信息平滑与重构方法对含噪声信号进行采集及重构的仿真结果图;图中,连续曲线为含噪声的原始信息,带点的曲线为经过本发明所述方法进行信息平滑以及重构后的信号,可以看出,重构信号对原始含噪信号进行了平滑,并且能较为精确的重构原始信号。
以上所述为本发明的较佳实施例,本发明不应该局限于该实施例和附图所公开的内容。凡是不脱离本发明所公开的精神下完成的等效或修改,都落入本发明保护的范围。