The radial Marchenko-Blažek integral equation for complete spectra

The same simple and standard way, by means of which an inverse problem for scattering of spinless particles by central potentials is solved in the Gelfand-Levitan method, is applied to the Marchenko method for general angular momenta including the bound states. We first derive an integral equation for the kernel with a triangularity property, which relates it to a potential and then the other one, which connects the kernel with spectral data. A solution corresponding to the general Yukawa potential is found and some formulae are checked by solving the problem of the restrained phase equivalent potentials.     

RIEMANN-HILBERT PROBLEM OF THE QUANTUM NONLINEAR SCHRÖDINGER MODEL

In this paper we present a representation of quantum Riemann-Hilbert problem of the nonlinear Schrödinger model, from which both the results of direct problem of the quantum inverse scattering method and the quantum Gelfand-Levitan equation can be given. It is pointed out that this method may be extended to any other quantum AKNS-type model.     

A new approach to the inverse scattering and spectral problems for the Sturm-Liouville equation

New proofs of the known uniqueness theorems for the one-dimensional inverse spectral and scattering problems are given. Proof of the invertibility of all of the steps in the inversion procedures of Gelfand-Levitan and Marchenko is given. The proposed method of investigation yields some new results, for example, a Marchenko-type equation at x = 0 which holds on the whole axis, rather than on a half-axis, as usual for the scattering theory on half-axis. It also yields a new method, shorter and simpler than earlier published, for proving that the potential in the class L_1,1, obtained by the Marchenko reconstruction procedure, generates the scattering data from which it was reconstructed.     

Exact solution of the inverse acoustic scattering problem for the one-dimensional density profile in cylindrical geometry

The inhomogeneous acoustic wave equation in cylindrically symmetric geometry, is reduced to a one-dimensional Schrödinger equation for time harmonic waves. The scattering potential so formed is uniquely reconstructed by the Gelfand-Levitan theory. In physical terms, the scattering potential is a function of the impedance profile, which is then recovered. If the velocity profile is known (or constant) the density profile, as a function of the radial coordinate, can be deduced. A practical experiment is envisaged whereby the back-scattered time sequence is measured in a pulse-echo mode, and the data are used to reconstruct the density profile.     

对称性及多群中子扩散方程数值解

在多群中子扩散方程解析解的基础上,利用方程及求解域的对称性建立了新的数值求解中子扩散方程的理论模型.该模型显著的优点是适用于各种对称区域(二维、三维区域)尤其是非正方形区域中子扩散方程的求解,它彻底避免了常规节块法应用于非正方形几何时所出现的奇异性问题,且所得的解在求解域内任意点上均满足扩散方程.以二、三维六角形几何为例建立了数学模型,并用基准问题校核了模型的正确性. 关键词： 中子扩散方程 对称群 数值解 解析     

基于约束优化的多光谱辐射真温反演算法

