• Journal of Institute of Control, Robotics and Systems
• /
• v.19 no.9
• /
• pp.827-833
• /
• 2013

In this paper, an efficient computational method is presented for state space analysis of bilinear systems via Legendre wavelets. The differential matrix equation is converted to a generalized Sylvester matrix equation by using Legendre wavelets as a basis. First, an explicit expression for the inverse of the integral operational matrix of the Legendre wavelets is presented. Then using it, we propose a preorder traversal algorithm to solve the generalized Sylvester matrix equation, which greatly reduces the computation time. Finally the efficiency of the proposed method is discussed using numerical examples.

• Journal of Institute of Control, Robotics and Systems
• /
• v.19 no.7
• /
• pp.577-582
• /
• 2013

In this paper Legendre wavelets are used to approximate the solutions of linear time-invariant system. The Legendre wavelet and its integral operational matrix are presented and an efficient algorithm to solve the Sylvester matrix equation is proposed. The algorithm is based on the decomposition of the Sylvester matrix equation and the preorder traversal algorithm. Using the special structure of the Legendre wavelet's integral operational matrix, the full order Sylvester matrix equation can be solved in terms of the solutions of pure algebraic matrix equations, which reduce the computation time remarkably. Finally a numerical example is illustrated to demonstrate the validity of the proposed algorithm.

• Journal of applied mathematics & informatics
• /
• v.28 no.5_6
• /
• pp.1055-1071
• /
• 2010

Recently, a number of algorithms have been proposed for numerical evaluation of Hankel transforms as these transforms arise naturally in many areas of science and technology. All these algorithms depend on separating the integrand $rf(r)J_{\upsilon}(pr)$ into two components; the slowly varying component rf(r) and the rapidly oscillating component $J_{\upsilon}(pr)$. Then the slowly varying component rf(r) is expanded either into a Fourier Bessel series or various wavelet series using different orthonormal bases like Haar wavelets, rationalized Haar wavelets, linear Legendre multiwavelets, Legendre wavelets and truncating the series at an optimal level; or approximating rf(r) by a quadratic over the subinterval using the Filon quadrature philosophy. The purpose of this communication is to take a different approach and replace rapidly oscillating component $J_{\upsilon}(pr)$ in the integrand by its Bernstein series approximation, thus avoiding the complexity of evaluating integrals involving Bessel functions. This leads to a very simple efficient and stable algorithm for numerical evaluation of Hankel transform.