一种分布式过程监控信源低计算复杂度高可靠编解码算法的制作方法

本发明涉及信号处理领域，具体涉及一种分布式过程监控信源低计算复杂度高可靠编解码算法。

背景技术：

分布式信源从监控需求上可以被分为分布式实时信源和分布式过程信源两类。分布式实时监控信源是指对信息实时性要求高的信源，例如瓦斯、风速、负压等，这类传感器节点需要在短时间内定时输出一个采样值并传输，要求做到实时编解码。分布式过程监控信源是指对信息实时性要求不高的信源，比如煤矿采空区温度、槽波、微震等，这类传感器节点不需要实时传输，可以多次采集后集中传输。本发明研究分布式过程监控信源的编解码算法，由于每个传感器传输的都是一个时间序列，因此分布式过程监控信源不仅具有信号间相关性还具有信号内相关性。分布式压缩感知是一种既能够同时利用信号内相关性又能利用信号间相关性的理论，所以特别适合这种一次通信过程中需要传输多个采样序列值的场合。

文章《distributedcompressivesensing》中drorbaron等人提出了分布式压缩感知理论，是一种结合了压缩感知和分布式信源编码优点的理论，不仅能够利用信号内相关性还能利用信号间相关性。在由任意数量传感器和一个中心节点组成的无线传感器网络中，每个传感器与其它传感器互不通信，独立压缩信号传输至中心节点，中心节点联合恢复每个传感器传过来的信号。分布式压缩感知自提出以来，受到了研究人员的大量关注，文章《distributedcompressivesensing》中drorbaron等人提出了三种联合稀疏表示模型(jsm)。λ:＝{1,2,……,j}表示全体传感器采集信号序号集合，jsm-1模型中，每一个信号由共同稀疏部分和独立稀疏部分构成。xj＝zc+zij,j∈λ，其中zc＝ψθc，zij＝ψθj，||θc||0＝kc，||θj||0＝kj，zc是信号的公共部分，在基ψ上稀疏度为kc，zij为每个信号的独立部分，在基ψ上稀疏度为kj。以三个信号为例，信号表示如公式(1)所示，信号间关系维恩图解如图1所示，信号间相关性结构如图2所示。

jsm-2模型中，公共信息部分为0，每个信号的独立信息部分可以稀疏表示并且有共同的稀疏支撑基，但是非零系数不同，如公式(2)所示。

xj＝ψθj,j∈{1,2,…,j}其中||θj||0＝k,j∈{1,2,…,j}(2)

在jsm-3模型中，公共信息部分在任何基下都不能稀疏表示，每个信号的独立部分可以稀疏表示，如公式(3)所示，其中θc没有非零值。

xj＝zc+zij,zc＝ψθc,zij＝ψθj,其中||θj||0＝kj,j∈{1,2,…,j}(3)

文章《jointrecoveryalgorithmsusingdifferenceofinnovationsfordistributedcompressedsensing》中diegovalsesia等人提出了基于单边信息信号差分的分布式压缩感知算法，在编码端使用随机观测矩阵，随机观测矩阵是在文献《compressedsensing》中提出的，已经被证明满足k阶rip条件。但是随机观测矩阵只在统计意义下以很高概率满足rip和弱相关性，不能保证每次随机观测的信号都能精确恢复原始信号。随机观测矩阵在应用时为乘法运算，计算复杂度高。稀疏二进制观测矩阵变乘法运算为加法运算，计算简单。本发明在编码端使用weizhilu等人在文献《sparsebinarymatricesofldpccodesforcompressedsensing》中提出的基于渐进边增长构造算法(peg)的稀疏二进制观测矩阵。

weizhilu等人提出的基于peg算法的确定性稀疏二进制观测矩阵构造流程如下：观测矩阵可以用tanner图表示，tanner图是一种双向图，由变量节点、校验节点以及这两类节点之间相连的边组成，变量节点对应于校验矩阵的列，校验节点对应于校验矩阵的行，tanner图中，与节点相连的边数目称为节点的度，它与校验矩阵的行重或列重一致。假设观测矩阵的变量节点数为n，校验节点数为m，且节点度分布已经给定，将某变量节点bj的度记为dbj，某校验节点的ci的度记为dci，与变量节点集合vb相连的的边的集合为与bj相连的第k条边记为将该树图中包含的bj深度为l的所有校验节点的集合记为表示，其补集记为

peg算法流程如下：