多光谱辐射测温是通过测量待测物某点的多个光谱辐射强度信息,通过普朗克公式反演获得真实温度。但是,通过普朗克公式获得的多光谱辐射测温方程组,是欠定方程组,即N个方程,N+1个未知数(N个未知的光谱发射率ε_λi和1个待求真温T)。目前,多采用事先假设一组发射率模型(发射率-波长或发射率-温度模型),假设模型与实际情况如果相符,则反演结果能够满足要求,如果假设模型与实际情况不符,则反演结果误差很大。但是,发射率模型受温度、表面状态、波长等诸多因素影响,难以事先确定发射率模型。因此受未知光谱发射率的制约一直是多光谱辐射测温理论面临的主要障碍,能否在无需任何光谱发射率假设模型的情况下,实现真温和光谱发射率的直接反演一直是多光谱辐射测温理论研究的热点和难点。通过对参考温度模型的分析表明,多光谱辐射测温反演过程的实质是寻找一组光谱发射率,使得每个通道方程解得的真温都相同,如不相同则继续寻找合适的光谱发射率,直到每个通道解得的真温都相等。为此,提出将多光谱辐射测温参考温度模型的求解过程转换为约束优化问题,即在光谱发射率0≤ε_λi≤1的约束条件下,通过梯度投影算法不断寻找光谱发射率,带入多光谱辐射测温参考温度模型方程组后,计算温度反演值的方差,直到每个光谱通道方程获得的温度值应该近似相等,此时各个光谱通道的温度反演值方差最小,这样就把多光谱辐射真温和发射率的反演问题转换为约束优化问题。约束优化算法是解决这一类问题的主要方法,但为了满足Ax≥b的约束条件,将0≤ε_λi≤1分解为ε_λi≥0和-ε_λi≥-1的两个约束条件,从而满足了约束优化问题Ax≥b的约束条件。这样就可以通过约束优化算法在无需任何光谱发射率假设模型的条件下,直接求解真温和光谱发射率。实验采用六种不同光谱发射率分布模式(随波长递增、递减、凸波动、凹波动、“M”型波动、“W”型波动)的材料为研究对象,以验证新算法对不同材料光谱发射率分布反演的适应性,利用Matlab的minRosen函数,选择光谱发射率的初始值均为0.5(取中间值,提高计算效率)。针对六种不同光谱发射率模型的仿真结果表明,新算法无需任何有关发射率的先验知识,对不同发射率模型反演结果均表现较好,在真温1 800 K的情况下,绝对误差均小于20 K,相对误差均小于1.2%,新算法具有无需考虑任何光谱发射率先验知识、反演精度较高及适合于各种发射率模型等优点,进一步完善了多光谱辐射测温理论,在高温测量领域具有良好的应用前景。     

基于逆散射理论的地震波速度正则化反演

本文基于逆散射理论利用正则化有限差分对比源反演算法对地震波传播速度进行反演, 该方法是基于波动方程的频率域波形反演算法, 利用非线性共轭梯度法, 通过最小化目标优化函数不断迭代更新速度模型. 由于地球物理反演问题的病态性和不稳定性, 通过基于反演参数总变差的正则化处理, 使反演问题变为良性问题且算法具有较强的抗噪声干扰能力. 反演过程中使用了频率-空间域9点差分正演算子以及PML吸收边界条件. 与其他反演算法相比, 由于背景模型在反演迭代过程中保持不变, 可以避免在每次迭代过程中重新构造正演算子及矩阵分解等相关计算过程, 使得该算法非常适合于大规模三维反演计算. 此外, 本文采用基于MPI的并行计算, 进一步提高了反演计算的效率. 二维CSEG模型反演结果表明该方法可以反演得到高分辨率的地震波速度重建结果, 为地震勘探数据处理及解释提供准确的速度信息.     

Refractive index imaging from radiative transfer equation-based reconstruction algorithm: Fundamentals

We present the first reconstruction algorithm for refractive index imaging, which is based on the radiative transfer equation (RTE). An objective function is iteratively minimized to find a solution to the problem of inversion of the refractive index field. The function describes the discrepancies of the emerging light measurements on the surface of the sample to be probed with predicted data from the corresponding numerical model. The unknown refractive index field is updated within each reconstruction iteration according to a search direction on the index distribution given by the adjoint model to the RTE. In this paper, emphasis is placed on the theoretical aspects. Preliminary tests are demonstrated on generic phantoms.     

限流扩散理论的局限性及其改进

对几种典型的限流扩散理论的局限性与含时间吸收项的渐近扩散理论的长处作了分析和数字计算比较。结论是,理论上所导出的限流扩散方程,虽然形式上解决了中子迁移的限流问题,但由于理论推导中所作的假定带有一定的局限性,往往使限流过度;而带有时间吸收项的渐近扩散理论避免了上述的局限性,因而它对中子迁移问题的描述优于限流扩散理论。     

Extension of the electrokinematics theorem to the electromagnetic field and quantum mechanics

