水平层状介质并矢Green函数的递推矩阵方法

采用递推矩阵方法计算任意数目水平层状介质的并矢Green函数.根据层界面处电场和磁场的连续性条件得到3个确定Sommerfeld积分待定系数的矩阵方程组,分别对应于垂向单位电偶极子产生的TM波、水平方向单位电偶极子产生的TE波和TM波,这些方程组均可通过递推方法快速求解.只需改变3个方程组中源项元素的位置,就可以方便地得到当源点和场点在任意层时的并矢Green函数.本文给出的并矢Green函数表达式形式简洁且不含指数增加项,计算时不会出现溢出现象.     

多体关联Green函数动力学Ⅲ截断近似

从多体关联Green函数的基本动力学方程组出发,通过对多体关联Green函数作低阶截断近似,得到了通常Green函数理论中所有主要的非微扰结果.而且二体关联Green函数动力学把梯形图对应的短程关联与环形图对应的长程关联统一在同一个动力学方程组中,这是通常的Green函数理论所没有得到的新结果.     

多体关联Green函数动力学——Ⅲ.截断近似

从多体关联Green函数的基本动力学方程组出发,通过对多体关联Green函数作低阶截断近似,得到了通常Green函数理论中所有主要的非微扰结果。而且二体关联Green函数动力学把梯形图对应的短程关联与环形图对应的长程关联统一在同一个动力学方程组中,这是通常的Green函数理论所没有得到的新结果。     

偶极子在径向非均匀介质中的电磁场分布

提出了柱状多层介质中的磁偶极子在任意介质层中的电磁场的递推算法，给出了磁偶极子产生的电磁波在柱状多层介质中的反射、透射系数和广义反射系数，从而得到了任意介质层中满足各层界面边界条件的电磁场解析表达式．     

多体关联Green函数动力学Ⅱ.等时极限

对多体关联Green函数的基本动力学方程组作等时极限,得到了密度矩阵形式的多体关联动力学方程组,并证明了各种不同时序的等时极限间的等价性.     

水平层状介质中任意方向磁偶极子的电磁场分布

将水平层状介质中的任意方向磁偶极子分解成水平磁偶极子(HMD)和垂直磁偶极子(VMD),给出了将VMD和HMD产生的波分解为电场矢量垂直于入射面的线形极化波(EV波)和平行于入射面的线形极化波(EP波)的方法.通过分别研究这两种波在各层介质中的反射和透射规律,导出确定各介质层中电磁波的递推公式,从而得到了任意介质层中电磁波的解析表达式.     

随机介质背散射二维Mueller矩阵的数值模拟与分析

推导了随机介质背散射Mueller矩阵的直接计算公式,并运用矢量Monte Carlo方法进行了数值模拟.结果表明随机介质背散射二维Mueller矩阵方位关系随散射系数的减小而增强,而与微粒大小关系不大;Mueller矩阵元素绝对值的空间分布随径向呈近似指数规律衰减,矩阵元素的方位变化具有周期性.对称系统的二维Mueller矩阵的花样图中仅有7幅独立,其余9幅可通过对称、旋转变换得到.     

用改进的传输线算法计算水平层状横向同性地层中海洋可控源电磁响应

利用二维Fourier变换与电磁场分解技术将层状横向同性地层中Maxwell方程转化成两个独立的关于横磁(TM)波和横电(TE)波的传输线方程; 借助传输线理论与叠加原理, 仅利用电流源传输线Green函数得到TM波和TE波的解, 改进传输线算法, 建立横向同性地层中频率-波数域电流源电场和磁场并矢Green函数的新算法与新的解析表达式, 提高海洋可控源电磁响应数值模拟效率. 在此基础上, 利用传输线Green函数的基本解以及边界条件, 推导出广义反射系数与振幅递推公式, 得到各个地层中传输线Green函数的解析解; 然后利用Fourier逆变换与Bessel公式将海洋可控源电磁响应表示为Sommerfeld形式的积分, 借助三次样条插值与Lommel积分公式快速计算其数值解. 通过数值模拟结果考察工作频率以及地层各向异性电阻率变化等对海洋电磁响应的影响. 关键词： 传输线法 横向同性地层 海洋可控源电磁 Sommerfeld积分     

Rssler超混沌系统多变量广义预测控制快速算法

提出了一种Rssler超混沌系统多变量广义预测控制快速算法.采用改进的时变遗忘因子递推最小二乘方法辨识混沌系统,通过在常规多变量广义预测控制性能指标函数中引入前馈增益矩阵与柔化矩阵,实现系统对参考信号的快速跟踪.该算法降低了求逆矩阵的维数,从而减少了计算量.仿真结果验证了该方法的有效性和实用性.     

R(o)ssler超混沌系统多变量广义预测控制快速算法

提出了一种R(o)ssler超混沌系统多变量广义预测控制快速算法.采用改进的时变遗忘因子递推最小二乘方法辨识混沌系统,通过在常规多变量广义预测控制性能指标函数中引入前馈增益矩阵与柔化矩阵,实现系统对参考信号的快速跟踪.该算法降低了求逆矩阵的维数,从而减少了计算量.仿真结果验证了该方法的有效性和实用性.     

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.     

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.     

