QRS complex detection based on multi wavelet packet decomposition

We present in this paper a wavelet packet based QRS complex detection algorithm. Our proposed algorithm consists of a particular combination of two vectors obtained by applying a designed routine of QRS detection process using ‘haar’ and ‘db10’ wavelet functions respectively. The QRS complex detection routine is based on the histogram approach where our key idea was to search for the node with highest number of histogram coefficients, at center, which we assume that they are related to the iso-electric baseline whereas the remaining least number coefficients reflect the R waves peaks. Following a classical approach based of a calculated fixed threshold, the possible QRS complexes will be determined. The QRS detection complex algorithm has been applied to the whole MIT-BIH arrhythmia Database to assess its robustness. The algorithm reported a global sensitivity of 98.68%, positive predictive value of 97.24% and a percentage error of 04.12%. Eventhough, the obtained global results are not as excellent as expected, we have demonstrate that our designed QRS detection algorithm performs good on a partial selected high percentage of the whole database, e.g., the partial results, obtained when applying the algorithm on 85.01% of the whole MIT-BIH arrhythmia Database, are 99.14% of sensitivity, 98.94% of positive predictive value and 01.92% of percentage error.     

Inversion of the wavelet transform using Riemannian sums

We study the approximation of the inverse wavelet transform using Riemannian sums. For a large class of wavelet functions, we show that the Riemannian sums converge to the original function as the sampling density tends to infinity. When the analysis and synthesis wavelets are the same, we also give some necessary conditions for the Riemannian sums to be convergent.     

A direct integral decomposition of the wavelet representation

In this paper we use the concept of wavelet sets, as introduced by X. Dai and D. Larson, to decompose the wavelet representation of the discrete group associated to an arbitrary integer dilation matrix as a direct integral of irreducible monomial representations. In so doing we generalize a result of F. Martin and A. Valette in which they show that the wavelet representation is weakly equivalent to the regular representation for the Baumslag-Solitar groups.

     

Numerical stability of biorthogonal wavelet transforms

Biorthogonal wavelets are essential tools for numerous practical applications. It is very important that wavelet transforms work numerically stable in floating point arithmetic. This paper presents new results on the worst-case analysis of roundoff errors occurring in floating point computation of periodic biorthogonal wavelet transforms, i.e. multilevel wavelet decompositions and reconstructions. Both of these wavelet algorithms can be realized by matrix–vector products with sparse structured matrices. It is shown that under certain conditions the wavelet algorithms can be remarkably stable. Numerous tests demonstrate the performance of the results.      

Structural optimization by wavelet transforms and neural networks

An efficient methodology is presented to achieve optimal design of structures for earthquake loading. In this methodology a combination of wavelet transforms, neural networks and evolutionary algorithms are employed. The stochastic nature of the evolutionary algorithms makes the slow convergence. Specially, when earthquake induced loads are taken into account. To reduce the computational burden, a discrete wavelet transform is used by means of which the number of points in the earthquake record is decreased. Then, by using a surrogate model, the dynamic responses of the structures are predicted. In order to investigate the efficiency of the proposed methodology, two structures are designed for optimal weight. The numerical results demonstrate the computational advantages of the proposed hybrid methodology to optimal dynamic design of structures.     

The high frequency behaviour of continuous wavelet transforms

The high frequency behaviour of continuous wavelet transforms is characterized by the number of vanishing moments of the corresponding basic wavelets. As a consequence we give satisfying answers to the following questions of theoretical as well as practical interest:

What are the differences or similarities between transforms to different wavelets?

Why do wavelet transforms react so sensitively to abrupt signal changes ?     

The wavelet transform method applied to a coupled chaotic system with asymmetric coupling schemes

The wavelet transform method originated by Wei et al. (2002) [19] is an effective tool for enhancing the transverse stability of the synchronous manifold of a coupled chaotic system. Much of the theoretical study on this matter is centered on networks that are symmetrically coupled. However, in real applications, the coupling topology of a network is often asymmetric; see Belykh et al. (2006)  [23], [24], Chavez et al. (2005)  [25], Hwang et al. (2005)  [26], Juang et al. (2007)  [17], and Wu (2003)  [13]. In this work, a certain type of asymmetric sparse connection topology for networks of coupled chaotic systems is presented. Moreover, our work here represents the first step in understanding how to actually control the stability of global synchronization from dynamical chaos for asymmetrically connected networks of coupled chaotic systems via the wavelet transform method. In particular, we obtain the following results. First, it is shown that the lower bound for achieving synchrony of the coupled chaotic system with the wavelet transform method is independent of the number of nodes. Second, we demonstrate that the wavelet transform method as applied to networks of coupled chaotic systems is even more effective and controllable for asymmetric coupling schemes as compared to the symmetric cases.     