Summary A recent electrokinematics theorem leads to a general equation that, through an arbitrary irrotational fieldF, connects the motion of the electric-charge carriers, the internal potential and the dielectric properties of a physical system with its external currents, voltages and powers. It has been proved for quasi-electrostatic fields,i.e. when the vector potential may be disregarded, and on the basis of classical mechanics. Here the theorem is extended to any type of electromagnetic field and to quasi-relativistic quantum mechanics, in the case of many-particle systems for which, moreover, the probability current density is suitably computed. The new equation so obtained, throughF, connects the external currents again with the internal electric permittivity and the scalar potential, in the same way as in the preceding approach, and with the carrier velocity that, however, has to be computed according to quantum mechanics. Moreover, it contains two new contributions, one deriving from the vector potential and the other from a current density arising from the electron spin. By means of proper choices ofF, new expressions of the external currents of the system are determined as functions of the motion of its internal carriers. In particular, the electrokinematics theorem is exploited to compute the output current in two-terminal nanoelectronic devices in which, owing to the small sizes, quantum effects cannot be disregarded. Finally, such results, when they are applied to the double-barrier tunnelling structures, allow us to show the splitting of the electron pulse into two uncorrelated pulses, and as a consequence, to obtain a possible shot noise suppression, up to fifty per cent of the full shot noise.     

严谨而精密的创新——瑞利进行研究的特点

着重介绍了瑞利在测量氮气密度过程中发现惰性气体氩的研究和创新历程,以及瑞利对光学、声学等物理学的各个领域的主要贡献,指出正是求实的科学态度、锲而不舍的钻研精神和以创新为乐趣的事业追求是瑞利一生中取得丰硕的研究成果的主要原因。     

Axisymmetric wave propagation of thermoelastic transversely isotropic half-space under buried loading using potential functions

By virtue of a new scalar potential function and Hankel integral transforms, the wave propagation analysis of a thermoelastic transversely isotropic half-space is presented under buried loading and heat flux. The governing equations of the problem are the differential equations of motion and the energy equation of the coupled thermoelasticity theory. Using a scalar potential function, these coupled equations have been uncoupled and a six-order partial differential equation governing the potential function is received. The displacements, temperature, and stress components are obtained in terms of this potential function in cylindrical coordinate system. Applying the Hankel integral transform to suppress the radial variable, the governing equation for potential function is reduced to a six-order ordinary differential equation with respect to z. Solving that equation, the potential function and therefore displacements, temperature, and stresses are derived in the Hankel transformed domain for two regions. Using inversion of Hankel transform, these functions can be obtained in the real domain. The integrals of inversion Hankel transform are calculated numerically via Mathematica software. Our numerical results for displacement and temperature are calculated for surface excitations and compared with the results reported in the literature and a very good agreement is achieved.     

扩展卡尔曼滤波和不敏卡尔曼滤波在实时雷达回波反演大气波导中的应用

为了改善雷达回波反演大气波导(RFC)方面存在的单时次、单方位角反演的问题,提出利用扩展卡尔曼滤波和不敏卡尔曼滤波的反演算法对大气波导结构的多方位角实时跟踪反演. 在卡尔曼滤波方法中分别给出大气波导结构的参数化方程、观测方程、滤波算法的状态转移方程,最后导出滤波反演算法的迭代求解流程. 在大气波导结构不随时间变化和随时间变化的两种条件下,对扩展卡尔曼滤波和不敏卡尔曼滤波算法进行数值实验. 实验结果表明,不敏卡尔曼滤波更适用于RFC这高度非线性反演问题,它可能今后为大气波导结构多方位角实时跟踪反演的业务化运行提供理论基础与技术保证. 关键词： 大气波导 雷达回波 扩展卡尔曼滤波 不敏卡尔曼滤波     

Quantum tunneling: A general study in multi-dimensional potential barriers with and without dissipative coupling

