小波分析在信号去噪声中的应用

小波分析是近年来发展起来的一种数学方法,在信号与图象处理中有重要的应用.中值滤波是信号处理中常用的一种非线性滤波器,它能够有效地消除瞬时脉冲干扰,并且能够很好地保持信号的边缘信息,在信号和图象处理中得到广泛应用.对中值滤波器与小波变换的结合进行了比较系统的研究.通过实例说明中值滤波器与小波变换相结合具有比单一滤波器更好的效果.     

Dual integral equations related to the kontorovich-lebedev transform : PMM vol. 38, n 6, 1974, pp. 1090–1097

We study the dual integral equations related to the Kontorovich-Lebedev integral transforms arising in the course of solution of the problems of mathematical physics, in particular of the mixed boundary value problems for the wedge-shaped regions. We show that the solutions of these equations can be expressed in quadratures, using the auxilliary functions satisfying the integral Fredholm equation of second kind with a symmetric kernel.At present, the dual equations investigated in most detail are those connected with the Fourier and Hankel integral transforms. The results obtained and their applications are given in [1–3]. A large number of papers also deal with the theory and applications of the dual integral equations connected with the Mehler-Fock integral transform and its generalizations [4–11]., The dual integral transforms considered in the present paper belong to a more complex class than those listed above, and so far, no effective solution has been obtained for them. The only relevant results known to the authors are those in [12, 13]. In [12] a method of solving the equations (1.2) is given for a single particular value of the parameter γ = π/2, while in [13] the dual equations of the type under consideration are reduced to a solution of an infinite system of linear algebraic equations.     

正交小波变换的向量空间算法

本文揭示了一个事实，小波不仅可构成Ｌ２空间中的正交基，小波分解与重构滤波还可产生Ｎ维空间中的正交基．在本文提出修改的小波变换算法之下，Ｎ点信号的小波变换等价于Ｎ维空间中的正交变换．用该算法进行信号或图象压缩，无需对信号或图象进行周期延拓，可严格地在Ｎ维空间中进行．     

Beyond sumudu transform and natural transform: J-transform properties and applications

In this paper, we introduce an efficient integral transform called the ${\mathbb J}$-transform which is a modification of the well-known Sumudu transform and the Natural transform for solving differential equations with real applications in applied physical sciences and engineering. The ${\mathbb J}$-transform is more advantageous than both the Sumudu transform and the Natural transform. Interestingly, our proposed ${\mathbb J}$-transform can be applied successfully to solve complex problems that are ordinarily beyond the scope of either Sumudu transform or Natural transform. As a proof of concept, we consider some classic examples and highlight the limitations of the previously reported integral transforms and lastly demonstrate the superiority of the proposed ${\mathbb J}$-transform in addressing those limitations.     

小波尺度函数计算的广义高斯积分法及其应用

对于小波尺度函数变换的分解系数的积分运算建立了以尺度函数为权的广义高斯积分方法的运算格式．借助于样条函数,证明了其广义高斯积分随小波分解水平（resolutionlevel）指标的上升而收敛．在此基础上给出了以小波尺度函数变换重构或逼近任一函数的显式解析式,并对具有函数算子、微分或积分算子的运算给出了变换规则．这对于求解复杂非线性方程（组）是一种强有力的工具．最后给出了用该文方法求解非线性二点边值问题的算例．     

A bilateral asymptotic method of solving the integral equation of the contact problem of the torsion of an elastic half-space inhomogeneous in depth

An approximate semi-analytical method for solving integral equations generated by mixed problems of the theory of elasticity for inhomogeneous media is developed. An effective algorithm for constructing approximations of transforms of the kernels of integral equations by analytical expressions of a special type is proposed, and closed analytical solutions are presented. A comparative analysis of the approximation algorithms is given. The accuracy of the method is analysed using the example of the contact problem of the torsion of a medium with a non-uniform coating by a stiff circular punch. The relation between the error of the approximation of the transform of a kernel by special analytical expressions, constructed using different algorithms and the error of approximate solutions of the corresponding contact problems is investigated using a numerical experiment.     

