A New Approach to the Numerical Solution of Fractional Order Optimal Control Problems

Similar documents
Some algorthim for solving system of linear volterra integral equation of second kind by using MATLAB 7 ALAN JALAL ABD ALKADER

Solvability of nonlinear Klein-Gordon equation by Laplace Decomposition Method

Spectral Galerkin Method for Optimal Control Problems Governed by Integral and Integro- Differential Equations

Contraction Mapping Principle Approach to Differential Equations

Technical Vibration - text 2 - forced vibration, rotational vibration

LECTURE 5. is defined by the position vectors r, 1. and. The displacement vector (from P 1 to P 2 ) is defined through r and 1.

Investigations of Electromagnetic Space-Charge Effects in Beam Physics

Science Advertisement Intergovernmental Panel on Climate Change: The Physical Science Basis 2/3/2007 Physics 253

Solutions for Nonlinear Partial Differential Equations By Tan-Cot Method

e t dt e t dt = lim e t dt T (1 e T ) = 1

Stress Analysis of Infinite Plate with Elliptical Hole

f(x) dx with An integral having either an infinite limit of integration or an unbounded integrand is called improper. Here are two examples dx x x 2

EFFECT OF TEMPERATURE ON NON-LINEAR DYNAMICAL PROPERTY OF STUFFER BOX CRIMPING AND BUBBLE ELECTROSPINNING

An Automatic Door Sensor Using Image Processing

() t. () t r () t or v. ( t) () () ( ) = ( ) or ( ) () () () t or dv () () Section 10.4 Motion in Space: Velocity and Acceleration

Theory of! Partial Differential Equations!

September 20 Homework Solutions

EXISTENCE AND UNIQUENESS OF SOLUTIONS FOR A SECOND-ORDER ITERATIVE BOUNDARY-VALUE PROBLEM

ÖRNEK 1: THE LINEAR IMPULSE-MOMENTUM RELATION Calculate the linear momentum of a particle of mass m=10 kg which has a. kg m s

P441 Analytical Mechanics - I. Coupled Oscillators. c Alex R. Dzierba

USING LOWER CLASS WEIGHTS TO CORRECT AND CHECK THE NONLINEARITY OF BALANCES

Theory of! Partial Differential Equations-I!

Solutions to Problems from Chapter 2

Fourier-Bessel Expansions with Arbitrary Radial Boundaries

Fractional operators with exponential kernels and a Lyapunov type inequality

Green s Functions and Comparison Theorems for Differential Equations on Measure Chains

(b) 10 yr. (b) 13 m. 1.6 m s, m s m s (c) 13.1 s. 32. (a) 20.0 s (b) No, the minimum distance to stop = 1.00 km. 1.

5.1-The Initial-Value Problems For Ordinary Differential Equations

Perturbation Analysis of the LMI-Based Continuous-time Linear Quadratic Regulator Problem for Descriptor Systems

ANALYSIS OF KINEMATICS AND KINETOSTATICS OF FOUR-BAR LINKAGE MECHANISM BASED ON GIVEN PROGRAM

ME 141. Engineering Mechanics

ON THE STABILITY OF DELAY POPULATION DYNAMICS RELATED WITH ALLEE EFFECTS. O. A. Gumus and H. Kose

Particle Horizons in Cosmology

). So the estimators mainly considered here are linear

Previously. Extensions to backstepping controller designs. Tracking using backstepping Suppose we consider the general system

Computer Aided Geometric Design

On the Pseudo-Spectral Method of Solving Linear Ordinary Differential Equations

Solution of fuzzy multi-objective nonlinear programming problem using interval arithmetic based alpha-cut

Convergence of Singular Integral Operators in Weighted Lebesgue Spaces

On Control Problem Described by Infinite System of First-Order Differential Equations

D zone schemes

On Some Hadamard-Type Inequalıtıes for Convex Functıons

Available online Journal of Scientific and Engineering Research, 2018, 5(10): Research Article

LAPLACE TRANSFORM OVERCOMING PRINCIPLE DRAWBACKS IN APPLICATION OF THE VARIATIONAL ITERATION METHOD TO FRACTIONAL HEAT EQUATIONS

Lecture 10. Solution of Nonlinear Equations - II

..,..,.,