在解码端，基于单边信息信号差分的分布式压缩感知算法中，边信息采用压缩感知方式采样，通过信号观测值与边信息观测值相减去除了信号与边信息之间的公共部分，获得对差分的观测值，然后通过恢复算法恢复出差分，差分与边信息相加得到待编码信号。由于不需要估计公共信息部分，该算法不存在任何误差。另外，边信息采用全采样的方式，信号就可以利用与边信息之间的公共部分少采样，相当于把工作量转移到了采集边信息的传感器，是一种不对称的分布式压缩感知算法。但是，该算法过度依赖一个边信息会造成解码可靠性和有效性低的问题，在该边信息传感器通信中断时会无边信息可用，在该边信息与待恢复信号相关性很小时会造成解码成功率很低。

技术实现要素：

本发明针对现有技术的不足提出了一种分布式过程监控信源低计算复杂度高可靠编解码算法。

本发明为实现以上目的，采用如下方案：

一种分布式过程监控信源低计算复杂度高可靠编解码算法，其中相关参数如下：边信息集合为s＝{s1,s2,…sq,…sq},q＝1,2,……,q，其中sq∈rⁿ；待编码信号集合为w＝{w1,w2,…wl,…wl},l＝1,2,……,l，其中wl∈rⁿ；信号wl的观测矩阵用φl表示，φl是一个ml×n大小的矩阵，ml＜＜n，φl为稀疏二进制观测矩阵；yl为使用观测矩阵对信号wl观测值，yl＝φlwl；△ylq为信号观测值与边信息观测值的差分值，△ylq＝yl-φlsq；为使用贪婪追踪算法从△ylq恢复得到的信号间差分的估计值；为使用的边信息sq得到的信号wl的估计值，rlq为恢复残差，sc为经过选择后使用的边信息，所述算法流程如下：

s1、编码端，采用稀疏二进制观测矩阵φl对信号进行编码得到yl＝φlwl，并将编码后的结果传输至解码端；

s2、在解码端，计算yl与每个边信息观测后差值，得到观测后差值△ylq＝yl-φlsq；

s3、用贪婪追踪算法从△ylq恢复得到差值估计

s4、依据非零元素个数和恢复残差rlq两个指标对边信息按照优先级排序，首先按照非零元素个数由小到大排序，对于非零元素个数相同的边信息按照恢复残差rlq由小到大排序；

s5、如果最优边信息通信中断，使用次优边信息，以此类推，直至得到使用的边信息sc；

s6、计算信号观测值与边信息sc观测值的差分△ylc＝yl-φlsc；

s7、从信号观测值差分△ylc恢复出信号间差分的估计

s8、计算信号估计值，

本发明和现有技术相比，具有如下优点和有益效果：本发明提出了分布式过程监控信源低计算复杂度高可靠编解码算法，编码端使用稀疏二进制观测矩阵，变乘法运算为加法运算，降低了计算复杂度；在解码端，使用信号差分估计稀疏度和恢复残差两个指标对边信息进行优先级排序，使用最优边信息提高恢复成功率，在最优边信息通信中断时，使用次优边信息，以此类推，保证使用的是能够得到的与待恢复信号相关性最大的边信息，从而提高解码可靠性。

附图说明

图1为分布式信息源相关性关系图一

图2为分布式信息源相关性关系图二；

图3为本发明算法结构图；

图4为本发明编解码流程图；

图5为本发明算法与基于单边信息信号差分分布式压缩感知算法恢复误差比较图；

图6为基于多边信息的煤矿物联网传感节点部署示意图。

具体实施方式

如图3和图4所示，一种分布式过程监控信源低计算复杂度高可靠编解码算法，其中相关参数如下：边信息集合为s＝{s1,s2,…sq,…sq},q＝1,2,……,q，其中sq∈rⁿ；待编码信号集合为w＝{w1,w2,…wl,…wl},l＝1,2,……,l，其中wl∈rⁿ；信号wl的观测矩阵用φl表示，φl是一个ml×n大小的矩阵，ml＜＜n，φl为稀疏二进制观测矩阵；yl为使用观测矩阵对信号wl观测值，yl＝φlwl；△ylq为信号观测值与边信息观测值的差分值，△ylq＝yl-φlsq；为使用贪婪追踪算法从△ylq恢复得到的信号间差分的估计值；为使用的边信息sq得到的信号wl的估计值，rlq为恢复残差，sc为经过选择后使用的边信息，所述算法流程如下：

s1、编码端，采用稀疏二进制观测矩阵φl对信号进行编码得到yl＝φlwl，并将编码后的结果传输至解码端；

s2、在解码端，计算yl与每个边信息观测后差值，得到观测后差值△ylq＝yl-φlsq；

s3、用贪婪追踪算法从△ylq恢复得到差值估计

s4、依据非零元素个数和恢复残差rlq两个指标对边信息按照优先级排序，首先按照非零元素个数由小到大排序，对于非零元素个数相同的边信息按照恢复残差rlq由小到大排序；

s5、如果最优边信息通信中断，使用次优边信息，以此类推，直至得到使用的边信息sc；

s6、计算信号观测值与边信息sc观测值的差分△ylc＝yl-φlsc；

