• Abstract

We consider the inverse source problems with multi-frequency sparse near field measurements. In contrast to the existing near field operator based on the integral over the space variable, a multi-frequency near field operator is introduced based on the integral over the frequency variable. A factorization of this multi-frequency near field operator is further given and analyzed. Based on such a factorization, we introduce a single-receiver multi-frequency sampling method to reconstruct a shell support of the source. Its theoretical foundation is derived from the properties of the factorized operators and a properly chosen point spread function. Numerical examples are provided to illustrate the multi-frequency sampling method with sparse near field measurements. Finally we briefly discuss how to extend the near field case to the far field case.

• General Summary

• Abstract

Matrix splitting iteration methods play a vital role in solving large sparse linear systems. Their performance heavily depends on the splitting parameters, however, the approach of selecting optimal splitting parameters has not been well developed. In this paper, we present a multitask kernel-learning parameter prediction method to automatically obtain relatively optimal splitting parameters, which contains simultaneous multiple parameters prediction and a data-driven kernel learning. For solving time-dependent linear systems, including linear differential systems and linear matrix systems, we give a new matrix splitting Kronecker product method, as well as its convergence analysis and preconditioning strategy. Numerical results illustrate our methods can save an enormous amount of time in selecting the relatively optimal splitting parameters compared with the exists methods. Moreover, our iteration method as a preconditioner can effectively accelerate GMRES. As the dimension of systems increases, all the advantages of our approaches become significantly. Especially, for solving the differential Sylvester matrix equation, the speedup ratio can reach tens to hundreds of times when the scale of the system is larger than one hundred thousand.

• General Summary

• Abstract

This paper investigates the elastic scattering by unbounded deterministic and random rough surfaces, both of which are assumed to be graphs of Lipschitz continuous functions. For the deterministic case, an a priori bound explicitly dependent on frequencies is derived by the variational approach. For the scattering by random rough surfaces with a random source, well-posedness of the corresponding variation problem is proved. Moreover, a similar bound with explicit dependence on frequencies for the random case is also established based upon the deterministic result, Pettis measurability theorem and Bochner’s integrability theorem.

• General Summary

• Abstract

In this paper, we develop an asymptotic expansion-regularization (AER) method for inverse source problems in two-dimensional nonlinear and nonstationary singularly perturbed partial differential equations (PDEs). The key idea of this approach is the use of the asymptotic-expansion theory, which allows us to determine the conditions for the existence and uniqueness of a solution to a given PDE with a sharp transition layer. As a by-product, we derive a simpler link equation between the source function and first-order asymptotic approximation of the measurable quantities, and based on that equation we propose an efficient inversion algorithm, AER, for inverse source problems. We prove that this simplification will not decrease the accuracy of the inversion result, especially for inverse problems with noisy data. Various numerical examples are provided to demonstrate the efficiency of our new approach.

• General Summary

• Abstract

Capacity constrained optimal transport is a variant of optimal transport, which adds extra constraints on the set of feasible couplings in the original optimal transport problem to limit the mass transported between each pair of source and sink. Based on this setting, constrained optimal transport has numerous applications, e.g., finance, network flow. However, due to the large number of constraints in this problem, existing algorithms are both time-consuming and space-consuming. In this paper, inspired by entropic regularization for the classical optimal transport problem, we introduce a novel regularization term for capacity constrained optimal transport. The regularized problem naturally satisfies the capacity constraints and consequently makes it possible to analyze the duality. Unlike the matrix-vector multiplication in the alternate iteration scheme for solving classical optimal transport, in our algorithm, each alternate iteration step is to solve several single-variable equations. Fortunately, we find that each of these equations corresponds to a single-variable monotonic function, and we convert solving these equations into finding the unique zero point of each single-variable monotonic function with Newton’s method. Theoretical analysis further provides a convergence guarantee to our algorithm. Extensive numerical experiments demonstrate that our proposed method has a significant advantage in terms of accuracy, efficiency, and memory consumption compared with existing methods.

• General Summary

• Abstract

In this paper, we propose a two-level simultaneous orthogonal matching pursuit (TLSOMP) algorithm for simultaneous sparse approximation (SSA) problems. Most existing algorithms for SSA problems are directly generalized from the ones for the sparse approximation (SA) problems, for example, the simultaneous orthogonal matching pursuit (SOMP) method is generalized from the orthogonal matching pursuit (OMP) method. Our newly proposed algorithm is designed from another viewpoint. We first analyze the noiseless case and propose a selection algorithm. Motivated by the analysis and presuming noise as a perturbation, we extend the selection algorithm into a TLSOMP algorithm. This novel algorithm mainly uses the information from the subspace spanned by the multiple signals, which is not available in SA problems. Numerical experiments show the superiority of our TLSOMP algorithm over other traditional SSA solvers.

• General Summary