CALDERON S REPRODUCING FORMULA FOR DUNKL CONVOLUTION

Forms of Energy. Mass = Energy. Page 1. SPH4U: Introduction to Work. Work & Energy. Particle Physics:

Reinforcement learning

Lecture-V Stochastic Processes and the Basic Term-Structure Equation 1 Stochastic Processes Any variable whose value changes over time in an uncertain

Application on Inner Product Space with. Fixed Point Theorem in Probabilistic

4.8 Improper Integrals

«ACTA»2002. Introduction. Theory. Dissipative Systems. I. D. Chueshov Introduction to the Theory of Infinite-Dimensional Dissipative Systems

7 Wave Equation in Higher Dimensions

Compressive modulus of adhesive bonded rubber block

Numerical solution of fuzzy differential equations by Milne s predictor-corrector method and the dependency problem

On the Semi-Discrete Davey-Stewartson System with Self-Consistent Sources

Boundary layer problem for system of first order of ordinary differential equations with linear non-local boundary conditions

The Finite Element Method for the Analysis of Non-Linear and Dynamic Systems

Average & instantaneous velocity and acceleration Motion with constant acceleration

Hermite-Hadamard-Fejér type inequalities for convex functions via fractional integrals

EE757 Numerical Techniques in Electromagnetics Lecture 9

Improved Approximate Solutions for Nonlinear Evolutions Equations in Mathematical Physics Using the Reduced Differential Transform Method

Procedia Computer Science

Optimization and Application of initial Value of Non-equidistant New Information GM(1,1) Model

Ans: In the rectangular loop with the assigned direction for i2: di L dt , (1) where (2) a) At t = 0, i1(t) = I1U(t) is applied and (1) becomes

Analytic solution of linear fractional differential equation with Jumarie derivative in term of Mittag-Leffler function

Physics 2A HW #3 Solutions

THE DESIGN METHOD OF ROBUST CONTROL BY FLEXIBLE SPACECRAFT. V.Yu. Rutkovsky, S.D. Zemlyakov, V.M. Sukhanov, V.M. Glumov

The Hermite-Hadamard's inequality for some convex functions via fractional integrals and related results

Mathematics 805 Final Examination Answers

0 for t < 0 1 for t > 0

Some Inequalities variations on a common theme Lecture I, UL 2007

v T Pressure Extra Molecular Stresses Constitutive equations for Stress v t Observation: the stress tensor is symmetric

Modeling of the thermodilution curve for coronary flow assessment. Carel Stasse TU/e BMTE07.22

Think of the Relationship Between Time and Space Again

1.0 Electrical Systems

Narrow-band Receiver Radio Architectures

_ J.. C C A 551NED. - n R ' ' t i :. t ; . b c c : : I I .., I AS IEC. r '2 5? 9

Physics 201, Lecture 5

FRACTIONAL EULER-LAGRANGE EQUATION OF CALDIROLA-KANAI OSCILLATOR

Extension of Hardy Inequality on Weighted Sequence Spaces

SOME USEFUL MATHEMATICS

NEW FRACTIONAL DERIVATIVES WITH NON-LOCAL AND NON-SINGULAR KERNEL Theory and Application to Heat Transfer Model

On The Hermite- Hadamard-Fejér Type Integral Inequality for Convex Function

Development of a New Scheme for the Solution of Initial Value Problems in Ordinary Differential Equations

Fractional Order Thermoelastic Deflection in a Thin Circular Plate

Computer Propagation Analysis Tools

On Hadamard and Fejér-Hadamard inequalities for Caputo k-fractional derivatives

On The Estimation of Two Missing Values in Randomized Complete Block Designs

Picture Array Generation Based on Membrane Systems and 2D Context-Free Grammars

Lyapunov-type inequalities for Laplacian systems and applications to boundary value problems

A new approach to Kudryashov s method for solving some nonlinear physical models

ENGR 1990 Engineering Mathematics The Integral of a Function as a Function

Research Article Wirtinger-Type Inequality and the Stability Analysis of Delayed Lur e System

ItsApplication To Derivative Schrödinger Equation

FRACTIONAL-order differential equations (FDEs) are

Chapter 2. First Order Scalar Equations

[5] Solving Multiple Linear Equations A system of m linear equations and n unknown variables:

