正弦波经验模式分解中极点拟合的抛物线插值法

本文针对正弦渡用经验模式分解（EMD）进行分析时存在的问题,从实验上结出了采样率和信号间相互作用对END的影响。为了提高对极值点采样的准确性,本文提出了利用抛物线插值来拟舍极值点的方法,使结果得到明显的改善。本文还对信号间相互作用的问题进行了讨论。     

Blow-up phenomena and local well-posedness for a generalized Camassa–Holm equation in the critical Besov space

In this paper we mainly study the Cauchy problem for a generalized Camassa–Holm equation in a critical Besov space. First, by using the Littlewood–Paley decomposition, transport equations theory, logarithmic interpolation inequalities and Osgood’s lemma, we establish the local well-posedness for the Cauchy problem of the equation in the critical Besov space $$B^{\frac{1}{2}}_{2,1}$$. Next we derive a new blow-up criterion for strong solutions to the equation. Then we give a global existence result for strong solutions to the equation. Finally, we present two new blow-up results and the exact blow-up rate for strong solutions to the equation by making use of the conservation law and the obtained blow-up criterion.     

Characterizations of local upper Lipschitz property of perturbed solutions to nonlinear second-order cone programs

We characterize the local upper Lipschitz property of the stationary point mapping and the Karush–Kuhn–Tucker (KKT) mapping for a nonlinear second-order cone programming problem using the graphical derivative criterion. We demonstrate that the second-order sufficient condition and the strict constraint qualification are sufficient for the local upper Lipschitz property of the stationary point mapping and are both sufficient and necessary for the local upper Lipschitz property of the KKT mapping.     

The new test criterion for the homogeneity of parameters of several populations

Summary The new test criterion for testing the homogeneity of parameters of several populations is proposed and the test properties of it is discussed. The asymptotic expansions of the distributions of test criterion are discussed under (i) null hypothesis, (ii) fixed alternative hypothesis and (iii) local alternative hypothesis converging to the null hypothesis with appropriate rate of convergence as the sample size increases. As a particular case the asymptotic theory of a statistic for a homogeneity of variances of normal populations is also discussed and the exact moments of it under a null hypothesis can be used to obtain a percentage point by a Pearsonian curve fitting. This Institute of statistical Mathematics     

A criterion on a repeller being a null set of any limit measure for stochastic differential equations

This paper studies limit behaviors of stationary measures for stochastic ordinary differential equations with nondegenerate noise and presents a criterion to guarantee that a repeller with zero Lebesgue measure is a null set of any limit measure. Using this criterion, we first provide a series of nontrivial concrete examples to show that their repelling limit cycles or quasi-periodic orbits are null sets for all limit measures, which deduces that all their limit measures are concentrated on stable equilibria and stable limit cycles or quasi-periodic orbits,and saddles. Interesting open questions on exact supports of limit measures are proposed.     

基于EMD方法的混沌信号中周期分量的提取

提出一种从Duffing振子产生的混沌信号中提取谐波分量的方法．依据任何信号由不同的固有简单振动模态组成的概念,利用经验模式分解（EMD）方法,将混沌信号分离为不同的内在模态函数(IMF),并在特定参数下从中分解出单一频率成分的谐波信号,从而成功地将混沌信号和谐波分量分离．仿真实验都表明该方法非常有效．     

Quantum exotic PDE’s

Following the previous works on the Prástaro’s formulation of algebraic topology of quantum (super) PDE’s, it is proved that a canonical Heyting algebra (integral Heyting algebra) can be associated to any quantum PDE. This is directly related to the structure of its global solutions. This allows us to recognize a new inside in the concept of quantum logic for microworlds. Furthermore, the Prástaro’s geometric theory of quantum PDE’s is applied to the new category of quantum hypercomplex manifolds, related to the well-known Cayley–Dickson construction for algebras. Theorems of existence for local and global solutions are obtained for (singular) PDE’s in this new category of noncommutative manifolds. Finally, the extension of the concept of exotic PDE’s, recently introduced by Prástaro, has been extended to quantum PDE’s. Then a smooth quantum version of the quantum (generalized) Poincaré conjecture is given too. These results extend ones for quantum (generalized) Poincaré conjecture, previously given by Prástaro.     

