Journal of Computational and Applied Mathematics
|
|
- Emma Hudson
- 5 years ago
- Views:
Transcription
1 Journal of Computational and Applied Mathematics 234 (21) Contents lists available at ScienceDirect Journal of Computational and Applied Mathematics journal homepage: A numerical solution of problems in calculus of variation using direct method and nonclassical parameterization Mohammad Maleki a,, Mahmoud Mashali-Firouzi b a Department of Mathematics, Payame oor University (PU), Isfahan, Iran b Department of Mathematics, School of Mohajer, Hezar Jerib Street, Isfahan, Iran a r t i c l e i n f o a b s t r a c t Article history: Received 28 May 28 Received in revised form 11 December 29 Keywords: onclassical parameterization Calculus of variation onlinear programming umerical methods In this paper a direct method for solving variational problems using nonclassical parameterization is presented. A nonclassical parameterization based on nonclassical orthogonal polynomials is first introduced to reduce a variational problem to a nonlinear mathematical programming problem. Then, using the Lagrange multiplier technique the problem is converted to that of solving a system of algebraic equations. Illustrative examples are given to demonstrate the validity and applicability of the technique. 21 Elsevier B.V. All rights reserved. 1. Introduction In a large number of problems arising in analysis, mechanics, geometry and so forth, it is necessary to determine the maximum or minimum of a certain functional. Because of the important role of this subject in science and engineering, considerable attention has been received by this kind of problems. Such problems are called variational problems. Some popular methods for solving variational problems are direct methods. In [1,2] the direct method of Ritz and Galerkin is investigated for solving variational problems. Chen and Hsiao [3] introduced the Walsh series method to variational problems. Due to the nature of the Walsh functions, the solution obtained was piecewise constant. Some orthogonal polynomials are applied on variational problems to find the continuous solutions for these problems [4 6]. Also the authors of [7 12] introduced the Legendre Wavelets method, rationalized Haar method, Adomian decomposition method, He s variational iteration method and Chebyshev finite difference method for solving variational problems, respectively. More historical comments about variational problems are found in [1,2]. A direct method converts the variational problem into a mathematical programming problem by using either the discretization technique [13,14] or the parameterization technique [15 17]. In classical parameterization, the variables of the problem are estimated using a finite length series of classical orthogonal polynomials such as Chebyshev, Legendre, etc. with unknown parameters [14 17]. The class of solution methods based on orthogonal polynomials has become known as spectral methods. Spectral methods are implemented in various ways. For example, the tau, Galerkin, and collocation methods have all been proposed as implementation strategies [18,19]. The collocation method has established itself as the one that permits the most convenient computer implementation. However, in pseudospectral methods the nodes must correspond to the zeros of the derivatives of classical orthogonal polynomials on the interval [ 1, 1], including the end points. The points are generally based on the Legendre or Chebyshev polynomials. The idea of employing nonclassical weight functions has been used in [2 26]. Corresponding author. address: mm_maleki25@yahoo.com (M. Maleki) /$ see front matter 21 Elsevier B.V. All rights reserved. doi:1.116/j.cam
2 M. Maleki, M. Mashali-Firouzi / Journal of Computational and Applied Mathematics 234 (21) In the present paper, we propose a direct method using nonclassical parameterization and nonclassical orthogonal polynomials, for finding the extremal of variational problems. The proposed method requires the definition of collocation points and nonclassical orthogonal polynomials. The application of the method to variational problem leads to a nonlinear mathematical programming problem and using Lagrange multiplier technique, to an algebraic system. In fact, our method is a combination of the pseudospectral and parameterization methods which improves the parameterization technique. The proposed method differs from the traditional parameterization and pseudospectral methods in several ways. First, by using the quasilinearization technique, the presented method could be employed for problems with non quadratic functional. Second, in the presented method the unknown parameters that would be determined are the values of variables x(t), in nodes of interpolation but in traditional parameterizations the unknown parameters are coefficients of a series of orthogonal polynomials. Third, in the proposed method a weighted interpolant is used to approximate the unknown functions as opposed to Lagrange interpolating polynomials. Fourth, the nodes are based on the roots of derivatives of the th order orthogonal polynomials based on arbitrary positive weight function across an arbitrary interval, whereas the traditional pseudospectral methods use nodes based on the roots of the derivatives of th order Jacobi polynomials across the interval [ 1, 1]. The proposed method has the following advantages: it can be employed for time-varying and nonlinear variational problems; the number of unknown parameters to be determined is less than the discretization technique [13,14]; this method converts the problems with non quadratic functional to a sequence of quadratic programming problems and gives the information of quadratic programming problem, the Hessian and the gradient, explicitly, and there is no need to use the differentiation or the finite difference method to construct these information; in this method by taking a fixed degree of weighted interpolation and changing the weight functions the obtained results can be improved. This paper is organized as follows: The following section is devoted to the explanation of the classical and nonclassical estimations of functions and generation of collocation points. In Section 3, we introduce the statement of problems in calculus of variations. In Section 4, we describe the basic formulation of the nonclassical parameterization method. umerical examples are then given in Section 5 to illustrate the applicability of the proposed method. 2. onclassical parameterization 2.1. Expansion in classical and nonclassical orthogonal polynomials The classical parameterization technique [15 17] approximates the function f L 2 [ 1, 1] by a finite length series of classical orthogonal polynomials such as Chebyshev or Legendre polynomials, with unknown parameters f (ˆt) = a j P j (ˆt), ˆt [ 1, 1], j= where P j (ˆt) is the jth order classical orthogonal polynomial and aj s are the unknown parameters. In this paper, we propose to expand the function using weighted interpolations of degree of the form [2,21] f (t) = F (t) = W(t) j= W(t j ) L j(t)f (t j ), t [a, b], (1) where t j, j =, 1,...,, are a set of distinct interpolation nodes in [a, b], W(t) is a positive weight function and L j (t), are a set of interpolating functions satisfies L j (t k ) = δ jk. Thus, F (t), defined by Eq. (1), is an interpolant to the function f (t) in the sense that f (t k ) = F (t k ), k =, 1,...,. The functions L j (t) are often chosen to be sets of interpolating polynomials of degree, in which case they can be represented explicitly by Lagrange s formula Generation of collocation points Here the main idea for the generation of collocation points (interpolation nodes), t j, is to form new sets of polynomials P n (t) orthogonal with respect to a weight function w(t) on the interval [a, b], that is b a w(t)p n (t)p m (t)dt = δ nm. The polynomials satisfy a three-term recurrence relation [22] P k+1 (t) = (t α k )P k (t) β k P k 1 (t), k =, 1, 2,... (3) P 1 (t) =, P (t) = 1. (2)
3 1366 M. Maleki, M. Mashali-Firouzi / Journal of Computational and Applied Mathematics 234 (21) It is important to note that in Eq. (1), W(t) is some weight function which adds additional flexibility because by taking a fixed and changing the weight functions, variate weighted interpolations can be employed as approximations of f (t). Also in Eq. (2), w(t) defines the orthogonal polynomials and the collocation points (interpolation nodes). The collocation points t j may be determined by the method outlined in [27]. The approach is based on determining the eigenvalues of a modified tridiagonal Jacobi matrix α β1 β1 α 1 β2 β2 α 2 β3 J = ,.. β 1 α 1 β β α where α, β are obtained from the solution of the linear system of equations ( ) ( ) ( ) P (a) P 1 (a) α ap (a) =. P (b) P 1 (b) bp (b) β The collocation points t j, including the end points, are determined as the eigenvalues of J, and can be determined if the coefficients α k and β k, k =, 1, 2,..., are known. The procedure for calculating the coefficients is discussed in the next subsection Generation of nonclassical orthogonal polynomials The generation of a set of polynomials orthogonal with respect to some weight function w(t) has been discussed in several texts [28 3]. For a nonclassical weight function, the usual Schmidt procedure [31] which involves orthogonalization of a given member of the set to all the functions of lower order, is highly unstable due to roundoff errors, and is not practical. The Schmidt procedure is analogous to the methods based on the moment of the weight function. The best approach is one based on the three term recurrence relation of the polynomials {P k (t)} given in Eq. (3). The recurrence coefficients in Eq. (3) are given in [23] by b a α k = tw(t)p 2 k (t)dt b β = b a a w(t)p 2 k (t)dt, k =, 1, 2,... w(t)p 2 (t)dt, β k = b a w(t)p 2 k (t)dt b. a w(t)p 2 k 1 (t)dt 3. Problem statement The simplest form of a variational problem can be considered as finding the extremum of the functional J[x(t)] = t1 t F[x(t), ẋ(t), t]dt. To find the extreme value of J, the boundary points of the admissible curves are known in the following form x(t ) = α, x(t 1 ) = β. (4) (5) The necessary condition for x(t) to extremize J[x(t)] is that it should satisfy the Euler Lagrange equation F x d ( ) F =, dt ẋ with boundary conditions given in Eq. (5). The above boundary value problem, does not always have a solution and if the solution exists, it may not be unique. ote that if the solution of Euler s equation satisfies the boundary conditions, it is unique. The general form of the variational problem in Eq. (4) is J[x 1, x 2,..., x n ] = t1 t F[x 1, x 2,..., x n, ẋ 1, ẋ 2,..., ẋ n, t]dt, (6)
4 M. Maleki, M. Mashali-Firouzi / Journal of Computational and Applied Mathematics 234 (21) Table 1 Different w(t) and W(t) on [, 1]. Case w(t) W(t) (1 + t) 1 3 e t2 1 4 e 2t 1 5 e 2t e t 6 e 2t (1 + t) t + t cos(t) cos(t) (1 + t) sin(t) 1 with the given boundary conditions for all functions x 1 (t ) = α 1, x 2 (t ) = α 2,..., x n (t ) = α n, (7) x 1 (t 1 ) = β 1, x 2 (t 1 ) = β 2,..., x n (t 1 ) = β n. (8) Here the necessary condition for the extremum of the functional in Eq. (6) is to satisfy the following system of second-order differential equations F d ( ) F =, i = 1, 2,..., n, (9) x i dt ẋ i with boundary conditions given in Eqs. (7) (8). However, the above system of differential equations can be solved easily only for simple cases. Here we consider a direct method for finding the extremum of Eq. (6) using a nonclassical parameterization technique. 4. Direct method using nonclassical parameterization Approximating variables x i (t), i = 1, 2,..., n, using Eq. (1), gives x i (t) = W(t) W(t j ) L j(t)a ij = a T i P(t), i = 1, 2,..., n, (1) j= where a i = [a i, a i1,..., a i ] T, i = 1, 2,..., n, are (+1) 1 vectors of unknown parameters with property that a ij = x i (t j ), and P(t) = [P (t), P 1 (t),..., P (t)] T is ( + 1) 1 vector with P j (t) = W(t) W(t j ) L j(t), j =, 1,...,. Then, ẋ i (t), can be represented as ẋ i (t) = a T i P (t), i = 1, 2,..., n, (11) where, P (t) is derivative vector of P(t). Substituting Eqs. (1) (11) in Eq. (6), the functional J becomes a nonlinear mathematical programming problem of unknown parameters a ij, where i = 1, 2,..., n, and j =, 1,...,. Hence, to find the extremum of J, we solve J a ij =, i = 1, 2,..., n, j =, 1,...,. (12) The above procedure is now used to solve the following variational problems. 5. Illustrative examples In this section, three examples are given to demonstrate the applicability and accuracy of our method. There are many combinations of orthogonal weight functions w(t) and weighted interpolating polynomials W(t) that could be selected [21, 24]. In this paper only a limited number of combinations are used for each example. The different cases are summarized in Table 1, the similar to these weights were also used in [24]. ote that we have computed the numerical results by the well-known symbolic software Mathematica 5.2.
5 1368 M. Maleki, M. Mashali-Firouzi / Journal of Computational and Applied Mathematics 234 (21) Example 1 Consider the problem of finding the minimum of the time-varying functional [6] J(x) = [ẋ 2 + tẋ + x 2 ]dt, (13) with boundary conditions x() =, x(1) = 1 4. (14) The exact solution of this problem is x(t) = c 1e t + c 2 e t where c 1 = Using Eqs. (1) (11) and (13), we get J(x) = Eq. (15) can be simplified to 2 e 4(e 2 1) and c 2 = e 2e2 4(e 2 1). [a T P (t)p T (t)a + a T tp (t) + a T P(t)P T (t)a]dt. (15) J(x) = 1 2 at Ha + c T a, (16) where H = 2 [P (t)p T (t) + P(t)P T (t)]dt, c T = tp T (t)dt. Eq. (1) and the boundary conditions in Eq. (14) imply x() = a T P() =, x(1) = a T P(1) = 1 4. (17) The quadratic programming problem in Eqs. (16) (17) can be rewritten as follows: subject to where J(x) = 1 2 at Ha + c T a, F 1 a b 1 =, F 1 = ( ) ( ) P T (), P T b (1) 1 = 1. 4 The optimal values of unknown parameters, a, can be obtained easily using Lagrange multiplier technique as a = H 1 c + H 1 F T 1 (F 1H 1 F T 1 ) 1 (F 1 H 1 c + b 1 ), by substituting these optimal parameters in Eq. (1) the values of x(t) can be calculated. ow we define the maximum errors for x (t) as E = x (t) x exact (t) = max{ x (t) x exact (t), t 1}, where x (t) is the approximation of x exact (t). In Table 2 we give the errors E for = 6 and = 1 and different cases of Table 1. In Table 3 a comparison is made between the approximate values of x(t) using the present method with = 6 and case 8, using the RH functions method [8] for k = 8 and the exact solution. In Table 4 a comparison is made between the approximate values of x(t) using the present method with = 6 and case 8, using the Legendre wavelets method [7] for M = 3 and k = 3 and the exact solution. Table 2 shows that in the present method by taking a fixed and changing the weight functions the obtained results improve.
6 M. Maleki, M. Mashali-Firouzi / Journal of Computational and Applied Mathematics 234 (21) Table 2 The maximum errors of E for = 6, 1 for example E (case 1) E (case 2) E (case 3) E (case 4) E (case 5) E (case 6) E (case 7) E (case 8) E (case 9) E (case 1) Table 3 Estimated and exact values of x(t) for example 1. t RH functions for k = 8 [8] t Present method with = 6 and case 8 Exact t < 1/ /8 t < 2/ /8 t < 3/ /8 t < 4/ /8 t < 5/ /8 t < 5/ /8 t < 6/ /8 t < 7/ /8 t < t = Table 4 Estimated and exact values of x(t) for example 1. t Legendre wavelets method for M = 3 and k = 3 [7] Present method with = 6 and case 8 Exact Example 2 In this example, we consider the following variational problem [2,1]. J = 1 + x 2 (t) dt, ẋ 2 (t) that satisfies the conditions x() =, x(1) =.5. The exact solution of this problem is x(t) = sinh( t). In this problem the functional has not quadratic form, hence, for solving this problem using the presented method, we first (using the quasilinearization method [32]) expand the performance index and boundary conditions around x k (t) and ẋ k (t) up to the second and first order respectively, to get the following sequence of linear-quadratic problems. [ 1 J k+1 = 2 C k (t)x k ] 2 Ek (t)(ẋ k+1 ) 2 + α k (t)x k+1 + β k (t)ẋ k+1 + D k (t)x k+1 ẋ k+1 + h k (t) dt, (2) subject to (18) (19) x k+1 () =, x k+1 (1) =.5, (21)
7 137 M. Maleki, M. Mashali-Firouzi / Journal of Computational and Applied Mathematics 234 (21) where h k (t) = f (x k, ẋ k, t) A k (t)x k B k (t)ẋ k C k (t)x k Ek (t)(ẋ k ) 2 + D k (t)x k ẋ k, A k f (x, ẋ, t) (t) =, x x k,ẋ k B k f (x, ẋ, t) (t) =, ẋ x k,ẋ k C k (t) = 2 f (x, ẋ, t), x 2 x k,ẋ k D k (t) = 2 f (x, ẋ, t), x ẋ x k,ẋ k E k (t) = 2 f (x, ẋ, t), ẋ 2 x k,ẋ k α k (t) = A k (t) x k C k (t) ẋ k D k (t), β k (t) = B k (t) x k D k (t) ẋ k E k (t), and in this problem we have f (x, ẋ, t) = (1 + x 2 (t))/ẋ 2 (t). The initial guesses of x (t) and ẋ (t) is requested. We suggest to start from the linear function between the initial value x() = and the final value x(1) =.5 for x (t) and its derivation for ẋ (t). Substituting Eqs. (1) (11) in sequence of problems in Eqs. (2) (21) gives a sequence of problems as follows: [ 1 J k+1 = 2 C k (t) ( ) a T PP T a Ek (t) ( a T P ) P T a + α k (t) ( ) a T P + β k (t) ( a T P ) + D k (t) ( ) ] a T PP T a + h k (t) dt, (22) subject to a T P() =, a T P(1) =.5, where J k+1 is the approximate value of J k+1. The sequence of problems in Eqs. (22) (23) can be rewritten as a sequence of quadratic programming problems as follows: (23) subject to where J k+1 = 1 2 at Ha + c T a + d, F 1 a b 1 =, H = c T = d = h k (t)dt, ( ) P F 1 = T (), P T (1) ( ) b 1 =..5 [ C k (t)pp T + E k (t)p P T + 2D k (t)pp T ] dt, [ α k (t)p T + β k (t)p T ] dt, (24) (25) These quadratic programming problems are solved similarly to Example 1. After obtaining the optimal solution of the unknown parameters a, we substitute these parameters in Eqs. (1) (11) to obtain the new nominal rate x k (t) and ẋ k (t) to be used in the next iteration. These new nominal trajectories have to be substituted in Eqs. (2) (21) to get the next constrained quadratic problem. This procedure has to be repeated until an acceptable convergence is achieved. In this paper the computations are terminated if (J k+1 J k k+1 )/J < 1 1. The convergence is achieved after three iterations for each.
8 M. Maleki, M. Mashali-Firouzi / Journal of Computational and Applied Mathematics 234 (21) Table 5 The maximum errors of E for example E (case 1) E (case 2) E (case 3) E (case 4) E (case 5) E (case 6) E (case 7) E (case 8) E (case 9) E (case 1) Method of Saadatmandi and Dehghan [12] ote that we can employ the method explained in this example for all nonlinear problems that their functional in not quadratic. In Table 5 a comparison is made between the errors E using the present method with = 3, 4, 5 and 6 and the method of Saadatmandi and Dehghan [12] by using Chebyshev finite difference method. ot that their method is not direct and it solves the two point boundary value problem arrives in Eq. (9) Example 3 Consider the extremization of functional [2,1] J(x 1 (t), x 2 (t)) = π 2 with the boundary conditions ( π ) x 1 () =, x 1 = 1, 2 ( π ) x 2 () =, x 2 = 1. 2 [ẋ ẋ x 1x 2 ]dt, (26) The exact solutions of the problem are x 1 (t) = sin(t) and x 2 (t) = sin(t). For this problem the transformation t = π τ is 2 used to change the variational problem in Eqs. (26) (28) to the following form J(x 1 (τ), x 2 (τ)) = π [ 4 2 π 2 ẋ2 + 4 ] 1 π 2 ẋ x 1x 2 dτ, (29) with the boundary conditions x 1 () =, x 1 (1) = 1, x 2 () =, x 2 (1) = 1. We approximate x 1 (τ) and x 2 (τ) and their derivations using Eqs. (1) (11) to get the following nonlinear mathematical programming as approximation to the original problem. subject to where J = a T 1 H 1a 1 + a T 2 H 1a 2 + a T 1 H 2a 2, F 1 a 1 b 1 =, F 1 a 2 b 2 =, H 1 = 2 π P (τ)p T (τ)dτ, H 2 = 2 P(τ)P T (τ)dτ, π ( ) P F 1 = T (), P T (1) (27) (28) (3) (31) (32) (33) (34)
9 1372 M. Maleki, M. Mashali-Firouzi / Journal of Computational and Applied Mathematics 234 (21) Table 6 The maximum errors of E and E for example E (case 1) E E (case 2) E E (case 3) E E (case 4) E E (case 5) E E (case 6) E E (case 7) E E (case 8) E E (case 9) E E (case 1) E Method of Saadatmandi and Dehghan [12] E E ( ) b 1 =, 1 ( ) b 2 =. 1 The nonlinear mathematical programming problem in Eqs. (32) (34) can be solved easily by Lagrange multiplier technique. In Table 6 a comparison is made between the maximum errors E and E for x 1(t) and x 2 (t), respectively, for = 6, 1, 12 and the method of Saadatmandi and Dehghan [12] by using Chebyshev finite difference method. 6. Conclusion This paper describes an efficient method for finding the minimum of a functional over the specified domain. The main objective is to estimate the rate variables using nonclassical parameterization based on nonclassical orthogonal polynomials, to arrive at a nonlinear mathematical programming problem. Applications of the presented method are demonstrated through illustrative examples. The obtained results show that this approach can solve the problem effectively. Acknowledgements The authors are very grateful to the anonymous referees and also Prof. T. Mitsui (editor) for carefully reading the paper and for their comments and suggestions. References [1] I.M. Gelfand, S.V. Fomin, Calculus of Variations, Prentice-Hall, J, 1963 (revised English edition translated and edited by R.A. Silverman). [2] L. Elsgolts, Differential Equations and Calculus of Variations, Mir, Moscow, 1977 (translated from the Russian by G. Yankovsky). [3] C.F. Chen, C.H. Hsiao, A walsh series direct method for solving variational problems, J. Franklin Inst. 3 (1975) [4] R.Y. Chang, M.L. Wang, Shifted Legendre direct method for variational problems, J. Optim. Theory Appl. 39 (1983) [5] I.R. Horng, J.H. Chou, Shifted Chebyshev direct method for solving variational problems, Internat. J. Systems Sci. 16 (1985) [6] C. Hwang, Y.P. Shih, Laguerre series direct method for variational problems, J. Optim. Theory Appl. 39 (1) (1983) [7] M. Razzaghi, S. Yousefi, Legendre wavelets direct method for variational problems, Math. Comput. Simulation 53 (2) [8] M. Razzaghi, Y. Ordokhani, An application of rationalized Haar functions for variational problems, Appl. Math. Comput. 122 (21) [9] M. Dehghan, M. Tatari, The use of Adomian decomposition method for solving problems in calculus of variations, Math. Probl. Eng. 26 (26) [1] M. Tatari, M. Dehghan, Solution of problems in calculus of variations via He s variational iteration method, Phys. Lett. A 362 (27) [11] S.A. Yousefi, M. Dehghan, The use of He s variational iteration method for solving variational problems, Int. J. Comput. Math. (29), doi:1.18/
10 M. Maleki, M. Mashali-Firouzi / Journal of Computational and Applied Mathematics 234 (21) [12] A. Saadatmandi, M. Dehghan, The numerical solution of problems in calculus of variation using Chebyshev finite difference method, Phys. Lett. A 372 (28) [13] O. Stryk, R. Bulirsch, Direct and indirect methods for trajectory optimization, Ann. Oper. Res. 37 (1992) [14] D. Kraft, On converting optimal control problems into nonlinear programming problems, in: K. Schittkowski (Ed.), Computational Mathematical Programming, vol. F15, Springer, Berlin, 1985, pp [15] C.J. Goh, K.L. Teo, Control parameterization: A unified approach to optimal control problems with general constraints, Automatica 24 (1988) [16] J. Vlassenbroeck, R. Van Dooren, A Chebyshev technique for solving nonlinear optimal control problems, IEEE Trans. Automat. Control 33 (1988) [17] H. Jaduu, Direct solution of nonlinear optimal control problems using quasilinearization and Chebyshev polynomials, J. Franklin Inst. 339 (22) [18] D. Gottlieb, A. Orszag, umerical Analysis of Spectral Methods: Theory and Applications, SIAM, Philadelphia, [19] B. Fornberg, A Practical Guide to Pseudospectral Methods, Cambridge University Press, [2] J.A.C. Weidman, Spectral methods based on nonclassical orthogonal polynomials, approximations and computation of orthogonal polynomials, in: International Series of umerical Mathematics, vol. 131, Birkhauser, Basel, 1999, pp [21] H. Chen, B.D. Shizgal, A spectral solution of the Sturm Liouville equation: Comparison of classical and nonclassical basis sets, J. Comput. Appl. Math. 136 (21) [22] B. Shizgal, A Gaussian quadrature procedure for use in the solution of the Boltzman equation and related problems, J. Comput. Phys. 41 (1981) [23] B. Shizgal, H. Chen, The quadrature discretization method (QDM) in the solution of the Schrödinger equation with nonclassical basis functions, J. Chem. Phys. 14 (1996) [24] P. Williams, A quadrature discretization for solving optimal control problems, in: Proceedings of the AAS/AIAA Space Flight Mechanic Meeting, February 8 12, 24, Hawaii, in: Advances in the Astronautical Sciences, Spaceflight Mechanics, vol. 119, 24, pp [25] A. Alipanah, M. Razzaghi, M. Dehghan, onclassical pseudospectral method for the solution of brachistochrone problem, Chaos Solitons Fractals 34 (27) [26] M. Maleki, A. Hadi-Vencheh, Combination of non-classical pseudospectral and direct methods for the solution of brachistichrone problem, Int. J. Comput. Math. (29), doi:1.18/ [27] G.H. Golub, J.H. Welsch, Calculation of Gauss quadrature rules, Math. Comp. 23 (1969) [28] T.S. Chihara, An Introduction to Orthogonl Polynomials, Gordon and Breach, ew York, [29] P. evai, Orthogonal Polynomials: Theory and Practice, Kluwer Academic, Cambridge University Press, orwell, MA, 199. [3] W. Gautschi, Orthogonal polynomials-constractive theory and applications, J. Comput. Appl. Math (1985) [31] J.D. Jakson, Mathematics for Quantum Mechanics, Academic Press, ew York, 1975, pp [32] R. Bellman, R. Kalaba, Quasilinearization and onlinear Boundary Value Problems, Elsevier, ew York, 1965.
Research Article The Numerical Solution of Problems in Calculus of Variation Using B-Spline Collocation Method
Applied Mathematics Volume 2012, Article ID 605741, 10 pages doi:10.1155/2012/605741 Research Article The Numerical Solution of Problems in Calculus of Variation Using B-Spline Collocation Method M. Zarebnia
More informationDirect method for variational problems by using hybrid of block-pulse and Bernoulli polynomials
Direct method for variational problems by using hybrid of block-pulse and Bernoulli polynomials M Razzaghi, Y Ordokhani, N Haddadi Abstract In this paper, a numerical method for solving variational problems
More informationTHE USE OF ADOMIAN DECOMPOSITION METHOD FOR SOLVING PROBLEMS IN CALCULUS OF VARIATIONS
THE USE OF ADOMIAN DECOMPOSITION METHOD FOR SOLVING PROBLEMS IN CALCULUS OF VARIATIONS MEHDI DEHGHAN AND MEHDI TATARI Received 16 March 25; Revised 9 August 25; Accepted 12 September 25 Dedication to Professor
More informationThe Series Solution of Problems in the Calculus of Variations via the Homotopy Analysis Method
The Series Solution of Problems in the Calculus of Variations via the Homotopy Analysis Method Saeid Abbasbandy and Ahmand Shirzadi Department of Mathematics, Imam Khomeini International University, Ghazvin,
More informationResearch Article A Nonclassical Radau Collocation Method for Nonlinear Initial-Value Problems with Applications to Lane-Emden Type Equations
Applied Mathematics Volume 2012, Article ID 103205, 13 pages doi:10.1155/2012/103205 Research Article A Nonclassical Radau Collocation Method for Nonlinear Initial-Value Problems with Applications to Lane-Emden
More informationNumerical solution of delay differential equations via operational matrices of hybrid of block-pulse functions and Bernstein polynomials
Computational Methods for Differential Equations http://cmdetabrizuacir Vol, No, 3, pp 78-95 Numerical solution of delay differential equations via operational matrices of hybrid of bloc-pulse functions
More informationA Gauss Lobatto quadrature method for solving optimal control problems
ANZIAM J. 47 (EMAC2005) pp.c101 C115, 2006 C101 A Gauss Lobatto quadrature method for solving optimal control problems P. Williams (Received 29 August 2005; revised 13 July 2006) Abstract This paper proposes
More informationState Analysis and Optimal Control of Linear. Time-Invariant Scaled Systems Using. the Chebyshev Wavelets
Contemporary Engineering Sciences, Vol. 5, 2012, no. 2, 91-105 State Analysis and Optimal Control of Linear Time-Invariant Scaled Systems Using the Chebyshev Wavelets M. R. Fatehi a, M. Samavat b *, M.
More informationNumerical solution of optimal control problems by using a new second kind Chebyshev wavelet
Computational Methods for Differential Equations http://cmde.tabrizu.ac.ir Vol. 4, No. 2, 2016, pp. 162-169 Numerical solution of optimal control problems by using a new second kind Chebyshev wavelet Mehdi
More informationA Legendre Computational Matrix Method for Solving High-Order Fractional Differential Equations
Mathematics A Legendre Computational Matrix Method for Solving High-Order Fractional Differential Equations Mohamed Meabed KHADER * and Ahmed Saied HENDY Department of Mathematics, Faculty of Science,
More informationKybernetika. Terms of use: Persistent URL: Institute of Information Theory and Automation AS CR, 2015
Kybernetika Yousef Edrisi Tabriz; Mehrdad Lakestani Direct solution of nonlinear constrained quadratic optimal control problems using B-spline functions Kybernetika, Vol. 5 (5), No., 8 98 Persistent URL:
More informationGeneralized B-spline functions method for solving optimal control problems
Computational Methods for Differential Equations http://cmde.tabrizu.ac.ir Vol. 2, No. 4, 24, pp. 243-255 Generalized B-spline functions method for solving optimal control problems Yousef Edrisi Tabriz
More informationA Numerical Solution of Volterra s Population Growth Model Based on Hybrid Function
IT. J. BIOAUTOMATIO, 27, 2(), 9-2 A umerical Solution of Volterra s Population Growth Model Based on Hybrid Function Saeid Jahangiri, Khosrow Maleknejad *, Majid Tavassoli Kajani 2 Department of Mathematics
More informationSOLUTION OF NONLINEAR VOLTERRA-HAMMERSTEIN INTEGRAL EQUATIONS VIA SINGLE-TERM WALSH SERIES METHOD
SOLUTION OF NONLINEAR VOLTERRA-HAMMERSTEIN INTEGRAL EQUATIONS VIA SINGLE-TERM WALSH SERIES METHOD B. SEPEHRIAN AND M. RAZZAGHI Received 5 November 24 Single-term Walsh series are developed to approximate
More informationSOLUTION FOR A CLASSICAL PROBLEM IN THE CALCULUS OF VARIATIONS VIA RATIONALIZED HAAR FUNCTIONS
KYBERNETIKA VOLUME 3 7 (2001), NUMBER 5, PAGES 575-583 SOLUTION FOR A CLASSICAL PROBLEM IN THE CALCULUS OF VARIATIONS VIA RATIONALIZED HAAR FUNCTIONS MOHSEN RAZZAGHI 1 AND YADOLLAH ORDOKHANI A numerical
More informationA Numerical Solution of the Lax s 7 th -order KdV Equation by Pseudospectral Method and Darvishi s Preconditioning
Int. J. Contemp. Math. Sciences, Vol. 2, 2007, no. 22, 1097-1106 A Numerical Solution of the Lax s 7 th -order KdV Equation by Pseudospectral Method and Darvishi s Preconditioning M. T. Darvishi a,, S.
More informationChebyshev finite difference method for a Two Point Boundary Value Problems with Applications to Chemical Reactor Theory
Iranian Journal of Mathematical Chemistry, Vol. 3, o., February 22, pp. 7 IJMC Chebyshev finite difference method for a Two Point Boundary Value Problems with Applications to Chemical Reactor Theory ABBAS
More informationRATIONAL CHEBYSHEV COLLOCATION METHOD FOR SOLVING NONLINEAR ORDINARY DIFFERENTIAL EQUATIONS OF LANE-EMDEN TYPE
INTERNATIONAL JOURNAL OF INFORMATON AND SYSTEMS SCIENCES Volume 6, Number 1, Pages 72 83 c 2010 Institute for Scientific Computing and Information RATIONAL CHEBYSHEV COLLOCATION METHOD FOR SOLVING NONLINEAR
More informationDirect method to solve Volterra integral equation of the first kind using operational matrix with block-pulse functions
Journal of Computational and Applied Mathematics 22 (28) 51 57 wwwelseviercom/locate/cam Direct method to solve Volterra integral equation of the first kind using operational matrix with block-pulse functions
More informationA Third-degree B-spline Collocation Scheme for Solving a Class of the Nonlinear Lane-Emden Type Equations
Iranian Journal of Mathematical Sciences and Informatics Vol. 12, No. 2 (2017), pp 15-34 DOI: 10.7508/ijmsi.2017.2.002 A Third-degree B-spline Collocation Scheme for Solving a Class of the Nonlinear Lane-Emden
More informationCOMPUTATION OF BESSEL AND AIRY FUNCTIONS AND OF RELATED GAUSSIAN QUADRATURE FORMULAE
BIT 6-85//41-11 $16., Vol. 4, No. 1, pp. 11 118 c Swets & Zeitlinger COMPUTATION OF BESSEL AND AIRY FUNCTIONS AND OF RELATED GAUSSIAN QUADRATURE FORMULAE WALTER GAUTSCHI Department of Computer Sciences,
More informationConvergence of a Gauss Pseudospectral Method for Optimal Control
Convergence of a Gauss Pseudospectral Method for Optimal Control Hongyan Hou William W. Hager Anil V. Rao A convergence theory is presented for approximations of continuous-time optimal control problems
More informationSOLUTION OF DIFFERENTIAL EQUATIONS BASED ON HAAR OPERATIONAL MATRIX
Palestine Journal of Mathematics Vol. 3(2) (214), 281 288 Palestine Polytechnic University-PPU 214 SOLUTION OF DIFFERENTIAL EQUATIONS BASED ON HAAR OPERATIONAL MATRIX Naresh Berwal, Dinesh Panchal and
More informationExact Solutions of Fractional-Order Biological Population Model
Commun. Theor. Phys. (Beijing China) 5 (009) pp. 99 996 c Chinese Physical Society and IOP Publishing Ltd Vol. 5 No. 6 December 15 009 Exact Solutions of Fractional-Order Biological Population Model A.M.A.
More informationA Modification in Successive Approximation Method for Solving Nonlinear Volterra Hammerstein Integral Equations of the Second Kind
Journal of Mathematical Extension Vol. 8, No. 1, (214), 69-86 A Modification in Successive Approximation Method for Solving Nonlinear Volterra Hammerstein Integral Equations of the Second Kind Sh. Javadi
More informationNumerical Solution of Two-Dimensional Volterra Integral Equations by Spectral Galerkin Method
Journal of Applied Mathematics & Bioinformatics, vol.1, no.2, 2011, 159-174 ISSN: 1792-6602 (print), 1792-6939 (online) International Scientific Press, 2011 Numerical Solution of Two-Dimensional Volterra
More informationRational Chebyshev pseudospectral method for long-short wave equations
Journal of Physics: Conference Series PAPER OPE ACCESS Rational Chebyshev pseudospectral method for long-short wave equations To cite this article: Zeting Liu and Shujuan Lv 07 J. Phys.: Conf. Ser. 84
More informationProperties of BPFs for Approximating the Solution of Nonlinear Fredholm Integro Differential Equation
Applied Mathematical Sciences, Vol. 6, 212, no. 32, 1563-1569 Properties of BPFs for Approximating the Solution of Nonlinear Fredholm Integro Differential Equation Ahmad Shahsavaran 1 and Abar Shahsavaran
More informationAn Efficient Numerical Method for Solving. the Fractional Diffusion Equation
Journal of Applied Mathematics & Bioinformatics, vol.1, no.2, 2011, 1-12 ISSN: 1792-6602 (print), 1792-6939 (online) International Scientific Press, 2011 An Efficient Numerical Method for Solving the Fractional
More informationAn Analytic Study of the (2 + 1)-Dimensional Potential Kadomtsev-Petviashvili Equation
Adv. Theor. Appl. Mech., Vol. 3, 21, no. 11, 513-52 An Analytic Study of the (2 + 1)-Dimensional Potential Kadomtsev-Petviashvili Equation B. Batiha and K. Batiha Department of Mathematics, Faculty of
More informationChebyshev finite difference method for solving a mathematical model arising in wastewater treatment plants
Computational Methods for Differential Equations http://cmde.tabrizu.ac.ir Vol. 6, No. 4, 2018, pp. 448-455 Chebyshev finite difference method for solving a mathematical model arising in wastewater treatment
More informationDecentralized control synthesis for bilinear systems using orthogonal functions
Cent Eur J Eng 41 214 47-53 DOI: 12478/s13531-13-146-1 Central European Journal of Engineering Decentralized control synthesis for bilinear systems using orthogonal functions Research Article Mohamed Sadok
More informationNUMERICAL SOLUTION OF DELAY DIFFERENTIAL EQUATIONS VIA HAAR WAVELETS
TWMS J Pure Appl Math V5, N2, 24, pp22-228 NUMERICAL SOLUTION OF DELAY DIFFERENTIAL EQUATIONS VIA HAAR WAVELETS S ASADI, AH BORZABADI Abstract In this paper, Haar wavelet benefits are applied to the delay
More informationA Numerical Approach for Solving Optimal Control Problems Using the Boubaker Polynomials Expansion Scheme
2014 (2014) 1-18 Available online at www.ispacs.com/jiasc Volume 2014, Year 2014 Article ID jiasc-000, 18 Pages doi:10.5899/2014/jiasc-000 Research Article A Numerical Approach for Solving Optimal Control
More informationAnalytical solution for determination the control parameter in the inverse parabolic equation using HAM
Available at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 1932-9466 Vol. 12, Issue 2 (December 2017, pp. 1072 1087 Applications and Applied Mathematics: An International Journal (AAM Analytical solution
More informationApplication of fractional-order Bernoulli functions for solving fractional Riccati differential equation
Int. J. Nonlinear Anal. Appl. 8 (2017) No. 2, 277-292 ISSN: 2008-6822 (electronic) http://dx.doi.org/10.22075/ijnaa.2017.1476.1379 Application of fractional-order Bernoulli functions for solving fractional
More informationExistence of solution and solving the integro-differential equations system by the multi-wavelet Petrov-Galerkin method
Int. J. Nonlinear Anal. Appl. 7 (6) No., 7-8 ISSN: 8-68 (electronic) http://dx.doi.org/.75/ijnaa.5.37 Existence of solution and solving the integro-differential equations system by the multi-wavelet Petrov-Galerkin
More informationOptimal Control of Nonlinear Systems Using the Shifted Legendre Polynomials
Optial Control of Nonlinear Systes Using Shifted Legendre Polynoi Rahan Hajohaadi 1, Mohaad Ali Vali 2, Mahoud Saavat 3 1- Departent of Electrical Engineering, Shahid Bahonar University of Keran, Keran,
More informationdifferentiable functions in all arguments. Our aim is to minimize the quadratic objective functional (x T (t)qx(t)+u T (t)ru(t))dt, (2)
SOLVING NON-LINEAR QUADRATIC OPTIMAL... 49 differentiable functions in all arguments. Our aim is to minimize the quadratic objective functional J[x, u] = 1 2 tf t 0 (x T (t)qx(t)+u T (t)ru(t))dt, (2) subject
More informationESTIMATES OF THE TRACE OF THE INVERSE OF A SYMMETRIC MATRIX USING THE MODIFIED CHEBYSHEV ALGORITHM
ESTIMATES OF THE TRACE OF THE INVERSE OF A SYMMETRIC MATRIX USING THE MODIFIED CHEBYSHEV ALGORITHM GÉRARD MEURANT In memory of Gene H. Golub Abstract. In this paper we study how to compute an estimate
More informationAn Elegant Perturbation Iteration Algorithm for the Lane-Emden Equation
Volume 32 - No.6, December 205 An Elegant Perturbation Iteration Algorithm for the Lane-Emden Equation M. Khalid Department of Mathematical Sciences Federal Urdu University of Arts, Sciences & Techonology
More informationAdomian Decomposition Method with Laguerre Polynomials for Solving Ordinary Differential Equation
J. Basic. Appl. Sci. Res., 2(12)12236-12241, 2012 2012, TextRoad Publication ISSN 2090-4304 Journal of Basic and Applied Scientific Research www.textroad.com Adomian Decomposition Method with Laguerre
More informationThe restarted QR-algorithm for eigenvalue computation of structured matrices
Journal of Computational and Applied Mathematics 149 (2002) 415 422 www.elsevier.com/locate/cam The restarted QR-algorithm for eigenvalue computation of structured matrices Daniela Calvetti a; 1, Sun-Mi
More informationNumerical Solution of Fredholm Integro-differential Equations By Using Hybrid Function Operational Matrix of Differentiation
Available online at http://ijim.srbiau.ac.ir/ Int. J. Industrial Mathematics (ISS 8-56) Vol. 9, o. 4, 7 Article ID IJIM-8, pages Research Article umerical Solution of Fredholm Integro-differential Equations
More informationApplied Numerical Analysis
Applied Numerical Analysis Using MATLAB Second Edition Laurene V. Fausett Texas A&M University-Commerce PEARSON Prentice Hall Upper Saddle River, NJ 07458 Contents Preface xi 1 Foundations 1 1.1 Introductory
More informationNumerical integration formulas of degree two
Applied Numerical Mathematics 58 (2008) 1515 1520 www.elsevier.com/locate/apnum Numerical integration formulas of degree two ongbin Xiu epartment of Mathematics, Purdue University, West Lafayette, IN 47907,
More informationIntroduction to Numerical Analysis
J. Stoer R. Bulirsch Introduction to Numerical Analysis Second Edition Translated by R. Bartels, W. Gautschi, and C. Witzgall With 35 Illustrations Springer Contents Preface to the Second Edition Preface
More informationSpecial Functions of Mathematical Physics
Arnold F. Nikiforov Vasilii B. Uvarov Special Functions of Mathematical Physics A Unified Introduction with Applications Translated from the Russian by Ralph P. Boas 1988 Birkhäuser Basel Boston Table
More informationThe Approximate Solution of Non Linear Fredholm Weakly Singular Integro-Differential equations by Using Chebyshev polynomials of the First kind
AUSTRALIA JOURAL OF BASIC AD APPLIED SCIECES ISS:1991-8178 EISS: 2309-8414 Journal home page: wwwajbaswebcom The Approximate Solution of on Linear Fredholm Weakly Singular Integro-Differential equations
More informationDifferential transformation method for solving one-space-dimensional telegraph equation
Volume 3, N 3, pp 639 653, 2 Copyright 2 SBMAC ISSN -825 wwwscielobr/cam Differential transformation method for solving one-space-dimensional telegraph equation B SOLTANALIZADEH Young Researchers Club,
More informationImproved Newton s method with exact line searches to solve quadratic matrix equation
Journal of Computational and Applied Mathematics 222 (2008) 645 654 wwwelseviercom/locate/cam Improved Newton s method with exact line searches to solve quadratic matrix equation Jian-hui Long, Xi-yan
More informationOrthogonal polynomials
Orthogonal polynomials Gérard MEURANT October, 2008 1 Definition 2 Moments 3 Existence 4 Three-term recurrences 5 Jacobi matrices 6 Christoffel-Darboux relation 7 Examples of orthogonal polynomials 8 Variable-signed
More informationEN Applied Optimal Control Lecture 8: Dynamic Programming October 10, 2018
EN530.603 Applied Optimal Control Lecture 8: Dynamic Programming October 0, 08 Lecturer: Marin Kobilarov Dynamic Programming (DP) is conerned with the computation of an optimal policy, i.e. an optimal
More informationIntroduction. Chapter One
Chapter One Introduction The aim of this book is to describe and explain the beautiful mathematical relationships between matrices, moments, orthogonal polynomials, quadrature rules and the Lanczos and
More informationA Comparison between Solving Two Dimensional Integral Equations by the Traditional Collocation Method and Radial Basis Functions
Applied Mathematical Sciences, Vol. 5, 2011, no. 23, 1145-1152 A Comparison between Solving Two Dimensional Integral Equations by the Traditional Collocation Method and Radial Basis Functions Z. Avazzadeh
More informationThe Chebyshev Collection Method for Solving Fractional Order Klein-Gordon Equation
The Chebyshev Collection Method for Solving Fractional Order Klein-Gordon Equation M. M. KHADER Faculty of Science, Benha University Department of Mathematics Benha EGYPT mohamedmbd@yahoo.com N. H. SWETLAM
More informationMODIFIED LAGUERRE WAVELET BASED GALERKIN METHOD FOR FRACTIONAL AND FRACTIONAL-ORDER DELAY DIFFERENTIAL EQUATIONS
MODIFIED LAGUERRE WAVELET BASED GALERKIN METHOD FOR FRACTIONAL AND FRACTIONAL-ORDER DELAY DIFFERENTIAL EQUATIONS Aydin SECER *,Neslihan OZDEMIR Yildiz Technical University, Department of Mathematical Engineering,
More informationMULTISTAGE HOMOTOPY ANALYSIS METHOD FOR SOLVING NON- LINEAR RICCATI DIFFERENTIAL EQUATIONS
MULTISTAGE HOMOTOPY ANALYSIS METHOD FOR SOLVING NON- LINEAR RICCATI DIFFERENTIAL EQUATIONS Hossein Jafari & M. A. Firoozjaee Young Researchers club, Islamic Azad University, Jouybar Branch, Jouybar, Iran
More informationComparing Numerical Methods for Solving Nonlinear Fractional Order Differential Equations
Comparing Numerical Methods for Solving Nonlinear Fractional Order Differential Equations Farhad Farokhi, Mohammad Haeri, and Mohammad Saleh Tavazoei Abstract This paper is a result of comparison of some
More informationResearch Article A Legendre Wavelet Spectral Collocation Method for Solving Oscillatory Initial Value Problems
Applied Mathematics Volume 013, Article ID 591636, 5 pages http://dx.doi.org/10.1155/013/591636 Research Article A Legendre Wavelet Spectral Collocation Method for Solving Oscillatory Initial Value Problems
More informationACTA UNIVERSITATIS APULENSIS No 20/2009 AN EFFECTIVE METHOD FOR SOLVING FRACTIONAL INTEGRO-DIFFERENTIAL EQUATIONS. Wen-Hua Wang
ACTA UNIVERSITATIS APULENSIS No 2/29 AN EFFECTIVE METHOD FOR SOLVING FRACTIONAL INTEGRO-DIFFERENTIAL EQUATIONS Wen-Hua Wang Abstract. In this paper, a modification of variational iteration method is applied
More informationHybrid Functions Approach for the Fractional Riccati Differential Equation
Filomat 30:9 (2016), 2453 2463 DOI 10.2298/FIL1609453M Published by Faculty of Sciences and Mathematics, University of Niš, Serbia Available at: http://www.pmf.ni.ac.rs/filomat Hybrid Functions Approach
More informationON THE NUMERICAL SOLUTION FOR THE FRACTIONAL WAVE EQUATION USING LEGENDRE PSEUDOSPECTRAL METHOD
International Journal of Pure and Applied Mathematics Volume 84 No. 4 2013, 307-319 ISSN: 1311-8080 (printed version); ISSN: 1314-3395 (on-line version) url: http://www.ijpam.eu doi: http://dx.doi.org/10.12732/ijpam.v84i4.1
More informationNumerical analysis of the one-dimensional Wave equation subject to a boundary integral specification
Numerical analysis of the one-dimensional Wave equation subject to a boundary integral specification B. Soltanalizadeh a, Reza Abazari b, S. Abbasbandy c, A. Alsaedi d a Department of Mathematics, University
More informationNumerical solution of linear time delay systems using Chebyshev-tau spectral method
Available at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 93-9466 Applications and Applied Mathematics: An International Journal (AAM) Vol., Issue (June 07), pp. 445 469 Numerical solution of linear time
More informationFractional Spectral and Spectral Element Methods
Fractional Calculus, Probability and Non-local Operators: Applications and Recent Developments Nov. 6th - 8th 2013, BCAM, Bilbao, Spain Fractional Spectral and Spectral Element Methods (Based on PhD thesis
More information(Received 10 December 2011, accepted 15 February 2012) x x 2 B(x) u, (1.2) A(x) +
ISSN 1749-3889 (print), 1749-3897 (online) International Journal of Nonlinear Science Vol13(212) No4,pp387-395 Numerical Solution of Fokker-Planck Equation Using the Flatlet Oblique Multiwavelets Mir Vahid
More informationVariational iteration method for solving multispecies Lotka Volterra equations
Computers and Mathematics with Applications 54 27 93 99 www.elsevier.com/locate/camwa Variational iteration method for solving multispecies Lotka Volterra equations B. Batiha, M.S.M. Noorani, I. Hashim
More informationAn efficient hybrid pseudo-spectral method for solving optimal control of Volterra integral systems
MATHEMATICAL COMMUNICATIONS 417 Math. Commun. 19(14), 417 435 An efficient hybrid pseudo-spectral method for solving optimal control of Volterra integral systems Khosrow Maleknejad 1, and Asyieh Ebrahimzadeh
More informationA numerical method for solving Linear Non-homogenous Fractional Ordinary Differential Equation
Appl. Math. Inf. Sci. 6, No. 3, 441-445 (2012) 441 Applied Mathematics & Information Sciences An International Journal A numerical method for solving Linear Non-homogenous Fractional Ordinary Differential
More informationA New Operational Matrix of Derivative for Orthonormal Bernstein Polynomial's
A New Operational Matrix of Derivative for Orthonormal Bernstein Polynomial's MayadaN.Mohammed Ali* Received 16, May, 2013 Accepted 26, September, 2013 Abstract: In this paper, an orthonormal family has
More informationRecurrence Relations and Fast Algorithms
Recurrence Relations and Fast Algorithms Mark Tygert Research Report YALEU/DCS/RR-343 December 29, 2005 Abstract We construct fast algorithms for decomposing into and reconstructing from linear combinations
More informationNumerical Optimal Control Overview. Moritz Diehl
Numerical Optimal Control Overview Moritz Diehl Simplified Optimal Control Problem in ODE path constraints h(x, u) 0 initial value x0 states x(t) terminal constraint r(x(t )) 0 controls u(t) 0 t T minimize
More informationNumerical Investigation of the Time Invariant Optimal Control of Singular Systems Using Adomian Decomposition Method
Applied Mathematical Sciences, Vol. 8, 24, no. 2, 6-68 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/.2988/ams.24.4863 Numerical Investigation of the Time Invariant Optimal Control of Singular Systems
More informationHAAR AND LEGENDRE WAVELETS COLLOCATION METHODS FOR THE NUMERICAL SOLUTION OF SCHRODINGER AND WAVE EQUATIONS. H. Kheiri and H.
Acta Universitatis Apulensis ISSN: 1582-5329 No. 37/214 pp. 1-14 HAAR AND LEGENDRE WAVELETS COLLOCATION METHODS FOR THE NUMERICAL SOLUTION OF SCHRODINGER AND WAVE EQUATIONS H. Kheiri and H. Ghafouri Abstract.
More informationChapter 2 Optimal Control Problem
Chapter 2 Optimal Control Problem Optimal control of any process can be achieved either in open or closed loop. In the following two chapters we concentrate mainly on the first class. The first chapter
More informationAn Accurate Fourier-Spectral Solver for Variable Coefficient Elliptic Equations
An Accurate Fourier-Spectral Solver for Variable Coefficient Elliptic Equations Moshe Israeli Computer Science Department, Technion-Israel Institute of Technology, Technion city, Haifa 32000, ISRAEL Alexander
More informationBLOCK-PULSE FUNCTIONS AND THEIR APPLICATIONS TO SOLVING SYSTEMS OF HIGHER-ORDER NONLINEAR VOLTERRA INTEGRO-DIFFERENTIAL EQUATIONS
Electronic Journal of Differential Equations, Vol. 214 (214), No. 54, pp. 1 9. ISSN: 172-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu BLOCK-PULSE FUNCTIONS
More informationDivya Garg Michael A. Patterson Camila Francolin Christopher L. Darby Geoffrey T. Huntington William W. Hager Anil V. Rao
Comput Optim Appl (2011) 49: 335 358 DOI 10.1007/s10589-009-9291-0 Direct trajectory optimization and costate estimation of finite-horizon and infinite-horizon optimal control problems using a Radau pseudospectral
More informationA Numerical Scheme for Generalized Fractional Optimal Control Problems
Available at http://pvamuedu/aam Appl Appl Math ISSN: 1932-9466 Vol 11, Issue 2 (December 216), pp 798 814 Applications and Applied Mathematics: An International Journal (AAM) A Numerical Scheme for Generalized
More informationHomotopy Analysis Method for Nonlinear Differential Equations with Fractional Orders
Homotopy Analysis Method for Nonlinear Differential Equations with Fractional Orders Yin-Ping Liu and Zhi-Bin Li Department of Computer Science, East China Normal University, Shanghai, 200062, China Reprint
More informationInstructions for Matlab Routines
Instructions for Matlab Routines August, 2011 2 A. Introduction This note is devoted to some instructions to the Matlab routines for the fundamental spectral algorithms presented in the book: Jie Shen,
More informationNUMERICAL SOLUTION OF FRACTIONAL ORDER DIFFERENTIAL EQUATIONS USING HAAR WAVELET OPERATIONAL MATRIX
Palestine Journal of Mathematics Vol. 6(2) (217), 515 523 Palestine Polytechnic University-PPU 217 NUMERICAL SOLUTION OF FRACTIONAL ORDER DIFFERENTIAL EQUATIONS USING HAAR WAVELET OPERATIONAL MATRIX Raghvendra
More informationEngineering Notes. NUMERICAL methods for solving optimal control problems fall
JOURNAL OF GUIANCE, CONTROL, AN YNAMICS Vol. 9, No. 6, November ecember 6 Engineering Notes ENGINEERING NOTES are short manuscripts describing ne developments or important results of a preliminary nature.
More informationComputers and Mathematics with Applications. A new application of He s variational iteration method for the solution of the one-phase Stefan problem
Computers and Mathematics with Applications 58 (29) 2489 2494 Contents lists available at ScienceDirect Computers and Mathematics with Applications journal homepage: www.elsevier.com/locate/camwa A new
More informationSpectral Collocation Method for Parabolic Partial Differential Equations with Neumann Boundary Conditions
Applied Mathematical Sciences, Vol. 1, 2007, no. 5, 211-218 Spectral Collocation Method for Parabolic Partial Differential Equations with Neumann Boundary Conditions M. Javidi a and A. Golbabai b a Department
More informationApplication of the Bernstein Polynomials. for Solving the Nonlinear Fredholm. Integro-Differential Equations
Journal of Applied Mathematics & Bioinformatics, vol., no.,, 3-3 ISSN: 79-66 (print), 79-6939 (online) International Scientific Press, Application of the Bernstein Polynomials for Solving the Nonlinear
More informationADOMIAN-TAU OPERATIONAL METHOD FOR SOLVING NON-LINEAR FREDHOLM INTEGRO-DIFFERENTIAL EQUATIONS WITH PADE APPROXIMANT. A. Khani
Acta Universitatis Apulensis ISSN: 1582-5329 No 38/214 pp 11-22 ADOMIAN-TAU OPERATIONAL METHOD FOR SOLVING NON-LINEAR FREDHOLM INTEGRO-DIFFERENTIAL EQUATIONS WITH PADE APPROXIMANT A Khani Abstract In this
More informationInterpolation and Cubature at Geronimus Nodes Generated by Different Geronimus Polynomials
Interpolation and Cubature at Geronimus Nodes Generated by Different Geronimus Polynomials Lawrence A. Harris Abstract. We extend the definition of Geronimus nodes to include pairs of real numbers where
More informationSolving a class of nonlinear two-dimensional Volterra integral equations by using two-dimensional triangular orthogonal functions.
Journal of Mathematical Modeling Vol 1, No 1, 213, pp 28-4 JMM Solving a class of nonlinear two-dimensional Volterra integral equations by using two-dimensional triangular orthogonal functions Farshid
More informationSolving a class of linear and non-linear optimal control problems by homotopy perturbation method
IMA Journal of Mathematical Control and Information (2011) 28, 539 553 doi:101093/imamci/dnr018 Solving a class of linear and non-linear optimal control problems by homotopy perturbation method S EFFATI
More informationInterval solutions for interval algebraic equations
Mathematics and Computers in Simulation 66 (2004) 207 217 Interval solutions for interval algebraic equations B.T. Polyak, S.A. Nazin Institute of Control Sciences, Russian Academy of Sciences, 65 Profsoyuznaya
More informationApplication of Variational Iteration Method to a General Riccati Equation
International Mathematical Forum,, 007, no. 56, 759-770 Application of Variational Iteration Method to a General Riccati Equation B. Batiha, M. S. M. Noorani and I. Hashim School of Mathematical Sciences
More informationFour Point Gauss Quadrature Runge Kuta Method Of Order 8 For Ordinary Differential Equations
International journal of scientific and technical research in engineering (IJSTRE) www.ijstre.com Volume Issue ǁ July 206. Four Point Gauss Quadrature Runge Kuta Method Of Order 8 For Ordinary Differential
More informationApproximate Solution of an Integro-Differential Equation Arising in Oscillating Magnetic Fields Using the Differential Transformation Method
IOSR Journal of Mathematics (IOSR-JM) e-issn: 2278-5728, p-issn: 2319-765X. Volume 13, Issue 5 Ver. I1 (Sep. - Oct. 2017), PP 90-97 www.iosrjournals.org Approximate Solution of an Integro-Differential
More informationA Differential Quadrature Algorithm for the Numerical Solution of the Second-Order One Dimensional Hyperbolic Telegraph Equation
ISSN 1749-3889 (print), 1749-3897 (online) International Journal of Nonlinear Science Vol.13(01) No.3,pp.59-66 A Differential Quadrature Algorithm for the Numerical Solution of the Second-Order One Dimensional
More informationA collocation method for solving the fractional calculus of variation problems
Bol. Soc. Paran. Mat. (3s.) v. 35 1 (2017): 163 172. c SPM ISSN-2175-1188 on line ISSN-00378712 in press SPM: www.spm.uem.br/bspm doi:10.5269/bspm.v35i1.26333 A collocation method for solving the fractional
More informationNumerical Analysis of Electromagnetic Fields
Pei-bai Zhou Numerical Analysis of Electromagnetic Fields With 157 Figures Springer-Verlag Berlin Heidelberg New York London Paris Tokyo Hong Kong Barcelona Budapest Contents Part 1 Universal Concepts
More informationSMOOTHNESS PROPERTIES OF SOLUTIONS OF CAPUTO- TYPE FRACTIONAL DIFFERENTIAL EQUATIONS. Kai Diethelm. Abstract
SMOOTHNESS PROPERTIES OF SOLUTIONS OF CAPUTO- TYPE FRACTIONAL DIFFERENTIAL EQUATIONS Kai Diethelm Abstract Dedicated to Prof. Michele Caputo on the occasion of his 8th birthday We consider ordinary fractional
More informationThe approximation of solutions for second order nonlinear oscillators using the polynomial least square method
Available online at www.isr-publications.com/jnsa J. Nonlinear Sci. Appl., 1 (217), 234 242 Research Article Journal Homepage: www.tjnsa.com - www.isr-publications.com/jnsa The approximation of solutions
More information