Short-run electricity load forecasting with combinations of stationary wavelet transforms

Short-term forecasting of electricity load is an essential issue for the management of power systems and for energy trading. Specific modeling approaches are needed given the strong seasonality and volatility in load data. In this paper, we investigate the benefit of combining stationary wavelet transforms to produce one day-ahead forecasts of half-hourly electric load in France. First, we assess the advantage of decomposing the aggregate load into several subseries with a wavelet transform. Each component is predicted separately and aggregated to get the final forecast. One innovation of this paper is to propose several approaches to deal with the boundary problem which is particularly detrimental in electricity load forecasting. Second, we examine the benefit of combining forecasts over individual models. An extensive out-of-sample evaluation shows that a careful treatment of the border effect is required in the multiresolution analysis. Combinations including the wavelet predictions provide the most accurate forecasts. This result is valid with several assumptions about the forecast error in temperature and for different types of hours (peak, normal, off-peak), different days of the week and various forecasting periods.     

Analysis of some injection bounds for Sobolev spaces by wavelet decomposition

We consider the Sobolev spaces $H^s(\Omega )$ and $H_0^s(\Omega )$ and the Besov spaces $B_{2,\infty }^{1/2}(\Omega )$, where Ω is a sufficiently regular (see Lemma 2) subdomain of ${\mathbb{R}}^2$. It is well known that for the values of $s\in [0,1/2)$ the two Sobolev spaces coincide, with equivalence of the norms, and that the inclusion $B_{2,\infty }^{1/2}(\Omega )?H^s(\Omega )$ holds. The Note is concerned with the explicit analysis of the constants appearing in the continuity bounds for the injections $H^s(\Omega )?H_0^s(\Omega )$ and $B_{2,\infty }^{1/2}(\Omega )?H^s(\Omega )$ and of their dependence on the regularity s of the spaces. The analysis is carried out by using the wavelet characterization of the corresponding norms.     

A novel approach to solve linear systems arising from BEM using fast wavelet transforms

This paper uses Daubechies orthogonal wavelets to change dense and fully populated matrices of boundary element method (BEM) systems into sparse and semi‐banded matrices. Then a novel algorithm based on hierarchical nature of multiresolution analysis is introduced to solving resultant sparse linear systems. This algorithm decomposes NS‐form of transformed parent matrix into descendant systems with reduced sizes and solves them iteratively using GMRES algorithm. Both parts, changing dense matrices to sparse systems and the novel solver, can be added as a black box to the existing BEM codes. Transforming matrices into wavelet space needs less time than saved by solving sparse large systems. Numerical results with a precise study on sensitivity of solution for physical variables to the thresholding parameter, and savings in computer time and memory are presented. Also, the suitable value for thresholding parameter is recommended for elasticity problems. The results indicate that the proposed method is efficient for large problems. Copyright © 2009 John Wiley & Sons, Ltd.     

Three-dimensional simulation of flat rolling through a combined finite element–boundary element approach

This paper presents an innovative approach for analysing three-dimensional flat rolling. The proposed approach is based on a solution resulting from the combination of the finite element method with the boundary element method. The finite element method is used to perform the rigid–plastic numerical modelling of the workpiece allowing the estimation of the roll separating force, rolling torque and contact pressure along the surface of the rolls. The boundary element method is applied for computing the elastic deformation of the rolls. The combination of the two numerical methods is made using the finite element solution of the contact pressure along the surface of the rolls to define the boundary conditions to be applied on the elastic analysis of the rolls. The validity of the proposed approach is discussed by comparing the theoretical predictions with experimental data found in the literature.     

Analysis of the stress singularity for a bi-material V-notch by the boundary element method

A new algorithm coupling the boundary element technique with the characteristic expansion method is proposed for the computation of the singular stress field in the V-notched bi-material structure. After the stress asymptotic expansions are introduced into the linear elasticity equilibrium equations, the governing equations at the small sector dug out from the bi-material V-notch tip region are transformed into the ordinary differential eigen-equations. All the parameters in the asymptotic expansions except the combination coefficients can be achieved by solving the established eigen-equations with the interpolating matrix method. Furthermore, the conventional boundary element method is applied to modeling the remaining structure without the notch tip region. The combination coefficients in the asymptotic expansion forms can be computed by the discretized boundary integral equations. Thus, the singular stress field at the V-notch tip and the generalized stress intensity factors of the bi-material notch are successfully calculated. The accurate singular stress field obtained here is very useful in the evaluation of the fracture property and the fatigue life of the V-notched bi-material structure.     

