Discontinuous Legendre wavelet element method for elliptic partial differential equations

By incorporating the Legendre multiwavelet into the discontinuous Galerkin (DG) method, this paper presents a novel approach for solving Poisson’s equation with Dirichlet boundary, which is known as the discontinuous Legendre multiwavelet element (DLWE) method, derive an adaptive algorithm for the method, and estimate the approximating error of its numerical fluxes. One striking advantage of our method is that the differential operator, boundary conditions and numerical fluxes involved in the elementwise computation can be done with lower time cost. Numerical experiments demonstrate the validity of this method. Furthermore, this paper generalizes the DLWE method to the general elliptic equations defined on a bounded domain and describes the possibilities of constructing optimal adaptive algorithm. The proposed method and its generalizations are also applicable to some other kinds of partial differential equations.     

L~2(R~n)上的半正交多小波框架

本文研究L2(Rn)上伸缩矩阵A满足|detA|1的半正交多小波框架.本文得到半正交和严格半正交框架的一系列性质及刻画.本文证明半正交Parseval多小波框架与广义多分辨分析(GMRA)Parseval多小波框架是等价的.特别地,本文利用最小频率支撑(MSF)多小波框架和小波集,构造若干半正交多小波框架的例子.     

Parametrization of balanced multiwavelet

This paper deals with the parametrization of balanced multiwavelets and different properties associated with them. We introduce the property balancing symmetry and orthogonal properties of multiwavelet and link these properties to the matrix of the low pass synthesis multifilter. Using these new results, we present the parametrization of orthogonal multiwavelets of flip-symmetry with length two and three. This is a direct construction method, making the construction of the balanced multiwavelet as easy as the scalar wavelet.     

对称-反对称多小波的构造

本文从多小波的多分辨分析出发 ,给出了一种由双正交单小波构造对称 反对称多小波的新方法 ,并以传统 9 7单小波为例构造了二重对称 反对称多小波 .     

多尺度函数与多小波的构造

在多小波和单小波的基础上利用矩阵卷积构造出了一类多尺度函数与多小波,并通过实例对构造算法加以说明.     

基于广义交叉认证的多小波阈值的图像降噪

提出一种新的小波收缩阈值降噪方法,该方法是通过对噪声图像进行多小波变换,然后用广义交叉认证的方法来确定小波阈值参数．由于本文采用的是多小波变换,而多小波一般同时具有正交性和线性相位,另外广义交叉认证方法不需要对噪声的强度进行估计,所以这种方法有比较好的降噪效果．实验结果表明该方法与基于小波变换的广义交叉认证的图像降噪方法相比较,其降噪效果有一定的提高;同时也表明在一定的条件下,其降噪效果要明显好于传统的Wiener滤波方法．     

正交共轭滤波器的构造

It is very importent for generating an orthonormal multiwavelet system to construct a conjugate quadrature filter(CQF). In this paper, a general method of deriving a length-L 1 conjugate quadrature filter from a length-L conjugate quadrature filter and vice versa is obtained. As a special case, we study generally the construction of any length-L 1 compactly supported symmetric-antisymmetric orthonormal multiwavelet system with multiplicity 2 from a length-L multiwavelet system and vice versa. Examples of conjugate quadrature filter are given.     

多小波框架的构造理论

本文给出了多小波框架的sub-QMF条件,提出了多小波框架低通滤波器的参数化设计,由正交分解和矩阵的酉扩张得到其相应的高通滤波器表示的整套多小波框架设计的参数化方法,同时针对多描述编码的需求,构造了两个长折叠对称带参数的多小波紧框架.     

求解第一类积分方程的正则化—小波方法及其数值试验

1 方法的描述 第一类(Fredholm)积分方程是指形如 (1.1)的积分方程,其中核k(x,y)和右端函数f(x)给定,u(x)是未知函数.许多物理、化学、力学和工程应用问题都能导致第一类积分方程.求解第一类积分方程的一个本质性困难是方程的不适定性,即解的存在性、唯一性和稳定性遭到破坏.常用的数值方法有奇异值分解(SVD)方法、Tikhonov正则化方法、投影方法、正则化-样条方法、再生核方法等.本文提出一种新的正则化-小波方法,在第一类积分方程有多个解时,可以求出具有最小范数的数值解;如果原积分方程有唯一解,则所得的数值解收敛于准确解.数值试验表明,该方法是可行的. 我们在L~2[a,b]中考虑第一类(Fredholm)积分方程,即假设方程(1.1)中积分算子K∈L~2([a,b]×[a,b])及右端f(x)∈L~2[a,b]给定.为保证数值求解算法的稳定性,我们先用正则化方法处理该方程,将不适定问题化为泛函极值问题来求解,然后利用多重正交样条小波基构造求解格式.由于我们给出了直接计算低阶的多重正交样条小波基函数的一般公式,使得解法可以在计算机迅速实现.     

The comparison of two reliable methods for accurate solution of time‐fractional Kaup–Kupershmidt equation arising in capillary gravity waves

In this paper, we compared two different methods, one numerical technique, viz Legendre multiwavelet method, and the other analytical technique, viz optimal homotopy asymptotic method (OHAM), for solving fractional‐order Kaup–Kupershmidt (KK) equation. Two‐dimensional Legendre multiwavelet expansion together with operational matrices of fractional integration and derivative of wavelet functions is used to compute the numerical solution of nonlinear time‐fractional KK equation. The approximate solutions of time fractional Kaup–Kupershmidt equation thus obtained by Legendre multiwavelet method are compared with the exact solutions as well as with OHAM. The present numerical scheme is quite simple, effective, and expedient for obtaining numerical solution of fractional KK equation in comparison to analytical approach of OHAM. Copyright © 2016 John Wiley & Sons, Ltd.     

