Time reversal and the inverse filter

To focus ultrasonic waves in an unknown inhomogeneous medium using a phased array, one has to calculate the optimal set of signals to be applied on the transducers of the array. In the case of time-reversal mirrors, one assumes that a source is available at the focus, providing the Green's function of this point. In this paper, the robustness of this time-reversal method is investigated when loss of information breaks the time-reversal invariance. It arises in dissipative media or when the field radiated by the source is not entirely measured by the limited aperture of a time-reversal mirror. However, in both cases, linearity and reciprocity relations ensure time reversal to achieve a spatiotemporal matched filtering. Nevertheless, though it provides robustness to this method, no constraints are imposed on the field out of the focus and sidelobes may appear. Another approach consists of measuring the Green's functions associated to the focus but also to neighboring points. Thus, the whole information characterizing the medium is known and the inverse source problem can be solved. A matrix formalism of the propagation operator is introduced to compare the time-reversal and inverse filter techniques. Moreover, experiments investigated in various media are presented to illustrate this comparison.     

Optimal adaptive focusing through heterogeneous media with the minimally invasive inverse filter

The inverse filter is a technique used to adaptively focus waves through heterogeneous media. It is based on the inversion of the Green's functions matrix between the M transducers of a focusing array and N control points in the focal area. The inverse filter minimizes the pressure deposited around the focal point. However it is highly invasive, requiring the presence of N transducers or hydrophones in the focal area at the control points' locations to measure the Green's functions. This paper presents a way of reaching the inverse filter's focusing quality with a minimally invasive setup: only one transducer (at the desired focal location) is needed. This minimally invasive inverse filter takes advantage of the fact all the information about the propagation medium can be retrieved from the signals backscattered by the medium towards the focusing array, if the propagation medium is lossless. A numerical simulation is performed to test this minimally invasive inverse filter through a scattering, lossless medium. The focusing quality equals the conventional, highly invasive inverse filter's. The average spatial and temporal contrast is increased by up to 10 dB compared to the time reversal focusing.     

AC shot noise through a vibrating quantum dot

In this paper we investigate the joint effects of the electron-phonon interaction and an external alternating (ac) gate voltage on the spectral density of shot noise through a vibrating quantum dot system, by applying the Lang-Firsov canonical transformation and the Keldysh nonequilibrium Green's function (NGF) technique. We find that the effects of the electron-phonon and electron-photon interaction on the differential shot noise are different. The main resonant peak of the differential shot noise is decreased only when a time-dependent gate voltage is imposed on quantum dot. With the ac field amplitude increasing, the main resonant peak of the differential shot noise decreases. The Fano factor of the system, which exhibits the deviation of shot noise from the Schottky formula, is also studied. Super-Poissonian shot noise appears due to the photon absorption (emission) processes in the low bias voltage region.     

基于时间反演技术的电磁器件端口场与内部场转换方法

随着电磁器件的集成化,器件搭载的模块、实现的功能愈发多样.各模块间的耦合难以忽略,设计难度陡然增加,传统设计方法逐渐力不从心,迫切需要寻找一种新的电磁综合设计方法.本文利用时间反演电磁波的时空同步聚焦特性,探索了将时间反演技术应用于器件设计的可能性.首先,基于通用的器件逆设计流程,利用时间反演技术、并矢格林函数及电磁学的基本原理,提出了将器件端口场分布转换为内部场分布的方法,并证明由端口期望场的时间反演场在空间某一位置获得的连续等效源的共轭分布可在端口处产生与期望场接近的场分布.且在单点频逆设计过程中,只需知道端口电场或磁场的切向分量即可完成端口场与内部场的转换.同时,借助格林函数的互易性对本文所提理论做适当变换后,进行数值仿真验证,分析讨论了不同初始信息条件下该方法的适用性.仿真结果与理论相符,证明了理论的正确性,为将时间反演技术应用于电磁器件的逆设计提供了可能.     

Real time inverse filter focusing through iterative time reversal

