In this paper,we propose a Quasi-Orthogonal Matching Pursuit(QOMP)algorithm for constructing a sparse approximation of functions in terms of expansion by orthonormal polynomials.For the two kinds of sampled data,data ...In this paper,we propose a Quasi-Orthogonal Matching Pursuit(QOMP)algorithm for constructing a sparse approximation of functions in terms of expansion by orthonormal polynomials.For the two kinds of sampled data,data with noises and without noises,we apply the mutual coherence of measurement matrix to establish the convergence of the QOMP algorithm which can reconstruct s-sparse Legendre polynomials,Chebyshev polynomials and trigonometric polynomials in s step iterations.The results are also extended to general bounded orthogonal system including tensor product of these three univariate orthogonal polynomials.Finally,numerical experiments will be presented to verify the effectiveness of the QOMP method.展开更多
Orthogonal matching pursuit(OMP for short)algorithm is a popular method of sparse signal recovery in compressed sensing.This paper applies OMP to the sparse polynomial reconstruction problem.Distinguishing from classi...Orthogonal matching pursuit(OMP for short)algorithm is a popular method of sparse signal recovery in compressed sensing.This paper applies OMP to the sparse polynomial reconstruction problem.Distinguishing from classical research methods using mutual coherence or restricted isometry property of the measurement matrix,the recovery guarantee and the success probability of OMP are obtained directly by the greedy selection ratio and the probability theory.The results show that the failure probability of OMP given in this paper is exponential small with respect to the number of sampling points.In addition,the recovery guarantee of OMP obtained through classical methods is lager than that of ℓ_(1)-minimization whatever the sparsity of sparse polynomials is,while the recovery guarantee given in this paper is roughly the same as that of ℓ_(1)-minimization when the sparsity is less than 93.Finally,the numerical experiments verify the availability of the theoretical results.展开更多
In this paper,we investigate the coupling of the Multi-dimensional Optimal Order Detection(MOOD)method and the Arbitrary high order DERivatives(ADER)approach in order to design a new high order accurate,robust and com...In this paper,we investigate the coupling of the Multi-dimensional Optimal Order Detection(MOOD)method and the Arbitrary high order DERivatives(ADER)approach in order to design a new high order accurate,robust and computationally efficient Finite Volume(FV)scheme dedicated to solve nonlinear systems of hyperbolic conservation laws on unstructured triangular and tetrahedral meshes in two and three space dimensions,respectively.The Multi-dimensional Optimal Order Detection(MOOD)method for 2D and 3D geometries has been introduced in a recent series of papers for mixed unstructured meshes.It is an arbitrary high-order accurate Finite Volume scheme in space,using polynomial reconstructions with a posteriori detection and polynomial degree decrementing processes to deal with shock waves and other discontinuities.In the following work,the time discretization is performed with an elegant and efficient one-step ADER procedure.Doing so,we retain the good properties of the MOOD scheme,that is to say the optimal high-order of accuracy is reached on smooth solutions,while spurious oscillations near singularities are prevented.The ADER technique permits not only to reduce the cost of the overall scheme as shown on a set of numerical tests in 2D and 3D,but it also increases the stability of the overall scheme.A systematic comparison between classical unstructured ADER-WENO schemes and the new ADER-MOOD approach has been carried out for high-order schemes in space and time in terms of cost,robustness,accuracy and efficiency.The main finding of this paper is that the combination of ADER with MOOD generally outperforms the one of ADER and WENO either because at given accuracy MOOD is less expensive(memory and/or CPU time),or because it is more accurate for a given grid resolution.A large suite of classical numerical test problems has been solved on unstructured meshes for three challenging multi-dimensional systems of conservation laws:the Euler equations of compressible gas dynamics,the classical equations of ideal magneto-Hydrodynamics(MHD)and finally the relativistic MHD equations(RMHD),which constitutes a particularly challenging nonlinear system of hyperbolic partial differential equation.All tests are run on genuinely unstructured grids composed of simplex elements.展开更多
A background removal method based on two-dimensional notch filtering in the frequency domain for polarization interference imaging spectrometers(PIISs) is implemented. According to the relationship between the spati...A background removal method based on two-dimensional notch filtering in the frequency domain for polarization interference imaging spectrometers(PIISs) is implemented. According to the relationship between the spatial domain and the frequency domain, the notch filter is designed with several parameters of PIISs, and the interferogram without a background is obtained. Both the simulated and the experimental results demonstrate that the background removal method is feasible and robust with a high processing speed. In addition, this method can reduce the noise level of the reconstructed spectrum, and it is insusceptible to a complicated background, compared with the polynomial fitting and empirical mode decomposition(EMD) methods.展开更多
基金supported by National Natural Science Foundation of China no.12071019.
文摘In this paper,we propose a Quasi-Orthogonal Matching Pursuit(QOMP)algorithm for constructing a sparse approximation of functions in terms of expansion by orthonormal polynomials.For the two kinds of sampled data,data with noises and without noises,we apply the mutual coherence of measurement matrix to establish the convergence of the QOMP algorithm which can reconstruct s-sparse Legendre polynomials,Chebyshev polynomials and trigonometric polynomials in s step iterations.The results are also extended to general bounded orthogonal system including tensor product of these three univariate orthogonal polynomials.Finally,numerical experiments will be presented to verify the effectiveness of the QOMP method.
基金supported by the National Natural Science Foundation of China(Grant No.12071019).
文摘Orthogonal matching pursuit(OMP for short)algorithm is a popular method of sparse signal recovery in compressed sensing.This paper applies OMP to the sparse polynomial reconstruction problem.Distinguishing from classical research methods using mutual coherence or restricted isometry property of the measurement matrix,the recovery guarantee and the success probability of OMP are obtained directly by the greedy selection ratio and the probability theory.The results show that the failure probability of OMP given in this paper is exponential small with respect to the number of sampling points.In addition,the recovery guarantee of OMP obtained through classical methods is lager than that of ℓ_(1)-minimization whatever the sparsity of sparse polynomials is,while the recovery guarantee given in this paper is roughly the same as that of ℓ_(1)-minimization when the sparsity is less than 93.Finally,the numerical experiments verify the availability of the theoretical results.
基金the European Research Council(ERC)under the European Union’s Seventh Framework Programme(FP7/2007-2013)the research project STiMulUs,ERC Grant agreement no.278267+1 种基金.R.L.has been partially funded by the ANR under the JCJC project“ALE INC(ubator)3D”the reference LA-UR-13-28795.The authors would like to acknowledge PRACE for awarding access to the SuperMUC supercomputer based in Munich,Germany at the Leibniz Rechenzentrum(LRZ)。
文摘In this paper,we investigate the coupling of the Multi-dimensional Optimal Order Detection(MOOD)method and the Arbitrary high order DERivatives(ADER)approach in order to design a new high order accurate,robust and computationally efficient Finite Volume(FV)scheme dedicated to solve nonlinear systems of hyperbolic conservation laws on unstructured triangular and tetrahedral meshes in two and three space dimensions,respectively.The Multi-dimensional Optimal Order Detection(MOOD)method for 2D and 3D geometries has been introduced in a recent series of papers for mixed unstructured meshes.It is an arbitrary high-order accurate Finite Volume scheme in space,using polynomial reconstructions with a posteriori detection and polynomial degree decrementing processes to deal with shock waves and other discontinuities.In the following work,the time discretization is performed with an elegant and efficient one-step ADER procedure.Doing so,we retain the good properties of the MOOD scheme,that is to say the optimal high-order of accuracy is reached on smooth solutions,while spurious oscillations near singularities are prevented.The ADER technique permits not only to reduce the cost of the overall scheme as shown on a set of numerical tests in 2D and 3D,but it also increases the stability of the overall scheme.A systematic comparison between classical unstructured ADER-WENO schemes and the new ADER-MOOD approach has been carried out for high-order schemes in space and time in terms of cost,robustness,accuracy and efficiency.The main finding of this paper is that the combination of ADER with MOOD generally outperforms the one of ADER and WENO either because at given accuracy MOOD is less expensive(memory and/or CPU time),or because it is more accurate for a given grid resolution.A large suite of classical numerical test problems has been solved on unstructured meshes for three challenging multi-dimensional systems of conservation laws:the Euler equations of compressible gas dynamics,the classical equations of ideal magneto-Hydrodynamics(MHD)and finally the relativistic MHD equations(RMHD),which constitutes a particularly challenging nonlinear system of hyperbolic partial differential equation.All tests are run on genuinely unstructured grids composed of simplex elements.
基金supported by the Major Program of the National Natural Science Foundation of China(No.41530422)the National Science and Technology Major Project of the Ministry of Science and Technology of China(No.32-Y30B08-9001-13/15)+1 种基金the National Natural Science Foundation of China(Nos.61275184,61540018,61405153,and 60278019)the National High Technology Research and Development Program of China(No.2012AA121101)
文摘A background removal method based on two-dimensional notch filtering in the frequency domain for polarization interference imaging spectrometers(PIISs) is implemented. According to the relationship between the spatial domain and the frequency domain, the notch filter is designed with several parameters of PIISs, and the interferogram without a background is obtained. Both the simulated and the experimental results demonstrate that the background removal method is feasible and robust with a high processing speed. In addition, this method can reduce the noise level of the reconstructed spectrum, and it is insusceptible to a complicated background, compared with the polynomial fitting and empirical mode decomposition(EMD) methods.