Quintic B-spline collocation method for numerical solution of the Kuramoto–Sivashinsky equation

In this paper, the quintic B-spline collocation scheme is implemented to find numerical solution of the Kuramoto–Sivashinsky equation. The scheme is based on the Crank–Nicolson formulation for time integration and quintic B-spline functions for space integration. The accuracy of the proposed method is demonstrated by four test problems. The numerical results are found to be in good agreement with the exact solutions. Results are also shown graphically and are compared with results given in the literature.     

Numerical solution of the coupled viscous Burgers’ equation

In the present paper, a numerical method is proposed for the numerical solution of a coupled system of viscous Burgers’ equation with appropriate initial and boundary conditions, by using the cubic B-spline collocation scheme on the uniform mesh points. The scheme is based on Crank–Nicolson formulation for time integration and cubic B-spline functions for space integration. The method is shown to be unconditionally stable using von-Neumann method. The accuracy of the proposed method is demonstrated by applying it on three test problems. Computed results are depicted graphically and are compared with those already available in the literature. The obtained numerical solutions indicate that the method is reliable and yields results compatible with the exact solutions.     

On the quartic curve of Han

The quartic curve of Han [X. Han, Piecewise quartic polynomial curves with shape parameter, Journal of Computational and Applied Mathematics 195 (2006) 34–45] can be considered as the generalization of the cubic B-spline curve incorporating shape parameters into the polynomial basis functions. We show that this curve can be considered as the linear blending of the original cubic B-spline curve and a fixed quartic curve. Moreover, we present the Bézier form of the curve, which is useful in terms of incorporating the curve into existing CAD systems. Geometric effects of the alteration of shape parameters is also discussed, including design oriented computational methods for constrained shape control of the curve.     

On integro quartic spline interpolation

In this paper, we use quartic B-spline to construct an approximating function to agree with the given integral values of a univariate real-valued function over the same intervals. It is called integro quartic spline interpolation. Our interpolation method is new and easy to implement. Moreover, it can work successfully even without any boundary conditions. The interpolation errors are studied. The super convergence (sixth order and fourth order, respectively) in approximating function values and second-order derivative values at the knots is proved. Numerical examples illustrate that our method is very effective and our integro-interpolating quartic spline has higher approximation ability than others.     

B-spline based finite element method in one-dimensional discontinuous elastic wave propagation

The B-spline variant of the finite element method (FEM) is tested in one-dimensional discontinuous elastic wave propagation. The B-spline based FEM (called Isogeometric analysis IGA) uses spline functions as testing and shape functions in the Galerkin continuous content. Here, the accuracy of stress distribution and spurious oscillations of the B-spline based FEM are studied in numerical modeling of one-dimensional propagation of stress discontinuities in a bar, where the analytical solution is known. For time integration, the Newmark method, implicit form of the generalized-α method, the central difference method and the predictor/multi-corrector method are tested and compared. The use of lumped and consistent mass matrices in explicit time integration is discussed. Due to accuracy, the consistent mass matrix is preferred in explicit time integration in IGA.     

A closed-form solution for nonlinear oscillation and stability analyses of the elevator cable in a drum drive elevator system experiencing free vibration

Responses of the dynamical systems to some extent are affected by the natural frequencies. In the present paper, the parameter expansion method (PEM) is employed to investigate nonlinear oscillation and stability of the elevator’s drum as a single-degree-of-freedom (SDOF) swing system. A sensitivity analysis to observe the influence of various parameters on the nonlinear dynamic response, stability and natural frequency is performed. Comparing the results of the proposed closed-form analytical solution, the traditional numerical iterative time integration solution, and the linearized governing equations confirm the accuracy and efficiency of the proposed approach. Based on the results of the proposed closed-form solution, the linearization errors in calculating the natural frequencies in different cases are discussed as well. In contrast to the available numerical methods, the proposed method is free from the numerical damping and the time integration accumulated errors. Moreover, in comparison with the traditional multistep numerical iterative time integration methods, a much less computational time is required for the method in this research. Results reveal that for nonlinear systems, the natural frequency is remarkably affected by the initial conditions. Furthermore, the stability decreases as the dimensions of the mechanism increase.     

An efficient cubic B-spline and bicubic B-spline collocation method for numerical solutions of multidimensional nonlinear stochastic quadratic integral equations

This paper aims to develop a novel numerical approach on the basis of B-spline collocation method to approximate the solution of one-dimensional and two-dimensional nonlinear stochastic quadratic integral equations. The proposed approach is based on the hybrid of collocation method, cubic B-spline, and bi-cubic B-spline interpolation and Itô approximation. Using this method, the problem solving turns into a nonlinear system solution of equations that is solved by a suitable numerical method. Also, the convergence analysis of this numerical approach has been discussed. In the end, examples are given to test the accuracy and the implementation of the method. The results are compared with the results obtained by other methods to verify that this method is accurate and efficient.     