s7、从信号观测值差分△ylc恢复出信号间差分的估计

s8、计算信号估计值，

本发明使用多个边信息，通过信号间差分估计稀疏度和恢复残差两个指标对边信息按照优先级排序，使用最优边信息提高解码准确率，在最优边信息无法获得时使用次优边信息。估计稀疏度和恢复残差两个指标是通过贪婪追踪算法得到的，本发明使用的是基于残差收敛的正交匹配追踪算法，流程如下：

输入：观测矩阵φl、观测后差值△ylq、迭代终止条件ε

输出：信号间差分的估计恢复残差rlq

以下流程中t表示迭代次数，表示空集，λt表示t次迭代的索引(列序号集合)，λt表示第t次找到索引(列序号)，aj表示矩阵φl的第j列，at表示按照索引λt选出的矩阵φl的列集合。

(1)初始化：r0＝△ylq，t＝1；

(2)找到索引λt，使得：

(3)令λt＝λt-1∪λt，at＝at-1∪aλ；

(4)求差分信号估计

(5)更新残差：rlq＝△ylq-at(at^tat)^-1at^t△ylq；

(6)令t＝t+1，如果rlq<e,返回差值信号估计和恢复残差rlq；否则，进入第二步。

本发明建立与方法研究思路是，为了解决基于单边信息信号差分分布式压缩感知算法编码计算复杂度高和解码可靠性低的问题，本发明提出了基于稀疏二进制观测矩阵和多边信息的信号差分分布式压缩感知算法。传统基于差分分布式压缩感知算法编码端使用随机观测矩阵，随机观测矩阵是稠密观测矩阵，编码运算时为乘法运算，计算复杂度高，这对计算能力弱的编码端传感器造成了很大的负担。本发明的解决方案为使用weizhilu等人提出的基于peg算法的稀疏二进制稀疏观测矩阵，将非相关线性测量中的乘法运算变为加法运算，降低了编码复杂度。

基于单边信息信号差分的分布式压缩感知算法解码时过度依赖一个边信息，在该边信息传感器通信中断时会无边信息可用，在该边信息与信号相关性很小时会解码成功率很低，本发明使用多个边信息，使用信号差分估计稀疏度和差分恢复残差两个指标对边信息进行优先级排序，使用最优边信息提高恢复成功率，在最优边信息通信中断时，使用次优边信息，以此类推，保证使用的是能够得到的与待恢复信号相关性最大的边信息，从而提高解码可靠性。

本发明算法适用于部分传感器有线供电的传感网络，有线供电的传感器不受能量限制，能够采用全采样的编码方法，提供更多的信息作为边信息，在有多个传感器采集的信息都可以作为边信息时，本发明算法能够对边信息优先级排序，始终利用最优边信息，从而提高解码可靠性。以煤矿监控传感网络为例，如图6所示，汇聚节点和有线传感器都是有线供电的，均可以采用全采样的编码方式作为边信息，是一种典型的存在多个边信息的分布式压缩感知算法应用场景，此时需要对边信息进行优先级排序，解码时利用与待恢复信号相关性最大的边信息，在最优边信息通信中断时使用次优边信息解码，以此类推，提高解码可靠性。

本节通过仿真验证本发明算法对边信息优先级排序和解码算法的有效性。以三个信号为例，x1和x2作为边信息，x3作为待编码信号。每个信号采用标准高斯分布，支撑集随机选择，信号长度n＝100，测量矩阵采用基于peg算法的稀疏二进制观测矩阵，默认在无噪声条件下进行仿真，信号稀疏度kc、kc{1，2}、kc{2,3}、kc{1,3}、ki1、ki2、ki3随机整数且和为n/2。观测值m从10变化至80，重复试验500次，每次信号和观测矩阵随机产生，对边信息x1、x2进行优先级排序，比较使用本发明算法和基于单边信息分布式压缩感知算法的恢复误差实验效果如图5所示。由图5可以看到，m<30基于多边信息分布式压缩感知恢复算法与使用基于单边信息分布式压缩感知算法差别不大，这主要是由于观测率低造成信号与边信息差分恢复误差大导致边信息优先级排序不准确，随着观测率的增加，40≤m≤80使用基于多边信息分布式压缩感知算法明显优于使用单边信息分布式压缩感知算法，并且随着观测率增加，这种优势越来越明显，m>70时使用基于多边信息分布式压缩感知算法误差趋近于0，使用基于单边信息分布式压缩感知算法要到m>78左右之后误差才趋近于零，说明使用本发明能够选择到与信号相关性最大的边信息，信号恢复误差小于基于单边信息信号差分分布式压缩感知算法。

以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明技术原理的前提下，还可以做出若干改进和变型，这些改进和变型也应视为本发明的保护范围。