起伏海面环境下水声信道特性及估计方法*

为了研究随机起伏海面的散射对水声信道特性的影响,以及不同信道估计算法在该信道条件下的性能问题,该文从利用Monte Carlo方法模拟的基于PM海浪谱的二维随机起伏海面出发,提取了与风向不同夹角的一维随机起伏海面,并利用Kirchhoff近似法计算了其散射强度。将得到的声散射特性和海面环境信息融入到信道模型中,建立了基于PM海浪谱的随机信道模型。通过仿真分析了不同风速、传播距离条件下,声波的传播损失和信道冲激响应变化,并利用l_0-最小均方误差法、匹配追踪算法和正交匹配追踪算法进行了信道估计。给出了不同环境条件下,三种算法估计性能的综合比较结果,验证了算法针对该信道模型的有效性。     

Effects of red blood cell aggregates dissociation on the estimation of ultrasound speckle image velocimetry

Ultrasound speckle image of blood is mainly attributed by red blood cells (RBCs) which tend to form RBC aggregates. RBC aggregates are separated into individual cells when the shear force is over a certain value. The dissociation of RBC aggregates has an influence on the performance of ultrasound speckle image velocimetry (SIV) technique in which a cross-correlation algorithm is applied to the speckle images to get the velocity field information. The present study aims to investigate the effect of the dissociation of RBC aggregates on the estimation quality of SIV technique. Ultrasound B-mode images were captured from the porcine blood circulating in a mock-up flow loop with varying flow rate. To verify the measurement performance of SIV technique, the centerline velocity measured by the SIV technique was compared with that measured by Doppler spectrograms. The dissociation of RBC aggregates was estimated by using decorrelation of speckle patterns in which the subsequent window was shifted as much as the speckle displacement to compensate decorrelation caused by in-plane loss of speckle patterns. The decorrelation of speckles is considerably increased according to shear rate. Its variations are different along the radial direction. Because the dissociation of RBC aggregates changes ultrasound speckles, the estimation quality of SIV technique is significantly correlated with the decorrelation of speckles. This degradation of measurement quality may be improved by increasing the data acquisition rate. This study would be useful for simultaneous measurement of hemodynamic and hemorheological information of blood flows using only speckle images.     

脉冲电子束对材料破坏效应的数值研究

利用ＭｏｎｔｅＣａｒｌｏ方法计算了电子束的能量沉积，用流体动力学方程计算了热击波的传播，所用物态方程为ＧＲＡＹ三相物态方程。全部计算均由程序ＤＲＡＭ１Ｄ完成。讨论了波的传播规律，并给出了铝材迎光面反冲速度峰值及比质量亏损与电子束的入射通量和能量之间的关系。     

Inverse problems in cancellous bone: estimation of the ultrasonic properties of fast and slow waves using Bayesian probability theory

Quantitative ultrasonic characterization of cancellous bone can be complicated by artifacts introduced by analyzing acquired data consisting of two propagating waves (a fast wave and a slow wave) as if only one wave were present. Recovering the ultrasonic properties of overlapping fast and slow waves could therefore lead to enhancement of bone quality assessment. The current study uses Bayesian probability theory to estimate phase velocity and normalized broadband ultrasonic attenuation (nBUA) parameters in a model of fast and slow wave propagation. Calculations are carried out using Markov chain Monte Carlo with simulated annealing to approximate the marginal posterior probability densities for parameters in the model. The technique is applied to simulated data, to data acquired on two phantoms capable of generating two waves in acquired signals, and to data acquired on a human femur condyle specimen. The models are in good agreement with both the simulated and experimental data, and the values of the estimated ultrasonic parameters fall within expected ranges.     

Theoretical Analysis of Shear Wave Interference Patterns by Means of Dynamic Acoustic Radiation Forces