Algorithms for numerical solution of the modified equal width wave equation using collocation method

Quintic B-spline collocation algorithms for numerical solution of the modified equal width wave (MEW) equation have been proposed. The algorithms are based on Crank–Nicolson formulation for time integration and quintic B-spline functions for space integration. Quintic B-spline collocation method over the finite intervals is also applied to the time split MEW equation and space split MEW equation. Results for the three algorithms are compared by studying the propagation of the solitary wave, interaction of the solitary waves, wave generation and birth of solitons.     

线性定常系统非齐次两点边值问题的扩展精细积分方法

提出了一种求解非齐次线性两点边值问题的高精度和高稳定的扩展精细积分方法（EPIM）.首先引入了区段量（即区段矩阵和区段向量）来离散非齐次线性微分方程,建立了非齐次两点边值问题基于区段量的求解框架.在该框架下,不同区段的区段量可以并行计算,整体代数方程组的集成不依赖于边界条件.然后引入区段响应矩阵来处理两点边值问题的非齐次项,导出了多项式函数、指数函数、正/余弦函数及其组合函数形式的非齐次项对应的区段响应矩阵的加法定理,结合增量存储技术提出了EPIM.对具有上述函数形式的非齐次项,该方法可以得到计算机上的精确解,一般形式的非齐次项则利用上述函数近似求解.最后通过两个具有刚性特征的数值算例验证了该方法的高精度和高稳定性.     

A new efficient DQ algorithm for the solution of elliptic problems in higher dimensions

Two new efficient algorithms are developed to approximate the derivatives of sufficiently smooth functions. The new techniques are based on differential quadrature method with quartic B-spline bases as test functions. To obtain the weighting coefficients of differential quadrature method (DQM), we use the midpoints of a uniform partition mixed with near-boundary grid points. This enables us to obtain the weighting coefficients without adding the new extra relations. By obtaining the error bounds, it is proved that the method in its classic form is non-optimal. Then, some new weighting coefficients are constructed to obtain higher accuracy. By obtaining the error bounds, it is proved that the new algorithm is superconvergent. Afterwards, by defining some new symbols, we find a way to approximate the partial derivatives of multivariate functions. Also, some approximations are constructed to the mixed derivatives of multivariate functions. Finally, the applicability of the methods is examined by solving some well-known problems of partial differential equations. Some examples of 2D and 3D biharmonic, Poisson, and convection-diffusion equations are solved and compared to the existing methods to show the efficiency of the proposed algorithms.     

High order approximations using spline-based differential quadrature method: Implementation to the multi-dimensional PDEs

A new differential quadrature method based on cubic B-spline is developed for the numerical solution of differential equations. In order to develop the new approach, the B-spline basis functions are used on the grid and midpoints of a uniform partition. Some error bounds are obtained by help of cubic spline collocation, which show that the method in its classic form is second order convergent. In order to derive higher accuracy, high order perturbations of the problem are generated and applied to construct the numerical algorithm. A new fourth order method is developed for the numerical solution of systems of second order ordinary differential equations. By solving some test problems, the performance of the proposed methods is examined. Also the implementation of the method for multi-dimensional time dependent partial differential equations is presented. The stability of the proposed methods is examined via matrix analysis. To demonstrate the applicability of the algorithms, we solve the 2D and 3D coupled Burgers’ equations and 2D sine-Gordon equation as test problems. Also the coefficient matrix of the methods for multi-dimensional problems is described to analyze the stability.     

Extended Entropy Functional for Nonlinear Systems in Stochastic Dynamics

The maximum entropy approach is utilized for deriving the stationary probability density function of nonlinear stochastic systems to white noise excitation. To this aim a variational formulation is proposed where by means of the Lagrange multiplier methods the entropy functional is maximised constrained to the Fokker Planck equation. Some exact solutions in terms of Lagrange function of MDOF linear systems and for a class of SDOF nonlinear systems, are obtained.     

A numerical technique based on B-spline for a class of time-fractional diffusion equation

This paper presents an efficient numerical technique for solving a class of time-fractional diffusion equation. The time-fractional derivative is described in the Caputo form. The L1 scheme is used for discretization of Caputo fractional derivative and a collocation approach based on sextic B-spline basis function is employed for discretization of space variable. The unconditional stability of the fully-discrete scheme is analyzed. Two numerical examples are considered to demonstrate the accuracy and applicability of our scheme. The proposed scheme is shown to be sixth order accuracy with respect to space variable and (2 − α)-th order accuracy with respect to time variable, where α is the order of temporal fractional derivative. The numerical results obtained are compared with other existing numerical methods to justify the advantage of present method. The CPU time for the proposed scheme is provided.     