In order to achieve an optimal focusing through heterogeneous media we need to build the inverse filter of the propagation operator. Time reversal is an easy and robust way to achieve such an inverse filter in nondissipative media. However, as soon as losses appear in the medium, time reversal is not equivalent to the inverse filter anymore. Consequently, it does not produce the optimal focusing and beam degradations may appear. In such cases, we showed in previous works that the optimal focusing can be recovered by using the so-called spatiotemporal inverse filter technique. This process requires the presence of a complete set of receivers inside the medium. It allows one to reach the optimal focusing even in extreme situations such as ultrasonic focusing through human skull or audible sound focusing in strongly reverberant rooms. But, this technique is time consuming and implied fastidious numerical calculations. In this paper we propose a new way to process this inverse filter focusing technique in real time and without any calculation. The new process is based on iterative time reversal process. Contrary to the classical inverse filter technique, this iteration does not require any computation and achieves the inverse filter in an experimental way using wave propagation instead of computational power. The convergence from time reversal to inverse filter during the iterative process is theoretically explained. Finally, the feasibility of this iterative technique is experimentally demonstrated for ultrasound applications.     

Equivalence of Euclidean and Wightman field theories

A new inversion formula for the Laplace transformation of tempered distributions with supports in the closed positive semiaxis is obtained. The inverse Laplace transform of a tempered distribution is defined by means of a limit of a special distribution constructed from this distribution. The weak spectral condition on the Euclidean Green's functions implies that some of the limits needed for the inversion formula exist for any Euclidean Green's function with an even number of variables. We then prove that the initial Osterwalder-Schrader axioms [1] and the weak spectral condition are equivalent with the Wightman axioms.The research described in this publication was made possible in part by Grant No. 93-011-147 from the Russian Foundation for Basic Research     

Trading detection for resolution in active sonar receivers

This paper proposes an active sonar receivers that offers a smooth trade-off between detection and resolution. A matched filter is the optimal detector of known signals in white Gaussian noise but may fail to resolve the targets if the time separation of targets is less than the mainlobe width of the autocorrelation function of the transmitted signal. An inverse filter achieves optimal resolution performance for multiple targets in the absence of noise, but amplifies the noise outside the signal bandwidth in a manner that makes it impractical in many realistic scenarios. The proposed active sonar receiver, the variable resolution and detection receiver (VRDR) combines the matched and inverse filter properties to achieve a smooth trade-off between detection and resolution. Simulated receiver operating characteristics demonstrate that for a range of dipole sonar targets, the performance of the VRDR is superior to the matched and inverse filter, as well as another previously proposed bandlimited inverse filter.     

激光相干场成像散斑噪声复合去噪方法

噪声是影响激光相干场高分辨成像系统像质的重要因素,激光相干场成像系统既受背景光加性噪声影响,又受激光乘性散斑噪声影响.为解决激光相干场成像系统受激光乘性散斑噪声和背景光加性噪声叠加引起的成像像质退化效应问题,从噪声抑制角度提高激光相干场系统高分辨成像像质,研究建立了激光散斑乘性噪声和背景光加性噪声对大气下行链路激光回波场信号影响干扰模型,并基于该模型提出了一种基于同态滤波和稀疏基追踪级联复合去噪算法.首先基于同态滤波理论将激光乘性散斑噪声转化为加性噪声,再由高通滤波器滤除散斑噪声,最后采用基追踪稀疏理论方法抑制背景光等加性噪声对像质的影响.研究表明,较现有单一去噪方法,该级联复合去噪方法可一次性消除激光乘性散斑噪声和背景加性噪声两种不同性质的噪声,有效改善了激光相干场成像质量.     

Optimisation using measured Green's function for improving spatial coherence in acoustic measurements

Aberrating materials can degrade acoustic measurements by distorting the acoustic wavefront and causing acoustic speckle (as opposed to speckle noise which is a manifestation of coherent backscatter). The amplitude and phase fluctuations associated with acoustic speckle can introduce considerable measurement uncertainty which is difficult to deal with. This paper demonstrates a new technique which optimises the spatial distribution of the generation of the ultrasound to compensate for the aberration. This technique uses experimentally measured Green's functions to allow the calculation of the field resulting from the generation wavefront during optimisation. The technique is used to improve the accuracy of velocity measurements in a steel sample using 82 MHz SAW waves. This is achieved by optimising for improved spatial coherence in the measurement region which suppresses the speckle noise. Experimental evidence of acoustic aberration arising from grain structure is shown for steel and aluminium and the measured Green's function optimisation technique is shown to overcome the resulting acoustic speckle. The technique was performed using the Adaptive Optical Scanning Acoustic Microscope (AOSAM) at Nottingham University, UK.     

Continuation of acoustic near-fields