Acoustic radiation forces associated with high intensity focused ultrasound stimulate shear wave propagation allowing shear wave speed and shear viscosity estimation of tissue structures. As wave speeds are meters per second, real time displacement tracking over an extend field-of-view using ultrasound is problematic due to very high frame rate requirements. However, two spatially separated dynamic external sources can stimulate shear wave motion leading to shear wave interference patterns. Advantages are shear waves can be imaged at lower frame rates and local interference pattern spatial properties reflect tissue's viscoelastic properties. Here a theoretical analysis of shear wave interference patterns by means of dynamic acoustic radiation forces is detailed. Using a viscoelastic Green's function analysis, tissue motion due to a pair of focused ultrasound beams and associated radiation forces are presented. Overall, this paper theoretically demonstrates shear wave interference patterns can be stimulated using dynamic acoustic radiation forces and tracked using conventional ultrasound imaging.     

Speckle intensity images of target based on Monte Carlo method

Speckle intensity in the detector plane is deduced in the free-space optical system and imaging system based on Van Cittert-Zemike theorem. The speckle intensity images of plane target and conical target are obtained by using the Monte Carlo method and measured experimentally. The results show that when the range extent of target is smaller, the speckle size along the same direction become longer, and the speckle size increase with increasing incident light wavelengths. The speckle size increases and the speckle intensity images of target is closer to the actual object when the aperture scale augments. These findings are useful to access the target information by speckle in laser radar systems.     

Deterministic phase retrieval from diffracted intensities speckle fields

We present a model where the phase of a wave front is obtained by processing the intensity of speckle field pattern. By using the light propagation equation and knowing the intensity of the wavefront it is possible to build a system of linear equations which can be solved to give the wanted phase information. Simulations are presented and the convergence of the solution is discussed.     

Beam and phase effects in angle-beam-through-transmission method of ultrasonic velocity measurement

The accuracy of a plane wave approximation for phase velocity measurements in isotropic and anisotropic material using the angle-beam-through-transmission method has been investigated numerically and experimentally. In this method the velocity is measured in different propagation directions as a function of incidence angle. The effect of two factors on the measurement accuracy have been discussed: intrinsic phase shift of the transmitted signal through a fluid-solid interface and beam diffraction due to the finite beam size of receiver and transmitter. It is shown that the interface-induced phase shift can introduce an error in time delay measurements of the shear wave after the first critical angle and that this time delay error can be accurately corrected for. Numerical results obtained by a time-domain beam model show that except at the critical angles, the finite width of the transmitter and receiver only affects the amplitudes of the transmitted signals and has almost no effect on the measured zero-cross time delay; therefore the plane wave approximation for obtaining phase velocity from the measured time delay data by this method and the plane wave interface-induced phase correction are fully applicable.     

Speckle interferometry and Doppler diagnostics of microflows of scattering liquids

A relation for the correlation function of intensity fluctuations of the speckle field formed upon scattering of focused light radiation in blood flows in microvessels is obtained. The width of the intensity-fluctuation spectrum is shown to depend not only on the velocity of the flow, but also on its scattering characteristics. Simulation of the processes of scattering of light beams in random flows is performed by the Monte Carlo method. The expression for the first spectral moment of the Doppler signal is derived both on the basis of the theory of radiative-energy transfer and within the framework of the speckle-interference approach; comparison of the obtained results is carried out.     

Diffusion in a 2D bond percolation model. A Monte Carlo simulation

A Monte Carlo computer simulation of a simple stochastic hopping model is presented. The velocity autocorrelation function and the diffusion coefficient are analyzed and compared with the results obtained by using kinetic theory. The theory and the simulation agree within the statistical error.     

受力腓肠肌剪切波速度的超声检测

肌肉组织受力时的弹性特征变化,是其功能特性、健康状态的重要评估参量.该文利用在离体肌肉组织表面施加低频振动引起组织发生形变,在组织内部产生剪切波的方法,经过快速超声成像,通过剪切波图像估计剪切波的传播速度.研究结果表明顺肌纤维方向传播的剪切波速度快于垂直肌纤维方向传播的剪切波速度,剪切波的传播速度随受力的增加而变大,且...     