A Partitioned ∈-Relaxation Algorithm for Separable Convex Network Flow Problems

A relaxation method for separable convex network flow problems is developed that is well-suited for problems with large variations in the magnitude of the nonlinear cost terms. The arcs are partitioned into two sets, one of which contains only arcs corresponding to strictly convex costs. The algorithm adjusts flows on the other arcs whenever possible, and terminates with primal-dual pairs that satisfy complementary slackness on the strictly convex arc set and -complementary slackness on the remaining arcs. An asynchronous parallel variant of the method is also developed. Computational results demonstrate that the method is significantly more efficient on ill-conditioned networks than existing methods, solving problems with several thousand nonlinear arcs in one second or less.     

Construction of Biorthogonal Discrete Wavelet Transforms Using Interpolatory Splines

We present a new family of biorthogonal wavelet and wavelet packet transforms for discrete periodic signals and a related library of biorthogonal periodic symmetric waveforms. The construction is based on the superconvergence property of the interpolatory polynomial splines of even degrees. The construction of the transforms is performed in a “lifting” manner that allows more efficient implementation and provides tools for custom design of the filters and wavelets. As is common in lifting schemes, the computations can be carried out “in place” and the inverse transform is performed in a reverse order. The difference with the conventional lifting scheme is that all the transforms are implemented in the frequency domain with the use of the fast Fourier transform. Our algorithm allows a stable construction of filters with many vanishing moments. The computational complexity of the algorithm is comparable with the complexity of the standard wavelet transform. Our scheme is based on interpolation and, as such, it involves only samples of signals and it does not require any use of quadrature formulas. In addition, these filters yield perfect frequency resolution.     

The Fourier-series method for inverting transforms of probability distributions

This paper reviews the Fourier-series method for calculating cumulative distribution functions (cdf's) and probability mass functions (pmf's) by numerically inverting characteristic functions, Laplace transforms and generating functions. Some variants of the Fourier-series method are remarkably easy to use, requiring programs of less than fifty lines. The Fourier-series method can be interpreted as numerically integrating a standard inversion integral by means of the trapezoidal rule. The same formula is obtained by using the Fourier series of an associated periodic function constructed by aliasing; this explains the name of the method. This Fourier analysis applies to the inversion problem because the Fourier coefficients are just values of the transform. The mathematical centerpiece of the Fourier-series method is the Poisson summation formula, which identifies the discretization error associated with the trapezoidal rule and thus helps bound it. The greatest difficulty is approximately calculating the infinite series obtained from the inversion integral. Within this framework, lattice cdf's can be calculated from generating functions by finite sums without truncation. For other cdf's, an appropriate truncation of the infinite series can be determined from the transform based on estimates or bounds. For Laplace transforms, the numerical integration can be made to produce a nearly alternating series, so that the convergence can be accelerated by techniques such as Euler summation. Alternatively, the cdf can be perturbed slightly by convolution smoothing or windowing to produce a truncation error bound independent of the original cdf. Although error bounds can be determined, an effective approach is to use two different methods without elaborate error analysis. For this purpose, we also describe two methods for inverting Laplace transforms based on the Post-Widder inversion formula. The overall procedure is illustrated by several queueing examples.     

The solution of infinitely ill-conditioned weakly-singular problems

It is demonstrated on examples that a weak singularity (i.e., with converging improper integral) may produce in computations (depending on the algorithm employed) an infinitely ill-conditioned situation when arbitrarily small imprécisions introduced by the algorithm or by a software create divergent approximations for mathematically convergent integrals. The possibility of hidden singularities is shown, and the double error phenomenon is identified and demonstrated in a simple example. Construction of test problems is proposed to check the applicability of existing software prior to its use for the solution of real life problems with weakly-singular equations. It is shown that the application of the integration by parts formula to weakly-singular integrals may create strong singularities (i.e., unbounded terms or divergent improper integrals). Methods of removal of singularities with and without compensation are studied for the numerical solution of infinitely ill-conditioned weakly-singular problems.     