Local mass conservative Hermite interpolation for the solution of flow problems by a multi-domain boundary element approach

This paper presents a local Hermite radial basis function interpolation scheme for the velocity and pressure fields. The interpolation for velocity satisfies the continuity equation (mass conservative interpolation) while the pressure interpolation obeys the pressure equation. Additionally, the Dual Reciprocity Boundary Element method (DRBEM) is applied to obtain an integral representation of the Navier-Stokes equations. Then, the proposed local interpolation is used to obtain the values of the field variables and their partial derivatives at the boundary of the sub-domains. This interpolation allows one to obtain the boundary values needed for the integral formulas for velocity and pressure at some nodes within the sub-domains. In the proposed approach the boundary elements are merely used to parameterize the geometry, but not for the evaluation of the integrals as it is usually done. The presented multi-domain approach is different from the traditional ones in boundary elements because the resulting integral equations are non singular and the boundary data needed for the boundary integrals are approximated using a local interpolation. Some accurate results for simple Stokes problems and for the Navier-Stokes equations at low Reynolds numbers up to Re = 400 were obtained.     

Null‐field approach for Laplace problems with circular boundaries using degenerate kernels

In this article, a semi‐analytical method for solving the Laplace problems with circular boundaries using the null‐field integral equation is proposed. The main gain of using the degenerate kernels is to avoid calculating the principal values. To fully utilize the geometry of circular boundary, degenerate kernels for the fundamental solution and Fourier series for boundary densities are incorporated into the null‐field integral equation. An adaptive observer system is considered to fully employ the property of degenerate kernels in the polar coordinates. A linear algebraic system is obtained without boundary discretization. By matching the boundary condition, the unknown coefficients can be determined. The present method can be seen as one kind of semianalytical approaches since error only attributes to the truncated Fourier series. For the eccentric case, vector decomposition technique for the normal and tangential directions is carefully considered in implementing the hypersingular equation in mathematical essence although we transform it to summability to divergent series. The five advantages, well‐posed linear algebraic system, principal value free, elimination of boundary‐layer effect, exponential convergence, and mesh free, are achieved. Several examples involving infinite, half‐plane, and bounded domains with circular boundaries are given to demonstrate the validity of the proposed method. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

Null space correction and adaptive model order reduction in multi-frequency Maxwell’s problem

A model order reduction method is developed for an operator with a non-empty null-space and applied to numerical solution of a forward multi-frequency eddy current problem using a rational interpolation of the transfer function in the complex plane. The equation is decomposed into the part in the null space of the operator, calculated exactly, and the part orthogonal to it which is approximated on a low-dimensional rational Krylov subspace. For the Maxwell’s equations the null space is related to the null space of the curl. The proposed null space correction is related to divergence correction and uses the Helmholtz decomposition. In the case of the finite element discretization with the edge elements, it is accomplished by solving the Poisson equation on the nodal elements of the same grid. To construct the low-dimensional approximation we adaptively choose the interpolating frequencies, defining the rational Krylov subspace, to reduce the maximal approximation error. We prove that in the case of an adaptive choice of shifts, the matrix spanning the approximation subspace can never become rank deficient. The efficiency of the developed approach is demonstrated by applying it to the magnetotelluric problem, which is a geophysical electromagnetic remote sensing method used in mineral, geothermal, and groundwater exploration. Numerical tests show an excellent performance of the proposed methods characterized by a significant reduction of the computational time without a loss of accuracy. The null space correction regularizes the otherwise ill-posed interpolation problem.     

Abstract robust coarse spaces for systems of PDEs via generalized eigenproblems in the overlaps