双横波声速法检测单向受压混凝土构件永存应力

该文力求寻找一种高效准确检测混凝土构件永存应力的方法.基于声弹性理论提出了一种双横波声速法检测单向受压混凝土构件永存应力的方法,该方法通过测试受力构件第一波速和第二波速,以第一波速和第二波速的平方差为基础构建综合声学参数来检测单向受压混凝土构件永存应力.所述第一波速为声波传播方向与应力方向垂直、质点振动方向与应力方向平...     

Interferometric evaluation of angular displacements using phase retrieval

Phase retrieval is carried out using sequential intensity measurements of a volume speckle field and a wave propagation equation. Retrieved phases and phase subtraction facilitate the analysis of wavefronts before and after undergoing a small rotation. Angular displacement between incident planar wavefronts is determined from the unwrapped phase difference, phase diffuser aperture diameter, and the light source wavelength. Numerical simulations confirm the experimental results.     

A Metropolis Monte Carlo Implementation of Bayesian Time-Domain Parameter Estimation: Application to Coupling Constant Estimation from Antiphase Multiplets

The Bayesian perspective on statistics asserts that it makes sense to speak of a probability of an unknown parameter having a particular value. Given a model for an observed, noise-corrupted signal, we may use Bayesian methods to estimate not only the most probable value for each parameter but also their distributions. We present an implementation of the Bayesian parameter estimation formalism developed by G. L. Bretthorst (1990,J. Magn. Reson.88, 533) using the Metropolis Monte Carlo sampling algorithm to perform the parameter and error estimation. This allows us to make very few assumptions about the shape of the posterior distribution, and allows the easy introduction of prior knowledge about constraints among the model parameters. We present evidence that the error estimates obtained in this manner are realistic, and that the Monte Carlo approach can be used to accurately estimate coupling constants from antiphase doublets in synthetic and experimental data.     

纺织纤维光散射特性的模拟与实验

纺织纤维的光散射特性在纺织材料微观结构、光学性质以及无损检测中至关重要。利用纺织纤维各向异性的结构特点改进现有的角谱展开法,得到纺织纤维对倾斜入射平面波的散射以及表征其光散射特性的Mueller矩阵。采用蒙特卡罗算法模拟了平行棉纤维束对倾斜入射偏振光的多次散射,并对其计算结果进行了实验验证。结果表明理论计算和实验测量的散射光斑在形状特征、光强分布上均一致,从而验证了蒙特卡洛方法模拟纺织纤维光散射的正确性。结果同时也表明蒙特卡罗方法在纺织材料光传播特性的理论研究中具有重要的指导意义。     

海底沉积物压缩波速度与切变波速度的关系

基于连续介质假设,根据无吸收各向同性弹性介质通用方程分析沉积物声波传播关系,提出应用弹性结构分布因子表达的声速通用模型(GMSS,General Model of Sound Speed)分析海底沉积物的声速特性;通过研究Willey时间平均模型、Wood方程、Gassmann方程、Buckingham模型、Biot-Stoll模型和EDFM模型,可以表述成GMSS通用模型中的弹性结构分布因子的具体表达形式,得出GMSS通用模型在解释压缩波速度和切变波速度特性上具有一致性的特点。GMSS通用模型具有弹性结构分布因子、孔隙度、孔隙海水的等效密度和等效弹性模量、固相颗粒的等效密度、固相颗粒的等效体积弹性模量和等效切变弹性模量共7个参数,为研究海底沉积物压缩波和切变波速度提供了一种模型简单、参数少、通用性强的方法。但也需要从物理结构上以及应力应变关系上开展更为深入的分析和探寻GMSS模型的物理意义和参数测量的方法。     

A ray-optic, Monte Carlo, description of a Gaussian beam waist – applied to reverse saturable absorption

