First, we build a computational procedure to reconstruct the convex body from a pre-given surface area measure. Nontrivially, we prove the convergence of this procedure. Then, the sufficient and necessary conditions of a Sobolev binary function to be a lightness function of a convex body are investigated. Finally, we present a computational procedure to compute the curvature function from a pre-given lightness function, and then we reconstruct the convex body from this curvature function by using the first procedure. The convergence is also proved. The main computations in both procedures are simply about the spherical harmonics.

The Orthogonal Least Squares (OLS) algorithm is an efficient sparse recovery algorithm that has received much attention in recent years. On one hand, this paper considers that the OLS algorithm recovers the supports of sparse signals in the noisy case. We show that the OLS algorithm exactly recovers the support of $K$-sparse signal $\boldsymbol{x}$ from $\boldsymbol{y}=\boldsymbol{\unicode[STIX]{x1D6F7}}\boldsymbol{x}+\boldsymbol{e}$ in $K$ iterations, provided that the sensing matrix $\boldsymbol{\unicode[STIX]{x1D6F7}}$ satisfies the restricted isometry property (RIP) with restricted isometry constant (RIC) $\unicode[STIX]{x1D6FF}_{K+1}<1/\sqrt{K+1}$, and the minimum magnitude of the nonzero elements of $\boldsymbol{x}$ satisfies some constraint. On the other hand, this paper demonstrates that the OLS algorithm exactly recovers the support of the best $K$-term approximation of an almost sparse signal $\boldsymbol{x}$ in the general perturbations case, which means both $\boldsymbol{y}$ and $\boldsymbol{\unicode[STIX]{x1D6F7}}$ are perturbed. We show that the support of the best $K$-term approximation of $\boldsymbol{x}$ can be recovered under reasonable conditions based on the restricted isometry property (RIP).

In the paper the correspondence between a formal multiple power series and a special type of branched continued fractions, the so-called ‘multidimensional regular C-fractions with independent variables’ is analysed providing with an algorithm based upon the classical algorithm and that enables us to compute from the coefficients of the given formal multiple power series, the coefficients of the corresponding multidimensional regular C-fraction with independent variables. A few numerical experiments show, on the one hand, the efficiency of the proposed algorithm and, on the other, the power and feasibility of the method in order to numerically approximate certain multivariable functions from their formal multiple power series.

The generalized orthogonal matching pursuit $\left( \text{gOMP} \right)$ algorithm has received much attention in recent years as a natural extension of the orthogonal matching pursuit $\left( \text{OMP} \right)$. It is used to recover sparse signals in compressive sensing. In this paper, a new bound is obtained for the exact reconstruction of every $K$-sparse signal via the $\text{gOMP}$ algorithm in the noiseless case. That is, if the restricted isometry constant $\left( \text{RIC} \right)$${{\delta }_{NK+1}}$ of the sensing matrix $A$ satisfies

${{\delta }_{NK+1}}\,<\frac{1}{\sqrt{\frac{K}{N}+\,1}},$

then the $\text{gOMP}$ can perfectly recover every $K$-sparse signal $x$ from $y\,=\,Ax$. Furthermore, the bound is proved to be sharp. In the noisy case, the above bound on $\text{RIC}$ combining with an extra condition on the minimum magnitude of the nonzero components of $K$-sparse signals can guarantee that the $\text{gOMP}$ selects all of the support indices of the $K$-sparse signals.

In this paper, a new stopping rule is proposed for orthogonal multi-matching pursuit (OMMP). We show that, for ℓ2 bounded noise case, OMMP with the new stopping rule can recover the true support of any K-sparse signal x from noisy measurements y = Фx + e in at most K iterations, provided that all the nonzero components of x and the elements of the matrix Ф satisfy certain requirements. The proposed method can improve the existing result. In particular, for the noiseless case, OMMP can exactly recover any K-sparse signal under the same RIP condition.