Coarse spaces are instrumental in obtaining scalability for domain decomposition methods for partial differential equations (PDEs). However, it is known that most popular choices of coarse spaces perform rather weakly in the presence of heterogeneities in the PDE coefficients, especially for systems of PDEs. Here, we introduce in a variational setting a new coarse space that is robust even when there are such heterogeneities. We achieve this by solving local generalized eigenvalue problems in the overlaps of subdomains that isolate the terms responsible for slow convergence. We prove a general theoretical result that rigorously establishes the robustness of the new coarse space and give some numerical examples on two and three dimensional heterogeneous PDEs and systems of PDEs that confirm this property.     

A local hermitian RBF meshless numerical method for the solution of multi‐zone problems

The local Hermitian interpolation (LHI) method is a strong‐form meshless numerical technique in which the solution domain is covered by a series of small and heavily overlapping radial basis function (RBF) interpolation systems. Aside from its meshless nature and the ability to work on very large scattered datasets, the main strength of the LHI method lies in the formation of local interpolations, which themselves satisfy both boundary and governing PDE operators, leading to an accurate and stable reconstruction of partial derivatives without the need for artificial upwinding or adaptive stencil selection. In this work, an extension is proposed to the LHI formulation which allows the accurate capture of solution profiles across discontinuities in governing equation parameters. Continuity of solution value and mass flux is enforced between otherwise disconnected interpolation systems, at the location of the discontinuity. In contrast to other local meshless methods, due to the robustness of the Hermite RBF formulation, it is possible to impose both matching conditions simultaneously at the interface nodes. The procedure is demonstrated for 1D and 3D convection–diffusion problems, both steady and unsteady, with discontinuities in various PDE properties. The analytical solution profiles for these problems, which experience discontinuities in their first derivatives, are replicated to a high degree of accuracy. The technique has been developed as a tool for solving flow and transport problems around geological layers, as experienced in groundwater flow problems. The accuracy of the captured solution profiles, in scenarios where the local convective velocities exceed those typically encountered in such Darcy flow problems, suggests that the technique is indeed suitable for modeling discontinuities in porous media properties. © 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 27: 1201–1230, 2011     

Are the eigenvalues of preconditioned banded symmetric Toeplitz matrices known in almost closed form?

The problem of computing oscillatory integrals with general oscillators is considered. We employ a Filon-type method, where the interpolation basis functions are chosen in such a way that the moments are in terms of elementary functions and the oscillator only. This allows us to evaluate the moments rapidly and easily without needing to engage hypergeometric functions. The proposed basis functions form a Chebyshev set for any oscillator function even if it has some stationary points in the integration interval. This property enables us to employ the Filon-type method without needing any information about the stationary points if any. Interpolation by the proposed basis functions at the Fekete points (which are known as nearly optimal interpolation points), when combined with the idea of splines, leads to a reliable convergent method for computing the oscillatory integrals. Our numerical experiments show that the proposed method is more efficient than the earlier ones with the same advantages.     

R-D optimized tree-structured compression algorithms with discrete directional wavelet transform

A new image coding method based on discrete directional wavelet transform (S-WT) and quad-tree decomposition is proposed here. The S-WT is a kind of transform proposed in [V. Velisavljevic, B. Beferull-Lozano, M. Vetterli, P.L. Dragotti, Directionlets: anisotropic multidirectional representation with separable filtering, IEEE Trans. Image Process. 15(7) (2006)], which is based on lattice theory, and with the difference with the standard wavelet transform is that the former allows more transform directions. Because the directional property in a small region is more regular than in a big block generally, in order to sufficiently make use of the multidirectionality and directional vanishing moment (DVM) of S-WT, the input image is divided into many small regions by means of the popular quad-tree segmentation, and the splitting criterion is on the rate-distortion sense. After the optimal quad-tree is obtained, by means of the embedded property of SPECK, a resource bit allocation algorithm is fast implemented utilizing the model proposed in [M. Rajpoot, Model based optimal bit allocation, in: IEEE Data Compression Conference, 2004, Proceedings, DCC 2004.19]. Experiment results indicate that our algorithms perform better compared to some state-of-the-art image coders.     

A new smoothed aggregation multigrid method for anisotropic problems

