(X) 1. Introduction Orthogonal Families of Polynomial Are Impo
Total Page:16
File Type:pdf, Size:1020Kb
Orthogonality for a class of generalised Jacobi polynomial $P^{\alpha,\beta}_{\nu}(x)$ Item Type Article Authors Ford, Neville J.; Moayyed, H.; Rodrigues, M. M. Citation Ford, N. J., Moayyed, H. & Rodrigues, M. M. (2018) Orthogonality for a class of generalised Jacobi polynomial $P^{\alpha, \beta}_{\nu}(x)$. Fractional Differential Calculus, 8(1), 95-110 DOI 10.7153/fdc-2018-08-06 Publisher Ele-Math Journal Fractional Differential Calculus Download date 02/10/2021 17:29:02 Item License http://creativecommons.org/licenses/by-nc-nd/4.0/ Link to Item http://hdl.handle.net/10034/620536 F ractional D ifferential C alculus www.ele-math.com ORTHOGONALITY FOR A CLASS OF a;b GENERALISED JACOBI POLYNOMIAL Pn (x) N. J. FORD, H. MOAYYED AND M. M. RODRIGUES Submitted to Fracttional Differ. Calc. Abstract. This work considers g-Jacobi polynomials, a fractional generalisation of the classical Jacobi polynomials. We discuss the polynomials and compare some of their properties to the classical case. The main result of the paper is to show that one can derive an orthogonality prop- (a;b) erty for a sub-class of g-Jacobi polynomials Pn . The paper concludes with an application in modelling of ophthalmic surfaces. 1. Introduction Orthogonal families of polynomial are important both from theoretical and practi- cal viewpoints. They are widely used when approximating functions from mathematics and other sciences and they allow us to solve and interpret solutions of certain differen- tial equations. (a;b) The Jacobi Polynomials Pn are a classical family of orthogonal polynomials (see [7], [23]) that has been used in many applications due to its ability to approximate general classes of function. The Gegenbauer polynomials (which are special cases of Jacobi polynomials) are orthogonal polynomials on the interval [−1;1] with respect to the weight function (1 − x2)a−1=2 . They are a generalisation of the Legendre and Chebyshev polynomi- l als. Further, the Gegenbauer polynomials fCn (t)g constitute a complete orthogonal − 1 system of the real (or complex) L2−space L2 [−1;1];(1 − x2)n 2 dx with respect to the (1 − x2)a−1=2 on the interval [−1;1] (see [24]). One can also consider the question of orthogonality over a semi-infinite interval. The authors of [6] discuss this property for the Jacobi polynomials and this was also considered by V. Romanovski [20] and (more recently) by S.D. Bajpai[4]. In[11] the authors introduced a family of generalised Jacobi polynomials/functions with indexes a;b 2 R which are mutually orthogonal with respect to the corresponding Mathematics subject classification (2010): 33F05, 33C90, 12E10. Keywords and phrases: Jacobi polynomials; Approximation; Optic modelling.. The work of second and third authors were supported by Portuguese funds through the CIDMA - Center for Research and Development in Mathematics and Applications, and the Portuguese Foundation for Science and Technology (“FCT– Fundac¸ao˜ para a Cienciaˆ e a Tecnologia”), within project UID/MAT/ 0416/2013. Additionally the research grant (BIP- LtPICS-2016A) of the second author was also supported by FCT / MEC through national funds and where appropriate co-financed by the FEDER, under the PT2020 Partnership Agreement.. 1 Jacobi weights and which inherit selected important properties of the classical Jacobi polynomials. a;b;M;N ¥ In [13] the authors deal with generalised Jacobi polynomials fPn (x)gn=0 , a;b but the generalisation we consider in this paper, to the form fPn (x)g with n non in- teger is much more recent (see [15]). We will show that these fractional generalisations can have important properties of use in applications. The classical Jacobi Polynomials are usually defined by a Rodrigues’ type formula dn h i P(a;b)(x) = (−2)n(n!)−1(1 − x)−a (1 + x)−b (1 − x)n+a (1 + x)n+b ; n dxn where a > −1;b > −1 (see [7]). In [15] the authors defined the generalised Jacobi (or g-Jacobi) functions as an extension of the Jacobi polynomials. In the previous formula they replaced the integer derivative by the Riemann-Liouville fractional of order n 2 R. They also derived some formal properties of this new family of polynomials. The main theme of this paper is to find a sub-class of the g-Jacobi Polynomials, (a;b) Pn , with n non-integer, where a orthogonality relation holds, and to show (through an application in optics) how this new class of orthogonal polynomials can have useful practical applications in function aproximation. The paper is organised as follows. In the second section, various definitions and properties are recalled. In the third section, we prove the orthogonality property for a certain class of generalised Jacobi polynomials and some of their properties will be presented in section four. Finally, we compare the behaviour of Zernike, classical and modified Jacobi polynomials in modelling of ophthalmic surfaces. 2. Preliminaries 2.1. Special functions Here we recall some aspects of some of the special functions that feature in this paper: the Gamma function and the Pochhammer symbol. We recall that the Gamma function G(s) is defined by the integral (see [3]) Z ¥ G(s) = e−tts−1dt: (1) 0 The integral (1) defines a holomorphic function in the half-plane Re(s) > 0. In addition G satisfies the functional equation G(s + 1) = sG(s); Re(s) > 0: (2) Hence, since G(1) = 1, we have G(n + 1) = n! for all n 2 N. The reflection formula for the Gamma function establishes a relation between the sine and gamma functions, i.e., p csc(pa) G(−a) = − ; a 2= : G(a + 1) Z 2 The Pochhammer symbol is defined as see [3] G(a + k) (a) := ; k G(a) − for all a 2 CnZ0 . The regularised hypergeometric function 2F˜1(a;b;c;z) is defined by means of a hypergeometric series by the expression ¥ a b ˜ ( )k ( )k k 2F1(a;b;c;z) = ∑ z ; jzj ≤ 1: (3) k=0 G(c + k) k! It is related to the hypergeometric function 2F1 by the formula 1 F˜ (a;b;c;z) = F (a;b;c;z); (4) 2 1 G(c) 2 1 and the hypergeometric function 2F1 has the following series representation ¥ (a)k (b)k k 2F1(a;b;c;z) = ∑ z ; jzj ≤ 1: (5) k=0 (c)k k! When z = 1 the hypergeometric function 2F1 assumes the particular form (see formula (7.3.5.2) in [18]) G(c) G(c − a − b) F (a;b;c;1) = : (6) 2 1 G(c − a) G(c − b) 2.2. Generalised Jacobi polynomials We start by reviewing and defining the Jacobi polynomials. (a;b) We recall that the classical Jacobi polynomials Pn (x) (n ≥ 0) are defined by dn h i P(a;b)(x) = (−2)n(n!)−1(1 − x)−a (1 + x)−b (1 − x)n+a (1 + x)n+b : (7) n dxn a b Let wa;b (x) := (1 − x) (1 + x) be the Jacobi weight function. For a;b > −1 the Jacobi polynomials are mutually orthogonal in L2 ([−1;1]), i.e., wa;b 1 a+b+1 Z (a;b) (a;b) 2 G(a + n + 1)G(b + n + 1) Pm (x)Pn (x)wa;b (x)dx = dmn ; (8) −1 (2n + a + b + 1)G(n + a + b + 1)n! where dmn denotes the Kronecker delta m;n 2 N0 . This class of Jacobi weight func- tions leads to Jacobi polynomials with many attractive properties that are not shared by general orthogonal polynomials. The generalised Jacobi polynomials may be considered as fractional index gener- alisations of the classical Jacobi polynomials. These polynomials appear as a special case of fractional Gauss functions, defined as solutions of the fractional generalisation of the Gauss hypergeometric equation (see [10]). The representation of the generalised Jacobi polynomial in terms of the hypergeo- metric function 2F1 is given by 3 THEOREM 1. (c.f. [15]) The g-Jacobi Polynomial can be represented by G(a + n + 1) (1 − x) Pa;b (x) = F −n;a + b + n + 1;a + 1; : (9) n G(n + 1))G(1 + a) 2 1 2 Moreover, (9) satisfies the following differential equation (1 − x2)Y 00(x) + (b − a − (a + b + 2)x)Y 0(x) + n(n + 2l)Y(x) = 0: (10) The g-Jacobi polynomials were introduced in [15]. There, the authors replaced the integer derivatives in (7) by the Riemann-Liouville fractional derivative Dn , obtaining (−1)n h i Pa;b (x) = Dn (1 − x)n+a (1 + x)n+b ; (11) n G(n + 1)2n (1 − x)a (1 + x)b where the Riemann-Liouville fractional derivative of order a is defined as 8 n 1 d R t f (t) dt; n − 1 < a < n 2 <> G(n−a) dtn a (t−t)a+1−n N Da f (t) = (12) :> dn dtn ; a = n 2 N (see [21]). In [15] the authors presented the following explicit formula for the g-Jacobi poly- nomials THEOREM 2. (cf. [15]) For the generalised Jacobi functions the following repre- sentation holds ¥ ( ; ) n + b n + a P a b (x) = 2−n (x + 1)n−k(x − 1)k (13) n ∑ k n − k k=0 with n > 0;a > −1;b > −1 and where a G(1 + a) = b G(1 + b)G(1 + a − b) is the binomial coefficient with real arguments. We know that the classical Jacobi polynomials satisfy an orthogonality property. In the following example we can see that the orthogonality property (in general) for the g-Jacobi polynomials is not applicable.