Matrix Thresholding for Multiwavelet Image Denoising

Vector thresholding is a recently proposed technique for the denoising of one-dimensional signals by means of multiwavelet shrinkage. It is more suited both to dealing with the multiwavelet vector coefficients and to taking into account the correlations which can be introduced among the starting vector coefficients by the use of a suitable prefilter. Motivated by the successful results of the multiwavelet transform when used in image processing, the aim of this paper is to extend vector thresholding to the two-dimensional case by introducing the notion of matrix thresholding. This new method allows us to easily exploit the matrix nature of the two-dimensional multiwavelet transform, and represents the natural extension of vector thresholding to the 2-D case. Afterwards, as the choice of the threshold level is very important in the practical application of thresholding methods, we propose a first attempt to extend the recently introduced method of H-curve to a multiple wavelet setting. The results of extensive numerical simulations confirm the effectiveness of our proposals and encourage us to keep going in this direction with further studies.     

求解间断系数椭圆型问题的一种改进的DG方法

本文考虑对间断系数椭圆型问题的普通DG方法进行改进,提出了一种综合了DG方法及区域分解方法的优点的新方法.对此法进行了先验误差分析并给出其残量型后验误差估计,且通过数值实验验证了该方法及其自适应方法的有效性.     

Extension Principles for Dual Multiwavelet Frames of $$L_2(\mathbb {R}^s)$$ constructed from Multirefinable Generators

In this work we prove that any pair of homogeneous dual multiwavelet frames of \(L_2(\mathbb {R}^s)\) constructed from a pair of refinable function vectors gives rise to a pair of nonhomogeneous dual multiwavelet frames and vice versa. We also prove that the Mixed Oblique Extension Principle characterizes dual multiwavelet frames. Our results extend recent characterizations of affine dual frames derived from scalar refinable functions obtained in [3].     

Legendre multiwavelet Galerkin method for solving the hyperbolic telegraph equation

Recently, it is found that telegraph equation is more suitable than ordinary diffusion equation in modeling reaction diffusion for such branches of sciences. In this article a numerical method for solving the one‐dimensional hyperbolic telegraph equation is presented. The method is based upon Legendre multiwavelet approximations. The properties of Legendre multiwavelet are first presented. These properties together with Galerkin method are then utilized to reduce the telegraph equation to the solution of algebraic equations. Illustrative examples are included to demonstrate the validity and applicability of the technique. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

Construction of two-direction tight wavelet frames

We investigate the construction of two-direction tight wavelet frames First, a sufficient condition for a two-direction refinable function generating two-direction tight wavelet frames is derived. Second, a simple constructive method of two-direction tight wavelet frames is given. Third, based on the obtained two-direction tight wavelet frames, one can construct a symmetric multiwavelet frame easily. Finally, some examples are given to illustrate the results.     

Construction of orthogonal multiwavelet bases

We propose a method for constructing orthogonal multiwavelet bases of the space L ²(?) for any known multiscaling functions that generate a multiresolution analysis of dimension greater than 1.     

Symmetric–Antisymmetric Orthonormal Multiwavelets and Related Scalar Wavelets

For compactly supported symmetric–antisymmetric orthonormal multiwavelet systems with multiplicity 2, we first show that any length-2Nmultiwavelet system can be constructed from a length-(2N+1) multiwavelet system and vice versa. Then we present two explicit formulations for the construction of multiwavelet functions directly from their associated multiscaling functions. This is followed by the relationship between these multiscaling functions and the scaling functions of related orthonormal scalar wavelets. Finally, we present two methods for constructing families of symmetric–antisymmetric orthonormal multiwavelet systems via the construction of the related scalar wavelets.     

Statistically based multiwavelet denoising

In this work, we consider a statistically based multiwavelet thresholding method which acts on the empirical wavelet coefficients in groups, rather than individually, in order to obtain an edge-preserving image denoising technique. Our strategy allows us to exploit the dependencies between neighboring coefficients to make a simultaneous thresholding decision, so that estimation accuracy is increased.

By interpreting the multiwavelet analysis in a statistical context, we propose a new weighted multiwavelet matrix thresholding rule, based on the statistical modeling of empirical coefficients. This allows the thresholding decision to be adapted to the local structure of the underlying image, hence producing edge-preserving denoising. Extensive numerical results are presented showing the performance of our denoising procedure.     

Analysis of the parareal approach based on discontinuous Galerkin method for time-dependent Stokes equations

This paper analyzes a parareal approach based on discontinuous Galerkin (DG) method for the time-dependent Stokes equations. A class of primal discontinuous Galerkin methods, namely variations of interior penalty methods, are adopted for the spatial discretization in the parareal algorithm (we call it parareal DG algorithm). We study three discontinuous Galerkin methods for the time-dependent Stokes equations, and the optimal continuous in time error estimates for the velocities and pressure are derived. Based on these error estimates, the proposed parareal DG algorithm is proved to be unconditionally stable and bounded by the error of discontinuous Galerkin discretization after a finite number of iterations. Finally, some numerical experiments are conducted which confirm our theoretical results, meanwhile, the efficiency of the parareal DG algorithm can be seen through a parallel experiment.     

Global Support of a Scaling Vector

Multiwavelet decompositions are based on scaling vectors satisfying matrix refinement equations. The support and linear independence of scaling vectors play an essential role in the study of multiwavelets. In this paper we relate these properties with the coefficients in the matrix refinement equation satisfied by the scaling vector.