Research Papers: Modeling

Fractional Optimal Control of a Distributed System Using Eigenfunctions

[+] Author and Article Information
Om P. Agrawal

Mechanical Engineering and Energy Processes,  Southern Illinois University, Carbondale, IL 62901om@engr.siu.edu

J. Comput. Nonlinear Dynam 3(2), 021204 (Feb 04, 2008) (6 pages) doi:10.1115/1.2833873 History: Received May 23, 2007; Revised September 28, 2007; Published February 04, 2008

This paper presents a formulation and a numerical scheme for fractional optimal control (FOC) for a class of distributed systems. The fractional derivative is defined in the Caputo sense. The performance index of an FOC problem (FOCP) is considered as a function of both the state and the control variables, and the dynamic constraints are expressed by partial fractional differential equations. Eigenfunctions are used to eliminate the space parameter and to define the problem in terms of a set of state and control variables. This leads to a multi-FOCP in which each FOCP could be solved independently. Several other strategies are pointed out to reduce the problem to a finite dimensional space, some of which may not provide a decoupled set of equations. The calculus of variations, the Lagrange multiplier, and the formula for fractional integration by parts are used to obtain Euler–Lagrange equations for the problem. In the proposed technique, the FOC equations are reduced to Volterra-type integral equations. The time domain is discretized into several segments and a time marching scheme is used to obtain the response at discrete time points. For a linear case, the numerical technique results into a set of algebraic equations, which can be solved using a direct or an iterative scheme. The problem is solved for different number of eigenfunctions and time discretizations. Numerical results show that only a few eigenfunctions are sufficient to obtain good results, and the solutions converge as the size of the time step is reduced. The formulation presented is simple and can be extended to FOC of other distributed systems.

Copyright © 2008 by American Society of Mechanical Engineers
Your Session has timed out. Please sign back in to continue.



Grahic Jump Location
Figure 1

Analytical and numerical results for state x0(t) and control u0(t) for α=1.0 (top lines: state x0(t), bottom lines: control u0(t). Δ,O: analytical, ×,+: numerical)

Grahic Jump Location
Figure 2

State function y(z,t) as a function of time (O:z=0.1, +:z=0.2)

Grahic Jump Location
Figure 3

Control function u(z,t) as a function of time (O:z=0.1, +:z=0.2)

Grahic Jump Location
Figure 4

State y(z,t) as a function of time for z=0.1 (X:N=10, O:N=20, Δ:N=40)

Grahic Jump Location
Figure 5

Control u(z,t) as a function of time for z=0.1 (X:N=10, O:N=20, Δ:N=40)



Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In