Generalisation on the Zeros of a Family of Complex Polynomials

Transcription:

Avilble hp://pvu.edu/ Appl. Appl. Mh. ISSN: 93-9466 Vol. 8 Issue Decebe 3 pp. 53 534 Applicions nd Applied Mheics: An Inenionl ounl AAM A Ne Appoch o he Nueicl Soluion o Fcionl Ode Opil Conol Pobles T. Abin nd M. Keynpou Depen o Applied Mheics Fculy o Mheicl Sciences Univesiy o Guiln Rsh In Abin.@gil.co;.eynpou@gil.co Received: uly 4 ; Acceped: Sepebe 9 3 Absc In his icle ne nueicl ehod is poposed o solving clss o cionl ode opil conol pobles. The cionl deivive is consideed in he Cpuo sense. This ppoch is bsed on cobinion o he peubion hooopy nd peeizion ehods. The conol uncion u is ppoied by polynoil uncions ih unnon coeiciens. This ehod conves he cionl ode opil conol poble o n opiizion poble. Nueicl esuls e included o deonse he vlidiy nd pplicbiliy o he ehod. Keyods: Fcionl ode opil conol Hooopy peubion ehod Cpuo cionl deivive. AMS-MSC No.: 49 6A33. Inoducion Fcionl opil conol pobles FOCPs e opil conol pobles ssocied ih cionl dynic syses. The cionl opil conol heoy is vey ne opic in heics. FOCPs y be deined in es o dieen ypes o cionl deivives. Bu he os ipon ypes o cionl deivives e he Rienn-Liouville nd he Cpuo cionl deivives. In Agl 4 Agl nd Blenu 7 he uhos obined 53

54 T. Abin nd M. Keynpou necessy condiions o FOCPs ih he Rienn-Liouville deivive nd ee ble o solve he poble nueiclly. Agl 8 pesened qudic nueicl schee o clss o cionl opil conol pobles FOCPs. In Agl 8c he FOCPs e ouled o clss o disibued syses hee he cionl deivive is deined in he Cpuo sense nd nueicl echnique o FOCPs pesened. Blenu e l. 9 used diec nueicl schee o ind soluion o he FOCPs. In Biss nd Sen FOCPs ih ied inl ie e consideed nd nsvesely condiion is obined. FOCPs ih dynic consin involving inege nd cionl deivives e lso consideed Biss nd Sen b. Bsed on he epnsion oul o cionl deivives ne soluion schee s poposed in elicic nd Peovci 9. Loi e l. used Legende ohonol polynoil bsis o solve he FOCPs. A diec ehod using Eigen uncions o solve he FOCPs o -diensionl syse s pesened in zdei e l. 9 hee he G nld-leniov ppoiion s used o ppoie he cionl deivives. Siil eps hve been de by sevel eseches o solving he FOCP o disibued syses Hsn e l. Rpic nd elicic. Ticud nd Chen pesened nueicl schee o FOCPs bsed on inege ode opil conols poble. In Youse e l. he usge o Legende ulivele bsis nd collocion ehod s poposed o obining he ppoie soluion o FOCPs. Ticud nd Chen b poposed ionl ppoiion bsed on he Hnel d i o he ipulse esponse o obin soluion o he genel ie-opil poble. The ineesed ede is eeed o Evigen nd zdei Wng nd Zhou d zdei 9 Agl 8b nd Fedeico e l. 8 o uhe inoion. The hooopy peubion ehod HPM s pplied o solve ODE nd PDE equions in He 3 He 5 Biz nd Ghzvini 9. The nueicl ehod o FOCPs pesened in his ppe ollos he ppoch pesened in Keynpou nd Azizse Keynpou nd Abin Bozbdi e l. Bozbdi nd Azizse. O couse in his ppe e develop hybid o peizion nd odiied hooopy peubion ehod o solve FOCPs. This ppe is ognized s ollos: In Secion e pesen soe bsic deiniions. In secion 3 e descibe ou ehod. In secion 4 e epo ou nueicl esuls. Finlly secion 5 e conclude he ppe.. Fcionl Opil Conol Poble Seen.. Bsic Deiniions We pesen soe bsic deiniions eled o cionl deivives. The Le Rienn-Liouville deivive o cionl ode o uncion is deined by:

AAM: Inen.. Vol. 8 Issue Decebe 3 55 d d d D hile he Righ Rienn-Liouville cionl deivive is given s:. d d d D b b Anohe cionl deivive he le Cpuo cionl deivive is deined s: * d d d D hile he igh Cpuo cionl deivive given by: * d d d D b b hee... Fcionl Opil Conol Poble Foulion The FOCPs in he sense o Cpuo e ouled s ollos: b d u u Miniize subjec o: * u G D he iniil condiions o he pobles e = -. s b c hee 3 4 5 6 7

56 T. Abin nd M. Keynpou G g g l u u C l b R R l u in hich nd u e he se nd conol vibles especively b snds o he ie nd nd G e given nonline uncions. Hee e ssue h FOCPs hve unique soluion. The bsic eisence nd uniqueness ollo o he Lipschiz condiion by using concion pping heoe nd eighed no ih Mig-Lele in Y e l. Lshinh nd Mohp Podlubny 999 So e l. 993 Sh nd Y 999. 3. Descipion o Mehod In his secion he poposed ehod is descibed nd n ssocied lgoih is pesened. The coninuous conol uncion u is ppoied ih inie cobinion o eleens o bsis Rudin 976 s ollo: u j j q j. 8 Since he FOCPs e solved by hooopy peubion ehod e consuc conve hooopy s ollos: d d d p d G u D* p 9 nd suppose he soluion o Equion6 hs he olloing o: p p hee j j e unnon uncions. Subsiuing Equion ino Equion 9 o nd equing he coeiciens o he es ih idenicl poes o p e deive: d p : c d d d p : G u D* d d d d p : G u D* d d.

AAM: Inen.. Vol. 8 Issue Decebe 3 57 As p Equion 9 ends o Equion 6 nd Equion in os cses conveges o n ppoie soluion o Equion 6 i.e.. By subsiuion o Equion 8 nd Equion ino Equion 5 nd Equion 7 e obin n ppoie soluion o FOCPs s ollos: Le in s. : b b c s. j q d * be he opil vlue o Equion 3. A sopping cieion is poposed s ollos: j j 3 * * 4 hee he sll posiive nube is chosen ccoding o he ccucy desied. We popose he olloing lgoih hich is pesened in o sges. Algoih: Iniilizion sep: Choose o he ccucy desied nd se nd go o he in sep. Min sep Sep. Se u by Equion 8 nd go o Sep. Sep. Copue by Equion nd go o Sep 3. Sep 3. Then copue Sep 4. * in by Equion 3. I go o sep 5. Oheise go o Sep 4. I he sopping cieion Equion 4 holds hen sop; else go o Sep 5. Sep 5. nd go sep. 3. Nueicl Resuls In his secion e pply he ehod pesened in Secion 3 o solve he olloing o es eples. All copuions cied ou by he pcge MAPLE 3.

58 T. Abin nd M. Keynpou Eple. Conside he olloing ie invin poble u d 5 subjec o: D* u ih he iniil condiion. The ec soluion o is hee sinh cosh cosh sinh u cosh sinh sinh cosh.98 * ih objecive vlue u.999. In Figue he se vible nd he conol vible u e ploed o. I is obvious h by pplying he lgoih pesened in secion 3 he ppoie vlues o nd u convege o he ec soluions. Figue shos he se nd he conol inpu u s uncions o ie o dieen vlues o. Choosing 5 he esuls o he pplying he given lgoih e pesened in Tble. Tble. Nueicl esuls in Eple n * 4 7.999.99 6 6.953.9 7 6.795.8 6 5.679

AAM: Inen.. Vol. 8 Issue Decebe 3 59 Figue. Appoie soluions nd ec soluion o u o n 7 b Appoie soluions nd ec soluion o o n 7 o Eple b Figue.: Se s uncion o o dieen vlues o b Conol u s uncion o o dieen vlues o o Eple Eple. In his eple ie vying FOCP is consideed o ind he conol u hich iniizes he peonce inde

53 T. Abin nd M. Keynpou u d 6 subjec o: D* u ih ee einl condiion nd he iniil condiion. Figue 3 deonses he ppoiion o nd u o dieen vlues o. The esuls o pplying he lgoih e pesened in Tble. Figue 4 shos he se nd he conol vibles especively s uncion o ie o. 8 o dieen vlues o n. I is obvious h he ppoie vlues nd u convege o he ec soluions by incesing he vlues o n. Tble. Nueicl esuls in Eple n * 3 5.4847.99 3 5.48347.9 3 5.4765.8 3 5.467 b Figue 3. Se s uncion o o dieen vlues o u s uncion o o dieen vlues o o Eple b Conol

AAM: Inen.. Vol. 8 Issue Decebe 3 53 b Figue : Convegence o he se vible o he ie-vying syse o. 8 b Convegence o he conol vible o he ie- vying syse o.8. Tes pobles nd ee solved in Agl 8 in dieen y. Ou esuls shon in Figues -4 e in good geeen ih he esuls deonsed in Agl 8. Bu e chieved siscoy nueicl esuls in only 5 ieions hile in Agl 8 he nube o ppoiions ss in nd inceses up o 3. So i is signiicn h e chieved ou nueicl esuls ih vey sll ode o ppoiions. Also e ind he ppoie opil vlue o he objecive uncion o ech. 4. Conclusion In his ppe e hve developed he hooopy peubion nd peeizion ehods o solving clss o cionl opil conol pobles. By he poposed ehod e e ble o educe he in poble o n opiizion poble. The nueicl esuls hve deonsed he high ccucy o he poposed ehod. Acnoledgens We e vey geul o o nonyous eeees o hei ceul eding nd vluble coens hich led o he ipoveen o his ppe.

53 T. Abin nd M. Keynpou REFERENCES Agl O Psh 8. A qudic nueicl schee o cionl opil conol pobles. ASME ounl o Dynic Syses Mesueen nd Conol Vol. 3 No. ---6. Agl O Psh 8b. A oulion nd nueicl schee o cionl opil conol pobles.. Vib. Conol. Vol. 4 9-99. Agl O Psh 4. A genel oulion nd soluion schee o cionl opil conol pobles Nonline Dynics Vol. 38 No. 4 33-337. Agl O Psh nd Blenu Duiu 7. A Hilonin oulion nd diec nueicl schee o cionl opil conol pobles ounl o Vibion nd Conol Vol. 3 No. 9-69-8. Agl O Psh 8c. Fcionl opil conol o disibued syse using eigenuncions ASME ounl o Copuionl nd Nonline Dynics Vol. 3 No. 4--4-6. Bozbdi Ab Hshei Azizse Mojb nd Fd Oid Solini. Hooopy peubion ehod o opil conol pobles govened by Vole inegl equions ARDCS Vol. No. 3-5. Blenu Duiu Deeli Ozle nd Agl O Psh 9. A Cenl Dieence Nueicl Schee o Fcionl Opil Conol Pobles ounl o Vibion nd Conol Vol. 5 No. 4 583-597. Biz nd Ghzvini Hosein 9. Convegence o he hooopy peubion ehod o pil dieenil equions Nonline AnlysisRel Wod Applicions Vol. 633--64. Bozbdi Ab Hshei nd Azizse Mojb. Nueicl soluion o opil conol pobles govened by inego-dieenil equions vi hybid ieive schee WAS Vol. IssN 88-495 538-543. Biss Rj Ku nd Sen Siddhh. Fcionl opil conol pobles: pseudo-se-spce ppoch. ounl o Vibion nd Conol Vol. 7 No.7 34 4. Biss Rj Ku nd Sen Siddhh b. Fcionl opil conol pobles ih speciied inl ie. ASME ounl o Copuionl nd Nonline Dynics Vol. 6 9. 9.6. Evigen Fi nd zdei Neci. A Fcionl Ode Dynicl Tjecoy Appoch o Opiizion Poble ih HPM. Spinge Eds. Blenu D. Mchdo.A.T. Luo A.C.. Fcionl Dynics nd Conol ISBN 978--464-456-9 45-55. Fedeico Gso S. F. nd Toes Deli F. M. 8. Fcionl consevion ls in opil conol heoy Nonline Dynics Vol. 53 No.3 5-. He i-hun 3. Hooopy peubion ehod A ne nonline nlyicl echnique Appl.Mhe. Copu Vol. 35 73-79. He i-hun 5. Applicion o hooopy peubion ehod o nonline ve equions Chos Solion. Fc Vol. 6 695-7. Hsn M. Mehedi Tngpong Xingqing W nd Agl O Psh. Fcionl opil conol o disibued syses in spheicl nd cylindicl coodines.. Vib. Conol doi:.77/7754634847. d Fhd Abdeljd Thbe nd Blenu Duiu. Fcionl viionl opil conol pobles ih delyed guens. Nonline Dynics Vol. 6 No. 3 69 64.