The intensity distribution in a Gaussian beam is represented by a bundle of rays in which each ray has a random pointing error. The ray bundle is allowed to pass through a simple lens using the usual geometrical optics formulae and it is shown that the position of the focused beam waist agrees exactly with the usual Gaussian mode analysis. The interaction of the beam waist with a reverse saturable absorber (silicon napthalocyanaine) is modelled using a Monte Carlo method and a null absorption coefficient is introduced to increase the speed of the calculation. An f/5 geometry is considered and it is shown that the results from the ray model agree well with those obtained from standard split-step beam propagation techniques. Received: 23 July 1999/Revised version: 28 September 1999/Published online: 21 January 2000     

基于EKF的水下LD通信精对准控制算法

高速激光通信中接收机与光斑中心很难处于精对准状态,导致水下光通信链路难以稳定建立.首先采用蒙特卡洛仿真统计法分析激光光子在海水中传输的接收光强分布规律,再通过实验对接收端的光斑图像进行采样分析,利用曲线拟合得到接收器位置与接收光强的关系.仿真与实验结果表明:光束经过25 m的水下传输,接收光强分布仍近似为高斯分布.采用非线性估计算法(扩展卡尔曼滤波)与基本状态控制反馈理论,根据接收光强度估计接收器当前位置与最大光强处的距离,通过反馈算法实现接收端与光斑中心的主动跟踪对准.算法仿真结果显示,接收端对准误差在2 mm以下,稳定后接收效率超过98%.     

The estimation of lower refractivity uncertainty from radar sea clutter using the Bayesian-MCMC method

The estimation of lower atmospheric refractivity from radar sea clutter(RFC) is a complicated nonlinear optimization problem.This paper deals with the RFC problem in a Bayesian framework.It uses the unbiased Markov Chain Monte Carlo(MCMC) sampling technique,which can provide accurate posterior probability distributions of the estimated refractivity parameters by using an electromagnetic split-step fast Fourier transform terrain parabolic equation propagation model within a Bayesian inversion framework.In contrast to the global optimization algorithm,the Bayesian-MCMC can obtain not only the approximate solutions,but also the probability distributions of the solutions,that is,uncertainty analyses of solutions.The Bayesian-MCMC algorithm is implemented on the simulation radar sea-clutter data and the real radar seaclutter data.Reference data are assumed to be simulation data and refractivity profiles are obtained using a helicopter.The inversion algorithm is assessed(i) by comparing the estimated refractivity profiles from the assumed simulation and the helicopter sounding data;(ii) the one-dimensional(1D) and two-dimensional(2D) posterior probability distribution of solutions.     

新型二次电子倍增阴极的蒙特卡罗模拟研究

提出了一种新型二次电子倍增阴极强流二极管,并对其进行了动力学理论简化模型和蒙特卡罗数值模拟的对比验证研究。首先,基于设计结构原型,根据二次电子发射特性进行合理简化,建立了动力学模型,获得了电子速度、位移以及渡越时间的解析结果,并结合Vaughan的二次电子产额模型,确定了该新型二次电子倍增阴极强流二极管的理论工作区间;其次,理论分析了施加径向电场的重要意义,并给出了二次电子运动特征参数(最大位移、渡越时间、碰撞能量等)的理论预估结果;最后,对该新型二次电子倍增阴极强流二极管进行了蒙特卡罗模拟研究,获得了电子的运动轨迹、碰撞能量以及二次电子倍增工作区间等物理图像,并将蒙特卡罗数值模拟结果与理论结果进行了比对,两者吻合程度较好,对可能的误差来源进行了分析讨论。理论和模拟结果表明：新型二次电子倍增阴极强流二极管概念可行,工作区间内通过调整施加电场与磁场幅值,可有效达到电子运动状态可控的目标。另外,理论粗估了二次电子倍增饱和条件下的阴极发射电流密度,结果表明：发射电流密度可达kA/cm²水平,具备强流发射特性;增加外加径向场强幅值可有效提升发射电流密度。最后,对该新型二次电子倍增阴极设计步骤和依据进行了讨论。