A new prolongator is proposed for smoothed aggregation (SA) multigrid. The proposed prolongator addresses a limitation of standard SA when it is applied to anisotropic problems. For anisotropic problems, it is fairly standard to generate small aggregates (used to mimic semi‐coarsening) in order to coarsen only in directions of strong coupling. Although beneficial to convergence, this can lead to a prohibitively large number of non‐zeros in the standard SA prolongator and the corresponding coarse discretization operator. To avoid this, the new prolongator modifies the standard prolongator by shifting support (non‐zeros within a prolongator column) from one aggregate to another to satisfy a specified non‐zero pattern. This leads to a sparser operator that can be used effectively within a multigrid V‐cycle. The key to this algorithm is that it preserves certain null space interpolation properties that are central to SA for both scalar and systems of partial differential equations (PDEs). We present two‐dimensional and three‐dimensional numerical experiments to demonstrate that the new method is competitive with standard SA for scalar problems, and significantly better for problems arising from PDE systems. Copyright © 2008 John Wiley & Sons, Ltd.     

Initial transient detection in simulations using the second-order cumulant spectrum

Procedures for detecting an initial transient in simulation output data are developed. The tests use the second-order cumulant spectrum which differs from the power spectrum in that the stationarity constraint is not required for the former. The second-order cumulant spectrum can be interpreted as the nonstationary power spectrum and is an orthogonal decomposition of the variance of a nonstationary process. The null hypothesis is that the simulation output data series is a covariance stationary process. Equivalently, all estimates of the second-order cumulant spectrum in the region which excludes the estimates of the power spectrum will have an expected value of zero. The test procedures are designed to detect initialization bias in the estimation of the mean and the variance. These procedures can be extended to detect bias in the moments of cumulants of ordern, wheren>2. Results are presented from the application of the test to simulated processes with superimposed mean and variance transients and anM/M/1 queue example.     

闭区间[0,1]上的分形插值函数的分数阶积分

介绍了分形插值函数和迭代函数系统以及v阶黎曼-刘准尔分数阶积分的概念及相关定理.利用这些概念及定理讨论了分形插值函数的分数阶积分在[0,1]上连续性及判定[0,1]上的分形插值函数的分数阶积分也是[0,1]上的分形插值函数,并给予了证明.     

A new compound faults detection method for rolling bearings based on empirical wavelet transform and chaotic oscillator

The rolling bearings often suffer from compound faults in practice. The concurrence of different faults increases the fault detection difficulty and the decoupling detection of compound faults is attracting considerable attentions. Recent publications report the application of the multiwavelets and empirical mode decomposition (EMD) for compound faults decoupling. However, due to limited adaptability they would induce mode mixing or/and overestimation problems in the signal processing. Particularly, the mode mixing would greatly degrade their performance on compound faults detection. To address this issue, this work presents a new method based on the empirical wavelet transform-duffing oscillator (EWTDO) for compound faults decoupling diagnosis of rolling bearings. The empirical wavelet transform (EWT) is able to extract intrinsic modes of a signal by fully adaptive wavelet basis. Hence, the mode mixing and overestimation can be resolved in decoupling processing and the compound faults can be correctly decomposed into different single faults in the form of empirical modes. Then, each single fault frequency was incorporated into a duffing oscillator to establish its corresponding fault isolator. By directly observing the chaotic motion from the Poincar mapping of the isolator outputs the single faults were identified one by one from the empirical modes. Experimental tests were carried out on a rolling bearing fault tester to examine the efficacy of the proposed EWTDO method on compound faults detection. The analysis results show attractive performance with respect to existing decoupling approaches based on the multiwavelets and EMD. In particular, our proposed method is much more reliable in decoupling the compound faults. Hence, the proposed method has practical importance in compound faults decoupling diagnosis for rolling bears.     

关于三次插值样条函数的局部性质

近年来,三次插值样条函数的局部性质,受到了人们的注意,翁祖荫考察了在一般边界条件下三次插值样条的局部逼近性质及插值样条导数的局部界值,不过他对于分