WEBVTT
Kind: captions
Language: en
00:00:02.210 --> 00:00:08.120
Now welcome back to one dimensional flows,
where we are still dealing with flows on the
00:00:08.120 --> 00:00:20.599
line. The focus of this lecture will be on
fixed points and stability. Now the geometric
00:00:20.599 --> 00:00:29.009
way of thinking can readably be applied to
any one dimensional system x dot is equal
00:00:29.009 --> 00:00:37.519
to f of x. Now essentially you just need to
draw the graph of f of x and sketch the vector
00:00:37.519 --> 00:00:41.489
field on the real line.
00:00:41.489 --> 00:00:51.089
So, just as an example, x dot verses x and
you plot any arbitrary function f of x may
00:00:51.089 --> 00:01:00.469
be, highlight the fixed points of the system
and draw the direction of the flows. Now what
00:01:00.469 --> 00:01:11.369
we do is you place an imaginary particle which
is called a phase point at x0 and see how
00:01:11.369 --> 00:01:22.049
it is actually carried by the flow. Now, as
time moves on the phase point moves along
00:01:22.049 --> 00:01:29.469
the x axis according to some function x(t).
00:01:29.469 --> 00:01:40.060
Now this function is called the trajectory
based at x0 and represents the solution of
00:01:40.060 --> 00:01:50.140
the ordinary differential equations starting
from the initial condition x0. The above picture
00:01:50.140 --> 00:02:01.890
which shows all the qualitatively different
trajectories of the system is called a phase
00:02:01.890 --> 00:02:05.820
portrait.
00:02:05.820 --> 00:02:18.470
The fixed points x* defined by f of x* is
equal to zero correspond to stagnation points
00:02:18.470 --> 00:02:31.590
of the flow. The solid dot is a stable fixed
point and the local flow is towards the stable
00:02:31.590 --> 00:02:46.450
fixed point. The open dot is an unstable fixed
point and the local flow is away from it.
00:02:46.450 --> 00:02:56.310
In terms of the ordinary differential equation
fixed points represents equilibrium solutions,
00:02:56.310 --> 00:03:08.150
if initially x is equal to x* then x(t) is
equal to x* for all t and equilibrium is stable
00:03:08.150 --> 00:03:19.580
if all sufficiently small disturbances away
from it actually dumped out in time. Unstable
00:03:19.580 --> 00:03:33.870
equilibrium in which disturbances grow in
time are represented by unstable fixed points.
00:03:33.870 --> 00:03:38.570
Now let us consider some examples, we talked
about some theory, let us look at some examples
00:03:38.570 --> 00:03:48.940
now. So find the fixed points for x dot is
equal to x square -1 and classify the stability
00:03:48.940 --> 00:03:56.120
of those fixed points. So x dot is equal to
f of x, which is equal to x square -1, which
00:03:56.120 --> 00:04:08.710
is equal to x -1 times x + 1. So to find all
the fixed points, we set f of x* is equal
00:04:08.710 --> 00:04:13.200
to zero and solve for x*.
00:04:13.200 --> 00:04:24.449
This gives x* is equal to + 1 and -1 the algebra
is straight forward. Now we need to determine
00:04:24.449 --> 00:04:35.430
if these fixed points are actually stable
or unstable. So we have to plot x square -1
00:04:35.430 --> 00:04:43.009
and sketch the resulting vector field. So
we go ahead and plot x dot versus x and make
00:04:43.009 --> 00:04:52.610
a plot of f of x is equal to x square -1,
highlight the direction of the flows and highlight
00:04:52.610 --> 00:05:03.199
the fixed points. The flow is to the right
where x square -1 is greater than zero and
00:05:03.199 --> 00:05:06.930
to the left where x square -1 is less than
zero.
00:05:06.930 --> 00:05:19.590
Thus x* is equal to -1 is a stable fixed point
and x* is equal to one is an unstable fixed
00:05:19.590 --> 00:05:28.889
point. Now we go ahead and highlight the stable
fixed point and the unstable fixed point in
00:05:28.889 --> 00:05:36.200
the diagram. Now here is a question for you
can you actually find an explicit analytical
00:05:36.200 --> 00:05:42.979
solution to the above ordinary differential
equation? Now note that we have we are able
00:05:42.979 --> 00:05:49.690
classify stability of fixed points without
actually knowing the analytical solution.
00:05:49.690 --> 00:05:58.849
But do go ahead find, try and see an explicit
analytic solution to the above ODEs.
00:05:58.849 --> 00:06:09.599
Now here are some notes the definition of
stable equilibrium is actually based on the
00:06:09.599 --> 00:06:21.100
notion of small disturbances. So small disturbances
to x* is equal to -1, will actually decay
00:06:21.100 --> 00:06:35.520
a large disturbance, where x finds itself
to the right of x = 1 will actually not decay.
00:06:35.520 --> 00:06:44.110
In fact the phase point will get repelled
towards plus infinity to get this point across
00:06:44.110 --> 00:06:56.030
x* is equal to -1 is locally stable, but not
globally stable. In the sense, that we will
00:06:56.030 --> 00:07:01.909
not have convergence from any arbitrary initial
conditions.
00:07:01.909 --> 00:07:12.620
Now let us consider another example, consider
the equation x dot is equal to x - cos x,
00:07:12.620 --> 00:07:22.770
sketch the phase portrait and determine the
stability of the fixed points. Now let us
00:07:22.770 --> 00:07:37.099
plot y = x and y = cos x and that is the curve
for y = x and that is simple minded curve
00:07:37.099 --> 00:07:51.719
for y = cos x. Now note that they intersect
only at one point. The intersection is a fixed
00:07:51.719 --> 00:08:02.120
point since x* is equal to cos x* and therefore
f of x* is equal to zero.
00:08:02.120 --> 00:08:09.430
When the line is above the cosecant curve,
when x is greater than cos x and x dot is
00:08:09.430 --> 00:08:20.930
greater than zero and so the flow is to the
right when the line is below the cosine curve
00:08:20.930 --> 00:08:28.750
and x dot is less than zero and the flow is
to the left. Now, we go ahead and actually
00:08:28.750 --> 00:08:49.830
identify the fixed point x*, now x* is the
only fixed point and it is unstable. We do
00:08:49.830 --> 00:09:01.330
not have a formula for x*, but we were still
able to classify its stability.
00:09:01.330 --> 00:09:09.140
Now let us look at our final example and consider
an electrical circuit, so let us go ahead
00:09:09.140 --> 00:09:24.260
and actually make a diagram of our electrical
circuit. We have resistor R and capacitor
00:09:24.260 --> 00:09:36.920
C which are in series with a battery of constant
dc voltage V not. Now suppose that switch
00:09:36.920 --> 00:09:50.070
is closed at t is equal to zero and that there
is no charge initially on the capacitor. Let
00:09:50.070 --> 00:09:59.470
Q of t denote the charge on the capacitor
at the t greater than or equal to zero.
00:09:59.470 --> 00:10:08.720
Now let us try and understand some aspects
of this particular system. Now, this type
00:10:08.720 --> 00:10:17.810
of circuit is actually governed by linear
differential equations and can actually be
00:10:17.810 --> 00:10:28.640
solved analytically. However, we work with
a geometric approach to try and get some intuition
00:10:28.640 --> 00:10:34.390
about the about the above circuit.
00:10:34.390 --> 00:10:44.061
So now, if you go around the circuit, the
total voltage drop must equal zero. Which
00:10:44.061 --> 00:11:02.420
gives us -V not +RI + Q/C = 0. Where, I is
the current flowing through the resistor.
00:11:02.420 --> 00:11:12.330
The current causes the charge to accumulate
on the capacitor at a rate Q dot is equal
00:11:12.330 --> 00:11:25.670
to I, this gives us -V not + RQ dot + Q/C
= 0. Rewriting this we get Q dot is equal
00:11:25.670 --> 00:11:39.370
to f is the function Q = Vnot / R - Q / RC.
Now the graph of the f of Q is just a straight
00:11:39.370 --> 00:11:43.110
line with negative slope.
00:11:43.110 --> 00:11:57.310
Now that is the plot of Q dot versus Q and
we plot f of Q on this plot the vector field
00:11:57.310 --> 00:12:14.920
as a fixed point where f of Q = 0 at Q star
= CVnot. So highlight the fixed point the
00:12:14.920 --> 00:12:22.300
flow is to the right where f of Q is greater
than zero and the flow is to the left where
00:12:22.300 --> 00:12:37.880
f of Q is less than zero. So the flow is towards
Q star and hence we have stable fixed point.
00:12:37.880 --> 00:12:48.690
Now let us plot the Q of t versus t, start
a phase point at the origin of the figure
00:12:48.690 --> 00:13:00.970
and let us imagine how it would move the flow
would carry the phase point monotonically
00:13:00.970 --> 00:13:12.430
towards the fixed point Q star. The speed
Q dot decreases linearly as it approaches
00:13:12.430 --> 00:13:22.610
Q star. So Q of t is increasing and concave
down. Now let us go ahead and plot Q of t
00:13:22.610 --> 00:13:36.100
versus t and highlight the equilibrium value
Q star and denote the stable fixed point and
00:13:36.100 --> 00:13:46.150
that is the curve for Q of t versus t. Now
in fact the fixed point is approached from
00:13:46.150 --> 00:13:53.740
all initial conditions and so in that sense
it is actually globally stable.
00:13:53.740 --> 00:14:01.500
Ok so let us have some final comments and
remarks about this lecture. Now this lecture
00:14:01.500 --> 00:14:08.180
was centred around fixed points and stability
and we started dealing with stability in little
00:14:08.180 --> 00:14:16.410
bit more detail because stability is key concept
in study of control and nonlinear dynamics.
00:14:16.410 --> 00:14:22.810
Now we have these special points referred
to as fixed points and we were able to classify
00:14:22.810 --> 00:14:28.570
whether fixed points was stable or they were
unstable.
00:14:28.570 --> 00:14:36.670
And we were able to do this without any formal
algebraic technique. But there was one point
00:14:36.670 --> 00:14:42.800
of subtlety that arouse and that was of local
versus global stability. Now let just take
00:14:42.800 --> 00:14:49.720
an example, for example we got this pen of
mine, and let assume that the tip of the pen
00:14:49.720 --> 00:14:58.690
is a stable fixed point. So if I found myself
to be a silent neighbourhood of this stable
00:14:58.690 --> 00:15:04.360
fixed point. Then there is decent chance that
I would actually get attracted towards that
00:15:04.360 --> 00:15:06.030
stable fixed point.
00:15:06.030 --> 00:15:12.940
On the other hand, if there was a large disturbance
and I found myself rather far away from this
00:15:12.940 --> 00:15:17.450
particular stable fixed point, then it is
not very clear that the trajectories would
00:15:17.450 --> 00:15:22.700
actually converge towards to the stable fixed
point. In fact one of our examples highlighted
00:15:22.700 --> 00:15:29.130
that, you could actually diverge away from
that stable fixed point, they could actually
00:15:29.130 --> 00:15:35.350
have large disturbance. So there is one subtle
point which came out in the lecture and that
00:15:35.350 --> 00:15:41.380
was of local versus global stability and that
should be one of the key take away from this
00:15:41.380 --> 00:15:41.710
lecture.