Sign truncated matching pursuit (STrMP) algorithm is presented in this paper. STrMP is a new greedy algorithm for the recovery of sparse signals from the sign measurement, which combines the principle of consistent reconstruction with orthogonal matching pursuit (OMP). The main part of STrMP is as concise as OMP and hence STrMP is simple to implement. In contrast to previous greedy algorithms for one-bit compressed sensing, STrMP only need to solve a convex and unconstrained subproblem at each iteration. Numerical experiments show that STrMP is fast and accurate for one-bit compressed sensing compared with other algorithms.

An iterative discontinuous Galerkin (DG) method is proposed to solve the nonlinear Poisson Boltzmann (PB) equation. We first identify a function space in which the solution of the nonlinear PB equation is iteratively approximated through a series of linear PB equations, while an appropriate initial guess and a suitable iterative parameter are selected so that the solutions of linear PB equations are monotone within the identified solution space. For the spatial discretization we apply the direct discontinuous Galerkin method to those linear PB equations. More precisely, we use one initial guess when the Debye parameter λ = (1), and a special initial guess for λ ≫1 to ensure convergence. The iterative parameter is carefully chosen to guarantee the existence, uniqueness, and convergence of the iteration. In particular, iteration steps can be reduced for a variable iterative parameter. Both one and two-dimensional numerical results are carried out to demonstrate both accuracy and capacity of the iterative DG method for both cases of λ = (1) and λ ≪ 1. The (m + 1)th order of accuracy for L2 and mth order of accuracy for H1 for Pm elements are numerically obtained.

A new representation of spatio-temporal random processes is proposed in this work. In practical applications, such processes are used to model velocity fields, temperature distributions, response of vibrating systems, to name a few. Finding an efficient representation for any random process leads to encapsulation of information which makes it more convenient for a practical implementations, for instance, in a computational mechanics problem. For a single-parameter process such as spatial or temporal process, the eigenvalue decomposition of the covariance matrix leads to the well-known Karhunen-Loève (KL) decomposition. However, for multiparameter processes such as a spatio-temporal process, the covariance function itself can be defined in multiple ways. Here the process is assumed to be measured at a finite set of spatial locations and a finite number of time instants. Then the spatial covariance matrix at different time instants are considered to define the covariance of the process. This set of square, symmetric, positive semi-definite matrices is then represented as a third-order tensor. A suitable decomposition of this tensor can identify the dominant components of the process, and these components are then used to define a closed-form representation of the process. The procedure is analogous to the KL decomposition for a single-parameter process, however, the decompositions and interpretations vary significantly. The tensor decompositions are successfully applied on (i) a heat conduction problem, (ii) a vibration problem, and (iii) a covariance function taken from the literature that was fitted to model a measured wind velocity data. It is observed that the proposed representation provides an efficient approximation to some processes. Furthermore, a comparison with KL decomposition showed that the proposed method is computationally cheaper than the KL, both in terms of computer memory and execution time.

A class of numerical methods is developed for second order Volterra integro-differential equations by using a Legendre spectral approach. We provide a rigorous error analysis for the proposed methods, which shows that the numerical errors decay exponentially in the L∞-norm and L2-norm. Numerical examples illustrate the convergence and effectiveness of the numerical methods.

Reflected Brownian motion is used in areas such as physiology, electrochemistry and nuclear magnetic resonance. We study the first-passage-time problem of this process which is relevant in applications; specifically, we find a Volterra integral equation for the distribution of the first time that a reflected Brownian motion reaches a nondecreasing barrier. Additionally, we note how a numerical procedure can be used to solve the integral equation.

In, Tory and Pickard show that a simple subclass of unilateral AR processes identifies with Gaussian Pickard random fields on Z2. First, we extend this result to the whole class of unilateral AR processes, by showing that they all satisfy a Pickard-type property, under which correlation matching and maximum entropy properties are assessed. Then, it is established that the Pickard property provides the ‘missing’ equations that complement the two-dimensional Yule-Walker equations, in the sense that the conjunction defines a one-to-one mapping between the set of AR parameters and a set of correlations. It also implies Markov chain conditions that allow exact evaluation of the likelihood and an exact sampling scheme on finite lattices.