A Fast, Direct Algorithm for the Lippmann–Schwinger Integral Equation in Two Dimensions

The state-of-the-art, large-scale numerical simulations of the scattering problem for the Helmholtz equation in two dimensions rely on iterative solvers for the Lippmann–Schwinger integral equation, with an optimal CPU time O(m ³ log(m)) for an m-by-m wavelength problem. We present a method to solve the same problem directly, as opposed to iteratively, with the obvious advantage in efficiency for multiple right-hand sides corresponding to distinct incident waves. Analytically, this direct method is a hierarchical, recursive scheme consisting of the so-called splitting and merging processes. Algebraically, it amounts to a recursive matrix decomposition, for a cost of O(m ³), of the discretized Lippmann–Schwinger operator. With this matrix decomposition, each back substitution requires only O(m ² log(m)); therefore, a scattering problem with m incident waves can be solved, altogether, in O(m ³ log(m)) flops.     

Jordan标准形定理的初等变换证明

本文对 Jordan标准形定理给出了一种使用初等变换的证明 ,直观意义明显、易于理解 ,可用于线性代数教学 .     

Block plan for a fractional 2 m factorial design derived from a 2 r factorial design

Summary For a given fractional 2 ^m factorial (2 ^m -FF) designT, the constitution of a block plan to divideT intok (2 ^r−1<k≦2 ^r ) blocks withr block factors each at two levels is proposed and investigated. The well-known three norms of the confounding matrix are used as measures for determining a “good” block plan. Some theorems concerning the constitution of a block plan are derived for a 2 ^m -FF design of odd or even resolution. Two norms which may be preferred over the other norm are slightly modified. For each value ofN assemblies with 11≦N≦26, optimum block plans fork=2 blocks with block sizes [N/2] andN−[N/2] minimizing the two norms are presented forA-optimal balanced 2⁴-FF designs of resolutionV given by Srivastava and Chopra (Technometrics,13, 257–269).     

Achieving Optical Tri-Stability Through Periodic Metal-Dielectric Structure

Abstract

Through the modified transfer matrix method, the transmission properties of a one-dimensional coupled-resonator optical waveguide structure composed of metal layers and non-linear material layers is studied. Given proper incident frequency and structure parameters, an optical tri-stability has been achieved. The effect of loss has been considered.     

Derived equivalences between triangular matrix algebras

In this paper, we study derived equivalences between triangular matrix algebras using certain classical recollements. We show that special properties of these recollements actually characterize triangular matrix algebras and describe methods to construct tilting modules and tilting complexes inducing derived equivalences between them.     

An iterative method for solving the stable subspace of a matrix pencil and its application

This work is to propose an iterative method of choice to compute a stable subspace of a regular matrix pencil. This approach is to define a sequence of matrix pencils via particular left null spaces. We show that this iteration preserves a semigroup property depending only on the initial matrix pencil. Via this recursion relationship, we propose an accelerated iterative method to compute the stable subspace and use it to provide a theoretical result to solve the principal square root of a given matrix, both nonsingular and singular. We show that this method can not only find out the matrix square root, but also construct an iterative approach which converges to the square root with any desired order.     

Difference operators,Green's matrix,sampling theory and applications in signal processing

This paper introduces sampling representations for discrete signals arising from self adjoint difference operators with mixed boundary conditions. The theory of linear operators on finite-dimensional inner product spaces is employed to study the second-order difference operators. We give necessary and sufficient conditions that make the operators self adjoint. The equivalence between the difference operator and a Hermitian Green's matrix is established. Sampling theorems are derived for discrete transforms associated with the difference operator. The results are exhibited via illustrative examples, involving sampling representations for the discrete Hartley transform. Families of discrete fractional Fourier-type transforms are introduced with an application to image encryption.     

Integrable Couplings of the Boiti-Pempinelli-Tu Hierarchy and Their Hamiltonian Structures

Integrable couplings of the Boiti-Pempinelli-Tu hierarchy are constructed by a class of non-semisimple block matrix loop algebras. Further, through using the variational identity theory, the Hamiltonian structures of those integrable couplings are obtained. The method can be applied to obtain other integrable hierarchies.     

Is breaking of ensemble equivalence monotone in the number of constraints?

Breaking of ensemble equivalence between the microcanonical ensemble and the canonical ensemble may occur for random graphs whose size tends to infinity, and is signaled by a non-zero specific relative entropy between the two ensembles. In Garlaschelli et al. (2017) and Garlaschelli et al. (0000) it was shown that breaking occurs when the constraint is put on the degree sequence (configuration model). It is not known what is the effect on the relative entropy when the number of constraints is reduced, i.e., when only part of the nodes are constrained in their degree (and the remaining nodes are left unconstrained). Intuitively, the relative entropy is expected to decrease. However, this is not a trivial issue because when constraints are removed both the microcanonical ensemble and the canonical ensemble change. In this paper a formula for the relative entropy valid for generic discrete random structures, recently formulated by Squartini and Garlaschelli, is used to prove that the relative entropy is monotone in the number of constraints when the constraint is on the degrees of the nodes. It is further shown that the expression for the relative entropy corresponds, in the dense regime, to the degrees in the microcanonical ensemble being asymptotically multivariate Dirac and in the canonical ensemble being asymptotically Gaussian.