ECE 664 COMPUTABILITY, COMPLEXITY, and LANGUAGES Fall ...

icon

27

pages

icon

English

icon

Documents

Lire un extrait
Lire un extrait

Obtenez un accès à la bibliothèque pour le consulter en ligne En savoir plus

Découvre YouScribe en t'inscrivant gratuitement

Je m'inscris

Découvre YouScribe en t'inscrivant gratuitement

Je m'inscris
icon

27

pages

icon

English

icon

Ebook

Lire un extrait
Lire un extrait

Obtenez un accès à la bibliothèque pour le consulter en ligne En savoir plus

  • cours - matière potentielle : web site
ECE 664 COMPUTABILITY, COMPLEXITY, and LANGUAGES Fall 2011 Instructor: Avi Kak Office: EE 340 Phone: 494-3551 Document Revised on December 20, 2011 1
  • minutes days years centuries
  • greedy algorithms
  • right choices for the various components
  • unit of compu- tation
  • theory of np
  • model primitive
  • problems
  • complexity
Voir icon arrow

Publié par

Nombre de lectures

16

Langue

English

1 Gradient and Hamiltonian systems
1.1t systems
These are quite special systems of ODEs, Hamiltonian ones arising in conserva-
tive classical mechanics, and gradient systems, in some ways related to them,
arise in a number of applications. They are certainly nongeneric, but in view of
their origin, they are common.
A system of the form
0X = r V (X) (1)
n 1whereV :R !R is, say,C , is called, for obvious reasons, a gradient system.
A critical point of V is a point whererV = 0.
These systems have special properties, easy to derive.
Theorem 1. For the system (1), if V is smooth, we have (i) If c is a regular
1point ofV , then the vector eld is perpendicular to the level hypersurface V (c)
1along V (c).
(ii) A point is critical for V i it is critical for (1).
(iii) At any equilibrium, the eigenvalues of the linearized system are real.
More properties, related to stability, will be discussed in that context.
Proof.
(i) It is known that the gradient is orthogonal to level surface.
(ii) This is clear essentially by de nition.
(iii) The linearization matrix elements are a = V (the subscript no-ij x ;xi j
tation of di erentiation is used). Since V is smooth, we have a =a , and allij ji
eigenvalues are real.
1.2 Hamiltonian systems
If F is a conservative eld, then F = r V and the Newtonian equations of
motion (the mass is normalized to one) are
0
q =p (2)
0
p = r V (3)
n nwhere q2R is the position and p2R is the momentum. That is
@H0q = (4)
@p
@H0p = (5)
@q
where
2p
H = +V (q) (6)
2
1is the Hamiltonian. In general, the motion can take place on a manifold, and
then, by coordinate changes,H becomes a more general function ofq andp. The
coordinates q are called generalized positions, and q are the called generalized
momenta; they are canonical coordinates on the phase on the cotangent manifold
of the given manifold.
An equation of the form (4) is called a Hamiltonian system.
0Exercise 1. Show that a system x =F (x) is at the same time a Hamiltonian
system and a gradient system i the Hamiltonian H is a harmonic function.
Proposition 1. (i) The Hamiltonian is a constant of motion, that is, for any
solution X(t) = (p(t);q(t)) we have
H(p(t);q(t)) =const (7)
where the constant depends on the solution.
(ii) The c level surfaces of a smooth function F (p;q) are solutions of
a Hamiltonian system
@F0q = (8)
@p
@F0p = (9)
@x
Proof. (i) We have
dH dp dq
=r H +r = r Hr +r Hr = 0 (10)p q p q q p
dt dt dt
(ii) This is obtained very similarly.
1.2.1 Integrability: a few rst remarks
Hamiltonian systems in one dimension are integrable: the solution can be writ-
ten in closed form, implicitly, asH(y(x);x) =c. Note that for an equation of the
0form y =G(y;x), this is equivalent to the system having a constant of motion.
The latter is de ned as a function K(x;y) de ned globally in the phase space,
(perhaps with the exception of some isolated points where it may have \simple"
singularities, such as poles), and with the property that K(y(x);x) =const for
any given trajectory (the constant can depend on the trajectory, but not on x).
Indeed, in this case we have
d @C @C0C(y(x);x) = y + = 0
dx @y @x
or
@C @C0
y = =
@x @y
which is equivalent to the system
@C @C
x_ = ; y_ = (11)
@y @x
which is a Hamiltonian system.
26
1.2.2 Local versus global
It is important to mention that a system is actually called Hamiltonian if the
functionH is de ned over a su ciently large region, preferably the whole phase
space.
0Indeed, take any smooth rst order ODE, y =f(y;x) and di erentiate with
respect to the initial condition (we know already that the dependence is smooth;
we let dy=d(y ) =y_):0
@f0
y_ = y_ (12)
@y
with the solution Z x @f(y(s);s)
y_ = exp ds (13)
@yx0
and thus, in the local solution y = G(x;x ) we have G (x;x ) = 0, if G is0 x 00
smooth {i.e. the eld is regular{, and the implicit function theorem provides
a local function K so that x =K(y(x);x), that is a constant of motion! The0
big di erence between integrable and nonintegrable systems comes from the
possibility to extend K globally.
1.3 Example
As an example for both systems, we study the following problem: draw the
contour plot (constant level curves) of
2 2 2F (x;y) =y +x (x 1) (14)
and draw the lines of steepest descent of F .
For the rst part we use Proposition 1 above and we write
@F0
x = = 2y (15)
@y
@F0y = = 2x(x 1)(2x 1) (16)
@x
The critical points are (0; 0); (1=2; 0); (1; 0). It is easier to analyze them using
2 2the Hamiltonian. Near (0; 0) H is essentially x +y , that is the origin is a
center, and the trajectories are near-circles. We can also note the symmetry
x! (1 x) so the same conclusion holds for x = 1, and the phase portrait is
symmetric about 1=2.
2 2 2Near x = 1=2 we write x = 1=2 +s, H = y + (1=4 s ) and the leading
2 2Taylor approximation givesHy 1=2s . Then, 1=2 is a saddle (check). Now
we can draw the phase portrait easily, noting that for large x the curves essen-
4 2tially become x +y =C \ attened circles". Clearly, from the interpretation
of the problem and the expression of H we see that all trajectories are closed.
3Figure 1:
Figure 2:
The perpendicular lines solve the equations
@F0x = = 2x(x 1)(2x 1) (17)
@x
@F0y = = 2y (18)4@yWe note that this equation is separated! In any case, the two equation obviously
share the critical points, and the sign diagram can be found immediately from
the rst gure.
Exercise 2. Find the phase portrait for this system, and justify rigorously its
qualitative features. Find the expression of the trajectories of (17). I found

1
y =C 4
2(x 1=2)
2 Flows, revisited
Often in nonlinear systems, equilibria are of higher order (the linearization has
zero eigenvalues). Clearly such points are not hyperbolic and the methods we
have seen so far do not apply.
There are no general methods to deal with all cases, but an important one
is based on Lyapunov (or Lyapounov,...) functions.
De nition. A ow is a smooth map
(X;t)! (X)t
A di erential system
x_ =F (x) (19)
generates a ow
(X;t)!x(t;X)
where x(t;X) is the solution at time t with initial condition X.
The derivative of a function G along a vector eld F is, as usual,
D (G) =rGFF
If we write the di erential equation associated to F , (19), then clearly
d
D G = G(x(t))F jt=0
dt
2.1 Lyapunov stability
Consider the system (19) and assume x = 0 is an equilibrium.
Then
1. x = 0 is Lyapunov stable (or simply stable) if starting with initial condi-e
tions near 0 the ow remains in a neighborhood of zero. More precisely,
the condition is: for every > 0 there is a > 0 so that ifjxj < then0
jx(t)j< for all t> 0.
2. x = 0 is asymptotically stable if furthermore, trajectories that start closee
to the equilibrium converge to the equilibrium. That is, the equilibrium
x is stable if it is Lyapunov stable and if there exists> 0e
so that ifjxj<, then lim x(t) = 0:0 t!1
56
2.2 Lyapunov functions
Let X be a xed point of (19). A Lyapunov function for (19) is a function
de ned in a neighborhood O of X with the following properties
(1) L is di erentiable in O.
(2) L(X ) = 0 (this can be arranged by subtracting a constant).
(3) L(x)> 0 inOnfX g.
(4) D L 0 inO.F
A strict Lyapunov function is a Lyapunov function for which
(4’) D L< 0 inO.F
Finding a Lyapunov function is often nontrivial. In systems coming from
physics, the energy is a good candidate. In general systems, one may try to nd
an exactly integrable equation which is a good approximation for the actual one
in a neighborhood of X and look at the various constants of motion of the
approximation as candidates for Lyapunov functions.
Theorem 2 (Lyapunov stability). Assume X is a xed point for which there
exists a Lyapunov function L. Then
(i) X is stable.
(ii) If L is a strict Lyapunov function then X is asymptotically stable.
Proof. (i) Consider a small ballB3X contained inO; we denote the boundary
ofB (a sphere) by@B. Let be the minimum ofL on the@B. By the denition
of a Lyapunov function, (3), > 0. Consider the following subset:
U =fxB :L(x)<g (20)
From the continuity of L, we see thatU is an open set. Clearly, X U. Let
X 2U. Then x(t;X) is a continuous curve, and it cannot have components
outside B without intersecting @B. But an intersection is impossible since by
monotonicity,L(x(t))L(X)< for allt. Thus, trajectories starting inU are
con ned to U, proving stability.
(ii)
d 1. Note rst that X is the only critical point inO since L(x(t;X )) = 01dt
for any xed point.
2. Note that trajectoriesx(t;X) withX2U are contained in a compact set,
and thus they contain limit points. Any limit point x is strictly insideU
since L(x )<L(x(t);X)<.
3. Letx be a limit point of a trajectoryx(t;X) whereX2U, i.e. x(t ;X)!n
x . Then, by 1 and 2, x 2U and x is a regular point of the eld.
4. We want to show thatx =X . We will do so by contradiction. Assuming
x =X we have L(x ) => 0, again by (3) of the de nition of L.
5. By 3 the trajectoryfx(t;x ) :t 0g is well de ned and is contained in B.
6. We then have L(x(t;

Voir icon more
Alternate Text