基于格林函数的多终端量子链状体系电子输运性质的研究

利用递推格林函数技术计算了多终端量子体系的电子输运特性，首先运用递归方法给出介观 或量子体系的格林函数. 然后利用散射矩阵和输运方程给出体系的电导方程，可以将多终端 的输运简化为双终端的输运方程，以便得到体系电子输运的谱结构. 计算结果表明，由于中 间节点的存在，使器件的传输谱偏离一维链的对称性，在低能量端出现一个新的电导峰值. 此外，本方法可以被应用到各种复杂的带有吸附结构量子体系输运的研究中. 关键词： 格林函数 散射矩阵 量子体系     

Modeling of Coulomb coupling and acoustic wave propagation in LiNbO3

The angular spectrum propagation technique is applied to modeling of wave propagation through piezoelectric media. The calculations of the angular spectrum propagator are based on the relevant equation for the slowness surface resulting from the solution of the Christoffel equation with piezoelectrically stiffened elastic constants. A two-dimensional FFT algorithm is applied in the final field superposition. We concentrate on the case of Coulomb coupling through local electrical point contacts on both the excitation and detection side. To model that case we superpose solutions for acoustic Green's functions of different propagation modes convoluted with equivalent distributed effective sources. Calculated results are in good agreement with the measured ones.     

Vector Diffraction Integrals for Solving Inverse Problems of Radio-Holographic Sensing of the Earth's Surface and Atmosphere

Vector relationships between the fields on a certain surface confining an inhomogeneous three-dimensional volume and the fields inside this volume are obtained by the Stratton–Chu method developed for the case of homogeneous media. The vector relationships allow us to solve the direct and inverse problems of determining the fields inside an inhomogeneous medium given the field on its boundary. The vector equations take into acount the polarization changes of direct and inverse waves propagated in an inhomogeneous medium. In the case of a two-dimensional homogeneous medium, the vector equations reduce to the previously obtained scalar equations used in the approximation of spherical symmetry to describe the process of backward wave propagation during the atmospheric and ionospheric radio-occultation monitoring. It is shown that the Green's function of the scalar wave equation in an inhomogeneous medium should be used as the reference signal for solving the inverse problem of radio-occultation monitoring. This validates the method of focused synthetic aperture previously used for high-accuracy retrieval of the vertical refractive-index profiles in the ionosphere and atmosphere. In this method, the reference-signal phase was determined from a model which describes with sufficient accuracy the radiophysical parameters of a refracting medium in the region of radio-occultation sensing. The obtained equations can be used for the high-accuracy solving of inverse problems of radio-holographic sensing of the Earth's atmosphere and surface by precision signals from radio-navigation satellites.     

Topological Invariants in Terms of Green's Function for the Interacting Kitaev Chain

A one-dimensional closed interacting Kitaev chain and the dimerized version are studied. The topological invariants in terms of Green's function are calculated by the density matrix renormalization group method and the exact diagonalization method. For the interacting Kitaev chain, we point out that the calculation of the topological invariant in the charge density wave phase must consider the dimerized configuration of the ground states. The variation of the topological invariant is attributed to the poles of eigenvalues of the zero-frequency Green functions. For the interacting dimerized Kitaev chain, we show that the topological invariant defined by Green's functions can distinguish more topological nonequivalent phases than the fermion parity.     

Solution of near-field thermal radiation in one-dimensional layered media using dyadic Green's functions and the scattering matrix method

A general algorithm is introduced for the analysis of near-field radiative heat transfer in one-dimensional multi-layered structures. The method is based on the solution of dyadic Green's functions, where the amplitude of the fields in each layer is calculated via a scattering matrix approach. Several tests are presented where cubic boron nitride is used in the simulations. It is shown that a film emitter thicker than 1 μm provides the same spectral distribution of near-field radiative flux as obtained from a bulk emitter. Further simulations have pointed out that the presence of a body in close proximity to an emitter can alter the near-field spectrum emitted. This algorithm can be employed to study thermal one-dimensional layered media and photonic crystals in the near-field in order to design radiators optimizing the performances of nanoscale-gap thermophotovoltaic power generators.     

Green's function for arbitrarily shaped dielectric bodies of revolution

In this paper, a solution is developed to calculate the electric field at one point in space due to an electric dipole exciting an arbitrarily shaped dielectric body of revolution (BOR). Specifically, the electric field is determined from the solution of coupled surface integral equations (SIE) for the induced surface electric and magnetic currents on the dielectric body excited by an elementary electric current dipole source. Both the interior and exterior fields to the dielectric BOR may be accurately evaluated via this approach. For a highly lossy dielectric body, the numerical Green's function is also obtainable from an approximate integral equation (AIE) based on a surface boundary condition. If this equation is solved by the method of moments, significant numerical efficiency over SIE is realized. Numerical results obtained by both SIE and AIE approaches agree with the exact solution for the special case of a dielectric sphere. With this numerical Green's function, the complicated radiation and scattering problems in the presence of an arbitrarily shaped dielectric BOR are readily solvable by the method of moments.