The Bernstein Operational Matrix of Integration

Similar documents
VARIATIONAL ITERATION METHOD FOR DELAY DIFFERENTIAL-ALGEBRAIC EQUATIONS. Hunan , China,

FORCED VIBRATION of MDOF SYSTEMS

QR factorization. Let P 1, P 2, P n-1, be matrices such that Pn 1Pn 2... PPA

Solution of Impulsive Differential Equations with Boundary Conditions in Terms of Integral Equations

Key words: Fractional difference equation, oscillatory solutions,

Mixed Integral Equation of Contact Problem in Position and Time

Quintic B-Spline Collocation for Solving Abel's Integral Equation

Fully Fuzzy Linear Systems Solving Using MOLP

A Second Kind Chebyshev Polynomial Approach for the Wave Equation Subject to an Integral Conservation Condition

AML710 CAD LECTURE 12 CUBIC SPLINE CURVES. Cubic Splines Matrix formulation Normalised cubic splines Alternate end conditions Parabolic blending

Numerical Solution for the Variable Order Fractional Partial Differential

Fourth Order Runge-Kutta Method Based On Geometric Mean for Hybrid Fuzzy Initial Value Problems

Real-Time Systems. Example: scheduling using EDF. Feasibility analysis for EDF. Example: scheduling using EDF

Continuous Time Markov Chains

The Poisson Process Properties of the Poisson Process

The Linear Regression Of Weighted Segments

Determination of Antoine Equation Parameters. December 4, 2012 PreFEED Corporation Yoshio Kumagae. Introduction

Moments of Order Statistics from Nonidentically Distributed Three Parameters Beta typei and Erlang Truncated Exponential Variables

14. Poisson Processes

Asymptotic Behavior of Solutions of Nonlinear Delay Differential Equations With Impulse

Available online Journal of Scientific and Engineering Research, 2014, 1(1): Research Article

The Mean Residual Lifetime of (n k + 1)-out-of-n Systems in Discrete Setting

On an algorithm of the dynamic reconstruction of inputs in systems with time-delay

COMPARISON OF ESTIMATORS OF PARAMETERS FOR THE RAYLEIGH DISTRIBUTION

Least squares and motion. Nuno Vasconcelos ECE Department, UCSD

FALL HOMEWORK NO. 6 - SOLUTION Problem 1.: Use the Storage-Indication Method to route the Input hydrograph tabulated below.

Solving Fuzzy Equations Using Neural Nets with a New Learning Algorithm

Pricing Asian Options with Fourier Convolution

Stability Criterion for BAM Neural Networks of Neutral- Type with Interval Time-Varying Delays

Least Squares Fitting (LSQF) with a complicated function Theexampleswehavelookedatsofarhavebeenlinearintheparameters

An Exact Solution for the Differential Equation. Governing the Lateral Motion of Thin Plates. Subjected to Lateral and In-Plane Loadings

The algebraic immunity of a class of correlation immune H Boolean functions

Convexity Preserving C 2 Rational Quadratic Trigonometric Spline

Quantum Mechanics II Lecture 11 Time-dependent perturbation theory. Time-dependent perturbation theory (degenerate or non-degenerate starting state)

The Properties of Probability of Normal Chain

Modified Integrated Multi-Point Approximation And GA Used In Truss Topology Optimization

Solving fuzzy linear programming problems with piecewise linear membership functions by the determination of a crisp maximizing decision

Real-time Classification of Large Data Sets using Binary Knapsack

General Complex Fuzzy Transformation Semigroups in Automata

Comparison of the Bayesian and Maximum Likelihood Estimation for Weibull Distribution

The Signal, Variable System, and Transformation: A Personal Perspective

International Journal Of Engineering And Computer Science ISSN: Volume 5 Issue 12 Dec. 2016, Page No.

Bianchi Type II Stiff Fluid Tilted Cosmological Model in General Relativity

Probability Bracket Notation and Probability Modeling. Xing M. Wang Sherman Visual Lab, Sunnyvale, CA 94087, USA. Abstract

Optimal Eye Movement Strategies in Visual Search (Supplement)

The Optimal Combination Forecasting Based on ARIMA,VAR and SSM

θ = θ Π Π Parametric counting process models θ θ θ Log-likelihood: Consider counting processes: Score functions:

Iterated Bernstein polynomial approximations

As evident from the full-sample-model, we continue to assume that individual errors are identically and

IMPROVED PORTFOLIO OPTIMIZATION MODEL WITH TRANSACTION COST AND MINIMAL TRANSACTION LOTS

Solving Non-Linear Rational Expectations Models: Approximations based on Taylor Expansions

Partial Molar Properties of solutions

New approach for numerical solution of Fredholm integral equations system of the second kind by using an expansion method

Some Probability Inequalities for Quadratic Forms of Negatively Dependent Subgaussian Random Variables

Chapter 8. Simple Linear Regression

Equivalent Finite Element Formulations for the Calculation of Eigenvalues Using Higher-Order Polynomials

Other Topics in Kernel Method Statistical Inference with Reproducing Kernel Hilbert Space

Solution. The straightforward approach is surprisingly difficult because one has to be careful about the limits.

Integral Φ0-Stability of Impulsive Differential Equations

Delay-Dependent Robust Asymptotically Stable for Linear Time Variant Systems

Continuous Indexed Variable Systems

Time-Dependent Perturbation Theory

Regression Approach to Parameter Estimation of an Exponential Software Reliability Model

Solution set Stat 471/Spring 06. Homework 2

Orbital Euclidean stability of the solutions of impulsive equations on the impulsive moments

Chebyshev Polynomials for Solving a Class of Singular Integral Equations

Newton-Product Integration for a Stefan Problem with Kinetics

Asymptotic Regional Boundary Observer in Distributed Parameter Systems via Sensors Structures

CONTROLLABILITY OF A CLASS OF SINGULAR SYSTEMS

4. THE DENSITY MATRIX

Supplement Material for Inverse Probability Weighted Estimation of Local Average Treatment Effects: A Higher Order MSE Expansion

Integral Equations and their Relationship to Differential Equations with Initial Conditions

STABILITY CRITERION FOR HYBRID SYSTEMS WITH DELAY. Laviniu Bejenaru

GENERALIZED METHOD OF LIE-ALGEBRAIC DISCRETE APPROXIMATIONS FOR SOLVING CAUCHY PROBLEMS WITH EVOLUTION EQUATION

Voltage Sensitivity Analysis in MV Distribution Networks

Stochastic Petri Nets with Low Variation Matrix Exponentially Distributed Firing Time

Model for Optimal Management of the Spare Parts Stock at an Irregular Distribution of Spare Parts

JORIND 9(2) December, ISSN

Reliability Analysis of Sparsely Connected Consecutive-k Systems: GERT Approach

A Comparison of AdomiansDecomposition Method and Picard Iterations Method in Solving Nonlinear Differential Equations

Existence and Uniqueness Theorems for Generalized Set Differential Equations

The ray paths and travel times for multiple layers can be computed using ray-tracing, as demonstrated in Lab 3.

Mathematical Formulation

Average Consensus in Networks of Multi-Agent with Multiple Time-Varying Delays

An efficient approach to reliability-based design optimization within the enhanced sequential optimization and reliability assessment framework

AN INCREMENTAL QUASI-NEWTON METHOD WITH A LOCAL SUPERLINEAR CONVERGENCE RATE. Aryan Mokhtari Mark Eisen Alejandro Ribeiro

Stabilization of LTI Switched Systems with Input Time Delay. Engineering Letters, 14:2, EL_14_2_14 (Advance online publication: 16 May 2007) Lin Lin

Quantitative Portfolio Theory & Performance Analysis

A Constitutive Model for Multi-Line Simulation of Granular Material Behavior Using Multi-Plane Pattern

A note on Turán number Tk ( 1, kn, )

Research on portfolio model based on information entropy theory

Cyclone. Anti-cyclone

Fundamentals of Speech Recognition Suggested Project The Hidden Markov Model

Abstract. Keywords: Mutation probability, evolutionary computation, optimization, sensitivity, variability. 1. Introduction. 2. Proposed Algorithm

-distributed random variables consisting of n samples each. Determine the asymptotic confidence intervals for

Inner-Outer Synchronization Analysis of Two Complex Networks with Delayed and Non-Delayed Coupling

Numerical Methods for a Class of Hybrid. Weakly Singular Integro-Differential Equations.

STOCHASTIC CALCULUS I STOCHASTIC DIFFERENTIAL EQUATION

(1) Cov(, ) E[( E( ))( E( ))]

International Journal of Scientific & Engineering Research, Volume 3, Issue 10, October ISSN

Transcription:

Appled Mahemacal Sceces, Vol. 3, 29, o. 49, 2427-2436 he Berse Operaoal Marx of Iegrao Am K. Sgh, Vee K. Sgh, Om P. Sgh Deparme of Appled Mahemacs Isue of echology, Baaras Hdu Uversy Varaas -225, Ida Asrac A accurae mehod s proposed o solve prolems such as defcao, aalyss ad opmal corol usg he Berse orhoormal polyomals operaoal marx of egrao. he Berse polyomals are frs orhogoalzed, ormalzed ad he her operaoal marx of egrao s oaed. A example s gve o llusrae he proposed mehod. Mahemacs Sujec Classfcao: 4A, 49J5 Keywords: he Berse polyomals; Operaoal marx; Dffereal equao. Iroduco Approxmaos y orhoormal famly of fucos have played a val role he developme of physcal sceces, egeerg ad echology geeral ad mahemacal aalyss parcular sce log. I he las hree decades, hey have ee playg a mpora par he evaluao of ew echques o solve prolems such as defcao, aalyss ad opmal corol. he am of hese echques has ee o oa effecve algorhms ha are suale for he dgal compuers. he movao ad phlosophy ehd hs approach s ha rasforms he uderlyg dffereal equao of he prolem o a algerac equao, hus smplfyg he soluo process of he prolem o a grea exe. he asc dea of hs echque s as follows: Correspodg Auhor. E-mal : sghom@gmal.com (O. P. Sgh), vkshu@gmal.com (V. K. Sgh).

2428 A. K. Sgh, V. K. Sgh, O. P. Sgh () he dffereal equao s covered o a egral equao va mulple egrao. () Susequely, he varous sgals volved he egral equao are approxmaed y represeg hem as lear comaos of he orhoormal ass fucos ad rucag hem a opmal levels. () Fally, he egral equao s covered o a algerac equao y roducg he operaoal marx of egrao of he ass fucos. he key dea of he echque depeds o he followg egral propery of he ass vecor ϕ () a ϕ m+ a k k L ( σ )( dσ ) P ϕ( ), () where ϕ( ) = [ ϕ ( ), ϕ( ), K, ϕ m ( )] whch he elemes ϕ ( ), ϕ( ), K, ϕ m ( ) are he ass fucos, orhogoal o a cera erval [ a, ] ad P m+ s he operaoal marx for egrao of ϕ (). Noe ha P m+ s a cosa marx of order ( m + ) ( m + ). Usg he operaoal marx of a orhoormal sysem of fucos o perform egrao for solvg, defyg ad opmzg a lear dyamc sysem has several advaages: () he mehod s compuer oreed, hus solvg hgher order dffereal equao ecomes a maer of dmeso creasg, () he soluo s a mul-resoluo ype ad () he soluo s coverge, eve hough he sze of creme may e large. Ul ow, he operaoal marx of egrao has ee deermed for several ypes of orhogoal ass fucos, such as he Walsh fuco [-2], lock- pulse fuco [3-4], Laguerre seres [5-7], Cheyshev polyomals [8-9], Legedre polyomals [-], Fourer seres [2-3] ad Bessel seres [4]. Laer Gu ad Jag [5] derved he Haar waveles operaoal marx of egrao followed y Razzagh ad Yousef [6] who gave he Legedre waveles operaoal marx of egrao. he am of prese paper s o derve he Berse orhoormal polyomals marx of egrao P m+. he marx P m+ may e used o solve prolems of sysem aalyss ad syhess a maer smlar o hose of he oher orhogoal fucos. he Berse polyomals are frs orhoormalzed ad he operaoal marx of egrao s he derved. A umercal example s gve o llusrae he effcecy of he proposed mehod.

