CN106596724B - 一种窄带Lamb波检测数据稀疏压缩优化重建方法 - Google Patents
一种窄带Lamb波检测数据稀疏压缩优化重建方法 Download PDFInfo
- Publication number
- CN106596724B CN106596724B CN201611121562.7A CN201611121562A CN106596724B CN 106596724 B CN106596724 B CN 106596724B CN 201611121562 A CN201611121562 A CN 201611121562A CN 106596724 B CN106596724 B CN 106596724B
- Authority
- CN
- China
- Prior art keywords
- lamb wave
- narrow
- detection data
- wave detection
- band
- 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.)
- Active
Links
- 235000019687 Lamb Nutrition 0.000 title claims abstract description 73
- 238000001514 detection method Methods 0.000 title claims abstract description 66
- 238000000034 method Methods 0.000 title claims abstract description 28
- 238000005457 optimization Methods 0.000 title claims abstract description 28
- 238000007906 compression Methods 0.000 title claims abstract description 23
- 230000006835 compression Effects 0.000 title claims abstract description 23
- 238000005259 measurement Methods 0.000 claims abstract description 47
- 239000011159 matrix material Substances 0.000 claims abstract description 32
- BLRBOMBBUUGKFU-SREVYHEPSA-N (z)-4-[[4-(4-chlorophenyl)-5-(2-methoxy-2-oxoethyl)-1,3-thiazol-2-yl]amino]-4-oxobut-2-enoic acid Chemical compound S1C(NC(=O)\C=C/C(O)=O)=NC(C=2C=CC(Cl)=CC=2)=C1CC(=O)OC BLRBOMBBUUGKFU-SREVYHEPSA-N 0.000 claims abstract description 24
- 238000001914 filtration Methods 0.000 claims abstract description 10
- 238000012545 processing Methods 0.000 claims abstract description 6
- 239000000463 material Substances 0.000 claims abstract description 3
- 230000006870 function Effects 0.000 claims description 12
- 238000011084 recovery Methods 0.000 claims description 12
- 230000009466 transformation Effects 0.000 claims description 4
- 230000005540 biological transmission Effects 0.000 abstract description 9
- 238000003325 tomography Methods 0.000 abstract description 9
- 238000013144 data compression Methods 0.000 abstract description 4
- 238000003384 imaging method Methods 0.000 description 6
- 229910000831 Steel Inorganic materials 0.000 description 3
- 230000007547 defect Effects 0.000 description 3
- 239000010959 steel Substances 0.000 description 3
- 238000011161 development Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 229910052782 aluminium Inorganic materials 0.000 description 1
- XAGFODPZIPBFFR-UHFFFAOYSA-N aluminium Chemical compound [Al] XAGFODPZIPBFFR-UHFFFAOYSA-N 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 229910052751 metal Inorganic materials 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 238000009659 non-destructive testing Methods 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
- 238000004904 shortening Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/04—Analysing solids
- G01N29/06—Visualisation of the interior, e.g. acoustic microscopy
- G01N29/0654—Imaging
- G01N29/0672—Imaging by acoustic tomography
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/22—Details, e.g. general constructional or apparatus details
- G01N29/24—Probes
- G01N29/2412—Probes using the magnetostrictive properties of the material to be examined, e.g. electromagnetic acoustic transducers [EMAT]
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/36—Detecting the response signal, e.g. electronic circuits specially adapted therefor
- G01N29/42—Detecting the response signal, e.g. electronic circuits specially adapted therefor by frequency filtering or by tuning to resonant frequency
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/44—Processing the detected response signal, e.g. electronic circuits specially adapted therefor
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/44—Processing the detected response signal, e.g. electronic circuits specially adapted therefor
- G01N29/4409—Processing the detected response signal, e.g. electronic circuits specially adapted therefor by comparison
- G01N29/4427—Processing the detected response signal, e.g. electronic circuits specially adapted therefor by comparison with stored values, e.g. threshold values
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/44—Processing the detected response signal, e.g. electronic circuits specially adapted therefor
- G01N29/4481—Neural networks
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/44—Processing the detected response signal, e.g. electronic circuits specially adapted therefor
- G01N29/46—Processing the detected response signal, e.g. electronic circuits specially adapted therefor by spectral analysis, e.g. Fourier analysis or wavelet analysis
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/02—Indexing codes associated with the analysed material
- G01N2291/023—Solids
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/02—Indexing codes associated with the analysed material
- G01N2291/028—Material parameters
- G01N2291/0289—Internal structure, e.g. defects, grain size, texture
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/04—Wave modes and trajectories
- G01N2291/042—Wave modes
- G01N2291/0427—Flexural waves, plate waves, e.g. Lamb waves, tuning fork, cantilever
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/10—Number of transducers
- G01N2291/106—Number of transducers one or more transducer arrays
Landscapes
- Physics & Mathematics (AREA)
- Immunology (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- Pathology (AREA)
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Signal Processing (AREA)
- Acoustics & Sound (AREA)
- Mathematical Physics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Evolutionary Computation (AREA)
- Electromagnetism (AREA)
- Artificial Intelligence (AREA)
- Radar Systems Or Details Thereof (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
本发明公开了一种窄带Lamb波检测数据稀疏压缩优化重建方法，在待测材料检测区域两侧分别设置发射EMAT阵列和接收EMAT阵列，令发射阵列中的各个EMAT按照预设的顺序分别以窄带频率激发所需模态Lamb波，对侧接收阵列接收Lamb波，并对接收到的Lamb波检测信号进行窄带滤波处理，形成窄带Lamb波检测数据；得到检测数据后，经离散傅里叶变换，再采用教与学优化算法TLBO重构原始数据，最优化测量向量和恢复矩阵，并将低维测量向量传输到上位机；本发明解决了大面积板件电磁超声Lamb波层析成像数据量大、传输时间长导致检测效率低的问题，实现了快速、大比例和高质量的数据压缩和传输，为大面积板件Lamb波实时层析成像打下坚实基础，具有广阔的应用前景。
Description
技术领域
本发明属于金属无损检测技术领域，涉及一种窄带Lamb波检测数据稀疏压缩优化重建方法。
背景技术
电磁超声Lamb波缺陷层析成像技术可以快速有效地获得缺陷的轮廓和尺寸等具体信息，为实施大面积板件的高质量层析成像，所使用EMAT的数量可达数十对，相应一发一收EMAT检测数据可达数百种组，且每帧波形需要进行多次Lamb波重复收发以进行降噪处理，因此，单次成像产生的Lamb波检测数据总量庞大，以每一帧波形含有500kB检测数据估算，则单次成像数据总量最高可达5GB以上。在对大面积板件实施扫查成像的过程中，为保证一定的扫查速度，单次成像通常需要在十几秒内完成，装置主机需要在该规定时间内将大量检测数据传输到上位机。然而，受到目前主流数采设备和上位机接口传输速度的限制，装置主机与上位机之间的实际数据传输速度通常低于50MB/s，这使得装置完成上述大量检测数据的传输所需的时间将大于102s，难以满足单次成像的时间要求。
上述问题是制约Lamb波层析成像技术发展及缺陷实时扫查成像的一个瓶颈问题。目前，缺少切实可行的用于层析成像的电磁超声Lamb波检测数据压缩重建方法，开发快速、大比例和高质量的电磁超声Lamb波检测数据压缩和重建方法，对于促进电磁超声Lamb波实时层析成像技术的发展具有重要意义。
发明内容
本发明解决的问题在于提供一种窄带Lamb波检测数据稀疏压缩优化重建方法，解决了大面积板件电磁超声Lamb波层析成像数据量大、传输时间长导致检测效率低的问题。
本发明是通过以下技术方案来实现：
一种窄带Lamb波检测数据稀疏压缩优化重建方法，包括以下操作：
1)在待测材料检测区域两侧分别设置发射EMAT阵列和接收EMAT阵列；
2)发射阵列中的各个EMAT按照预设的顺序分别以窄带频率fc激发所需模态Lamb波；对侧的接收阵列接收Lamb波，并对接收到的Lamb波检测信号进行窄带滤波处理，滤波中心频率同样为fc，形成窄带Lamb波检测数据x(n)，其中n＝1,2，…，N，N为正整数；
3)利用离散傅里叶变换对窄带Lamb波检测数据进行分析，得到原始检测数据的稀疏表示及其稀疏度K，根据稀疏度K和原始数据长度N确定低维测量数的初始值M；其中，K、M为正整数；
4)采用高斯随机矩阵Φ测量原始数据向量x(n)，得到测量向量y(m)，高斯随机矩阵Φ的维度为M×N，测量向量的维度为低维测量数的初始值M，且M<<N，并得到恢复矩阵A；
5)根据测量向量和恢复矩阵，利用教与学优化算法TLBO重构窄带Lamb波原始检测数据向量，并计算出重构误差；
6)判断重构误差是否小于设定值，若是则进行步骤7；若不是判断该方法目前的迭代次数是否大于设定值，若是则进行步骤7，若不是对测量数M加1，并返回步骤4；
7)重构完成，提取此时的低维测量向量y(m)，并将其传输到上位机进行存储和处理。
所述EMAT阵列采用圆形回折线圈设计，通过工作频率和导线间距的匹配关系，激发和接收所需模态的窄带Lamb波；
用高阶滤波电路对接收EMAT阵列所接收的Lamb波信号进行窄带滤波，得到窄带Lamb波检测数据x(n)。
所述利用离散傅里叶变换对窄带Lamb波检测数据进行分析为：窄带Lamb波检测数据为x(n)，其中n＝1,2，…，N，N为正整数；则其离散傅里叶变换为：
其中，e是自然底数，j是虚数单位，k是计算傅里叶变换时用来对角度计数的变量；
采用矩阵变换形式表示为：
x＝Ψs
其中，Ψ为离散傅里叶变换正交基，矩阵维数为N×N；
s为原始数据x在频域下的稀疏表示，向量维数为N，稀疏度为K。
所述根据稀疏度K和原始数据长度N确定低维测量数的初始值M，确定方法为：
其中，min为取两者最小值的函数；
Floor为取整函数，取向负无穷方向最靠近的整数。
所述采用高斯随机矩阵Φ测量原始数据向量x(n)，得到测量向量y(m)为：
y＝Φx＝ΦΨs＝As
A＝ΦΨ
其中，恢复矩阵A的维度为M×N，M<<N。
所述采用教与学优化算法TLBO的重构为：
采用教与学优化算法TLBO求解欠定方程组：
y＝As1
求解出向量s1后，采用离散逆傅里叶变换得到窄带Lamb波检测重构数据x1，则重构误差e为：
其中，Norm为计算向量2范数的函数。
与现有技术相比，本发明具有以下有益的技术效果：
本发明提供的窄带Lamb波检测数据稀疏压缩优化重建方法，在深入分析窄带电磁超声Lamb波检测数据的稀疏特性之后，对系列窄带Lamb波检测数据先进行离散傅里叶变换，得到原始数据的稀疏表示并确定低维测量数初始值。再根据低维测量数，采用高斯随机矩阵测量原始数据并获得测量向量。根据测量向量和恢复矩阵，采用教与学优化算法TLBO重构原始数据。建立重构误差和测量数的迭代配合关系，最优化测量向量和恢复矩阵，最后将低维测量向量传输到上位机，实现了对窄带Lamb波检测数据的大比例压缩和高精度重构，压缩比例大于10，重构误差小于10％，从而大大缩短了数据传输时间。
本发明提供的窄带Lamb波检测数据稀疏压缩优化重建方法，计算准确、高效、快速，解决了大面积板件电磁超声Lamb波层析成像数据量大、传输时间长导致检测效率低的问题，对窄带Lamb波检测数据实现了快速、大比例和高质量的数据压缩和传输，为大面积板件Lamb波实时层析成像打下坚实基础，具有广阔的应用前景。
附图说明
图1为窄带Lamb波检测数据稀疏压缩优化重建流程图；
图2为窄带Lamb波检测数据重建结果图(2000点)；
图3为窄带Lamb波检测数据重建结果图(600点)。
具体实施方式
本发明提出的窄带Lamb波检测数据稀疏压缩优化重建方法，是建立一种全新的窄带Lamb波检测数据稀疏压缩优化重建模型上的。下面结合实施例对本发明做进一步说明：
一种窄带Lamb波检测数据稀疏压缩优化重建方法，包括以下操作：
步骤1：取一件厚度为5mm的待测钢板，在钢板两侧分别布置发射EMAT阵列(含12个发射EMAT)和接收EMAT阵列(含12个接收EMAT)，EMAT直径为40mm。
步骤2：用射频功率放大器分别激励发射阵列中的各个EMAT，用接收EMAT阵列接收Lamb波，激发频率fc为125kHz；采用四阶窄带滤波电路对接收到的Lamb波检测信号进行滤波，滤波中心频率fc为125kHz，取其中一帧原始窄带Lamb波检测数据x(n)，其中n＝1,2，…，N，数据点数N为2000。
步骤3：利用离散傅里叶变换(Discrete Fourier transform，DFT)对窄带Lamb波检测数据x(n)进行分析：
采用矩阵变换形式表示为：
x＝Ψs
其中，Ψ为离散傅里叶变换正交基，矩阵维数为N×N。得到向量s的稀疏度K＝142。
然后根据稀疏度K和原始数据长度N确定低维测量数的初始值M，确定方法为：
其中，min为取两者最小值的函数；
Floor为取整函数，取向负无穷方向最靠近的整数。
得到测量数M的初始值为M＝142。
步骤4：采用高斯随机矩阵Φ测量原始数据向量x(n)，得到测量向量y(m)，高斯随机矩阵Φ的维度为M×N，测量向量的维度为M，且M<<N，并得到恢复矩阵A：
y＝Φx＝ΦΨs＝As
A＝ΦΨ
步骤5：根据测量向量和恢复矩阵，利用教与学优化算法TLBO重构窄带Lamb波原始检测数据向量x1：
y＝As1
x1＝Ψs1
计算重构误差e：
其中，Norm为计算向量2范数的函数。
步骤6：判断重构误差是否小于设定值8％，若是，则进行步骤7；若不是，判断该方法目前的迭代次数是否大于设定值50，若是，则进行步骤7，若不是，对测量数M加1，并返回步骤4。
步骤7：重构完成，提取此时的低维测量向量y(m)，并将其传输到上位机进行存储和处理。
本实例最后得到的重构误差为7.6％，低维测量数为167，压缩比为12。重构结果如图2所示。
实施例2
步骤1：取一件厚度为3mm的待测铝板，在钢板两侧分别布置发射EMAT阵列(含14个发射EMAT)和接收EMAT阵列(含14个接收EMAT)，EMAT直径为50mm。
步骤2：用射频功率放大器分别激励发射阵列中的各个EMAT，用接收EMAT阵列接收Lamb波，激发频率为80kHz；采用四阶窄带滤波电路对接收到的Lamb波检测信号进行滤波，滤波中心频率为80kHz，取其中一帧原始窄带Lamb波检测数据x(n)，其中n＝1,2，…，N，数据点数N为600。
步骤3：利用离散傅里叶变换(Discrete Fourier transform，DFT)对窄带Lamb波检测数据x(n)进行分析：
采用矩阵变换形式表示为：
x＝Ψs
其中，
Ψ为离散傅里叶变换正交基，矩阵维数为N×N。得到向量s的稀疏度K＝62。根据稀疏度K和原始数据长度N确定低维测量数的初始值M，确定方法为：
其中，
min为取两者最小值的函数；
Floor为取整函数，取向负无穷方向最靠近的整数。
得到测量数M的初始值为M＝61。
步骤4：采用高斯随机矩阵Φ测量原始数据向量x(n)，得到测量向量y(m)，高斯随机矩阵Φ的维度为M×N，测量向量的维度为M，且M<<N，并得到恢复矩阵A：
y＝Φx＝ΦΨs＝As
A＝ΦΨ
步骤5：根据测量向量和恢复矩阵，利用教与学优化算法TLBO重构窄带Lamb波原始检测数据向量x1：
y＝As1
x1＝Ψs1
计算重构误差e：
其中，
Norm为计算向量2范数的函数。
步骤6：判断重构误差是否小于设定值7％，若是，则进行步骤7；若不是，判断该方法目前的迭代次数是否大于设定值60，若是，则进行步骤7，若不是，对测量数M加1，并返回步骤4。
步骤7：重构完成，提取此时的低维测量向量y(m)，并将其传输到上位机进行存储和处理。
本实例最后得到的重构误差为6.8％，低维测量数为63，压缩比为9.5。重构结果如图3所示。
图2和图3中重构信号的轮廓、幅值大小以及各信号出现的时间均与原始信号保持了极高的一致性，在大比例压缩的前提下，包含了原始信号足够丰富和准确的信息。
以上所述，仅为本发明较佳的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到的变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应该以权利要求的保护范围为准。
Claims (6)
1.一种窄带Lamb波检测数据稀疏压缩优化重建方法，其特征在于，包括以下操作：
1)在待测材料检测区域两侧分别设置发射EMAT阵列和接收EMAT阵列；其中EMAT阵列采用圆形回折线圈设计，通过工作频率和导线间距的匹配关系，激发和接收所需模态的窄带Lamb波；
2)发射阵列中的各个EMAT按照预设的顺序分别以窄带频率fc激发所需模态Lamb波；对侧的接收阵列接收Lamb波，并对接收到的Lamb波检测信号进行窄带滤波处理，滤波中心频率同样为fc，形成窄带Lamb波检测数据x(n)，其中n＝1,2，…，N，N为正整数；
3)利用离散傅里叶变换对窄带Lamb波检测数据进行分析，得到原始检测数据的稀疏表示及其稀疏度K，根据稀疏度K和原始数据长度N确定低维测量数的初始值M；其中，K、M为正整数；
4)采用高斯随机矩阵Φ测量原始数据向量x(n)，得到测量向量y(m)；其中高斯随机矩阵Φ的维度为M×N，测量向量的维度为低维测量数的初始值M，且M<<N，并得到恢复矩阵A；
5)根据测量向量和恢复矩阵，利用教与学优化算法TLBO重构窄带Lamb波原始检测数据向量，并计算出重构误差；
6)判断重构误差是否小于设定值，若是则进行步骤7；若不是判断该方法目前的迭代次数是否大于设定值，若是则进行步骤7，若不是对测量数M加1，并返回步骤4；
7)重构完成，提取此时的低维测量向量y(m)，并将其传输到上位机进行存储和处理。
2.根据权利要求1所述的窄带Lamb波检测数据稀疏压缩优化重建方法，其特征在于，用高阶窄带滤波电路对接收EMAT阵列所接收的Lamb波信号进行窄带滤波，得到窄带Lamb波检测数据x(n)。
3.根据权利要求1所述的窄带Lamb波检测数据稀疏压缩优化重建方法，其特征在于，所述利用离散傅里叶变换对窄带Lamb波检测数据进行分析为：窄带Lamb波检测数据为x(n)，其中n＝1,2，…，N，N为正整数；则其离散傅里叶变换为：
其中，e是自然底数，j是虚数单位，k是计算傅里叶变换时用来对角度计数的变量；
采用矩阵变换形式表示为：
x＝Ψs
其中，Ψ为离散傅里叶变换正交基，矩阵维数为N×N；
向量s为原始数据x在频域下的稀疏表示，向量维数为N，其稀疏度为K。
4.根据权利要求3所述的窄带Lamb波检测数据稀疏压缩优化重建方法，其特征在于，根据稀疏度K和原始数据长度N确定低维测量数的初始值M，确定方法为：
其中，min为取两者最小值的函数；
Floor为取整函数，取向负无穷方向最靠近的整数。
5.根据权利要求4所述的窄带Lamb波检测数据稀疏压缩优化重建方法，其特征在于，采用高斯随机矩阵Φ测量原始数据向量x(n)，得到测量向量y(m)为：
y＝Φx＝ΦΨs＝As
A＝ΦΨ
其中，恢复矩阵A的维度为M×N，M<<N。
6.根据权利要求1或5所述的窄带Lamb波检测数据稀疏压缩优化重建方法，其特征在于，采用教与学优化算法TLBO的重构为：
采用教与学优化算法TLBO求解欠定方程组：
y＝As1
求解出向量s1后，采用离散逆傅里叶变换得到窄带Lamb波检测重构数据x1，x1＝Ψs1，则重构误差e为：
其中，Norm为计算向量2范数的函数。
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611121562.7A CN106596724B (zh) | 2016-12-08 | 2016-12-08 | 一种窄带Lamb波检测数据稀疏压缩优化重建方法 |
| US15/829,586 US10962504B2 (en) | 2016-12-08 | 2017-12-01 | Method for compressing and reconstructing data |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611121562.7A CN106596724B (zh) | 2016-12-08 | 2016-12-08 | 一种窄带Lamb波检测数据稀疏压缩优化重建方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106596724A CN106596724A (zh) | 2017-04-26 |
| CN106596724B true CN106596724B (zh) | 2019-05-07 |
Family
ID=58598506
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201611121562.7A Active CN106596724B (zh) | 2016-12-08 | 2016-12-08 | 一种窄带Lamb波检测数据稀疏压缩优化重建方法 |
Country Status (2)
| Country | Link |
|---|---|
| US (1) | US10962504B2 (zh) |
| CN (1) | CN106596724B (zh) |
Families Citing this family (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP3309689B1 (en) * | 2016-10-17 | 2024-11-06 | Tata Consultancy Services Limited | Systems and methods for selecting optimal variables using modified teaching learning based search optimization technique |
| CN107607628B (zh) * | 2017-08-11 | 2019-05-24 | 西安交通大学 | 一种基于稀疏表示的超声导波频散补偿方法及其应用 |
| CN109506763B (zh) * | 2018-11-30 | 2020-09-25 | 中国人民解放军军事科学院国防科技创新研究院 | 一种基于教与学优化的矢量水听器稀疏布阵方法 |
| CN112748296B (zh) * | 2019-10-31 | 2023-05-23 | 青岛海尔智能技术研发有限公司 | 用于直流电的电参数监测方法、装置及直流家电 |
| CN111093166B (zh) * | 2019-12-06 | 2022-07-19 | 北京京航计算通讯研究所 | 在物联网中使用稀疏测量矩阵的压缩数据收集系统 |
| CN111130557B (zh) * | 2019-12-31 | 2023-11-17 | 秦皇岛职业技术学院 | 基于分布式拟牛顿投影追踪的数据重构方法 |
| CN112260693A (zh) * | 2020-09-14 | 2021-01-22 | 中国工程物理研究院电子工程研究所 | 一种数字调制信号重构方法 |
| CN113405689B (zh) * | 2021-06-22 | 2024-01-26 | 沈阳工业大学 | 一种基于压缩感知的声学ct温度场重建方法 |
| CN114324591B (zh) * | 2021-12-31 | 2023-09-01 | 上海理工大学 | 基于Radon变换的Lamb波波包飞行时间检测方法 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2136207A2 (en) * | 2008-06-20 | 2009-12-23 | The Boeing Company | Sub-wavelength ultrasound characterization of composite material |
| CN102279394A (zh) * | 2011-05-17 | 2011-12-14 | 西安电子科技大学 | 低频窄带高分辨超声探测成像方法 |
| CN104391040A (zh) * | 2014-11-25 | 2015-03-04 | 南京信息工程大学 | 一种相控阵结构健康监测Lamb波压缩及重构方法 |
| CN105486759A (zh) * | 2016-01-12 | 2016-04-13 | 清华大学 | 一种Lamb波时频能量密度析出走时提取方法 |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5629485A (en) * | 1994-12-13 | 1997-05-13 | The B.F. Goodrich Company | Contaminant detection sytem |
| US20120103097A1 (en) * | 2010-10-29 | 2012-05-03 | Borja Lopez Jauregui | Flexible EMAT Arrays for Monitoring Corrosion and Defect Propagation in Metal Components and Structures |
| US20140063314A1 (en) * | 2012-02-28 | 2014-03-06 | Aswin C Sankaranarayanan | System And Method Of Video Compressive Sensing For Spatial-Multiplexing Cameras |
| CN103020909B (zh) * | 2012-12-06 | 2015-02-18 | 清华大学 | 基于多尺度结构自相似与压缩感知的单图像超分辨率方法 |
| US8863044B1 (en) * | 2013-09-06 | 2014-10-14 | International Business Machines Corporation | Layout assessment method and system |
| US10240910B2 (en) * | 2016-09-07 | 2019-03-26 | Board Of Trustees Of Southern Illinois University On Behalf Of Southern Illinois University Carbondale | Systems and methods for compressive image sensor techniques utilizing sparse measurement matrices |
-
2016
- 2016-12-08 CN CN201611121562.7A patent/CN106596724B/zh active Active
-
2017
- 2017-12-01 US US15/829,586 patent/US10962504B2/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2136207A2 (en) * | 2008-06-20 | 2009-12-23 | The Boeing Company | Sub-wavelength ultrasound characterization of composite material |
| CN102279394A (zh) * | 2011-05-17 | 2011-12-14 | 西安电子科技大学 | 低频窄带高分辨超声探测成像方法 |
| CN104391040A (zh) * | 2014-11-25 | 2015-03-04 | 南京信息工程大学 | 一种相控阵结构健康监测Lamb波压缩及重构方法 |
| CN105486759A (zh) * | 2016-01-12 | 2016-04-13 | 清华大学 | 一种Lamb波时频能量密度析出走时提取方法 |
Non-Patent Citations (1)
| Title |
|---|
| 基于Lamb波结构健康监测的压缩感知方法研究;王强 等;《仪表技术与传感器》;20161031(第10期);摘要，第124页第左栏2-10段、右栏第1-17段，图1 |
Also Published As
| Publication number | Publication date |
|---|---|
| US20180164253A1 (en) | 2018-06-14 |
| US10962504B2 (en) | 2021-03-30 |
| CN106596724A (zh) | 2017-04-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106596724B (zh) | 一种窄带Lamb波检测数据稀疏压缩优化重建方法 | |
| CN106841402A (zh) | 一种基于贪婪算法的相控阵超声信号重构优化方法 | |
| CN103575808B (zh) | 基于多角度立体匹配的高实时定量超声检测方法 | |
| CN109649432B (zh) | 基于导波技术的云端平台钢轨完整性监测系统及方法 | |
| CN105004792A (zh) | 一种用于微裂纹检测的非线性超声相控阵成像方法 | |
| CN109340586B (zh) | 一种供水管道泄露的检测方法及系统 | |
| CN109085244A (zh) | 一种基于压电阵列的非线性Lamb波结构疲劳损伤层析成像方法 | |
| CN106814141A (zh) | 一种基于正交匹配追踪的相控阵超声信号压缩方法 | |
| CN113607817B (zh) | 一种管道环焊缝检测方法、系统、电子设备及介质 | |
| CN106644859B (zh) | 一种扫频式超声波聚焦式泥沙粒径分布在线测量方法 | |
| Li et al. | An Orthogonal Matching Pursuit based signal compression and reconstruction approach for electromechanical admittance based structural health monitoring | |
| Chang et al. | Corrosion monitoring using a new compressed sensing-based tomographic method | |
| Wang et al. | High resolution tomography of pipeline using multi-helical Lamb wave based on compressed sensing | |
| Kalimullah et al. | Multiresolution dynamic mode decomposition (mrDMD) of elastic waves for damage localisation in piezoelectric ceramic | |
| CN110057918B (zh) | 强噪声背景下的复合材料损伤定量识别方法及系统 | |
| CN108508093A (zh) | 一种工件缺陷高度的检测方法及系统 | |
| Jeon et al. | Compressive laser scanning with full steady state wavefield for structural damage detection | |
| Rosalie et al. | Structural health monitoring of composite structures using stress wave methods | |
| CN107727749A (zh) | 一种基于小波包融合特征提取算法的超声定量检测方法 | |
| CN112818762B (zh) | 一种大尺寸复合材料及其夹层结构快速无损检测方法 | |
| Xu et al. | Compressive sensing of ultrasonic array data with full matrix capture in nozzle welds inspection | |
| CN107045003A (zh) | 一种基于同步迭代重建算法的铝合金结构腐蚀损伤监测方法 | |
| CN113536894A (zh) | 缺陷检测方法、缺陷检测装置和计算机可读存储介质 | |
| CN117214294A (zh) | 管道损伤超声导波检测装置及其损伤信号三维重构方法 | |
| Yang et al. | Parametric identification of ultrasonic guided wave aliasing modes based on dispersion effect |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |