This article aims at studying two-direction refinable functions and two-direction wavelets in the setting R^s, s 〉 1. We give a sufficient condition for a two-direction refinable function belonging to L^2(R^s). The...This article aims at studying two-direction refinable functions and two-direction wavelets in the setting R^s, s 〉 1. We give a sufficient condition for a two-direction refinable function belonging to L^2(R^s). Then, two theorems are given for constructing biorthogonal (orthogonal) two-direction refinable functions in L^2(R^s) and their biorthogonal (orthogonal) two-direction wavelets, respectively. From the constructed biorthogonal (orthogonal) two-direction wavelets, symmetric biorthogonal (orthogonal) multiwaveles in L^2(R^s) can be obtained easily. Applying the projection method to biorthogonal (orthogonal) two-direction wavelets in L^2(R^s), we can get dual (tight) two-direction wavelet frames in L^2(R^m), where m ≤ s. From the projected dual (tight) two-direction wavelet frames in L^2(R^m), symmetric dual (tight) frames in L^2(R^m) can be obtained easily. In the end, an example is given to illustrate theoretical results.展开更多
Based on the discussion about working mechanism of horizontal reinforcement and that of vertical reinforcement,respectively,the working mechanism of two-direction reinforced composite foundation was studied.The enhanc...Based on the discussion about working mechanism of horizontal reinforcement and that of vertical reinforcement,respectively,the working mechanism of two-direction reinforced composite foundation was studied.The enhancing effect of horizontal reinforcement on vertical reinforced composite foundation was analyzed.A simplified calculation method for such two-direction reinforced working system was presented.A model experiment was carried out to validate the proposed method.In the experiment,geocell reinforcement worked as the horizontal reinforcement,while gravel pile composite foundation worked as the vertical reinforcement.The results show that the calculated curve is close to the measured one.The installation of geosynthetic reinforcement can increase the bearing capacity of composite foundation by nearly 68% at normal foundation settlement,which suggests that the enhancing effect by geosynthetic reinforcement should be taken into account in current design/analysis methods.展开更多
A new method for iris identification based on multiwavelets is proposed. By means of the properties of multiwavelets, such as orthogonality, symmetry, vanishing moments and approximation order, the iris texture can be...A new method for iris identification based on multiwavelets is proposed. By means of the properties of multiwavelets, such as orthogonality, symmetry, vanishing moments and approximation order, the iris texture can be simply presented. A brief overview of muhiwavelets is presented at first. Iris identification system and iris texture feature presentation and recognition based on multiwavelets a,e introduced subsequently. And the experiment indicates the validity of this method finally.展开更多
Assuming the material properties varying with an exponential law both in the thick- ness and radial directions, axisymmetric bending of two-directional functionally graded circular and annular plates is studied using ...Assuming the material properties varying with an exponential law both in the thick- ness and radial directions, axisymmetric bending of two-directional functionally graded circular and annular plates is studied using the semi-analytical numerical method in this paper. The deflections and stresses of the plates are presented. Numerical results show the well accuracy and convergence of the method. Compared with the finite element method, the semi-analytical nu- merical method is with great advantage in the computational efficiency. Moreover, study on ax- isymmetric bending of two-directional functionally graded annular plate shows that such plates have better performance than those made of isotropic homogeneous materials or one-directional functionally graded materials. Two-directional functionally graded material is a potential alternative to the one-directional functionally graded material. And the integrated design of materials and structures can really be achieved in two-directional functionally graded materials.展开更多
A new approach based on multiwavelets transformation and singular value decomposition (SVD) is proposed for the classification of image textures. Lower singular values are truncated based on its energy distribution to...A new approach based on multiwavelets transformation and singular value decomposition (SVD) is proposed for the classification of image textures. Lower singular values are truncated based on its energy distribution to classify the textures in the presence of additive white Gaussian noise (AWGN). The proposed approach extracts features such as energy, entropy, local homogeneity and max-min ratio from the selected singular values of multiwavelets transformation coefficients of image textures. The classification was carried out using probabilistic neural network (PNN). Performance of the proposed approach was compared with conventional wavelet domain gray level co-occurrence matrix (GLCM) based features, discrete multiwavelets transformation energy based approach, and HMM based approach. Experimental results showed the superiority of the proposed algorithms when compared with existing algorithms.展开更多
In this paper,the AUSMPW scheme based on adaptive algorithm of multi-wavelets is presented to solve two dimensional Euler equations.This scheme based on the original AUSMPW scheme uses the multiwavelets for multi-leve...In this paper,the AUSMPW scheme based on adaptive algorithm of multi-wavelets is presented to solve two dimensional Euler equations.This scheme based on the original AUSMPW scheme uses the multiwavelets for multi-level decomposition of the function and uses the method of the valve's value to construct adaptive grid to improve AUSMPW scheme.The obtained press and density have beed compared with those of results calculated by original AUSMPW scheme and WENO scheme.The numerical results demonstrate that this method has higher resolution.展开更多
A general method for constructing bidimensional orthogonal nonseparable mul- tiwavelets is presented. Moreover, this construction method can be extended to the con- struction of n-dimensional multiwavelets. In additio...A general method for constructing bidimensional orthogonal nonseparable mul- tiwavelets is presented. Moreover, this construction method can be extended to the con- struction of n-dimensional multiwavelets. In addition, we also study some properties of the multiwavelets such as balancing. Finally, we give an-example to illustrate our method to construct bidimensional nonseparable compactly supported orthogonal multiwavelets.展开更多
In this paper, we construct Chebyshev biorthogonal multiwavelets, and use this multiwavelets to approximate signals (functions). The convergence rate for signal approximation is derived. The fast signal decomposition ...In this paper, we construct Chebyshev biorthogonal multiwavelets, and use this multiwavelets to approximate signals (functions). The convergence rate for signal approximation is derived. The fast signal decomposition and reconstruction algorithms are presented. The numerical examples validate the theoretical analysis.展开更多
Based on Evans’spatial smoothing preprocessing scheme,a new approach calledtwo-direction spatial smoothing preprocessing method is presented.It is proved that the decorre-lation,the effective aperture and the maximum...Based on Evans’spatial smoothing preprocessing scheme,a new approach calledtwo-direction spatial smoothing preprocessing method is presented.It is proved that the decorre-lation,the effective aperture and the maximum number of distinguishable coherent signals(whenarray size is given)of the new method are better than those of the Evans’method.Simulationresults give a comparison between the eigenvector spectrums produced by the two methods.展开更多
Biorthogonal multiwavelets are generated from related scaling function vectors via multiresolution analysis. In this paper, we first show how to derive even-length biorthogonal lowpass filter pair from odd-length bior...Biorthogonal multiwavelets are generated from related scaling function vectors via multiresolution analysis. In this paper, we first show how to derive even-length biorthogonal lowpass filter pair from odd-length biorthogonal multiwavelet system with such properties as symmetry-antisymmetry and compactly support. So based on this, odd-length biorthogonal multiwavelet system can be constructed.展开更多
Deviation is essential to classic soft threshold denoising in wavelet domain. Texture features ofnoised image denoised by wavelet transform were weakened. Gibbs effect is distinct at edges of image.Image blurs compari...Deviation is essential to classic soft threshold denoising in wavelet domain. Texture features ofnoised image denoised by wavelet transform were weakened. Gibbs effect is distinct at edges of image.Image blurs comparing with original noised image. To solve the questions, a blind denoising method basedon single-wavelet transform and multiwavelets transform was proposed. The method doesn’t depend onsize of image and deviation to determine threshold of wavelet coefficients, which is different from classicalsoft-threshold denoising in wavelet domain. Moreover, the method is good for many types of noise. Gibbseffect disappeared with this method, edges of image are preserved well, and noise is smoothed andrestrained effectively.展开更多
In order to extract fault features of a weak signal from the strong noise and maintain signal smoothness, a new method of denoising based on the algorithm of balanced orthogonal multiwavelets is proposed. Multiwavelet...In order to extract fault features of a weak signal from the strong noise and maintain signal smoothness, a new method of denoising based on the algorithm of balanced orthogonal multiwavelets is proposed. Multiwavelets have several scaling functions and wavelet functions, and possess excellent properties that a scalar wavelet cannot satisfy simultaneously, and match the different characteristics of signals. Moreover, the balanced orthogonal multiwavelets can avoid the Gibbs phenomena and their processes have the advantages in denoising. Therefore, the denoising based on the algorithm of balanced orthogonal multiwavelets is introduced into the signal process. The algorithm of bal- anced orthogonal multiwavelet and the implementation steps of this denoising are described. The experimental compar- ison of the denoising effect between this algorithm and the traditional multiwavelet algorithm was done. The experi- ments indieate that this method is effective and feasible to extract the fault feature submerged in heavy noise.展开更多
A novel blind digital watermarking algorithm based on neural networks and multiwavelet transform is presented. The host image is decomposed through multiwavelet transform. There are four subblocks in the LL- level of ...A novel blind digital watermarking algorithm based on neural networks and multiwavelet transform is presented. The host image is decomposed through multiwavelet transform. There are four subblocks in the LL- level of the multiwavelet domain and these subblocks have many similarities. Watermark bits are added to low- frequency coefficients. Because of the learning and adaptive capabilities of neural networks, the trained neural networks almost exactly recover the watermark from the watermarked image. Experimental results demonstrate that the new algorithm is robust against a variety of attacks, especially, the watermark extraction does not require the original image.展开更多
Unlike scalar wavelets, multiscaling functions can be orthogonal, regular and symmetrical, and have compact support and high order of approximation simultaneously. For this reason, even if multiscaling functions are n...Unlike scalar wavelets, multiscaling functions can be orthogonal, regular and symmetrical, and have compact support and high order of approximation simultaneously. For this reason, even if multiscaling functions are not cardinal, they still hold for perfect A/D and D/A. We generalize the Walter's sampling theorem to multiwavelet subspaces based on reproducing kernel Hilbert space. The reconstruction function can be expressed by multiwavelet function using the Zak transform. The general case of irregular sampling is also discussed and the irregular sampling theorem for multiwavelet subspaces established. Examples are presented.展开更多
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 ...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 lowpass synthesis rnultifilter. Using these new results, we present the parametrization of orthogohal 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.展开更多
Synthetic aperture radar (SAR) images are corrupted by multiplicative speckle noise which limits the performance of the classical coder/decoder algorithm in spatial domain. The relatively new transform of multiwavel...Synthetic aperture radar (SAR) images are corrupted by multiplicative speckle noise which limits the performance of the classical coder/decoder algorithm in spatial domain. The relatively new transform of multiwavelets can possess desirable features simultaneously, such as orthogonality and symmetry, while scalar wavelets cannot. In this paper we propose a compression scheme combining with speckle noise reduction within the multiwavelet framework. Compared with classical set partitioning in hierarchical trees (SPIHT) algorithm, our method achieves favorable peak signal to noise ratio (PSNR) and superior speckle noise reduction performances.展开更多
The problem of estimating an image corrupted by additive white Gaussian noise has been of interest for practical reasons. Non-linear denoising methods based on wavelets, have become popular but Multiwavelets outperfor...The problem of estimating an image corrupted by additive white Gaussian noise has been of interest for practical reasons. Non-linear denoising methods based on wavelets, have become popular but Multiwavelets outperform wavelets in image denoising. Multiwavelets are wavelets with several scaling and wavelet functions, offer simultaneously Orthogonality, Symmetry, Short support and Vanishing moments, which is not possible with ordinary (scalar) wavelets. These properties make Multiwavelets promising for image processing applications, such as image denoising. The aim of this paper is to apply various non-linear thresholding techniques such as hard, soft, universal, modified universal, fixed and multivariate thresholding in Multiwavelet transform domain such as Discrete Multiwavelet Transform, Symmetric Asymmetric (SA4), Chui Lian (CL), and Bi-Hermite (Bih52S) for different Multiwavelets at different levels, to denoise an image and determine the best one out of it. The performance of denoising algorithms and various thresholding are measured using quantitative performance measures such as, Mean Square Error (MSE), and Root Mean Square Error (RMSE), Signal-to-Noise Ratio (SNR), Peak Signal-to-Noise Ratio (PSNR). It is found that CL Multiwavelet transform in combination with modified universal thresholding has given best results.展开更多
Based on the orthogonal multiwavelet model of 1/f signals, smoothing fractal signals from white Gaussian noise with multiwavelet filter is proposed. The proposed multiwavelet method is very simple and easy to realize....Based on the orthogonal multiwavelet model of 1/f signals, smoothing fractal signals from white Gaussian noise with multiwavelet filter is proposed. The proposed multiwavelet method is very simple and easy to realize. Compared with Wornell's single wavelet method, the new method has r filtering factors at each scale and has higher filtering speed, where r is the multiplicity of multiwavelet. Also due to the advantages of multiwavelet, the multiwavelet method performs better than that of Wornell's. Simulation results verify the analysis, and Wornell's method is the special case of our method when r = 1.展开更多
The element of pesedospectral-multiwavelet-Galerkin method, and how tocombine it with penalty method for treating boundary conditions are given. Multiwavelet bases don'toverlap on the given scale, and possess the ...The element of pesedospectral-multiwavelet-Galerkin method, and how tocombine it with penalty method for treating boundary conditions are given. Multiwavelet bases don'toverlap on the given scale, and possess the same compact set in a group of several functions, sothey can be directly used to the numerical discretion on the finite interval. Numerical tests showthat general boundary conditions can be enforced with the penalty method, and thatpesedospectral-multiwavelet-Galerkin method can well track the solutions' development. This alsoproves that pesedospectral-multiwavelet-Galerkin method is effective.展开更多
基金supported by the Natural Science Foundation China(11126343)Guangxi Natural Science Foundation(2013GXNSFBA019010)+1 种基金supported by Natural Science Foundation China(11071152)Natural Science Foundation of Guangdong Province(10151503101000025,S2011010004511)
文摘This article aims at studying two-direction refinable functions and two-direction wavelets in the setting R^s, s 〉 1. We give a sufficient condition for a two-direction refinable function belonging to L^2(R^s). Then, two theorems are given for constructing biorthogonal (orthogonal) two-direction refinable functions in L^2(R^s) and their biorthogonal (orthogonal) two-direction wavelets, respectively. From the constructed biorthogonal (orthogonal) two-direction wavelets, symmetric biorthogonal (orthogonal) multiwaveles in L^2(R^s) can be obtained easily. Applying the projection method to biorthogonal (orthogonal) two-direction wavelets in L^2(R^s), we can get dual (tight) two-direction wavelet frames in L^2(R^m), where m ≤ s. From the projected dual (tight) two-direction wavelet frames in L^2(R^m), symmetric dual (tight) frames in L^2(R^m) can be obtained easily. In the end, an example is given to illustrate theoretical results.
基金Project (2006AA11Z104) supported by the National High-Tech Research and Development Program("863" Program)
文摘Based on the discussion about working mechanism of horizontal reinforcement and that of vertical reinforcement,respectively,the working mechanism of two-direction reinforced composite foundation was studied.The enhancing effect of horizontal reinforcement on vertical reinforced composite foundation was analyzed.A simplified calculation method for such two-direction reinforced working system was presented.A model experiment was carried out to validate the proposed method.In the experiment,geocell reinforcement worked as the horizontal reinforcement,while gravel pile composite foundation worked as the vertical reinforcement.The results show that the calculated curve is close to the measured one.The installation of geosynthetic reinforcement can increase the bearing capacity of composite foundation by nearly 68% at normal foundation settlement,which suggests that the enhancing effect by geosynthetic reinforcement should be taken into account in current design/analysis methods.
文摘A new method for iris identification based on multiwavelets is proposed. By means of the properties of multiwavelets, such as orthogonality, symmetry, vanishing moments and approximation order, the iris texture can be simply presented. A brief overview of muhiwavelets is presented at first. Iris identification system and iris texture feature presentation and recognition based on multiwavelets a,e introduced subsequently. And the experiment indicates the validity of this method finally.
基金Project supported by the National Natural Science Foundation of China (No.10432030).
文摘Assuming the material properties varying with an exponential law both in the thick- ness and radial directions, axisymmetric bending of two-directional functionally graded circular and annular plates is studied using the semi-analytical numerical method in this paper. The deflections and stresses of the plates are presented. Numerical results show the well accuracy and convergence of the method. Compared with the finite element method, the semi-analytical nu- merical method is with great advantage in the computational efficiency. Moreover, study on ax- isymmetric bending of two-directional functionally graded annular plate shows that such plates have better performance than those made of isotropic homogeneous materials or one-directional functionally graded materials. Two-directional functionally graded material is a potential alternative to the one-directional functionally graded material. And the integrated design of materials and structures can really be achieved in two-directional functionally graded materials.
文摘A new approach based on multiwavelets transformation and singular value decomposition (SVD) is proposed for the classification of image textures. Lower singular values are truncated based on its energy distribution to classify the textures in the presence of additive white Gaussian noise (AWGN). The proposed approach extracts features such as energy, entropy, local homogeneity and max-min ratio from the selected singular values of multiwavelets transformation coefficients of image textures. The classification was carried out using probabilistic neural network (PNN). Performance of the proposed approach was compared with conventional wavelet domain gray level co-occurrence matrix (GLCM) based features, discrete multiwavelets transformation energy based approach, and HMM based approach. Experimental results showed the superiority of the proposed algorithms when compared with existing algorithms.
基金Sponsored by the National Natural Science Foundation of China(Grant No.50476028)
文摘In this paper,the AUSMPW scheme based on adaptive algorithm of multi-wavelets is presented to solve two dimensional Euler equations.This scheme based on the original AUSMPW scheme uses the multiwavelets for multi-level decomposition of the function and uses the method of the valve's value to construct adaptive grid to improve AUSMPW scheme.The obtained press and density have beed compared with those of results calculated by original AUSMPW scheme and WENO scheme.The numerical results demonstrate that this method has higher resolution.
文摘A general method for constructing bidimensional orthogonal nonseparable mul- tiwavelets is presented. Moreover, this construction method can be extended to the con- struction of n-dimensional multiwavelets. In addition, we also study some properties of the multiwavelets such as balancing. Finally, we give an-example to illustrate our method to construct bidimensional nonseparable compactly supported orthogonal multiwavelets.
文摘In this paper, we construct Chebyshev biorthogonal multiwavelets, and use this multiwavelets to approximate signals (functions). The convergence rate for signal approximation is derived. The fast signal decomposition and reconstruction algorithms are presented. The numerical examples validate the theoretical analysis.
文摘Based on Evans’spatial smoothing preprocessing scheme,a new approach calledtwo-direction spatial smoothing preprocessing method is presented.It is proved that the decorre-lation,the effective aperture and the maximum number of distinguishable coherent signals(whenarray size is given)of the new method are better than those of the Evans’method.Simulationresults give a comparison between the eigenvector spectrums produced by the two methods.
文摘Biorthogonal multiwavelets are generated from related scaling function vectors via multiresolution analysis. In this paper, we first show how to derive even-length biorthogonal lowpass filter pair from odd-length biorthogonal multiwavelet system with such properties as symmetry-antisymmetry and compactly support. So based on this, odd-length biorthogonal multiwavelet system can be constructed.
文摘Deviation is essential to classic soft threshold denoising in wavelet domain. Texture features ofnoised image denoised by wavelet transform were weakened. Gibbs effect is distinct at edges of image.Image blurs comparing with original noised image. To solve the questions, a blind denoising method basedon single-wavelet transform and multiwavelets transform was proposed. The method doesn’t depend onsize of image and deviation to determine threshold of wavelet coefficients, which is different from classicalsoft-threshold denoising in wavelet domain. Moreover, the method is good for many types of noise. Gibbseffect disappeared with this method, edges of image are preserved well, and noise is smoothed andrestrained effectively.
基金supported by Scientific and Technological Foundation of Henan Province under Grant No.112102210128Science Research Project of Educational Department of Henan Province under Grant No.2011C510005
文摘In order to extract fault features of a weak signal from the strong noise and maintain signal smoothness, a new method of denoising based on the algorithm of balanced orthogonal multiwavelets is proposed. Multiwavelets have several scaling functions and wavelet functions, and possess excellent properties that a scalar wavelet cannot satisfy simultaneously, and match the different characteristics of signals. Moreover, the balanced orthogonal multiwavelets can avoid the Gibbs phenomena and their processes have the advantages in denoising. Therefore, the denoising based on the algorithm of balanced orthogonal multiwavelets is introduced into the signal process. The algorithm of bal- anced orthogonal multiwavelet and the implementation steps of this denoising are described. The experimental compar- ison of the denoising effect between this algorithm and the traditional multiwavelet algorithm was done. The experi- ments indieate that this method is effective and feasible to extract the fault feature submerged in heavy noise.
基金The National Natural Science Foundation of China(No60473015)
文摘A novel blind digital watermarking algorithm based on neural networks and multiwavelet transform is presented. The host image is decomposed through multiwavelet transform. There are four subblocks in the LL- level of the multiwavelet domain and these subblocks have many similarities. Watermark bits are added to low- frequency coefficients. Because of the learning and adaptive capabilities of neural networks, the trained neural networks almost exactly recover the watermark from the watermarked image. Experimental results demonstrate that the new algorithm is robust against a variety of attacks, especially, the watermark extraction does not require the original image.
基金Project supported by the National Natural Science Foundation of China(Grant No.60672160)the Development Foundation of Shanghai Municipal Commission of Education(Grant No.05AZ42)
文摘Unlike scalar wavelets, multiscaling functions can be orthogonal, regular and symmetrical, and have compact support and high order of approximation simultaneously. For this reason, even if multiscaling functions are not cardinal, they still hold for perfect A/D and D/A. We generalize the Walter's sampling theorem to multiwavelet subspaces based on reproducing kernel Hilbert space. The reconstruction function can be expressed by multiwavelet function using the Zak transform. The general case of irregular sampling is also discussed and the irregular sampling theorem for multiwavelet subspaces established. Examples are presented.
文摘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 lowpass synthesis rnultifilter. Using these new results, we present the parametrization of orthogohal 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.
基金This work was supported by the National Natural Science Foundation of China under Grant No. 60472048.
文摘Synthetic aperture radar (SAR) images are corrupted by multiplicative speckle noise which limits the performance of the classical coder/decoder algorithm in spatial domain. The relatively new transform of multiwavelets can possess desirable features simultaneously, such as orthogonality and symmetry, while scalar wavelets cannot. In this paper we propose a compression scheme combining with speckle noise reduction within the multiwavelet framework. Compared with classical set partitioning in hierarchical trees (SPIHT) algorithm, our method achieves favorable peak signal to noise ratio (PSNR) and superior speckle noise reduction performances.
文摘The problem of estimating an image corrupted by additive white Gaussian noise has been of interest for practical reasons. Non-linear denoising methods based on wavelets, have become popular but Multiwavelets outperform wavelets in image denoising. Multiwavelets are wavelets with several scaling and wavelet functions, offer simultaneously Orthogonality, Symmetry, Short support and Vanishing moments, which is not possible with ordinary (scalar) wavelets. These properties make Multiwavelets promising for image processing applications, such as image denoising. The aim of this paper is to apply various non-linear thresholding techniques such as hard, soft, universal, modified universal, fixed and multivariate thresholding in Multiwavelet transform domain such as Discrete Multiwavelet Transform, Symmetric Asymmetric (SA4), Chui Lian (CL), and Bi-Hermite (Bih52S) for different Multiwavelets at different levels, to denoise an image and determine the best one out of it. The performance of denoising algorithms and various thresholding are measured using quantitative performance measures such as, Mean Square Error (MSE), and Root Mean Square Error (RMSE), Signal-to-Noise Ratio (SNR), Peak Signal-to-Noise Ratio (PSNR). It is found that CL Multiwavelet transform in combination with modified universal thresholding has given best results.
基金Supported by the National Laboratory of Space Microwave Technology Foundation(No.51473030105JB3201).
文摘Based on the orthogonal multiwavelet model of 1/f signals, smoothing fractal signals from white Gaussian noise with multiwavelet filter is proposed. The proposed multiwavelet method is very simple and easy to realize. Compared with Wornell's single wavelet method, the new method has r filtering factors at each scale and has higher filtering speed, where r is the multiplicity of multiwavelet. Also due to the advantages of multiwavelet, the multiwavelet method performs better than that of Wornell's. Simulation results verify the analysis, and Wornell's method is the special case of our method when r = 1.
基金This project is supported by National Natural Science Foundation of China(No. 19971020) Multidiseipline Scientific Research Foundation of Harbin Institute of Technology, China(No.HIT.MD2001.26).
文摘The element of pesedospectral-multiwavelet-Galerkin method, and how tocombine it with penalty method for treating boundary conditions are given. Multiwavelet bases don'toverlap on the given scale, and possess the same compact set in a group of several functions, sothey can be directly used to the numerical discretion on the finite interval. Numerical tests showthat general boundary conditions can be enforced with the penalty method, and thatpesedospectral-multiwavelet-Galerkin method can well track the solutions' development. This alsoproves that pesedospectral-multiwavelet-Galerkin method is effective.