Skip to content

✍内容

In this section we state the fundamental thm on rectification of a direction field and deduce from it thms on existence, uniqueness, and differentiable dependence of the solution on parameters and initial conditions, thms on extension, and thms on local phase flows.

Rectification of a Direction Field

Consider a smooth direction field in a domain U of n -dimensional space.

[!def] a rectification of a direction field is a diffeomorphism mapping it into a field of parallel directions. A field is said to be rectifiable if there exists a rectification of it.

![[Ordinary Differential Equations - V.I.Arnold(英文版).pdf#page=93&rect=112,126,316,229|Ordinary Differential Equations - V.I.Arnold(英文版), p.89]]

[!thm] (Fundamental). Every smooth direction field is rectifiable in a neighborhood of each point. If the field is r times continuously differentiable (of class Cr,1r), then the rectifying diffeomorphism can also be taken from the class Cr.

Example: The direction field of the equation x˙=x can be rectified by the diffeomorphism (t,x)(t,y=xet). Indeed this diffeomorphism maps the integral curves x=Cet in the (t,x) -plane to the parallel lines y=C in the (t,y) -plane.

Problem 1. Rectify the direction fields of the equations x˙=t and x˙=x2 in a neighborhood of the origin.

The first rectification is (x,t)(x,xt22). The second rectification is not (t,x)(t,1x+t) arround the origin.

NOTE

It is impossible to rectify the direction field of the equation x˙=x2 on the entire plane by a time-preserving diffeomorphism.

Problem 3. Suppose a (smooth) field of two-dimensional planes is given in R3 (a plane is attached at each point). Is it always possible to rectify this field (transform it into a field of parallel planes by a suitable diffeomorphism)? Hint: A rectifiable field is a field of planes tangent to a family of surfaces. Answer. No. Consider, for example, the field of planes given by the equation ydx+dz=0 (a vector belongs to a plane of this field if this 1-form vanishes on it). There is no surface tangent to the planes of the field.

[!thm] All smooth direction fields in domains of the same dimension are locally diffeomorphic (can be mapped into each other by a diffeomorphism).

12: By the fundamental thm all fields are locally diffeomorphic to one standard field. 21: The property of being locally diffeomorphic to any field implies, in particular the property of being locally diffeomorphic to a standard field, i.e., local rectifiability.

[!thm] The differential equation x˙=v(t,x) with smooth right-hand side v is locally equivalent to the very simple equation dy/dτ=0.

NOTE

In other words, in a neighborhood of each point of the extended phase space of (t,x) there exists an admissible coordinate system (τ,y) (transition to which is a diffeomorphic change of variables) in which the equation can be written in the very simple form dy/dτ=0.

13: We first rectify the direcion field v and then consider Cartesian coordinates in which the time axis τ is parallel to the lines of the rectified direction field. 31: Every direction field can be written loaclly as the direction field of a suitable field of a suitable differential equation. Passing to a local coordinate system in which the equation has the form dy/dτ=0 rectifies the given field.

Problem 5. Prove that the coordinate system in Thm 3 can be chosen so that time is left fixed (τt). The equivalence implies τ(t) is a smooth function of t, then dydt=dydτdτdt=0.

Problem 6. Rectify the direction field of the equation x˙=x+t on the entire plane by a time -preservng diffeomorphism (t,x)(t,y(t,x)). Put y(t,x)=(x+t+1)et, just by solving the ode x˙=x+t and let y be the parameter C.

The fundamental theorem on rectification was essentially discovered by Newton. In modern terms Newton's method consists of the following. Suppose given an equationn x˙=v(t,x). We shall seek a rectifying diffeomorphism y=h(t,x) for which y=x when t=0 (time is not changed). From the condition y˙=0 we obtain for h the equation yt+(yx)v0. We expand v nd h in series of powers of t :

h=h0+th1+,v=v0+tv1+

Then h0(x)x, and so hx=E+th1this might be wrong+. We then substitute the series for h and v into the equation for h. We then expand the left-hand side in a series in t. We then set the coefficients of t0,t1, equal to zero in this series (based on the uniqueness of the coefficients of a Taylor series). We obtain successively

h1+v0=0,2h2+h1v0+v1=0

Problem 8. Solve the equation x˙=x by Newton's method with the initial condition φ(0)=1. Solution. φ=1+tφ1+t2φ2+φ1+2φ2+3φ3t2+=1+φ1t+φ2t2+, and consequently φ1=1,φ2=φ1/2,φ3=φ2/3,, whence φk=1/k!. This is the way in which the series for the exponential was originallly discovered.

TIP

The proof of the convergence of the series constructed by Newton was much studied in the 19 th century. The convergence of the series for h in the analytic case was proved by Cauchy. Cauchy's thm was extended to the case of finite smoothness by Picard. The proof is discussed in §32.

Similar to linear algebra

The fundamental Theorem 1 is an assertion of the same character as the theorems of linear algebra on the reduction of quadratic forms or the matrices of linear operators to normal form. It give an exhaustive description of the possible local behaviors of direction fields, reducing all questions to the trivial case of a parallel field. A related theorem of analysis is the implicit function theorem. A smooth mapping f:RmRn is called nondegenerate at the point 0 if the rank of the derivative at this point is as large as possible (i.e., the smaller of m and n). Suppose f(0)=0. Two such mappings f and g are called locally equivalent at the point 0 if one of them maps into the other under the action of diffeomorphisms of the domain and target spaces that leave 0 fixed: h:RmRm,k:RnRn,fh=kg.

RmgRnhkRmfRn

In other words, two mappings are locally equivalent if under suitable choices of admissible local coordinate systems in the domain and target space (with origin at 0) they can be written by the same formulas.

[!thm] Implicit Function Theorem In some neighborhood of a nondegerate point any two smooth mappings (of spaces of fixed dimensions m and n ) are equivalent.

In particular every mapping is equivalent to its linear part at a nondegenerate point. Therefore the theorem just stated is one of numerous theorems on linearization. As a local normal form to which the mapping f reduces by diffeomorphisms h and k it is natural to take the following simplest one:

yi=xi for ir,yi=0 for i>r

where r=min(m,n) is the rank of the derivative of f at zero, xi are the coordinates of a point in the domain space, and yi are coordinates of a point in the target space. In other words, f is an imbedding if the dimension of the domain is less than that of the target space and a fibration otherwise. The reader accustomed to more complocated statements of the implicit unction theorem will easily verify that these more complicated statements are equivalent to the simple geometric statement given here. All the theorems just listed on normal forms describe the orbits of the actions of various groups ("changes of variable") on sets (matrices, forms, fields, and mappings, respectively).

Existence and Uniqueness Theorems

The following corollary is a consequence of the fundamental Theorem 1 on rectification.

[!cor] Through each point of a domain in which a smooth direction field is defined there passes an integral curve.

Consider a diffeomorphism that rectifies the given field. The rectified field consists of parallel directions. In that field an integral curve passes through each point (a straight line, to be specific). The diffeomorphism inverse to the rectifying diffeomorphism maps this line into the desired integral curve.

[!cor] Two integral curves of a smooth direction field having a point in common coincide in a neighborhood of that point.

For a rectified field this is obvious, but a rectifying diffeomorphism maps integral curves of the original field into integral curves of the rectified field.

[!cor] A solution of the differential equation x˙=v(t,x) with the initial condition (t0,x0) in the domain of smoothness of the right-hand side exists and is unique (in the sense that any two solutions with a common initial condition coincide in some neighborhood of the point t0)

We apply the above two corollaries to the direction field of the quation in the extended phase space. The result is Corollary 3.

In Corollary 3 and in what follows x is a point of a phase space of any (finite) dimension m. This corollary is called an existence and uniqueness theorem for solutions of a system of m first-order equations.

Theorems on Continuoes and Differentiable Dependence of the Solutions on the Initial Condition

Consider the value of the solution φ of the differential equation x˙=v(t,x) with initial condition φ(t0)=x0 at the instant of time t as a function Φ of (t0,x0;t) with values in the phase space. The following result is a consequence of the fundamental Theorem 1 on rectification.

[!cor] The solution of an equation with smooth right-hand side depends smoothly on the initial conditions.

NOTE

This means that the function Φ shown above is defined, continuous, and smooth in a neighborhood of each point (t0,x0;t0) (of class Cr if v is of class Cr).

For the simplest equation (v0) this is obvious (Φx0). The general equation reduces to this one by a diffeomorphism (the details are left to the reader).

The theorem on differentiablility with respect to the initial condition provides a quite efficient method of studying the influence exerted on the solution by a small ==perturbation== of the initial condition. If the solution is known for some initial condition, then to determine the deviation of the solution with a nearby initial condition ==from the given "unperturbed" solution== a linear-homogeneous equation is obtained in first approximation (the equation of variations). The "perturbation theory" that arises in this way is but one of the variants of Newton's series method.

Problem 1. Find the derivative of the solution φ of the equation x˙=x2+xsint with respect to the initial condition φ(0)=a for a=0. Solution. By Corollary 4 the solution can be expanded in a Taylor's series in a: φ=φ0+aφ1+ (the dots stand for an infinitesimal of higher order than first in a ). Here φ0 is the unperturbed solution (with zero as initial condition), and φ1 is the unknown derivative. For our equation φ00. Substituting the series into the equation and equating the terms of like degree in a on the left- and right-hand sides (by the uniqueness of Taylor series), we obtain for φ1 the equation of variations φ˙1=φ1sint with initial condition φ1(0)=1 (why?). Hence the answer is e1cost.

x˙=x2+xsintφ0+aφ1+(φ02+φ0sint)+(2φ0φ1+φ1sint)a+0+(φ1sint)a+

于是 φ1=φ1sint,由初值 φ(0)=a,可知 (φ0+φ1a+)(0)=a,于是 φ1(0)=1.

Problem 2. Find the segment of a phase curve of the genrealized Lotka-Volterra system x˙=x(1ya(x,y)), y˙=y(x1) passing through the point x=1,y=ε (with error of order ε2) Solution. The equation of the phase curves is dy/dx=y(x1)/(x(1ya)). The unperturbed solution is y0. The equation of variations is dy/dx=y(x1)/x. The answer is y=εex1/x independently of the form of the function a.

Problem 3. Find the derivative of the solution of the pendulum equation θ¨=sinθ with initial condition θ(0)=a,θ˙(0)=0 with respect to a at a=0. Solution. To apply Corollary 4 one must write the equation as a system of equations. The resulting system of equations of variations can be written in the form of a single second-order equation. It is convenient to write out only the second-order equations equivalent to these systems and their solutions rather than the systems themselves. The unperturbed solution is θ=0. The equation of variations is the equation for small oscillations of the pendulum: θ¨=θ. The answer is cost.

CAUTION

In using the approximate formulas for the perturbed solutions obtained through the equation of variations, one must not forget that they give a good approximation for fixed t and small deviation ε of the initial condition from the unperturbed condition: the error for fixed t is O(ε2), ut nonuniformly as t (the constant in O increases with t )

For example, the formula obtained in Problem 2 would give an incorrect picture of the form of the phase curves of the usual Lotka-Volterra model if we had tried to apply it to describe these curves in the large (as we know from §2 , these curves are closed; the part of the curve far from the x -axis is by no mains described by the solution to Problem 2 ).

In exactly the same way the solution of the full pendulum equation with initial condition (a,0) is near the solution of the equation of small oscillations (with the same initial condition ) for fixed t : the difference between them is of order O(a3) (why?).

Problem 4. Find the first term (linear in a) of the Taylor series expansion of the solution of the equation of the soft pendulum x¨=x3 with initial condition x(0)=0,x˙(0)=a. Solution. The unperturbed solution is x0. The quation of variations is φ¨+1=0. The initial condition is φ1(0)=0,φ˙(0)=1 (why?). The answer is xat.

利用 φ=φ0+φ1a+,φ00 简单验证。

It follows from the differentiability theorem that the error in this approximate formula is at most O(a2) for each fixed t. However for any fixed a0 the approximation becomes completely unsatisfactory for sufficiently large t. This can be seen, for example, from the fact that the approximate solution grows without bound, and the actual solution describes periodic oscillations of an amplitude that is small together with a (the size of the amplitude is of order a by similarity considerations). To estimate the range of applicability of the approximate formula we can compute the following approximations: x=at+a2φ2+a3φ3+. Substituting into the equation, we obtain a2φ¨2+a3φ¨3+=a3t3+. Hence φ20,φ¨3=t3, φ˙3=t4/4,φ3=t5/20, and xata3t5/20+. The second term is small in comparison with the first if a2t4/201, i.e., ta1/2. In other words the value of the approximate solution must be small in comparison with the amplitude of the true oscillation, at a.

Problem 5. Prove that under this condition the relative error of the approximate solution is indeed small. Solution. This follows from similarity conditions. The quasi-homogeneous dilations X=esx,T=est transform the equation x¨=x3 into itself. The solution with initial condition (0,a) becomes the solution with initial condition (0,A=e2sa). The approximate solution x at becomes XAT. We choose s so that A=1. For A=1 the solution XT has small relative error when T1. But dilations do not change the relative errors. Hence the relative error of the approximate solution xat is also small for T1. But T=est,a=e2s. Hence T1 for ta1/2. Thus for small a the approximation gives a small relative error, even for very large t, provided t is small in comparison with the large number 1/a.

In applications of the theory of differential equations it is always necessary to deal with a large number of quantities, some of which are "very small" and others "very large." To discern what is large in comparison with what (i.e., in what order to perform the passages to the limit) is not always easy; the study of this question is sometimes half of the job.

Transformation over the Time Interval from t0 to t

![[Ordinary Differential Equations - V.I.Arnold(英文版).pdf#page=100&rect=43,320,388,538|Ordinary Differential Equations - V.I.Arnold(英文版), p.96]] The following corollary is a consequence of the fundamental theorem on rectification.

[!cor] The transformations over the time interval from t0 to t for an equation with smoot right-hand side

  1. Are defined in a neighborhood of each phase point x0 for t sufficiently close to t0;
  2. Are local diffeomorphisms (of class Cr if the right-hand side is of class Cr) and depend smoothly on t and t0;
  3. Satisfy the identity gt0tx=gstgt0sx for s and t sufficiently closse to t0 (for all x in a sufficiently small neighborhood of the point x0);
  4. Are such that for fixed ξ the function φ(t)=gt0tξ is a solution of the equation x˙=v(t,x) satisfying the initial condition φ(x0)=ξ.

[!thm] The transformation over the time interval from t0 to t for an autonomous equation depends only on the length tt0 of the time interval and not on the initial instant t0.

A translation of the extended phase space of an sutonomous equation along the t -axis maps the direction field onto itself and hence maps integral curves onto one another. Under a translation by s the solution φ with initial condition φ(t0)=x0 maps into the solution ψ with the initial condition ψ(t0+s)=x0. For any t we have ψ(t+s)=φ(t). Consequently gt0tgt0+st+s, as asserted.

Problem 1. Assume that the equation x˙=v(t,x) has T -periodic coefficients (v(t+ T,x)v(t,x)) and that all the mappings over the time interval from t0 to t for it are defined everywhere. Prove that the transformations over the times that are multiples of T form a group: g0kT=Ak for any integer k. Which of the two following relations is true: g0kT+s=Akg0s,g0kT+s=g0sAk ? Answer. The second.

Theorems on Continuous and Differentiable Dependence on a Parameter

Assume that the right-hand side of a given equation x˙=v(t,x;α) depends smoothly on a parameter α ranging over some domain A of the space Ra. The following corollary is a consequence of the fundamental Theorem 1 on rectification.

[!cor] The value of the solution with initial condition φ(t0)=x0 are an instant t depends smoothly on the initial condition, the time, and the parmeter α.

We denote this value by Φ(t0,x0;α,t). The corollary asserts that the function Φ (with values in the phase space) is defined, continuous, and smooth in a neighborhood of each point (t0,x0;α0;t0) of the product of the extended phase space, the time axis, and the domain of variation of the parameter (it is of class Cr if the right-hand side is of class Cr).

A small trick is useful here. Consider the "extended equation" x˙=v(t,x;α),α˙=0 with phase space of dimension m+a (where m=dim{x}). The soluion of this equation with initial condition (t0,x0;α0) is a pair (x=φ(t),α=α0) whose first component φ is the solution of the original equation for α=α0 satisfying the initial condition φ(t0)=x0. By Corollary 4 this pair depends smoothly on (t0,x0;t;α0). Consequently the first component also depends smoothly on these arguments, which was to be proved.

This extension trick reduces the theorem on smooth dependence on the parameter to smooth dependence on the initial conditions. Conversely, given smooth dependence on the parameter (for a fixed initial condition), it is easy to deduce smooth dependence on the initial condition. It suffices to translate the equation so that the initial condition becomes the parameter vα(t,x)=v(t,xα).

Extension Theorems

The following corollary is an obvious consequence of the fundamental theorem on rectification.

[!cor] Corollary 7. A solution with initial condition in a compact set in the extended phase space can be extended forward and backward to the boundary of the compact set.

IMPORTANT

It is essential for the validity of the theorem that the direction field in the extended phase space be "nonvertical."

Proof of Corollary 7.

这就是蔓延的思想.

We begin by proving uniqueness. Consider the least upper bound of the time values at which two solutions with the same initial condition coincide. The solutions coincide to the left of this point. If both are defined at this point, then they coincide there also, since they are continuous. But then they must also coincide to the right of the point (by the local uniqueness theorem). Thus the point in question must be an endpoint of one of the intervals of definition. This proves uniqueness of the forward extension of the solution (for the backward extension the reasoning is similar). We now construct the extension.

The following result is an obvious consequence of Corollary 7 .

[!cor] Corollary 8. A solution with initial condition in a given compact set K in the phase space can be extended forward (resp. Backward) either indefinitely or to the boundary of the compact set K.

[!cor] Corollary 9. A solution of the autonomous equation x˙=v(x) with initial value in any compact set of the phase space can be continued forward (resp. Backward) either infinitely far or to the boundary of the compact set.