The Cubic B-Spline Method for a Class of Caputo-Fabrizio Fractional Differential Equations北大核心CSCD

基于分数阶微积分基本定理和三次B样条理论,构造了求解线性Caputo-Fabrizio型分数阶微分方程数值解的三次B样条方法,利用分数阶微积分基本定理将初值问题转化为关于解函数的表达式,再使用三次B样条函数逼近表达式中积分项的被积函数,进而计算了一类Caputo-Fabrizio型分数阶微分方程的数值解.给出了所构造的三次B样条方法的误差估计、收敛性和稳定性的理论证明.数值实验表明,该文数值方法在求解一类Caputo-Fabrizio型分数阶微分方程数值解时具有一定的可行性和有效性,且计算精度和计算效率优于现有的两种数值方法.     

二进小波的构造方法研究

提出两种二进小波的构造方法.首先,将Mallat构造的B-样条二进小波推广得到一种构造B-样条二进小波的新方法;其次,基于二进提升方案提出构造二进小波的另一种新方法—–构造定理,并通过调整定理中提升参数的形式、以新的B-样条二进小波作为初始二进小波,具体构造了具有有限长单位脉冲响应、高阶消失矩、线性相位的提升二进小波,这些提升二进小波不能由Sweldens提升方案得到.     

A Hilbert-Huang transform based identification method for general linear time-varying systems and weakly nonlinear systems

This paper proposes an identification method for general linear time-varying (LTV) MDOF systems and weakly nonlinear systems based on the Hilbert-Huang Transform (HHT)[1]. The proposed method uses Empirical mode decomposition (EMD) to decompose the response signals of systems into intrinsic mode functions (IMFs) and residues, and then analyzes the IMFs and the residues by Hilbert transform (HT) to obtain the analytical IMFs and analytical residues. After that, the above signals are synthesized to form new response signals and new analytical response signals. Finally, the new synthesized signals are used to identify the stiffness and damping coefficients of the systems. Three types of variation: smooth, abrupt and periodical variations are considered in the numerical simulations of LTV chainlike[2] and nonchainlike systems as well as weakly nonlinear systems such as Duffing oscillators and Van der Pol oscillators with white noise added in the system responses to demonstrate the effectiveness, accuracy and robustness of the proposed method. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Free vibration and buckling analysis of plates using B-spline wavelet on the interval Mindlin element

A finite element method (FEM) of B-spline wavelet on the interval (BSWI) is used in this paper to solve the free vibration and buckling problems of plates based on Reissner–Mindlin theory. By aid of the high accuracy of B-spline functions approximation for structural analysis, the proposed method could obtain a fast convergence and a satisfying numerical accuracy with fewer degrees of freedoms (DOF). The numerical examples demonstrate that the present BSWI method achieves the high accuracy compared to the exact solution and others existing approaches in the literatures. The BSWI finite element has potential to be used as a numerical method in analysis and design.     

Stabilization of the Lax-Wendroff method and a generalized one-step Runge-Kutta method for hyperbolic initial-value problems

In order ot integrate hyperbolic systems we distinguish explicit and implicit time integrators. Implicit methods allow large integration steps, but require more storage and are more difficult to implement than explicit methods. However explicit methods are subject to a restriction on the integration step. This restriction is a drawback if the variation of the solution in time is so small that accuracy considerations would allow a larger integration step. In this report we apply a smoothing technique in order to stabilize the Lax-Wendroff method and a generalized one-step Runge-Kutta method. Using this technique, the integration step is not limited by stability considerations.     

����ȱʧ���ݵ�B������ָ��ģ�͹���

??In this paper, we studied the inverse probability weighted least squares estimation of single-index model with response variable missing at random. Firstly, the B-spline technique is used to approximate the unknown single-index function, and then the objective function is established based on the inverse probability weighted least squares method. By the two-stage Newton iterative algorithm, the estimation of index parameters and the B-spline coefficients can be obtained. Finally, through many simulation examples and a real data application, it can be concluded that the method proposed in this paper performs very well for moderate sample     

Taylor–Galerkin and Taylor-collocation methods for the numerical solutions of Burgers’ equation using B-splines

In this paper, for the numerical solution of Burgers’ equation, we give two B-spline finite element algorithms which involve a collocation method with cubic B-splines and a Galerkin method with quadratic B-splines. In time discretization of the equation, Taylor series expansion is used. In order to verify the stabilities of the purposed methods, von-Neumann stability analysis is employed. To see the accuracy of the methods, L₂ and L_∞ error norms are calculated and obtained results are compared with some earlier studies.