This paper deals with the analytic continuation of a coherent pressure field specified on a finite sheet located close to and conformal to the surface of a vibrator. This analytic continuation is an extension or extrapolation of the given (measured) field into a region outside and tangential to the original finite sheet, and is based on the Green's function (the transfer function) relating acoustic quantities on the two conformal surfaces. The continuation of the measured pressure field is an inverse problem that requires the use or regularization theory, especially when noise is present in the data. An iteration algorithm is presented that is successful in continuing the pressure field into the tangential sheet. The results are accurate close to the original boundary and taper (decay) toward zero with distance away from it. The algorithm is tested on numerical and experimental data from a point-driven rectangular plate. Results show the successful extrapolation (continuation) of this data into an area nearly double that of the original pressure field. This algorithm is not limited to planar surfaces and can be applied to arbitrarily shaped surfaces.     

On the two-point cross-correlation function of anisotropic, spatially homogeneous ambient noise in the ocean and its relationship to the Green's function

It is well established that the free-space Green's function can be recovered from the two-point cross-correlation function of a random noise field if the noise is white and isotropic. Ambient noise in the ocean rarely satisfies either of these conditions. However, a non-uniform spectrum could be pre-whitened by the application of a suitable filter but anisotropy cannot be so readily eliminated. To investigate the effects of vertical anisotropy, three azimuthally uniform, spatially homogeneous noise fields are analyzed, two of which are idealized, while the third is representative of ambient noise in the deep ocean. In each case, the coherence function, the cross-correlation function, and the derivative of the latter with respect to the correlation delay, are derived for vertical and horizontal alignments of the sensor pair. With vertical sensors, any step-function discontinuity in the directional density function is mapped into a delta function at an appropriate time delay in the derivative (with respect to time delay) of the cross-correlation function. No such mapping occurs with horizontal sensors. In this case, only horizontally traveling noise can generate delta functions in the derivative of the cross-correlation function, and these always appear at the retarded time on either side of the origin.     

Acoustical inverse problems regularization: Direct definition of filter factors using Signal-to-Noise Ratio

Acoustic imaging aims at localization and characterization of sound sources using microphone arrays. In this paper a new regularization method for acoustic imaging by inverse approach is proposed. The method first relies on the singular value decomposition of the plant matrix and on the projection of the measured data on the corresponding singular vectors. In place of regularization using classical methods such as truncated singular value decomposition and Tikhonov regularization, the proposed method involves the direct definition of the filter factors on the basis of a thresholding operation, defined from the estimated measurement noise. The thresholding operation is achieved using modified filter functions. The originality of the approach is to propose the definition of a filter factor which provides more damping to the singular components dominated by noise than that given by the Tikhonov filter. This has the advantage of potentially simplifying the selection of the best regularization amount in inverse problems. Theoretical results show that this method is comparatively more accurate than Tikhonov regularization and truncated singular value decomposition.     

格林互易定理在求解均匀带电圆环的静电势分布中的应用

对于静电场,由普遍的格林互易定理得出特殊情形下的格林互易定理.应用该定理推导出均匀带电圆环的数种在形式上互不相同的静电势分布解式.     

Spatial coherence and cross correlation of three-dimensional ambient noise fields in the ocean

Ambient acoustic noise fields in the ocean are generally three dimensional in that they exhibit vertical and horizontal directivity. A model of spatially homogeneous noise is introduced in which the directionality is treated as separable, that is, the overall directionality of the field is the product of the individual directivities in the horizontal and vertical. A uni-modal von Mises circular distribution from directional statistics is taken to represent the noise in the horizontal, whilst the vertical component is consistent with a surface distribution of vertical dipoles. An analysis of the coherence and cross correlation of the noise at two horizontally aligned sensors is developed. The coherence function involves a single integral over finite limits, whilst the cross-correlation function, derived on the assumption that the noise has been pre-whitened, is given by an integral with limits that depend on the correlation delay time. Although the cross-correlation function does not exhibit delta functions that could be identified with the Green's function for propagation between the two sensors in the field, it does drop abruptly to zero at numerical time delays equal to the travel time between the sensors. Hence the noise could be used to recover the sound speed in the medium.     

数字全息周期像的产生机理及在抑制零级衍射上的应用