Reconstruction from Line Integrals in Spaces of Constant Curvature

New reconstruction formula for the line integral transformation in Euclidean spaces is found. The general k-plane integral transform in Euclidean space is related to a totally geodesic integral transform for an arbitrary Riemannian space of constant curvature by means of a factorization property. Duality theorems for the totally geodesic transforms are stated.     

On the solution of mixed problems for a fiber with a vertical circular cylindrical cavity

A method of solving paired integral equations that appear in considerations of mixed problems of elasticity and thermoelasticity theory is given, with the help of generalized integral Weber transforms. The paired equations are reduced to an integral Fredholm equation of the second kind on the semiaxis, which have a discontinuous kernel, or to Fredholm equations of the second kind on a finite interval and infinite systems of linear algebraic equations, which are normal in the sense of Poincare-Koch. As an example, contact problems for an inhomogeneous fiber with a cavity are considered. If the fiber is bonded with the elastic half-space, then a second appproach is realized, which is based on a reduction to an equation with a self-adjoint operator, for which some method of sequential iteractions and the Bubnov-Galerkin method are justified.Translated from Dinamicheskie Sistemy, No. 7, pp. 95–102, 1988.     

On the existence of L1(R) solutions to wiener-Hopf type equations

By means of the Fourier transforms of distributions we find necessary and sufficient conditions for the existence of L¹(R) solutions to Wiener-Hopf type integral equations. Thus we establish general criteria for the existence of L¹(R) filters operating on the observed signal to best approximate the true signal. The theorems apply to wide sense stationary stochastic processes.     

Noise reduction method for vibration signals 2D time‐frequency distribution using anisotropic diffusion equation

Threshold noise reduction methods of vibration signals have been widely researched and used. However, these methods are less efficient in such situation, including requiring a time‐consuming and subjective to manual editing because different degree of noise signal requires selecting different characterization for filtering. In this paper, an efficient denoising method based on PDE for mechanical vibration signals time‐frequency distribution is investigated, in which, a one‐dimensional vibration signal is transformed into 2D time‐frequency domain by using Gabor transform. This enables (i) simultaneously utilize both time and frequency characteristic for effectively multiple dimension signal denosing and (ii) isotropic and anisotropic characteristics to be imposed by employing PDE, which explicitly fit with the local structure of time‐frequency signal. This paper analyzes the basic methods of isotropic and anisotropic diffusion filtering, investigates the anisotropic diffusion method based on local feature structure of 2D information, and conducts a set of comparative tests. Experiments show that this proposed method has a better performance of denoising than that of thresholding. At the same time, it is more handy than that of other methods, such as independent component analysis. Finally, problems and ways of improving the PDE‐based filter method are analyzed in this paper. Copyright © 2014 John Wiley & Sons, Ltd.     

Maslov-Poisson measure and Feynman formulas for the solution of the Dirac equation

As the main step, the method used by V. P. Maslov for representing a solution of the initial-value problem for the classical Schrödinger equation and admitting an application to the Dirac equation includes the construction of a cylindrical countably additive measure (that is an analog of the Poisson distribution) on a certain space of functions (= trajectories in the impulse space) whose Fourier transform coincides with the factor in the formula for representation of the solution of the Schrödinger equation by the integral in the so-called cylindrical Feynman (pseudo) measure (in the trajectory space of the configurational space for the classical system). On the other hand, in the Maslov formula for the solution of the Schrödinger equation, the exponential factor is (with accuracy up to a shift) the Fourier transform of the Feynman pseudomeasure. In the case of the Dirac equation, historically, for the first time, there arose formulas for the impulse representation that use countably additive functional distributions of the Poisson-Maslov measure type but with noncommuting (matrix) values. The paper finds generalized measures whose Fourier transforms coincide with an analog of the exponential factor under the integral sign in the Maslov-type formula for the Dirac equation and integrals with respect to which yield solutions of the Cauchy problem for this equation in the configurational space.     

The Class of Clifford-Fourier Transforms

