CN105181122A - 机械振动信号数据压缩采集方法 - Google Patents
机械振动信号数据压缩采集方法 Download PDFInfo
- Publication number
- CN105181122A CN105181122A CN201510477122.4A CN201510477122A CN105181122A CN 105181122 A CN105181122 A CN 105181122A CN 201510477122 A CN201510477122 A CN 201510477122A CN 105181122 A CN105181122 A CN 105181122A
- Authority
- CN
- China
- Prior art keywords
- matrix
- vibration signal
- calculation matrix
- signal
- mechanical oscillation
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 42
- 238000013144 data compression Methods 0.000 title claims abstract description 18
- 239000011159 matrix material Substances 0.000 claims abstract description 184
- 238000007906 compression Methods 0.000 claims abstract description 46
- 230000006835 compression Effects 0.000 claims abstract description 46
- 238000005070 sampling Methods 0.000 claims abstract description 16
- 238000005457 optimization Methods 0.000 claims abstract description 13
- 238000005259 measurement Methods 0.000 claims abstract description 11
- 238000004364 calculation method Methods 0.000 claims description 118
- 230000010358 mechanical oscillation Effects 0.000 claims description 59
- 238000009826 distribution Methods 0.000 claims description 25
- 238000000354 decomposition reaction Methods 0.000 claims description 18
- 230000008447 perception Effects 0.000 claims description 16
- 238000010276 construction Methods 0.000 claims description 12
- 230000009466 transformation Effects 0.000 claims description 12
- 230000008569 process Effects 0.000 claims description 10
- 238000012545 processing Methods 0.000 claims description 7
- 230000005540 biological transmission Effects 0.000 claims description 6
- 230000008859 change Effects 0.000 claims description 4
- 238000003860 storage Methods 0.000 claims description 4
- 238000012938 design process Methods 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 6
- 238000013461 design Methods 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 3
- 238000012544 monitoring process Methods 0.000 description 3
- 238000011161 development Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000006855 networking Effects 0.000 description 2
- 206010068052 Mosaicism Diseases 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000006880 cross-coupling reaction Methods 0.000 description 1
- 238000013501 data transformation Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000009776 industrial production Methods 0.000 description 1
- 230000036244 malformation Effects 0.000 description 1
- 238000011089 mechanical engineering Methods 0.000 description 1
- 210000003765 sex chromosome Anatomy 0.000 description 1
- 230000035939 shock Effects 0.000 description 1
- 230000008054 signal transmission Effects 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Landscapes
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明公开了一种机械振动信号数据压缩采集方法,属于一种基于压缩采样理论的机械振动信号数据采集方法的技术领域。通过测量矩阵使得机械振动信号投影为低维的压缩采集值,并在振动信号重构端求解稀疏最优化问题,高概率地用低维的压缩采集值重构出高维机械振动信号。本发明很好地克服了依据传统采样定理(香农-奈奎斯特)得到海量机械振动数据的问题,测量矩阵对模数转换后的机械振动信号进行压缩采样,传输、存储和处理压缩采样值,这样大大减少了原始机械振动数据量,具有既可以得到较高的信号压缩比又有着精确的信号重构性能的优点。
Description
技术领域
本发明涉及一种机械振动信号数据压缩采集方法,属于振动信号数据采集和压缩方法的领域,具体说属于一种基于压缩采样理论的机械振动信号数据采集方法的技术领域。
背景技术
机械振动信号传递与承载着机械设备工作过程中所蕴含的大量重要信息,在线监测与采集机械振动信号是机械工程领域,尤其是故障诊断技术中的关键技术之一。传统的且目前被广泛采用的振动信号检测与采样用的各种传感器一直是以香农-奈奎斯特(Shannon-Nyquist)采样定理为基础的,如图1所示,经典的Shannon-Nyquist振动信号采集理论,要求采集频率高于振动信号最高频率的两倍,对采集后的振动信号采用线性加权能实现原始振动信号的精确重构。随着振动信号频带的加宽,基于经典的信号采集理论会得到巨量的振动数据,为了传输、存储和后期处理的方便,将采集的振动信号数据进行传统的数据变换压缩,这种采集后再压缩信号的方法效率较低,经典的振动信号采集理论是基于信号带宽的,在获取与采样机械振动信号时,为了不丢失振动信号中携带的机械设备相关信息,采样频率必须大于信号中最高频率的两倍,才能由采样信号精确重构原始振动信号。但是,随着实际工业生产要求的不断提高与机械工业的不断发展,机械装备愈趋大型化、成套化、高速化、网络化、集成化等;由于机械系统在工作过程中会产生撞击、速度突变、结构变形、摩擦变化等,且不同组成部件间相互交叉耦合,振动更加复杂,随机性振动频率越来越高且呈现非线性、非平稳性。而且随着大型机械振动监测向综合、高速、连续和网络化趋势的发展,基于香农-奈奎斯特的振动信号采样导致了巨量的采样数据,这些数据的实时传输与同步存储已成为亟待解决的成本与工程技术瓶颈问题,尤其在Internet远程设备状态监测领域,大量的采集数据对传输造成了巨大的压力。
现有的采集方法存在两方面的缺陷:其一,先采样再压缩的方式严重浪费了采样资源;其二,变换压缩后的数据在传输、存储过程中如果发生丢失,就会造成原始振动信号的不可恢复性问题。
发明内容
本发明针对现有技术的不足提供了一种机械振动信号数据压缩采集方法,以实现既可以得到较高的信号压缩比又有着精确的信号重构性能,在不丢失振动信息的情况下,可大大减少原始振动数据量的目的。
为达到所述的目的本发明的技术方案或方法是:
一种机械振动信号数据压缩采集方法,其具体步骤如下:
(1)振动信号的稀疏分解;假设机械振动信号为f∈RN,其中N为振动信号维数,采用傅里叶正交基Ψ对振动信号进行稀疏分解,得到振动信号的稀疏系数α=ΨTf,其中ΨT表示Ψ的转置矩阵;
(2)采集得到全局非自适应压缩测量值;确定一个测量矩阵Φ∈RM×N(M<<N)对振动信号的有用信息进行高效获取,其中M为压缩采集信号的维数;
(3)传输、存储和处理压缩测量值y;根据测量矩阵对机械振动信号在傅里叶正交字典下的稀疏系数进行测量,得到压缩测量值y=ΦΨα,并将压缩测量值存储在采样存储器中;
(4)基于振动信号的稀疏性,采用最优化算法,重构出原始振动信号;
在机械振动信号重构处理端,求解1-范数凸优化问题,
min||α||s.t.y=ΦΨα=Φf
其中Ψ为傅里叶基,具体算法采用正交匹配追踪算法,实现原始机械振动信号的重构;
形成通过测量矩阵使得高维机械振动信号投影为低维的压缩采集值,并在信号重构端求解最优化问题,高概率地用低维压缩采集值重构出高维的原始机械振动信号的结构。
所述的步骤(1)采用傅里叶正交基Ψ对振动信号进行稀疏分解为选取零范数最最小的字典作为稀疏分解基。
所述的步骤(1)经过稀疏变换后得到振动信号在傅里叶正交基上的稀疏系数,同时根据振动信号的特征,估计振动信号在稀疏域上的稀疏度k,稀疏度k是振动信号中所含频率成分的数量。
根据振动信号在傅里叶基上的稀疏度k,确定测量矩阵的行数M。
所述的稀疏度k确定后,测量矩阵的行数M必须满足Spark理论不等式M≥2k,以确定保证信息不丢失的最少测量次数;其中k为振动信号在傅里叶基上的稀疏度。
所述的步骤(2)中,测量矩阵Φ为第一测量矩阵和第二测量矩阵。
该第一测量矩阵为机械振动信号的最优测量矩阵,该最优测量矩阵的构造算法为:根据最优测量理论,其设计过程是:已知傅里叶正交基Ψ、初始测量矩阵Φ0、阈值t和尺度下降因子γ,得到使得μt{D=ΦΨ}最小的矩阵Φ∈RM×N作为测量矩阵,称此时的Φ为傅里叶稀疏变换基对应的最优测量矩阵。
该最优测量矩阵具体构造过程如下:
基于最优测量矩阵构造算法,得到振动信号的最优测量矩阵如下:
输入参数:傅里叶正交稀疏变换基ΨDFT∈CN×N,测量矩阵Φ∈RM×N的初始值Φ0为一个高斯随机矩阵或托普利兹矩阵,令压缩感知矩阵D=ΦΨ,t表示阈值,γ表示尺度下降因子(0<γ<1),M表示测量次数,Iter表示迭代次数,循环变量k,初始值k=0;
具体步骤:
(1)根据振动信号在傅里叶基上的稀疏性,确定测量矩阵Φ的行数,M≥2k,k表示信号稀疏度;
(2)计算感知矩阵D=ΦΨ;
(3)对感知矩阵D进行列单位化,得到矩阵
(4)计算感知矩阵的Gram矩阵,其中表示的转置矩阵;
(5)根据阈值t更新Gk,按如下关系得到矩阵
(6)减小的秩,对角化对角矩阵Λk的主对角线上的元素,保留绝对值前M大的元素,其他都设定为0,改变Λk后,重新计算
(7)求解Dk,Dk是M×N矩阵,通过感知矩阵得到经过一次迭代的振动信号优化测量矩阵;
(8)更新Φ,找到Φk+1,使得误差最小,即
(9)k=k+1;
直到k=Iter,循环结束。
输出:ΦDFT对应的最优测量矩阵ΦDFT。
该第二测量矩阵为高斯分布循环测量矩阵:
假设构造的高斯分布循环测量矩阵为Φ,Φ∈RM×N(其中M为矩阵的行数,N为列数,且M<<N)
高斯分布循环测量矩阵Φ具体形式为:
高斯分布循环测量矩阵构造的具体步骤为:
(1)生成N个服从高斯分布的随机数作为测量矩阵的第一行,即:
(2)将φ1通过循环移位生成φ2,即:
(3)在向量中,随机取出n个元素(1≤n≤N),记为[a1a2…an];
(4)其中αi~N(0,1),将放在[a1a2…an]在φ2中对应位置上,得到
(5)根据以上第(2)、(3)、(4)步的向量生成规律,重复第(2)、(3)、(4)步,得到M-1个行向量;
(6)将以上得到的M个行向量按生成的顺序组成矩阵。
高斯分布循环测量矩阵的构造过程中,在情况下,当n=N时,构造得到的高斯分布循环测量矩阵本质上是高斯随机测量矩阵;当n=1时,构造得到的高斯分布循环测量矩阵本质上是循环测量矩阵。
所述的步骤(4)机械振动信号的重构中的机械振动信号是可压缩信号或近似稀疏信号,根据凸优化理论,对振动信号的重构就是对振动信号L0-范数的求解,在信号具有稀疏性的情况下,L0-范数和L1-范数等价,所以求解更简洁的L1-范数会得到同样的解;在实际应用中,采用正交匹配追踪OMP完成对机械振动信号的重构。
采用本发明的方法由于针对目前机械振动信号频带越来越宽,依据传统香农-内奎斯特采样定理进行数据采集时,将会得到巨量振动数据,对存储、传输和处理带来困难的问题,提出了一种振动信号的数据压缩采集方法,通过测量矩阵使得机械振动信号投影为低维的压缩采集值,并在振动信号重构端求解稀疏最优化问题,高概率地用低维的压缩采集值重构出高维机械振动信号。通过试验分析了振动信号在傅里叶基上的稀疏性,融入振动信号在傅里叶基上稀疏性的结构信息到其测量矩阵的设计中,得到适合振动信号的测量矩阵并进行压缩采集,在不丢失振动信号信息的情况下,大大减少了原始振动数据量。具有压缩测量值可以重构原始振动信号,进行传统的处理分析,也可以直接从压缩测量值中提取信息的效果。
附图说明
图1为经典的机械振动信号压缩采集原理示意图;
图2为本发明基于稀疏性的机械振动信号压缩采集原理示意图;
图3为本发明振动信号的最优测量矩阵构造流程图;
图4为本发明高斯分布循环测量矩阵构造流程图。
具体实施方式
下面结合附图对本发明的技术方案或方法详细说明如下,以利全面的了解。
如图2,图3和图4所示;压缩采集机械振动信号方法,其特殊之处在于:通过测量矩阵使得机械振动信号投影为低维的压缩测量信号,并在信号重构端求解稀疏最优化问题,高概率地用低维的压缩测量信号重构出原始机械振动信号,基本原理如图2所示。
一种机械振动信号数据压缩采集方法,其具体步骤如下:
(1)振动信号的稀疏分解;假设机械振动信号为f∈RN(R在数学里表示全体实数;这里的意思是f属于R),其中N为振动信号维数,采用傅里叶正交基Ψ对振动信号进行稀疏分解,得到振动信号的稀疏系数α=ΨTf,其中ΨT表示Ψ的转置矩阵;
(2)采集得到全局非自适应压缩测量值;确定一个测量矩阵Φ∈RM×N(M<<N)对振动信号的有用信息进行高效获取,其中M为压缩采集信号的维数;
(3)传输、存储和处理压缩测量值y;根据测量矩阵对机械振动信号在傅里叶正交字典下的稀疏系数进行测量,得到压缩测量值y=ΦΨα,并将压缩测量值存储在采样存储器中;
(4)基于振动信号的稀疏性,采用最优化算法,重构出原始振动信号;
在机械振动信号重构处理端,求解1-范数凸优化问题,
min||α||s.t.y=ΦΨα=Φf
其中Ψ为傅里叶基,具体算法采用正交匹配追踪算法,实现原始机械振动信号的重构;
形成通过测量矩阵使得高维机械振动信号投影为低维的压缩采集值,并在信号重构端求解最优化问题,高概率地用低维压缩采集值重构出高维的原始机械振动信号的结构。
所述的步骤(1)采用傅里叶正交基Ψ对振动信号进行稀疏分解为选取零范数最小的字典作为稀疏分解基。
振动信号的稀疏分解可理解为首先对振动信号进行稀疏表示,其二利用完备字典(傅里叶基)和过完备字典进行稀疏分解,其三比较各种字典下的稀疏分解效果,选取零范数最小的字典作为稀疏分解基,最后经过分析,选取傅里叶基作为振动信号的稀疏分解基。
傅里叶基为:
其中W是1的n次方根的主值,大小为exp(-2πi/N),为正则化因子。
所述的步骤(1)经过稀疏变换后得到振动信号在傅里叶正交基上的稀疏系数,同时根据振动信号的特征,估计振动信号在稀疏域上的稀疏度k,稀疏度k是振动信号中所含频率成分的数量。
根据振动信号在傅里叶基上的稀疏度k,确定测量矩阵的行数M。
所述的稀疏度k确定后,测量矩阵的行数M必须满足Spark理论不等式M≥2k以确定保证信息不丢失的最少测量次数;其中k为振动信号在傅里叶基上的稀疏度。
所述的步骤(2)中,测量矩阵Φ为第一测量矩阵和第二测量矩阵。测量矩阵是得到原始振动信号测量值及其重构的关键,其性能直接影响着振动信号的压缩测量能否成功实现。融入振动信号在正交字典傅里叶基上稀疏性的结构信息到测量矩阵的设计中,得到其最优的测量矩阵能够显著提升压缩测量性能。
该第一测量矩阵为机械振动信号的最优测量矩阵,该最优测量矩阵的构造算法为:根据最优测量理论,其设计过程是:已知傅里叶正交基Ψ、初始测量矩阵Φ0、阈值t和尺度下降因子γ,得到使得μt{D=ΦΨ}最小的矩阵Φ∈RM×N作为测量矩阵,称此时的Φ为傅里叶稀疏变换基对应的最优测量矩阵。
如图3,该最优测量矩阵具体构造过程如下:
基于最优测量矩阵构造算法,得到振动信号的最优测量矩阵如下:
输入参数:傅里叶正交稀疏变换基ΨDFT∈CN×N,测量矩阵Φ∈RM×N的初始值Φ0为一个高斯随机矩阵或托普利兹矩阵,令压缩感知矩阵D=ΦΨ,t表示阈值,γ表示尺度下降因子(0<γ<1),M表示测量次数,Iter表示迭代次数,循环变量k,初始值k=0;离散傅里叶变换(DFT)(DiscreteFourierTransform)
具体步骤:
(1)根据振动信号在傅里叶基上的稀疏性,确定测量矩阵Φ的行数,M≥2k,k表示信号稀疏度;
(2)计算感知矩阵D=ΦΨ;
D=ΦΨ记为感知矩阵,D必须满足零空间特性(NSP),即D的零空间中不能包含稀疏度为2k的振动信号,才能重构稀疏度为k的两个不同的振动信号。但是,要验证一个矩阵是否满足NSP是一个NP难问题,为了寻求更容易操作的条件,因此,产生了很多等价形式,其中最著名的理论之一是Spark理论,当且仅当Spark(D)>2k时,能从某一测量值中最多恢复一个与其对应的某一原始振动信号,然而NSP和Spark都没有考虑测量值中含有噪声的情况,Tao等人提出了约束等距特性(RIP)。要验证和设计出的矩阵D满足NSP、Spark和RIP其中的一个,都是一个NP难问题,为了避开NP难问题,提出了矩阵D的相干性。
(3)对感知矩阵D进行列单位化,得到矩阵
(4)计算感知矩阵的Gram矩阵,其中表示的转置矩阵;
(5)根据阈值t更新Gk,按如下关系得到矩阵
(6)减小的秩,对角化对角矩阵Λk的主对角线上的元素,保留绝对值前M大的元素,其他都设定为0,改变Λk后,重新计算
(7)求解Dk,Dk是M×N矩阵,通过感知矩阵得到经过一次迭代的振动信号优化测量矩阵;
(8)更新Φ,找到Φk+1,使得误差最小,即
(9)k=k+1;
直到k=Iter,循环结束。
输出:ΦDFT对应的最优测量矩阵ΦDFT。
在最优测量矩阵构造的具体步骤及图3构造流程图中,初始测量矩阵与傅里叶稀疏变换基作乘积得到感知矩阵,对感知矩阵进行处理,最后再通过感知矩阵求得优化的测量矩阵,这样就把振动信号在傅里叶基上的稀疏性结构信息融入到了测量矩阵的设计中,并最终得到振动信号在傅里叶基上对应的最优测量矩阵。
如图4,该第二测量矩阵为高斯分布循环测量矩阵:
假设构造的高斯分布循环测量矩阵为Φ,Φ∈RM×N(其中M为矩阵的行数,N为列数,且M<<N)
高斯分布循环测量矩阵Φ具体形式为:
高斯分布循环测量矩阵构造的具体步骤为:
(1)生成N个服从高斯分布的随机数作为测量矩阵的第一行,即:
(2)将φ1通过循环移位生成φ2,即:
(3)在向量中,随机取出n个元素(1≤n≤N),记为[a1a2…an];
(4)其中αi~N(0,1),将放在[a1a2…an]在φ2中对应位置上,得到
(5)根据以上第(2)、(3)、(4)步的向量生成规律,重复第(2)、(3)、(4)步,得到M-1个行向量;
(6)将以上得到的M个行向量按生成的顺序组成矩阵。
高斯分布循环测量矩阵的构造过程中,在情况下,当n=N时,构造得到的高斯分布循环测量矩阵本质上是高斯随机测量矩阵;当n=1时,构造得到的高斯分布循环测量矩阵本质上是循环测量矩阵。
上述的测量矩阵构造与经典的循环矩阵构造方式不同,没有生成一个方阵再随机选取,而是直接生成了需要规模的测量矩阵。
所述的步骤(4)机械振动信号的重构中的机械振动信号是可压缩信号或近似稀疏信号,根据凸优化理论,对振动信号的重构就是对振动信号L0-范数的求解,在信号具有稀疏性的情况下,L0-范数和L1-范数等价,所以求解更简洁的L1-范数会得到同样的解;在实际应用中,采用正交匹配追踪OMP完成对机械振动信号的重构。
本发明专利在压缩感知理论基础上基于机械振动信号的稀疏性实现压缩采集机械振动信号。首先,新型的基于稀疏性机械振动信号压缩采集方法通过图2所示的原理框图实现;其次,通过机械振动信号的最优测量矩阵或者高斯分布循环测量矩阵对机械振动信号进行压缩采样;最后,采用正交匹配追踪算法(OMP)实现原始振动信号的重构。
原始机械振动信号为f∈RN,在傅里叶正交变换基Ψ上是稀疏信号,得到变换系数α=ΨTf,其中ΨT表示Ψ的转置矩阵;然后利用设计的测量矩阵对机械振动信号在傅里叶正交字典下的稀疏系数进行测量,得到压缩测量值y=ΦΨTf;最后,采用L1-范数意义下的优化问题求解原始机械振动信号:
其中表示f的重构信号。
本发明压缩采集机械振动信号方法的具体步骤如下:
1.机械振动信号的稀疏表示和稀疏度估计
假设机械振动信号为f∈RN,其中N为振动信号维数,采用傅里叶正交基Ψ对振动信号进行稀疏分解,得到振动信号的稀疏系数α=ΨTf,稀疏系数的绝对值大于等于3的数量确定为机械振动信号的稀疏度。
2.测量矩阵的构造
所述的稀疏度k确定后,根据测量次数M必须满足Spark理论不等式M≥2k确定保证信息不丢失的最少测量次数,即确定了测量矩阵的大小M×N。测量矩阵的具体构造形式是权利要求书中的机械振动信号的最优测量矩阵或者高斯分布循环测量矩阵。
3.压缩测量值的获得
y=Φf=ΦΨα,用机械振动信号的最优测量矩阵或者高斯分布循环测量矩阵采集原始机械振动信号,得到压缩测量值。
4.原始机械振动信号的重构
机械振动信号是可压缩信号(近似稀疏信号),根据凸优化理论,对振动信号的重构就是对振动信号进行L1-范数的求解,采用正交匹配追踪(OMP)完成对机械振动信号的重构。
Claims (10)
1.一种机械振动信号数据压缩采集方法,其特征在于具体步骤如下:
(1)振动信号的稀疏分解;假设机械振动信号为f∈RN,其中N为振动信号维数,采用傅里叶正交基Ψ对振动信号进行稀疏分解,得到振动信号的稀疏系数α=ΨTf,其中ΨT表示Ψ的转置矩阵;
(2)采集得到全局非自适应压缩测量值;确定一个测量矩阵Φ∈RM×N(M<<N)对振动信号的有用信息进行高效获取,其中M为压缩采集信号的维数;
(3)传输、存储和处理压缩测量值y;根据测量矩阵对机械振动信号在傅里叶正交字典下的稀疏系数进行测量,得到压缩测量值y=ΦΨα,并将压缩测量值存储在采样存储器中;
(4)基于振动信号的稀疏性,采用最优化算法,重构出原始振动信号;
在机械振动信号重构处理端,求解1-范数凸优化问题,
min‖α‖s.t.y=ΦΨα=Φf
其中Ψ为傅里叶基,具体算法采用正交匹配追踪算法,实现原始机械振动信号的重构;
形成通过测量矩阵使得高维机械振动信号投影为低维的压缩采集值,并在信号重构端求解最优化问题,高概率地用低维压缩采集值重构出高维的原始机械振动信号的结构。
2.根据权利要求1所述的机械振动信号数据压缩采集方法,其特征在于:所述的步骤(1)采用傅里叶正交基Ψ对振动信号进行稀疏分解为选取零范数最最小的字典作为稀疏分解基。
3.根据权利要求1所述的机械振动信号数据压缩采集方法,其特征在于:所述的步骤(1)经过稀疏变换后得到振动信号在傅里叶正交基上的稀疏系数,同时根据振动信号的特征,估计振动信号在稀疏域上的稀疏度k,稀疏度k是振动信号中所含频率成分的数量。
4.根据权利要求3所述的机械振动信号数据压缩采集方法,其特征在于:根据振动信号在傅里叶基上的稀疏度k,确定测量矩阵的行数M。
5.根据权利要求4所述的机械振动信号数据压缩采集方法,其特征在于:所述的稀疏度k确定后,测量矩阵的行数M必须满足Spark理论不等式M≥2k,以确定保证信息不丢失的最少测量次数;其中k为振动信号在傅里叶基上的稀疏度。
6.根据权利要求1所述的机械振动信号数据压缩采集方法,其特征在于:所述的步骤(2)中,测量矩阵Φ为第一测量矩阵和第二测量矩阵。
7.根据权利要求6所述的机械振动信号数据压缩采集方法,其特征在于:该第一测量矩阵为机械振动信号的最优测量矩阵,该最优测量矩阵的构造算法为:根据最优测量理论,其设计过程是:已知傅里叶正交基Ψ、初始测量矩阵Φ0、阈值t和尺度下降因子γ,得到使得μt{D=ΦΨ}最小的矩阵Φ∈RM×N作为测量矩阵,称此时的Φ为傅里叶稀疏变换基对应的最优测量矩阵。
8.根据权利要求7所述的机械振动信号数据压缩采集方法,其特征在于:该最优测量矩阵具体构造过程如下:
基于最优测量矩阵构造算法,得到振动信号的最优测量矩阵如下:
输入参数:傅里叶正交稀疏变换基ΨDFT∈CN×N,测量矩阵Φ∈RM×N的初始值Φ0为一个高斯随机矩阵或托普利兹矩阵,令压缩感知矩阵D=ΦΨ,t表示阈值,γ表示尺度下降因子(0<γ<1),M表示测量次数,Iter表示迭代次数,循环变量k,初始值k=0;
具体步骤:
(1)根据振动信号在傅里叶基上的稀疏性,确定测量矩阵Φ的行数,M≥2k,k表示信号稀疏度;
(2)计算感知矩阵D=ΦΨ;
(3)对感知矩阵D进行列单位化,得到矩阵
(4)计算感知矩阵的Gram矩阵,其中表示的转置矩阵;
(5)根据阈值t更新Gk,按如下关系得到矩阵
(6)减小的秩,对角化对角矩阵Λk的主对角线上的元素,保留绝对值前M大的元素,其他都设定为0,改变Λk后,重新计算
(7)求解Dk,Dk是M×N矩阵,通过感知矩阵得到经过一次迭代的振动信号优化测量矩阵;
(8)更新Φ,找到Φk+1,使得误差最小,即
(9)k=k+1;
直到k=Iter,循环结束。
输出:ΦDFT对应的最优测量矩阵ΦDFT。
9.根据权利要求6所述的机械振动信号数据压缩采集方法,其特征在于:该第二测量矩阵为高斯分布循环测量矩阵:
假设构造的高斯分布循环测量矩阵为Φ,Φ∈RM×N(其中M为矩阵的行数,N为列数,且M<<N)
高斯分布循环测量矩阵Φ具体形式为:
高斯分布循环测量矩阵构造的具体步骤为:
(1)生成N个服从高斯分布的随机数作为测量矩阵的第一行,即:
(2)将φ1通过循环移位生成φ2,即:
(3)在向量中,随机取出n个元素(1≤n≤N),记为[a1a2…an];
其中αi~N(0,1),将放在[a1a2…an]在φ2中对应位置上,得到
(5)根据以上第(2)、(3)、(4)步的向量生成规律,重复第(2)、(3)、(4)步,得到M-1个行向量;
(6)将以上得到的M个行向量按生成的顺序组成矩阵。
高斯分布循环测量矩阵的构造过程中,在情况下,当n=N时,构造得到的高斯分布循环测量矩阵本质上是高斯随机测量矩阵;当n=1时,构造得到的高斯分布循环测量矩阵本质上是循环测量矩阵。
10.根据权利要求1所述的机械振动信号数据压缩采集方法法,其特征在于:所述的步骤(4)机械振动信号的重构中的机械振动信号是可压缩信号或近似稀疏信号,根据凸优化理论,对振动信号的重构就是对振动信号L0-范数的求解,在信号具有稀疏性的情况下,L0-范数和L1-范数等价,所以求解更简洁的L1-范数会得到同样的解;在实际应用中,采用正交匹配追踪OMP完成对机械振动信号的重构。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510477122.4A CN105181122B (zh) | 2015-08-06 | 2015-08-06 | 机械振动信号数据压缩采集方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510477122.4A CN105181122B (zh) | 2015-08-06 | 2015-08-06 | 机械振动信号数据压缩采集方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105181122A true CN105181122A (zh) | 2015-12-23 |
CN105181122B CN105181122B (zh) | 2018-06-15 |
Family
ID=54903367
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510477122.4A Expired - Fee Related CN105181122B (zh) | 2015-08-06 | 2015-08-06 | 机械振动信号数据压缩采集方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105181122B (zh) |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106454892A (zh) * | 2016-08-30 | 2017-02-22 | 衢州学院 | 一种基于量化和事件的集值卡尔曼滤波方法 |
US10082520B2 (en) | 2016-03-18 | 2018-09-25 | Simmonds Precision Products, Inc. | Rotational frequency estimation from sensed vibrations based on a supervised learning method |
CN108599773A (zh) * | 2018-04-16 | 2018-09-28 | 兰州理工大学 | 一种基于确定性测量矩阵的振动信号数据压缩采集方法 |
CN108711151A (zh) * | 2018-05-22 | 2018-10-26 | 广东工业大学 | 一种焊接缺陷检测方法、装置、设备、存储介质及系统 |
CN109039341A (zh) * | 2018-07-26 | 2018-12-18 | 深圳大学 | 多量测压缩感知的感知矩阵构建方法、系统及存储介质 |
CN109540279A (zh) * | 2018-11-30 | 2019-03-29 | 厦门大学 | 机床高速铣削过程欠采样动态信号的逆压缩感知恢复方法 |
CN110044608A (zh) * | 2019-04-19 | 2019-07-23 | 石家庄铁道大学 | 基于mckd的振动信号压缩感知方法 |
CN110849461A (zh) * | 2019-12-04 | 2020-02-28 | 江苏方天电力技术有限公司 | 一种调相机振动信号采集与存储方法、系统 |
CN110987435A (zh) * | 2019-12-18 | 2020-04-10 | 北京化工大学 | 一种基于改进k-svd字典学习和压缩感知的轴承故障诊断方法 |
CN111010191A (zh) * | 2019-12-12 | 2020-04-14 | 国家电网有限公司大数据中心 | 数据采集方法、系统、设备及存储介质 |
CN112415439A (zh) * | 2020-11-04 | 2021-02-26 | 华北电力大学 | 一种利用稀疏投影的振声检测信号滤波方法和系统 |
CN112629651A (zh) * | 2020-10-16 | 2021-04-09 | 国网江苏省电力有限公司盐城供电分公司 | 一种基于压缩感知的输电线路舞动信息重构方法 |
CN112927710A (zh) * | 2021-01-21 | 2021-06-08 | 安徽南瑞继远电网技术有限公司 | 一种基于无监督方式的电力变压器工况噪声分离方法 |
CN113670435A (zh) * | 2021-08-20 | 2021-11-19 | 西安石油大学 | 一种基于压缩感知技术的井下振动测量装置及其测量方法 |
CN113776653A (zh) * | 2021-09-13 | 2021-12-10 | 中国人民解放军陆军工程大学 | 振动信号的采样重构系统及方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102879782A (zh) * | 2012-09-25 | 2013-01-16 | 北京理工大学 | 基于分数阶傅里叶变换的压缩感知sar成像方法 |
CN104458170A (zh) * | 2014-11-07 | 2015-03-25 | 桂林电子科技大学 | 机械装备监测振动信号的时频图处理方法及系统 |
CN104935349A (zh) * | 2015-06-04 | 2015-09-23 | 西南交通大学 | 一种振动信号压缩采样方法 |
-
2015
- 2015-08-06 CN CN201510477122.4A patent/CN105181122B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102879782A (zh) * | 2012-09-25 | 2013-01-16 | 北京理工大学 | 基于分数阶傅里叶变换的压缩感知sar成像方法 |
CN104458170A (zh) * | 2014-11-07 | 2015-03-25 | 桂林电子科技大学 | 机械装备监测振动信号的时频图处理方法及系统 |
CN104935349A (zh) * | 2015-06-04 | 2015-09-23 | 西南交通大学 | 一种振动信号压缩采样方法 |
Non-Patent Citations (3)
Title |
---|
刘晓静: "压缩感知中的测量矩阵设计与重构算法研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
王帅: "基于压缩感知的语音信号压缩重构算法研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
秦晓伟: "基于压缩感知超声图像重构研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
Cited By (23)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10082520B2 (en) | 2016-03-18 | 2018-09-25 | Simmonds Precision Products, Inc. | Rotational frequency estimation from sensed vibrations based on a supervised learning method |
CN106454892A (zh) * | 2016-08-30 | 2017-02-22 | 衢州学院 | 一种基于量化和事件的集值卡尔曼滤波方法 |
CN106454892B (zh) * | 2016-08-30 | 2020-02-11 | 衢州学院 | 一种基于量化和事件的集值卡尔曼滤波方法 |
CN108599773B (zh) * | 2018-04-16 | 2020-04-10 | 兰州理工大学 | 一种基于确定性测量矩阵的振动信号数据压缩采集方法 |
CN108599773A (zh) * | 2018-04-16 | 2018-09-28 | 兰州理工大学 | 一种基于确定性测量矩阵的振动信号数据压缩采集方法 |
CN108711151A (zh) * | 2018-05-22 | 2018-10-26 | 广东工业大学 | 一种焊接缺陷检测方法、装置、设备、存储介质及系统 |
CN109039341A (zh) * | 2018-07-26 | 2018-12-18 | 深圳大学 | 多量测压缩感知的感知矩阵构建方法、系统及存储介质 |
CN109039341B (zh) * | 2018-07-26 | 2020-08-11 | 深圳大学 | 多量测压缩感知的感知矩阵构建方法、系统及存储介质 |
CN109540279A (zh) * | 2018-11-30 | 2019-03-29 | 厦门大学 | 机床高速铣削过程欠采样动态信号的逆压缩感知恢复方法 |
CN109540279B (zh) * | 2018-11-30 | 2020-03-24 | 厦门大学 | 机床高速铣削过程欠采样动态信号的逆压缩感知恢复方法 |
CN110044608A (zh) * | 2019-04-19 | 2019-07-23 | 石家庄铁道大学 | 基于mckd的振动信号压缩感知方法 |
CN110849461B (zh) * | 2019-12-04 | 2021-08-31 | 江苏方天电力技术有限公司 | 一种调相机振动信号采集与存储方法、系统 |
CN110849461A (zh) * | 2019-12-04 | 2020-02-28 | 江苏方天电力技术有限公司 | 一种调相机振动信号采集与存储方法、系统 |
CN111010191A (zh) * | 2019-12-12 | 2020-04-14 | 国家电网有限公司大数据中心 | 数据采集方法、系统、设备及存储介质 |
CN111010191B (zh) * | 2019-12-12 | 2023-11-07 | 国家电网有限公司大数据中心 | 数据采集方法、系统、设备及存储介质 |
CN110987435A (zh) * | 2019-12-18 | 2020-04-10 | 北京化工大学 | 一种基于改进k-svd字典学习和压缩感知的轴承故障诊断方法 |
CN112629651A (zh) * | 2020-10-16 | 2021-04-09 | 国网江苏省电力有限公司盐城供电分公司 | 一种基于压缩感知的输电线路舞动信息重构方法 |
CN112415439B (zh) * | 2020-11-04 | 2021-11-19 | 华北电力大学 | 一种利用稀疏投影的振声检测信号滤波方法和系统 |
CN112415439A (zh) * | 2020-11-04 | 2021-02-26 | 华北电力大学 | 一种利用稀疏投影的振声检测信号滤波方法和系统 |
CN112927710A (zh) * | 2021-01-21 | 2021-06-08 | 安徽南瑞继远电网技术有限公司 | 一种基于无监督方式的电力变压器工况噪声分离方法 |
CN113670435A (zh) * | 2021-08-20 | 2021-11-19 | 西安石油大学 | 一种基于压缩感知技术的井下振动测量装置及其测量方法 |
CN113776653A (zh) * | 2021-09-13 | 2021-12-10 | 中国人民解放军陆军工程大学 | 振动信号的采样重构系统及方法 |
CN113776653B (zh) * | 2021-09-13 | 2023-08-11 | 中国人民解放军陆军工程大学 | 振动信号的采样重构系统及方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105181122B (zh) | 2018-06-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105181122A (zh) | 机械振动信号数据压缩采集方法 | |
CN107192878A (zh) | 一种基于压缩感知的电力系统谐波检测方法及装置 | |
CN103368578A (zh) | 基于压缩感知的分布式无线传感器网络节点信号采样方法 | |
CN106055919A (zh) | 一种基于遥测数据小波变换的卫星异常检测方法 | |
CN103280084A (zh) | 一种多参数实时监测的数据采集方法 | |
CN104935349A (zh) | 一种振动信号压缩采样方法 | |
CN107741579B (zh) | 基于压缩感知子空间重构的toa矿井目标定位方法 | |
CN107942377A (zh) | 一种地震数据压缩与重构方法 | |
CN103823133A (zh) | 一种基于压缩感知的在线电能质量监测系统 | |
CN102930172A (zh) | 一种基于emd的海浪多尺度特征及波动参数的提取方法 | |
CN102970707A (zh) | 基于压缩感知的无线传输数据丢失恢复方法 | |
CN102622331B (zh) | 一种基于压缩感知的高斯矩阵优化方法 | |
CN103558498B (zh) | 基于小波分析的绝缘子污闪泄漏电流信号稀疏表示方法 | |
CN110311685B (zh) | 时序贝叶斯压缩采样及信号解压缩重构方法及数据丢失恢复方法 | |
CN103269525A (zh) | 物联网系统的多传感器数据分布式压缩方法及其装置 | |
CN113225089B (zh) | 一种工业机器人多通道监测数据压缩方法及系统 | |
CN102355752A (zh) | 基于膨胀图的无线传感器网络压缩感知测量矩阵和重构方法 | |
CN103258134B (zh) | 一种高维的振动信号的降维处理方法 | |
CN103744828B (zh) | 一种测量值缺失的压缩感知重构方法 | |
CN106019369B (zh) | 一种改进的seg‑y文件中地震数据无损压缩算法 | |
CN105682141B (zh) | 基于邻居辅助的无线传感器网络数据采集方法 | |
CN104125459B (zh) | 基于支撑集和信号值检测的视频压缩感知重构方法 | |
CN104242947A (zh) | 一种基于Haar小波树的SAMP重构方法 | |
CN102802199B (zh) | 基于压缩感知的无线传感器实时监测系统的数据处理方法 | |
CN104867166A (zh) | 一种基于稀疏字典学习的油井示功图压缩存储方法 |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180615 |