在数字全息成像中, 利用CCD的RGB模式采样全息图时, 全息重构像会出现特定的周期性分布. 本文从理论和实验上详细研究了这种周期像产生的机理、分布特性和应用. 研究结果显示, 由于CCD的光谱滤镜会使全息图的RGB三个单色采样阵列出现部分像素信号的缺失, 因此, 需要通过特定的demosaicing数学算法对缺失的像素信号进行重建以形成完整的单色采样阵列, 这是数字全息再现像周期分布产生的根源. 而基于demosaicing算法的采样阵列重建会在全息图频谱中引入调制函数, 导致物体再现像和零级衍射斑的周期分布差异. 本文揭示了全息图的RGB采样、demosaicing算法与全息重构像周期性之间的内在关联. 最后, 讨论了结合空间移位和图像形态学技术, 利用重构像的周期性抑制零级衍射斑的应用. 所有理论与实验研究结果完全一致. 关键词： 数字全息 图像周期性 零级斑抑制     

Stabilization of Pattern in Complex Ginzburg-Landau Equation with Spatial Perturbation Scheme

In this paper, a new scheme of spatial perturbation is proposed to stabilize the pattern in the oscillatory media, which could be described with the complex Ginzburg-Landau equation. The numerical simulation results confirm that the spiral wave, antispiral wave and spatiotemporal chaos in the complex Ginzburg-Landau equation could be suppressed, and the scheme is also discussed with the conservative field theory. Furthermore, the spatiotemporal noise is also introduced into the whole media, it just confirms that it is robust to the spatiotemporal noise.     

Optical responses of dilute anisotropic composites：numerical calculations via Green‘s function formalism

We investigate the linear and nonlinear optical responses of dilute anisotropic networks using Green's function formalism [Gu Y et al 1999 Phys. Rev. B 59 12847]. For different applied fields, numerical calculations indicate that a large third-order nonlinear enhancement and a broad infrared absorption arise from the geometric anisotropy. We also show the overlap and separation between the absorption peak and nonlinear enhancement peak when the applied field is parallel and perpendicular to the anisotropy, respectively. The results can be understood in terms of the inverse participation ratios with q=2 and the spectral distribution of optical responses.     

Reconstruction of dynamic displacement and velocity from measured accelerations using the variational statement of an inverse problem

This paper presents two types of finite impulse response (FIR) filters to reconstruct dynamic displacement induced by structural vibration from measured acceleration. The governing equation for the reconstruction is derived by taking the variation of a minimization problem, which defines an inverse problem on displacement. A regularization function for overcoming the ill-posedness of the inverse problem is included in the minimization problem. The governing equation of the inverse problem becomes the same type of differential equation as that of a beam on an elastic foundation. The conventional FIR (CFIR) filter directly approximates the transfer function of the governing equation, while the FEM-based FIR (FFIR) filter is formulated by the discretization of the minimization problem with the finite element method. For the finite element discretization, the Hermitian shape function is utilized. The proposed FFIR filter is capable of reconstructing displacement and velocity simultaneously. The fundamental characteristics of the proposed filters are investigated in the frequency domain using the transfer and accuracy functions. It is shown that the proposed FIR filters suppress low frequency noise components in measured accelerations effectively, and reconstruct physically meaningful displacement accurately. The validity of the proposed filters is demonstrated through a numerical simulation study, a field experiment and an evaluation of flutter derivatives using measurements taken from a wind tunnel test.     

利用环境噪声互相关实现散射体无源成像

最近研究表明利用环境噪声的互相关可以恢复两点之间的时域格林函数(声脉冲响应),这一原理在文献中被称为格林函数恢复。基于此原理,通过对多个传声器所接收的环境噪声进行互相关处理,获取与散射体相关联的散射波的到达时延信息,结合最小二乘反演算法和改进克希霍夫移位算法,分别获得道路交通噪声场中石柱以及海浪噪声场中塑料桶的空间位置,且其估计结果与实际测量相一致。实验结果表明将环境噪声作为探测信号进行散射体无源成像是可行的。这为设计室内无源声监测系统以及通过海洋环境噪声实现对水中静默目标成像提供了新思路和有益参考。      

Detection of the spatiotemporal field of a single-shot terahertz pulse based on spectral holography

According to electro-optical sampling theory, we propose a new method to detect the spatiotemporal field of a single- shot terahertz pulse by spectral holography for the first time. The single-shot terahertz pulse is coupled into a broadened chirped femtosecond pulse according to electro-optical sampling theory in the detecting system. Then the reference wave and the signal wave are split by Dammann grating and spread into the interference band-pass filter. The filtered sub-waves are at different central-frequencies because of the different incident angles. These sub-waves at different central-frequencies interfere to form sub-holograms, which are recorded in a single frame of a charge coupled device （CCD）. The sub-holograms are numerically processed, and the spatiotemporal field distribution of the original terahertz pulse is reconstructed. The computer simulations verify the feasibility of the proposed method.