Quantum tunneling is formulated in terms of the time evolution of a localized state and thus shown to be dependent upon the eigenspectrum of the system Hamiltonian. A number of exactly solvable models with local and non-local double-well potentials are discussed, and it is shown how, for local potentials, other solvable models can be generated by using Gelfand-Levitan and Darboux transformations. Tunneling in multi-dimensional potential barriers is investigated under semi-classical approximation by developing the method of asymptotic expansion of the wave function for large quantum numbers and the WKB approximation for separable systems. General expressions for the imaginary-time tunneling trajectory are obtained in both methods and specific applications are discussed. Approximation schemes for non-separable systems are also presented. A general study of dissipative multi-dimensional tunneling is carried out by using the Gisin equation, the Schrödinger-Langevin equation and the complex potential model. It is shown that, in general, different models of dissipation are not equivalent in the tunneling context. Using these models one can show (a) the existence of critical damping beyond which no tunneling can occur, (b) that tunneling trajectories are dependent on the damping constant and (c) that dissipation may stabilize the excited state rather than the ground state. Finally the tunneling time delay in one-dimensional systems for undamped and for dissipative systems is formulated in terms of the phase shift, and this has been used to show that the effect of damping on the time delay is ignorable.     

Exact class of inverse scattering solutions

After presenting a procedure for finding exact solutions to the Gelfand-Levitan equation for inverse scattering, we present explicit solutions for 3- to 6-pole reflection coefficients. The method is applicable to an arbitrarily large number of poles. Here we give the first explicit treatment of 4, 5, and 6 poles.     

Generalized source method: New possibilities for waveguide and grating problems

A new, general and exact method for resolving waveguide and grating problems is presented whereby the unknown electric field in a given dielectric distribution in the presence of an arbitrary source distribution can be obtained exactly from the known complete electric field solution in a simpler dielectric structure in the presence of the same source distribution. The method can be used in the form of an iterative scheme. The solution can also be obtained through an implicit equation which can be solved numerically without convergence problem by matrix inversion. The application of the method to the case of abnormal reflection under normal incidence from a segmented waveguide is given as an example.     

A stable complex Jacobi iterative solution of 3D semivectorial wide-angle beam propagation using the iterated Crank–Nicholson method

An extension of the recently proposed three-dimensional (3D) wide-angle (WA) beam propagation method (BPM) using complex Jacobi iteration (CJI) taken into account polarization effects is presented. The resulting iterative BPM is faster than BPMs based on the traditional direct matrix inversion for waveguides with unchanging refractive index profiles during propagation direction. However, for varying refractive index waveguides the iterative method suffered from the fact that the iteration count between two successive cross-sections increases dramatically during the propagation direction. To overcome this problem, we propose the utility of the iterated Crank–Nicholson method. At each propagation step, the propagation equation is divided in multiple stages by the iterated Crank–Nicholson method and then each stage is recast in terms of a Helmholtz equation with source term, which is solved effectively by the complex Jacobi iterative method. The resulting approach is stable and well-suited for large structures with long propagation paths.     

A Weak Form of the Conjugate Gradient FFT Method for Two-Dimensional Elastodynamics

The problem of two-dimensional scattering of elastic waves by an elastic inclusion can be formulated in terms of a domain integral equation, in which the grad-div operator acts on a vector potential. The vector potential is the spatial convolution of a Green's function with the product of the density and the displacement over the domain of interest. A weak form of the integral equation for the unknown displacement is obtained by testing it with rooftop functions. This method shows excellent numerical performance.     

Plane-wave approximation in theS-matrix and the construction of bound-state wave functions

Three inversion problem approaches — byGelfand-Levitan, Marchenko andPetrá? [5] —in both non-relativistic and relativistic (Klein-Gordon) variants are used in an approximation scheme selected to construct bound-state wave functions which are advantageous for purposes of model hadron physics. This family of wave functions is created exclusively by theS-matrix quantities and derived in the approximation which requires the Jost function to be equal to the unity throughout the continuous spectrum (the plane-wave approximation). As a consequence of the difference in boundary conditions of the mentioned approaches, the resulting approximate wave functions are not identical, but it is shown that there exists a parallelism as to the form among them. This parallelism is explained more extensively in the non-relativistic case, where the transformation properties of alternative sets of functions are treated. In the present paper it is demonstrated that in the relativistic variants of the above approaches the non-relativistic plane-wave-approximation form of the constructed wave functions for a given bound state is preserved.