In this study,a wavelet multi-resolution interpolation Galerkin method(WMIGM)is proposed to solve linear singularly perturbed boundary value problems.Unlike conventional wavelet schemes,the proposed algorithm can be r...In this study,a wavelet multi-resolution interpolation Galerkin method(WMIGM)is proposed to solve linear singularly perturbed boundary value problems.Unlike conventional wavelet schemes,the proposed algorithm can be readily extended to special node generation techniques,such as the Shishkin node.Such a wavelet method allows a high degree of local refinement of the nodal distribution to efficiently capture localized steep gradients.All the shape functions possess the Kronecker delta property,making the imposition of boundary conditions as easy as that in the finite element method.Four numerical examples are studied to demonstrate the validity and accuracy of the proposedwavelet method.The results showthat the use ofmodified Shishkin nodes can significantly reduce numerical oscillation near the boundary layer.Compared with many other methods,the proposed method possesses satisfactory accuracy and efficiency.The theoretical and numerical results demonstrate that the order of theε-uniform convergence of this wavelet method can reach 5.展开更多
This paper presents discrete wavelet transform (DWT) and its inverse (IDWT) with Haar wavelets as tools to compute the variable size interpolated versions of an image at optimum computational load. As a human observer...This paper presents discrete wavelet transform (DWT) and its inverse (IDWT) with Haar wavelets as tools to compute the variable size interpolated versions of an image at optimum computational load. As a human observer moves closer to or farther from a scene, the retinal image of the scene zooms in or out, respectively. This zooming in or out can be modeled using variable scale interpolation. The paper proposes a novel way of applying DWT and IDWT in a piecewise manner by non-uniform down- or up-sampling of the images to achieve partially sampled versions of the images. The partially sampled versions are then aggregated to achieve the final variable scale interpolated images. The non-uniform down-or up-sampling here is a function of the required scale of interpolation. Appropriate zero padding is used to make the images suitable for the required non-uniform sampling and the subsequent interpolation to the required scale. The concept of zeroeth level DWT is introduced here, which works as the basis for interpolating the images to achieve bigger size than the original one. The main emphasis here is on the computation of variable size images at less computational load, without compromise of quality of images. The interpolated images to different sizes and the reconstructed images are benchmarked using the statistical parameters and visual comparison. It has been found that the proposed approach performs better as compared to bilinear and bicubic interpolation techniques.展开更多
With the increasing availability of precipitation radar data from space,enhancement of the resolution of spaceborne precipitation observations is important,particularly for hazard prediction and climate modeling at lo...With the increasing availability of precipitation radar data from space,enhancement of the resolution of spaceborne precipitation observations is important,particularly for hazard prediction and climate modeling at local scales relevant to extreme precipitation intensities and gradients.In this paper,the statistical characteristics of radar precipitation reflectivity data are studied and modeled using a hidden Markov tree(HMT)in the wavelet domain.Then,a high-resolution interpolation algorithm is proposed for spaceborne radar reflectivity using the HMT model as prior information.Owing to the small and transient storm elements embedded in the larger and slowly varying elements,the radar precipitation data exhibit distinct multiscale statistical properties,including a non-Gaussian structure and scale-to-scale dependency.An HMT model can capture well the statistical properties of radar precipitation,where the wavelet coefficients in each sub-band are characterized as a Gaussian mixture model(GMM),and the wavelet coefficients from the coarse scale to fine scale are described using a multiscale Markov process.The state probabilities of the GMM are determined using the expectation maximization method,and other parameters,for instance,the variance decay parameters in the HMT model are learned and estimated from high-resolution ground radar reflectivity images.Using the prior model,the wavelet coefficients at finer scales are estimated using local Wiener filtering.The interpolation algorithm is validated using data from the precipitation radar onboard the Tropical Rainfall Measurement Mission satellite,and the reconstructed results are found to be able to enhance the spatial resolution while optimally reproducing the local extremes and gradients.展开更多
A floating-point wavelet-based and an integer wavelet-based image interpolations in lifting structures and polynomial curve fitting for image resolution enhancement are proposed in this paper. The proposed prediction ...A floating-point wavelet-based and an integer wavelet-based image interpolations in lifting structures and polynomial curve fitting for image resolution enhancement are proposed in this paper. The proposed prediction methods estimate high-frequency wavelet coefficients of the original image based on the available low-frequency wavelet coefficients, so that the original image can be reconstructed by using the proposed prediction method. To further improve the reconstruction performance, we use polynomial curve fitting to build relationships between actual high-frequency wavelet coefficients and estimated high-frequency wavelet coefficients. Results of the proposed prediction algorithm for different wavelet transforms are compared to show the proposed prediction algorithm outperforms other methods.展开更多
The wavelet multiresolution interpolation for continuous functions defined on a finite interval is developed in this study by using a simple alternative of transformation matrix.The wavelet multiresolution interpolati...The wavelet multiresolution interpolation for continuous functions defined on a finite interval is developed in this study by using a simple alternative of transformation matrix.The wavelet multiresolution interpolation Galerkin method that applies this interpolation to represent the unknown function and nonlinear terms independently is proposed to solve the boundary value problems with the mixed Dirichlet-Robin boundary conditions and various nonlinearities,including transcendental ones,in which the discretization process is as simple as that in solving linear problems,and only common two-term connection coefficients are needed.All matrices are independent of unknown node values and lead to high efficiency in the calculation of the residual and Jacobian matrices needed in Newton’s method,which does not require numerical integration in the resulting nonlinear discrete system.The validity of the proposed method is examined through several nonlinear problems with interior or boundary layers.The results demonstrate that the proposed wavelet method shows excellent accuracy and stability against nonuniform grids,and high resolution of localized steep gradients can be achieved by using local refined multiresolution grids.In addition,Newton’s method converges rapidly in solving the nonlinear discrete system created by the proposed wavelet method,including the initial guess far from real solutions.展开更多
Aiming at harmonic detection, fast Fourier transform can only detect integer harmonics precisely, short time Fourier transform can detect non-integer harmonics with low resolution, and some former wavelet based method...Aiming at harmonic detection, fast Fourier transform can only detect integer harmonics precisely, short time Fourier transform can detect non-integer harmonics with low resolution, and some former wavelet based methods have no aliasing-reduction scheme which result in low measurement precision and poor robustness. A frequency-domain interpolation algorithm to detect harmonics is proposed by choosing Shannon wavelet. Shannon wavelet is an orthogonal wavelet possessing best ideal frequency domain localization ability, it can restrict wavelet aliasing but bring about Gibbs oscillation phenomenon simultaneously. An interpolation algorithm is developed to overcome this problem. Simulation reveals that the proposed method can effectively cancel aliasing, spectral leakage and Gibbs phenomenon, so it provides an effective means for power system harmonic analysis.展开更多
In the research of path planning for manipulators with many DOF, generally there is a problem in most traditional methods, which is that their computational cost (time and memory space) increases exponentially as DOF ...In the research of path planning for manipulators with many DOF, generally there is a problem in most traditional methods, which is that their computational cost (time and memory space) increases exponentially as DOF or resolution of the discrete configuration space increases. So this paper presents the collision-free trajectory planning for the space robot to capture a target based on the wavelet interpolation algorithm. We made wavelet sample on the desired trajectory of the manipulator’s end-effector to do trajectory planning by use of the proposed wavelet interpolation formula, and then derived joint vectors from the trajectory information of the end-effector based on the fixed-attitude-restrained generalized Jacobian matrix of multi-arm coordinated motion, so as to control the manipulator to capture a static body along the desired collision-free trajectory. The method overcomes the shortcomings of the typical methods, and the desired trajectory of the end-effector can be any kind of complex nonlinear curve. The algorithm is simple and highly effective and the real trajectory is close to the desired trajectory. In simulation, the planar dual-arm three DOF space robot is used to demonstrate the proposed method, and it shows that the algorithm is feasible.展开更多
The authors introduce nonseparable scaling function interpolation and show that its approximation can provide similar convergence properties as scalar wavelet system. Several equivalent statements of accuracy of nonse...The authors introduce nonseparable scaling function interpolation and show that its approximation can provide similar convergence properties as scalar wavelet system. Several equivalent statements of accuracy of nonseparable scaling function are also given. In the numerical experiments, it appears that nonseparable scaling function interpolation has better convergence results than scalar wavelet systems in some cases.展开更多
In this paper,we study a special class of fractal interpolation functions,and give theirHaar-wavelet expansions.On the basis of the expansions,we investigate the H(o|¨)ldersmoothness of such functions and their l...In this paper,we study a special class of fractal interpolation functions,and give theirHaar-wavelet expansions.On the basis of the expansions,we investigate the H(o|¨)ldersmoothness of such functions and their logical derivatives of order α.展开更多
The amount of image data generated in multimedia applications is ever increasing. The image compression plays vital role in multimedia applications. The ultimate aim of image compression is to reduce storage space wit...The amount of image data generated in multimedia applications is ever increasing. The image compression plays vital role in multimedia applications. The ultimate aim of image compression is to reduce storage space without degrading image quality. Compression is required whenever the data handled is huge they may be required to sent or transmitted and also stored. The New Edge Directed Interpolation (NEDI)-based lifting Discrete Wavelet Transfrom (DWT) scheme with modified Set Partitioning In Hierarchical Trees (MSPIHT) algorithm is proposed in this paper. The NEDI algorithm gives good visual quality image particularly at edges. The main objective of this paper is to be preserving the edges while performing image compression which is a challenging task. The NEDI with lifting DWT has achieved 99.18% energy level in the low frequency ranges which has 1.07% higher than 5/3 Wavelet decomposition and 0.94% higher than traditional DWT. To implement this NEDI with Lifting DWT along with MSPIHT algorithm which gives higher Peak Signal to Noise Ratio (PSNR) value and minimum Mean Square Error (MSE) and hence better image quality. The experimental results proved that the proposed method gives better PSNR value (39.40 dB for rate 0.9 bpp without arithmetic coding) and minimum MSE value is 7.4.展开更多
Wavelet methods are a very useful tool in solving integral equations. Both scaling functions and wavelet functions are the key elements of wavelet methods. In this article, we use scaling function interpolation method...Wavelet methods are a very useful tool in solving integral equations. Both scaling functions and wavelet functions are the key elements of wavelet methods. In this article, we use scaling function interpolation method to solve Volterra integral equations of the first kind, and Fredholm-Volterra integral equations. Moreover, we prove convergence theorem for the numerical solution of Volterra integral equations and Freholm-Volterra integral equations. We also present three examples of solving Volterra integral equation and one example of solving Fredholm-Volterra integral equation. Comparisons of the results with other methods are included in the examples.展开更多
A novel anti-aliasing wavelet packet transform method for harmonic detection is proposed. Aiming at the low measurement precision and poor robustness which exists in the former traditional wavelet methods for lack of ...A novel anti-aliasing wavelet packet transform method for harmonic detection is proposed. Aiming at the low measurement precision and poor robustness which exists in the former traditional wavelet methods for lack of the aliasing reduction scheme, an optimal interpolation wavelet packet filter is designed according to new optimal criteria. First, the limitation of anti-aliasing on the traditional wavelet filter bank is analyzed. Second, the designed optimal interpolation filters are denoted, and then the solution algorithm is given. This devised wavelet packet filter can seek a reasonable balance between signal preservation and aliasing reduction; it overcomes the inherent bug of traditional wavelet transforms, which rooted from just only concerning total aliasing cancellation but not aliasing-reduction in decomposition. Simulation and several comparative results indicate that the proposed method can effectively eliminate aliasing and precisely extract harmonic information.展开更多
There are many problems in science and engineering where the solution shows a boundary layer character. Near the boundary the gradient is large in contrast with the smooth behaviour in the central core. A uniform grid...There are many problems in science and engineering where the solution shows a boundary layer character. Near the boundary the gradient is large in contrast with the smooth behaviour in the central core. A uniform grid is, therefore, not suitable for a numerical solution. MHD flow problems belong to this category where a velocity and induced magnetic field profiles get flattened in a transverse flow. In the present paper an optimized grid has been generated using interpo-lating wavelets. The results are compared with those obtained using uniform grid, the finite element method and also from the analytical solution.展开更多
Some theorems of compactly supported non-tensor product form two-dimension Daubechies wavelet were analysed carefully. Compactly supported non-tensor product form two-dimension wavelet was constructed, then non-tensor...Some theorems of compactly supported non-tensor product form two-dimension Daubechies wavelet were analysed carefully. Compactly supported non-tensor product form two-dimension wavelet was constructed, then non-tensor product form two dimension wavelet finite element was used to solve the deflection problem of elastic thin plate. The error order was researched. A numerical example was given at last.展开更多
In this article, we use scaling function interpolation method to solve linear Fredholm integral equations, and we prove a convergence theorem for the solution of Fredholm integral equations. We present two examples wh...In this article, we use scaling function interpolation method to solve linear Fredholm integral equations, and we prove a convergence theorem for the solution of Fredholm integral equations. We present two examples which have better results than others.展开更多
A supported framework of a gyroscope′s rotor is designed and the B-Spline wavelet finite element model of nonlinear supported magnetic field is worked out. A new finite element space is studied in which the scaling f...A supported framework of a gyroscope′s rotor is designed and the B-Spline wavelet finite element model of nonlinear supported magnetic field is worked out. A new finite element space is studied in which the scaling function of the B-spline wavelet is considered as the shape function of a tetrahedron. The magnetic field is spited by an artificial absorbing body which used the condition of field radiating, so the solution is unique. The resolution is improved via the varying gradient of the B-spline function under the condition of unchanging gridding. So there are some advantages in dealing with the focus flux and a high varying gradient result from a nonlinear magnetic field. The result is more practical. Plots of flux and in the space is studied via simulating the supported system model. The results of the study are useful in the research of the supported magnetic system for the gyroscope rotor.展开更多
基金supported by the National Natural Science Foundation of China (No.12172154)the 111 Project (No.B14044)+1 种基金the Natural Science Foundation of Gansu Province (No.23JRRA1035)the Natural Science Foundation of Anhui University of Finance and Economics (No.ACKYC20043).
文摘In this study,a wavelet multi-resolution interpolation Galerkin method(WMIGM)is proposed to solve linear singularly perturbed boundary value problems.Unlike conventional wavelet schemes,the proposed algorithm can be readily extended to special node generation techniques,such as the Shishkin node.Such a wavelet method allows a high degree of local refinement of the nodal distribution to efficiently capture localized steep gradients.All the shape functions possess the Kronecker delta property,making the imposition of boundary conditions as easy as that in the finite element method.Four numerical examples are studied to demonstrate the validity and accuracy of the proposedwavelet method.The results showthat the use ofmodified Shishkin nodes can significantly reduce numerical oscillation near the boundary layer.Compared with many other methods,the proposed method possesses satisfactory accuracy and efficiency.The theoretical and numerical results demonstrate that the order of theε-uniform convergence of this wavelet method can reach 5.
文摘This paper presents discrete wavelet transform (DWT) and its inverse (IDWT) with Haar wavelets as tools to compute the variable size interpolated versions of an image at optimum computational load. As a human observer moves closer to or farther from a scene, the retinal image of the scene zooms in or out, respectively. This zooming in or out can be modeled using variable scale interpolation. The paper proposes a novel way of applying DWT and IDWT in a piecewise manner by non-uniform down- or up-sampling of the images to achieve partially sampled versions of the images. The partially sampled versions are then aggregated to achieve the final variable scale interpolated images. The non-uniform down-or up-sampling here is a function of the required scale of interpolation. Appropriate zero padding is used to make the images suitable for the required non-uniform sampling and the subsequent interpolation to the required scale. The concept of zeroeth level DWT is introduced here, which works as the basis for interpolating the images to achieve bigger size than the original one. The main emphasis here is on the computation of variable size images at less computational load, without compromise of quality of images. The interpolated images to different sizes and the reconstructed images are benchmarked using the statistical parameters and visual comparison. It has been found that the proposed approach performs better as compared to bilinear and bicubic interpolation techniques.
基金This study was funded by the National Natural Science Foundation of China(Grant No.41975027)the Natural Science Foundation of Jiangsu Province(Grant No.BK20171457)the National Key R&D Program on Monitoring,Early Warning and Prevention of Major Natural Disasters(Grant No.2017YFC1501401).
文摘With the increasing availability of precipitation radar data from space,enhancement of the resolution of spaceborne precipitation observations is important,particularly for hazard prediction and climate modeling at local scales relevant to extreme precipitation intensities and gradients.In this paper,the statistical characteristics of radar precipitation reflectivity data are studied and modeled using a hidden Markov tree(HMT)in the wavelet domain.Then,a high-resolution interpolation algorithm is proposed for spaceborne radar reflectivity using the HMT model as prior information.Owing to the small and transient storm elements embedded in the larger and slowly varying elements,the radar precipitation data exhibit distinct multiscale statistical properties,including a non-Gaussian structure and scale-to-scale dependency.An HMT model can capture well the statistical properties of radar precipitation,where the wavelet coefficients in each sub-band are characterized as a Gaussian mixture model(GMM),and the wavelet coefficients from the coarse scale to fine scale are described using a multiscale Markov process.The state probabilities of the GMM are determined using the expectation maximization method,and other parameters,for instance,the variance decay parameters in the HMT model are learned and estimated from high-resolution ground radar reflectivity images.Using the prior model,the wavelet coefficients at finer scales are estimated using local Wiener filtering.The interpolation algorithm is validated using data from the precipitation radar onboard the Tropical Rainfall Measurement Mission satellite,and the reconstructed results are found to be able to enhance the spatial resolution while optimally reproducing the local extremes and gradients.
文摘A floating-point wavelet-based and an integer wavelet-based image interpolations in lifting structures and polynomial curve fitting for image resolution enhancement are proposed in this paper. The proposed prediction methods estimate high-frequency wavelet coefficients of the original image based on the available low-frequency wavelet coefficients, so that the original image can be reconstructed by using the proposed prediction method. To further improve the reconstruction performance, we use polynomial curve fitting to build relationships between actual high-frequency wavelet coefficients and estimated high-frequency wavelet coefficients. Results of the proposed prediction algorithm for different wavelet transforms are compared to show the proposed prediction algorithm outperforms other methods.
基金supported by the National Natural Science Foundation of China(Nos.12172154 and 11925204)the 111 Project of China(No.B14044)the National Key Project of China(No.GJXM92579)。
文摘The wavelet multiresolution interpolation for continuous functions defined on a finite interval is developed in this study by using a simple alternative of transformation matrix.The wavelet multiresolution interpolation Galerkin method that applies this interpolation to represent the unknown function and nonlinear terms independently is proposed to solve the boundary value problems with the mixed Dirichlet-Robin boundary conditions and various nonlinearities,including transcendental ones,in which the discretization process is as simple as that in solving linear problems,and only common two-term connection coefficients are needed.All matrices are independent of unknown node values and lead to high efficiency in the calculation of the residual and Jacobian matrices needed in Newton’s method,which does not require numerical integration in the resulting nonlinear discrete system.The validity of the proposed method is examined through several nonlinear problems with interior or boundary layers.The results demonstrate that the proposed wavelet method shows excellent accuracy and stability against nonuniform grids,and high resolution of localized steep gradients can be achieved by using local refined multiresolution grids.In addition,Newton’s method converges rapidly in solving the nonlinear discrete system created by the proposed wavelet method,including the initial guess far from real solutions.
文摘Aiming at harmonic detection, fast Fourier transform can only detect integer harmonics precisely, short time Fourier transform can detect non-integer harmonics with low resolution, and some former wavelet based methods have no aliasing-reduction scheme which result in low measurement precision and poor robustness. A frequency-domain interpolation algorithm to detect harmonics is proposed by choosing Shannon wavelet. Shannon wavelet is an orthogonal wavelet possessing best ideal frequency domain localization ability, it can restrict wavelet aliasing but bring about Gibbs oscillation phenomenon simultaneously. An interpolation algorithm is developed to overcome this problem. Simulation reveals that the proposed method can effectively cancel aliasing, spectral leakage and Gibbs phenomenon, so it provides an effective means for power system harmonic analysis.
文摘In the research of path planning for manipulators with many DOF, generally there is a problem in most traditional methods, which is that their computational cost (time and memory space) increases exponentially as DOF or resolution of the discrete configuration space increases. So this paper presents the collision-free trajectory planning for the space robot to capture a target based on the wavelet interpolation algorithm. We made wavelet sample on the desired trajectory of the manipulator’s end-effector to do trajectory planning by use of the proposed wavelet interpolation formula, and then derived joint vectors from the trajectory information of the end-effector based on the fixed-attitude-restrained generalized Jacobian matrix of multi-arm coordinated motion, so as to control the manipulator to capture a static body along the desired collision-free trajectory. The method overcomes the shortcomings of the typical methods, and the desired trajectory of the end-effector can be any kind of complex nonlinear curve. The algorithm is simple and highly effective and the real trajectory is close to the desired trajectory. In simulation, the planar dual-arm three DOF space robot is used to demonstrate the proposed method, and it shows that the algorithm is feasible.
文摘The authors introduce nonseparable scaling function interpolation and show that its approximation can provide similar convergence properties as scalar wavelet system. Several equivalent statements of accuracy of nonseparable scaling function are also given. In the numerical experiments, it appears that nonseparable scaling function interpolation has better convergence results than scalar wavelet systems in some cases.
文摘In this paper,we study a special class of fractal interpolation functions,and give theirHaar-wavelet expansions.On the basis of the expansions,we investigate the H(o|¨)ldersmoothness of such functions and their logical derivatives of order α.
文摘The amount of image data generated in multimedia applications is ever increasing. The image compression plays vital role in multimedia applications. The ultimate aim of image compression is to reduce storage space without degrading image quality. Compression is required whenever the data handled is huge they may be required to sent or transmitted and also stored. The New Edge Directed Interpolation (NEDI)-based lifting Discrete Wavelet Transfrom (DWT) scheme with modified Set Partitioning In Hierarchical Trees (MSPIHT) algorithm is proposed in this paper. The NEDI algorithm gives good visual quality image particularly at edges. The main objective of this paper is to be preserving the edges while performing image compression which is a challenging task. The NEDI with lifting DWT has achieved 99.18% energy level in the low frequency ranges which has 1.07% higher than 5/3 Wavelet decomposition and 0.94% higher than traditional DWT. To implement this NEDI with Lifting DWT along with MSPIHT algorithm which gives higher Peak Signal to Noise Ratio (PSNR) value and minimum Mean Square Error (MSE) and hence better image quality. The experimental results proved that the proposed method gives better PSNR value (39.40 dB for rate 0.9 bpp without arithmetic coding) and minimum MSE value is 7.4.
文摘Wavelet methods are a very useful tool in solving integral equations. Both scaling functions and wavelet functions are the key elements of wavelet methods. In this article, we use scaling function interpolation method to solve Volterra integral equations of the first kind, and Fredholm-Volterra integral equations. Moreover, we prove convergence theorem for the numerical solution of Volterra integral equations and Freholm-Volterra integral equations. We also present three examples of solving Volterra integral equation and one example of solving Fredholm-Volterra integral equation. Comparisons of the results with other methods are included in the examples.
文摘A novel anti-aliasing wavelet packet transform method for harmonic detection is proposed. Aiming at the low measurement precision and poor robustness which exists in the former traditional wavelet methods for lack of the aliasing reduction scheme, an optimal interpolation wavelet packet filter is designed according to new optimal criteria. First, the limitation of anti-aliasing on the traditional wavelet filter bank is analyzed. Second, the designed optimal interpolation filters are denoted, and then the solution algorithm is given. This devised wavelet packet filter can seek a reasonable balance between signal preservation and aliasing reduction; it overcomes the inherent bug of traditional wavelet transforms, which rooted from just only concerning total aliasing cancellation but not aliasing-reduction in decomposition. Simulation and several comparative results indicate that the proposed method can effectively eliminate aliasing and precisely extract harmonic information.
文摘There are many problems in science and engineering where the solution shows a boundary layer character. Near the boundary the gradient is large in contrast with the smooth behaviour in the central core. A uniform grid is, therefore, not suitable for a numerical solution. MHD flow problems belong to this category where a velocity and induced magnetic field profiles get flattened in a transverse flow. In the present paper an optimized grid has been generated using interpo-lating wavelets. The results are compared with those obtained using uniform grid, the finite element method and also from the analytical solution.
文摘Some theorems of compactly supported non-tensor product form two-dimension Daubechies wavelet were analysed carefully. Compactly supported non-tensor product form two-dimension wavelet was constructed, then non-tensor product form two dimension wavelet finite element was used to solve the deflection problem of elastic thin plate. The error order was researched. A numerical example was given at last.
文摘In this article, we use scaling function interpolation method to solve linear Fredholm integral equations, and we prove a convergence theorem for the solution of Fredholm integral equations. We present two examples which have better results than others.
文摘A supported framework of a gyroscope′s rotor is designed and the B-Spline wavelet finite element model of nonlinear supported magnetic field is worked out. A new finite element space is studied in which the scaling function of the B-spline wavelet is considered as the shape function of a tetrahedron. The magnetic field is spited by an artificial absorbing body which used the condition of field radiating, so the solution is unique. The resolution is improved via the varying gradient of the B-spline function under the condition of unchanging gridding. So there are some advantages in dealing with the focus flux and a high varying gradient result from a nonlinear magnetic field. The result is more practical. Plots of flux and in the space is studied via simulating the supported system model. The results of the study are useful in the research of the supported magnetic system for the gyroscope rotor.