AAM: Inen.. Vol. 8 Issue Decebe 3 533 elicic Zon. D nd Peovci Nebojs 9. Opiliy condiions nd soluion schee o cionl opil conol pobles Sucul nd Mulidiscipliny Opiizion Vol. 38 No. 6 57-58. Keynpou Mohd nd Azizse Mojb. Nueicl soluion o opil conol pobles by n ieive schee AMO- Advnced Modeling nd Opiizion Vol. 3 No. 5-37. Keynpou Mohd nd Abin Theeh. Opil Conol o Se-Dely Syses vi Hybid o Peubion nd Peizion Mehod ou. Advnc. Resech. Dynic. Con. Sys Vol. 3 No. 45-58. Loi Ai Dehghnb Mehdi nd Yousei Sohb Ali. A nueicl echnique o solving cionl opil conol poblescopues nd Mheics ih Applicionsicle in pess doi:.6/j.c..3.44. Lshinh V. nd Mohp RN.. Sic sbiliy o dieen il equions Nonline Anl Vol. 46 No. 7 95 9. zdei Neci Agl O Psh sende Beyz Billu nd Kdeniz Dey 9. Fcionl opil conol o -diensionl disibued syse using eigenuncions Nonline Dyn Vol. 55 No. 3 5 6. zdei Neci Agl O Psh Kdeniz Dey nd sende Beyz Billu 9. Fcionl opil conol o n is-syeic diusion-ve popgion. Physic Scip T36 Vol. 34 44-49. zdei Neci Kdeniz Dey nd sende Beyz Billu 9. Fcionl opil conol poble o disibued syse in cylindicl coodines Physics Lees A Vol. 373 Issue: -6. Podlubny Igo 999. Fcionl dieenil equions. Acdeic Ne Yo. Rudin Wle 976. Pinciples o heicl nlysis 3d edn McG-Hill. Rpic Miln R nd elicic Zon D.. Opil Conol o Clss o Fcionl He Diusion Nonline Dyn Vol. 6-39-5. So Sen Kilbs Anolii Alesndovich Michev nd Oleg Igoevic 993. Fcionl inegls nd deivives: Theoy nd pplicions Godon nd Bech 6 p ISBN 884864. Sh Michel D. nd Y Coşun 999. Genelized viion o pees ih iniil ie dieence nd copison esul in e Lypunov-lie uncions In Non-line Di Equions Theoy Mehods Appl Vol. 5 86 8. Ticud Chisophe nd Chen Qun Yng. An ppoie ehod o nueiclly solving cionl ode opil conol pobles o genel o Copues nd Mheics ih Applicions Vol. 59 No. 5 644 655. Tngpong Xingqing W. nd Agl O Psh 9. Fcionl opil conol o coninuu syses. ASME ounl o Vibion nd Acousics Vol. 3 No.... Ticud Chisophe nd Chen Qun Yng b. Tie-Opil Conol o Syses ih Fcionl Dynics Inenionl ounl o Dieenil Equions Aicle ID 4648 6 pges. Yousei Sohb Ali Loi Ai nd Dehghnb Mehdi. The use o Legende ulivele collocion ehod o solving he cionl opil conol pobles ounl o Vibion nd Conol Vol. 7 No. 3 59 65.

534 T. Abin nd M. Keynpou Y Cosun Gucen Mus By nd Cice Muhed. Sic Sbiliy o Fcionl Peubed Syses in Tes o To Mesues Spinge Eds. Blenu D.Mchdo.A.T. Luo A.C.. Fcionl Dynics nd Conol ISBN 978--464-456- 9 59-57. Wng inrong nd Zhou Yong. A clss o cionl evoluion equions nd opil conols Nonline Anlysis: Rel Wold Applicions Vol. 6 7.