A parallel version of the preconditioned conjugate gradient method for boundary element equations

The parallel version of precondition techniques is developed for matrices arising from the Galerkin boundary element method for two-dimensional domains with Dirichlet boundary conditions. Results were obtained for implementations on a transputer network as well as on an nCUBE-2 parallel computer showing that iterative solution methods are very well suited for a MIMD computer. A comparison of numerical results for iterative and direct solution methods is presented and underlines the superiority of iterative methods for large systems.     

Some remarks on “On the windowed Fourier transform and wavelet transform of almost periodic functions,” by J.R. Partington and B. Ünalmı

J.R. Partington and B. Ünalmı consider in their 2001 paper [J.R. Partington, B. Ünalmı , Appl. Comput. Harmon. Anal. 10 (1) (2001) 45–60] the windowed Fourier transform and wavelet transform as tools for analyzing almost periodic signals. They establish Parseval-type identities and consider discretized versions of these transforms in order to construct generalized frame decompositions. We have found a gap in the construction of generalized frames in the windowed Fourier transform case; we comment on this gap and give an alternative proof. As for the wavelet transform case, in [J.R. Partington, B. Ünalmı , Appl. Comput. Harmon. Anal. 10 (1) (2001) 45–60] the generalized frame decomposition is done only for the simplest wavelet, the Haar wavelet; we show how to construct generalized frame decompositions for a wide family of wavelets.     

Efficiency of a posteriori BEM--error estimates for first-kind integral equations on quasi--uniform meshes

In the numerical treatment of integral equations of the first kind using boundary element methods (BEM), the author and E. P. Stephan have derived a posteriori error estimates as tools for both reliable computation and self-adaptive mesh refinement. So far, efficiency of those a posteriori error estimates has been indicated by numerical examples in model situations only. This work affirms efficiency by proving the reverse inequality. Based on best approximation, on inverse inequalities and on stability of the discretization, and complementary to our previous work, an abstract approach yields a converse estimate. This estimate proves efficiency of an a posteriori error estimate in the BEM on quasi--uniform meshes for Symm's integral equation, for a hypersingular equation, and for a transmission problem.

     

Boundary element bending analysis of heated elastic plates

Generally, the boundary integral equations for nonlinear and inhomogeneous problems need to be formulated in terms of boundary and domain integrals. The present paper declares the possibility of the transformation of the domain integral term into the boundary integral for the thermal bending analysis of a thin elastic plate,provided that the temperature can be expressed by a linear variation over the thickness.     

The wavelet multiscale method for inversion of porosity in the fluid-saturated porous media

In this paper, the wavelet multiscale method is applied to the inversion of porosity in the fluid-saturated porous media. The inverse problem is decomposed to multiple scales with wavelet transform and hence the original inverse problem is re-formulated to be a set of sub-inverse problem corresponding to different scales and is solved successively according to the size of scale from the smallest to the largest. On each scale, regularization Gauss–Newton method is carried out, which is stable and fast, until the optimum solution of original inverse problem is found. The results of numerical simulations demonstrate that the method is a widely convergent optimization method and exhibits the advantages of conventional regularization Gauss–Newton method methods on computational efficiency and precision.     

Numerical solution of the elastic body-plate problem by nonoverlapping domain decomposition type techniques

The purpose of this paper is to provide two numerical methods for solving the elastic body-plate problem by nonoverlapping domain decomposition type techniques, based on the discretization method by Wang. The first one is similar to an older method, but here the corresponding Schur complement matrix is preconditioned by a specific preconditioner associated with the plate problem. The second one is a ``displacement-force' type Schwarz alternating method. At each iteration step of the two methods, either a pure body or a pure plate problem needs to be solved. It is shown that both methods have a convergence rate independent of the size of the finite element mesh.

     

Interaction of multiple cracks in a rectangular plate

This paper is concerned with interaction of multiple cracks in a finite plate by using the hybrid displacement discontinuity method (a boundary element method). Detail solutions of the stress intensity factors (SIFs) of the multiple-crack problems in a rectangular plate are given, which can reveal the effect of geometric parameters of the cracked body on the SIFs. The numerical results reported here illustrate that the boundary element method is simple, yet accurate for calculating the SIFs of multiple crack problems in a finite plate.