Berse operaoal marx of egrao 2429 2. he Berse polyomals: A Berse polyomal, amed afer Serge Naaovch Berse, s a polyomal he Berse form, ha s a lear comao of Berse ass polyomals. he Berse ass polyomals of degree are defed y B, ( ) = ( ), for =,, 2, L,. (2) h here are ( + ) degree Berse ass polyomals formg a ass for he lear space V cossg of all polyomals of degree less ha or equal o R[x]-he rg of polyomals over he feld R. For mahemacal coveece, we usually se B, = f < or >. Ay polyomal B (x) R[x] may e wre as B x) = β B ( x). (3) (, = he B (x) s called a polyomal Berse form or Berse polyomal of degree. he coeffces β are called Berse or Bezer coeffces.bu several mahemacas call Berse ass polyomals B, ( x) as he Berse polyomals. We wll follow hs coveo as well. hese polyomals have he followg properes: ( ) B, () = δ ad B, () = δ, where δ s he Kroecker dela fuco. ( ) B, ( ) has oe roo, each of mulplcy ad, a = ad = respecvely. ( ) B, ( ) for [,] ad B, ( ) = B, ( ). (v) For, B, has a uque local maxmum [,] a = / ad he maxmum value ( ). (v) he Berse polyomals form a paro of uy.e. B ( ) =. (v) I has a degree rasg propery he sese ha ay of he lower-degree polyomals (degree < ) ca e expressed as a lear comaos of polyomals of degree. We have, + B, ( ) = B, ( ) + B +, ( ). (v) Le f ( x) C [,] (he class of couous fucos o [,] ), he =,

243 A. K. Sgh, V. K. Sgh, O. P. Sgh B f )( x) = f B ( x) coverges o f (x) uformly o [,] as. = (, (k ) (k ) (v) Le f ( x) C [,] (he class of k mes dffereale fuco wh f couous), he ( k ) ( ) k ( k ) ( k ) ( k ) B ( f ) f ad f B ( ) k f as k, where. s he ( ) k k sup. orm ad = L k s a ege value of B ; he correspodg ege fuco s a polyomal of degree k. 3. he orhoormal polyomals: Usg Gram- Schmd orhoormalzao process o B, ad ormalzg, we oa a class of orhoormal polyomals from Berse polyomals. We call hem orhoormal Berse polyomals of order ad deoe hem y,, L,. o For = 5 he fve orhoormal polyomals are gve y 5 5 ( ) = ( ) 4 5 5 ( ) = 6 5( ) ( ) 2 (4) 8 7 3 2 5 ( ) = ( ) 5 5 8 4 5 ( ) + ( ) 25 28 2 3 3 2 3 4 5 ( ) = ( ) 5( ) + ( ) ( 5 7 28 5 35 ) 4 2 3 3 2 4 ( ) = 7 3 5( ) 2( ) + 8( ) 4( ) + ( 7 5 45 ) 5 25 4 2 3 3 2 4 ( ) = 6 ( ) + ( ) 25( ) + 5( ) ( 2 3 6 5 55 )

Berse operaoal marx of egrao 243 A fuco f L 2 [,] may e wre as f ( ) = lm c ( ), (5) = where, c = f, ad, s he sadard er produc o L 2 [,]. If he seres (5) s rucaed a = m, he m f c = m m = C B( ), (6) where, ad = [ c m, cm, L, cm m (7) C ] B )] ( ) = [ m ( ), m ( ), L, m m (. (8) 4. he operaoal marx of egrao. he orhoormal Berse polyomals operaoal marx of egrao of order ( m + ) ( m + ) wll e derved ow. o acheve hs, cosder he followg egral m ( x) dx = ϕ ( ), <, =,, L, m. = m j= c j m j m ( ), = c, c, L, c ] B( ), for m. (9) [ m m m m Usg equaos (8) ad (9), we oa B x) dx = Pm + B( ) (, ()

2432 A. K. Sgh, V. K. Sgh, O. P. Sgh where he operaoal marx P m+ of egrao assocaed wh orhoormal Berse polyomals s gve y m Pm ( c j m ), j= + = () ad c = ϕ,. (2) j m j m For m = 5, he marx P 6 s deoed y P ad s gve as follows: P :=.52778.288948.24533.2663.59329.92228.2563.226.624.242.25.2248.62.246.242527.97222.27389.68.828.9725.69444.26896.46743.6686.34479.4667.8234.77868.52.65296.99.992.4375.92.783.3889 (3) 5. Numercal example he followg example shows he compuaoal power of he Berse polyomal operaoal marx of egrao. Cosder a lear me-varyg sysem a y& ( ) + y( ) = u( ), wh y() =, (4) / a where u() s he u sep fuco. he aalyc soluo of (4) s y( ) = e. Gu ad Jag [5] cosdered hs prolem wh a =. 25 ad gave a approxmae soluo y usg Haar waveles wh four, sx ad e ass fucos. Paraskevopoulos e all. [2] cosdered he same prolem wh a = ad used Fourer seres operaoal marx of egrao of orders ( ) ad ( 2 2) o oa approxmae soluos. I 2, Razzagh ad Yousef [6] used Legedre waveles operaoal marx of order ( 6 6) o solve hs prolem. We oa approxmae soluo of (4) usg he Berse operaoal marx of egrao P m+ y akg m = 4, 5 ad compare he soluos. Iegrag (4) from o, we ge a y( ) + y( x) dx = u( x) dx. (5)

Berse operaoal marx of egrao 2433 Usg (6), he ukow fuco y() ad u sep fuco u () are approxmaed as y( ) = C B( ) ad u( ) = d B( ), (6) where = [ c m, cm, K, cm m s o e deermed. Susug (6) (5), we oa C ] Usg (), we ge a C B( ) + C B( x) dx = d B( x) dx. (7) ac B ) + C P B( ) d P B( ). (8) ( m+ = m+ As (8) holds for all [,), reduces o a C. (9) + C Pm + = d Pm + akg raspose of (9), oe oas ( + a I + Pm + ) C = Pm d. Wrg Q a I + P m + =, ad E = P d, we ge m+ Q C = E, (2) where I s he ( m + ) ( m + ) u marx. Eq. (2) s a se of algerac equaos whose soluo gves c, m. Solvg (2), oe ges m C = Q E. (2) Fally, he soluo y () s oaed y susug (2) o (6).For m = 5, he gve y [.55277,.5,.44959,.372678,.288675,.66667] d s d =. (22) I fgures ad 2, graphs of he exac soluo as well as hose of he approxmae soluo y( ) = C B( ) for m = 4 ad 5 are gve akg a =. 25, respecvely. Fgure 3 depcs he correspodg errors ewee he approxmae soluos. From Fg.

2434 A. K. Sgh, V. K. Sgh, O. P. Sgh 3, ca e see ha he accuracy creases que fas as we go from level m = 4 o m = 5..75 y() Y().5.25.7.33.5.67.83 Fg.. he exac soluo y () Y ( ) (doed le) rucaed a level m = 4. (sold le) ad he approxmae soluo deoed y.75 y() Y ().5.25.7.33.5.67.83 Fg.2. he exac soluo y () (sold le) ad he approxmae soluo deoed y Y () (doed le) rucaed a level m = 5.

Berse operaoal marx of egrao 2435..5 y() Y () y( ) Y().5..7.33.5.67.83 Fg.3. Comparso ewee he errors a m = 5 (sold le) ad m = 4 (doed le). 6. Cocluso. he uform approxmao capales of Berse polyomals coupled wh he fac ha oly a small umer of polyomals are eeded o oa a sasfacory resul makes our mehod very aracve. I gves eer approxmao compared o ha of paraskevopoulos e all. [2], Gu ad Jag [5], ad Razzagh ad Yousef [6]. Refereces []. C. F. Che, C. H. Hsao, Walsh seres aalyss opmal corol, Ieraoal Joural of Corol, 2(975) 88 898. [2]. C. F. Che, C. H. Hsao, A Walsh seres drec mehod for solvg varaoal prolems, Joural of he Frakl Isue, 3 ( 975) 265 28. [3]. C.F. Che, Y.. say,.. Wu, Walsh operaoal marces for fracoal calculus ad her applcao o dsrued- parameer sysems, Joural of he Frakl Isue, 53(977), 267 284.

2436 A. K. Sgh, V. K. Sgh, O. P. Sgh [4]. P. Sau, Aalyss ad syhess of dyamc sysems va lock pulse fucos, Proceedgs of he Isuo of Elecrcal Egeers, 24(977) 569 -- 57. [5]. R-E. Kg, P. N. Paraskevopoulos, Paramerc defcao of dscree me SISO sysems, Ieraoal Joural of Corol, 3(979) 23 29. [6]. C. Hwag, Y. P. Shh, Parameer defcao va Laguerre polyomals, Ieraoal Joural of Sysems Scece, 3(982) 29 27. [7]. C. Hwag, Y. P. Shh, Soluo of egral equaos va Laguerre polyomals, Joural of Compuer ad Elecrcal Egeerg, 9(982) 23 29. [8]. P. N. Paraskevopoulos, Cheyshev seres approach o sysem defcao, aalyss ad corol, Joural of he Frakl Isue, 36(983) 35 57. [9]. I. R. Horg, J. H. Chou, Shfed Cheyshev drec mehod for solvg varaoal prolems, Ieraoal Joural of Sysems Scece, 6(985) 855 86. []. R. Y. Chag, M. L. Wag, Parameer defcao va shfed Legedre polyomals, Ieraoal Joural of Sysems Scece, 3(982) 25 35. []. P. N. Paraskevopoulos, Legedre seres approach o defcao ad aalyss of lear sysems, IEEE rasacos Auomac Corol, 3 (6) ( 985) 585 589. [2]. P. N. Paraskevopoulos, P. D. Spars, S. G. Mourousos, he Fourer seres operaoal marx of egrao, Ieraoal Joural of Sysems Scece, 6(985) 7 76. [3]. M. Razzagh, M. Razzagh, Fourer seres drec mehod for varaoal prolems, Ieraoal Joural of Corol, 48(988) 887 895. [4]. P. N. Paraskevopoulos, P. Sklavouos, G. Ch. Georgou, he operaoal marx of egrao for Bessel fucos, Joural of he Frakl Isue, 327(99) 329 34. [5]. J. S. Gu, W. S. Jag, he Haar waveles operaoal marx of egrao, Ieraoal Joural of Sysems Scece, 27(996) 623 628. [6]. M. Razzagh, S. Yousef, he legedre waveles operaoal marx of egrao, Ieraoal Joural of Sysems Scece, 32 (4) (2) 495 -- 52. Receved: Feruary, 29