基于时频分布的盲信号分离方法研究长骨中的多模式导波信号

由于超声导波对长骨皮质骨的材料特性及厚度变化非常敏感,因此用导波评价长骨状况的研究已成为近年来的一个研究热点.在接收到的导波中,往往多种模式混叠在一起,使得进一步求取每种导波模式的群速度等与骨质状况密切相关的参数时遇到困难.本文采用基于时频分布的盲信号分离(TFBSS)方法对长骨中混叠的多模式进行分离,然后对分离后的单...     

High-resolution bone microstructure imaging based on ultrasonic frequency-domain full-waveform inversion

The main challenge in bone ultrasound imaging is the large acoustic impedance contrast and sound velocity differences between the bone and surrounding soft tissue. It is difficult for conventional pulse-echo modalities to give accurate ultrasound images for irregular bone boundaries and microstructures using uniform sound velocity assumption rather than getting a prior knowledge of sound speed. To overcome these limitations, this paper proposed a frequency-domain fullwaveform inversion(FDFWI) algorithm for bone quantitative imaging utilizing ultrasonic computed tomography(USCT).The forward model was calculated in the frequency domain by solving the full-wave equation. The inverse problem was solved iteratively from low to high discrete frequency components via minimizing a cost function between the modeled and measured data. A quasi-Newton method called the limited-memory Broyden–Fletcher–Goldfarb–Shanno algorithm(L-BFGS) was utilized in the optimization process. Then, bone images were obtained based on the estimation of the velocity and density. The performance of the proposed method was verified by numerical examples, from tubular bone phantom to single distal fibula model, and finally with a distal tibia-fibula pair model. Compared with the high-resolution peripheral quantitative computed tomography(HR-p QCT), the proposed FDFWI can also clearly and accurately presented the wavelength scaled pores and trabeculae in bone images. The results proved that the FDFWI is capable of reconstructing high-resolution ultrasound bone images with sub-millimeter resolution. The parametric bone images may have the potential for the diagnosis of bone disease.     

带非局部阻尼项的粘弹性方程解的衰减估计

In this paper,we consider the following viscoelastic equation u tt- △u +∫t 0 g(t-s)△u(s)ds + a(x)u t + u |u|r = 0 with initial condition and Dirichlet boundary condition.The decay property of the energy function closely depends on the properties of the relaxation function g(t) at infinity.In the previous works of [3,7,11],it was required that the relaxation function g(t) decay exponentially or polynomially as t → +∞.In the recent work of Messaoudi [12,13],it was shown that the energy decays at a similar rate of decay of the relaxation function,which is not necessarily dacaying in a polynomial or exponential fashion.Motivated by [12,13],under some assumptions on g(x),a(x) and r,and by introducing a new perturbed energy,we also prove the similar results for the above equation.     

Barker码激励超声导波在长骨检测中的应用

单脉冲激励的超声导波在长骨中传播时,信号的衰减大,导致接收信号的幅度很小,且噪声严重。Barker码激励能有效增大接收信号幅度,提高信噪比(SNR)。将其应用到超声导波长骨检测中,进行仿真和长骨实验,得到的信号分别用加权匹配滤波器和有限冲激响应-最小均方误差(FIR-LS)逆滤波器进行压缩,并与单脉冲激励的结果进行了对比。结果表明,对于13位的Barker码,采用加权匹配滤波器进行解码时,压缩信号幅度是单脉冲激励接收信号的13倍;而FIR-LS逆滤波器则达到-63.59 dB的峰值旁瓣水平(PSL),更好地抑制噪声。说明可以将Barker码激励超声导波应用于长骨的检测中。      

采用超声导波特征参量的骨质疏松识别方法

骨质的定期检测对骨质疏松的防治至关重要。本文研究了骨质疏松对超声导波在人体长骨中传播的影响。提出采用多尺度小波变换方法对接收到的导波信号进行处理,通过分析在不同传播距离下高阶小波细节分量所占信号总能量的变化,来判断是否患有骨质疏松症。在13位志愿者的小腿胫骨上进行超声测量,得到导波信号。经多尺度小波变换方法的分析处理结果显示在13位志愿者中,有7位志愿者的超声导波信号随着传播距离的改变,其主要频率成分发生了明显的变化,显示这7位志愿者患有骨质疏松症。这一诊断结果与X射线技术诊断结果相比,准确率可以达到92.3%,表明本文所提出的利用小波多尺度变换方法对长骨进行超声诊断具有较好的潜力。      

一类双重退化抛物方程组解的整体存在与爆破

研究一类具有非线性边界流的双重退化抛物方程组解的整体存在与爆破,通过构造自相似的上下解,得到了整体存在曲线.借助一些新的结果,获得了Fujita临界指数.其中一个有趣的现象是:整体存在曲线和Fujita临界曲线分别是由一个矩阵和线性方程组来决定.     

Nonlinear inversion of ultrasonic guided waves for in vivo evaluation of cortical bone properties

Ultrasonic guided waves (UGWs), which propagate throughout the entire thickness of cortical bone, are attractive for the early diagnosis of osteoporosis. However, this is challenging due to the impact of soft tissue and the inherent difficulties related to multiparametric inversion of cortical bone quality factors, such as cortical thickness and bulk wave velocity. Therefore, in this research, a UGW-based multi-parameter inversion algorithm is developed to predict strength-related factors. In simulation, a free plate (cortical bone) and a bilayer plate (soft tissue and cortical bone) are used to validate the proposed method. The inversed cortical thickness (CTh), longitudinal velocity (V_L) and transverse velocity (V_T) are in accordance with the true values. Then four bovine cortical bone plates were used in in vitro experiments. Compared with the reference values, the relative errors for cortical thickness were 3.96%, 0.83%, 2.87%, and 4.25%, respectively. In the in vivo measurements, UGWs are collected from the tibias of 10 volunteers. The theoretical dispersion curves depicted by the estimated parameters (V_T, V_L, CTh) match well with the extracted experimental ones. In comparison with dual-energy x-ray absorptiometry, our results show that the estimated transverse velocity and cortical thickness are highly sensitive to osteoporosis. Therefore, these two parameters (CTh and V_T) of long bones have potential to be used for diagnosis of bone status in clinical applications.