Recently, there has been an increasing interest in the study of hypercomplex signals and their Fourier transforms. This paper aims to study such integral transforms from general principles, using 4 different yet equivalent definitions of the classical Fourier transform. This is applied to the so-called Clifford-Fourier transform (see Brackx et al., J. Fourier Anal. Appl. 11:669–681, 2005). The integral kernel of this transform is a particular solution of a system of PDEs in a Clifford algebra, but is, contrary to the classical Fourier transform, not the unique solution. Here we determine an entire class of solutions of this system of PDEs, under certain constraints. For each solution, series expressions in terms of Gegenbauer polynomials and Bessel functions are obtained. This allows to compute explicitly the eigenvalues of the associated integral transforms. In the even-dimensional case, this also yields the inverse transform for each of the solutions. Finally, several properties of the entire class of solutions are proven.     

Boundary value problems of electroelasticity for a plate with an inclusion and a half-space with a slit

Problems of determining the mechanical and electrical fields in a piezoelectric plate reinforced with an inclusion or in a half-space weakened by a cut are considered. Using the methods of the theory of analytic functions these problems are reduced to a system of singular integro-differential equations (for a plate) or to a singular integral equation with a fixed singularity (for a half-space). Approximate and exact solutions of the problems are obtained by the method of orthogonal polynomials and integral transforms.     

Plancherel Inversion as Unified Approach to Wavelet Transforms and Wigner Functions

We demonstrate that the Plancherel transform for Type-I groups provides one with a natural, unified perspective for the generalized continuous wavelet transform, on the one hand, and for a class of Wigner functions, on the other. We first prove that a Plancherel inversion formula, well known for Bruhat functions on the group, holds for a much larger class of functions. This result allows us to view the wavelet transform as essentially the inverse Plancherel transform. The wavelet transform of a signal is an L²-function on an appropriately chosen group while the Wigner function is defined on a coadjoint orbit of the group and serves as an alternative characterization of the signal, which is often used in practical applications. The Plancherel transform maps L²-functions on a group unitarily to fields of Hilbert-Schmidt operators, indexed by unitary irreducible representations of the group. The wavelet transform can essentially be looked upon as a restricted inverse Plancherel transform, while Wigner functions are modified Fourier transforms of inverse Plancherel transforms, usually restricted to a subset of the unitary dual of the group. Some known results on both Wigner functions and wavelet transforms, appearing in the literature from very different perspectives, are naturally unified within our approach. Explicit computations on a number of groups illustrate the theory. Communicated by Gian Michele Graf submitted 05/06/01, accepted: 19/09/02     

The Solution of Volterra Integral Equations of the Convolution Type Using Two-Point Rational Approximants

This paper applies the method of Laplace transform inversionusing two-point rational approximants, introduced by Grundy(1977), to the solution of Volterra integral equations of theconvolution type. The transforms which occur in solving theseequations using Laplace transform methods are, in many cases,of the type which can be inverted by the method introduced inGrundy (1977). There are however additional difficulties whichhave to be met. Examples are given to illustrate their resolutionand at the same time indicate the scope of the method.     

Axisymmetric Stokes' flow in a spherical shell revisited via the Fokas method. Part I: Irrotational flow

The axisymmetric irrotational Stokes' flow for a spherical shell is analysed by means of the recently developed Fokas method via the use of global relations. Alternative series and new integral representations concerning a system of concentric spheres, yielding, by a limiting procedure, the Dirichlet or Neumann problems for the interior and the exterior of a sphere, are presented. The boundary value problems considered can be classically solved using either the finite Gegenbauer transform or the Mellin transform. Application of the Gegenbauer transform yields a series representation which is uniformly convergent at the boundary, but not convenient for many applications. The Mellin transform, on the other hand, furnishes an integral representation which is not uniformly convergent at the boundary. Here, by algebraic manipulations of the global relation: (i) a Gegenbauer series representation is derived in a simpler manner, instead of solving ODEs and (ii) an alternative integral representation, different from the Mellin transform representation is derived which is uniformly convergent at the boundary. Copyright © 2011 John Wiley & Sons, Ltd.