5.3 Systems, Stability, Fourier, PDEs

Systems of DEs


First solution:

x(1)=eλ⁢t⁢ξλ x raised to the 1 power equals e raised to the lamda t power xi sub lamda

Second solution:

x(2)=t⁢ℇλ⁢t⁢ξ+ℇλ⁢t⁢η x raised to the 2 power equals t ℇ raised to the lamda t power xi plus ℇ raised to the lamda t power eta

where ξ=(A-λ⁢𝕀)⁢η xi equals open paren A minus lamda double struck I close paren times eta. Third solution:

x(3)=t22⁢eλ⁢t⁢ξλ+t⁢ℇλ⁢t⁢η+eλ⁢t⁢ζ x raised to the 3 power equals the fraction with numerator t squared and denominator 2 e raised to the lamda t power xi sub lamda plus t ℇ raised to the lamda t power eta plus e raised to the lamda t power zeta

where η=(A-λ⁢𝕀)⁢ζ eta equals open paren A minus lamda double struck I close paren times zeta.

Critical Points


Jacobian Matrix:

J=[∂x(F⁡(x0,y0))∂x(F⁡(x0,y0))∂x(G⁡(x0,y0))∂y(G⁡(x0,y0))] J equals the 2 by 2 matrix Row 1: Column 1, partial differential sub x open paren F of open paren x sub 0 comma y sub 0 close paren close paren Column 2, partial differential sub x open paren F of open paren x sub 0 comma y sub 0 close paren close paren Row 2: Column 1, partial differential sub x open paren G of open paren x sub 0 comma y sub 0 close paren close paren Column 2, partial differential sub y open paren G of open paren x sub 0 comma y sub 0 close paren close paren

Sometimes a nonlinear ODE has an exact phase portrait given by:

d(y)d(x)=G⁡(x,y)F⁡(x,y)→ d of y over d of x equals the fraction with numerator G of open paren x comma y close paren and denominator F of open paren x comma y close paren right arrow integrate H⁡(x,y)=c H of open paren x comma y close paren equals c

Eigenvalues

Critical Points

Stability

r1>r2>0 r sub 1 is greater than r sub 2 is greater than 0

Node (source)

Unstable

r1<r2<0 r sub 1 is less than r sub 2 is less than 0

Node (sink)

Asymp. Stable

r2<0<r1 r sub 2 is less than 0 is less than r sub 1

Saddle

Unstable

r1=r2>0 r sub 1 equals r sub 2 is greater than 0

Proper/Improper Node

Unstable

r1=r2<0 r sub 1 equals r sub 2 is less than 0

Proper/Improper Node

Asymp. Stable

r1,r2=λ±i⁢μ r sub 1 comma r sub 2 equals lamda plus or minus i mu (λ>0) open paren lamda is greater than 0 close paren

Spiral (Focus)

Unstable

r1,r2=λ±i⁢μ r sub 1 comma r sub 2 equals lamda plus or minus i mu (λ<0) open paren lamda is less than 0 close paren

Spiral (Focus)

Asymp. Stable

r1=i⁢μ, r sub 1 equals i mu comma r2=-i⁢μ r sub 2 equals negative i mu

Centre

Stable

Almost Linear Systems: The Proper/Improper bits become Node/Spiral Point, and the Centre becomes Centre or Spiral (Indeterminate).

Lyapunov Functions


Let V⁡(x,y) V of open paren x comma y close paren be defined on a domain D D containing (0,0) open paren 0 comma 0 close paren.

  • V⁡(x,y) V of open paren x comma y close parenis positive (negative) definite ifV⁡(0,0)=0 V of open paren 0 comma 0 close paren equals 0andV⁡(x,y)>0⁢(∀(x⁢y)∈D)⁢(V⁡(x,y)<0∀(x,y)∈D) V of open paren x comma y close paren is greater than 0 times open paren for all x y is a member of D close paren times open paren V of open paren x comma y close paren is less than 0 for all open paren x comma y close paren is a member of D close paren

  • V⁡(x,y) V of open paren x comma y close parenis positive (negative) semi-definite ifV⁡(0,0)=0 V of open paren 0 comma 0 close paren equals 0andV⁡(x,y)≥0⁢(∀(x⁢y)∈D)⁢(V⁡(x,y)≤0∀(x,y)∈D) V of open paren x comma y close paren is greater than or equal to 0 times open paren for all x y is a member of D close paren times open paren V of open paren x comma y close paren is less than or equal to 0 for all open paren x comma y close paren is a member of D close paren

Theorem: Given an autonomous system with critical point (0,0) open paren 0 comma 0 close paren, if ∃V⁡(x,y) there exists V of open paren x comma y close paren continuous with continuous first partial derivatives, is positive definite, then (0,0) open paren 0 comma 0 close paren is:

  • Asymptotically Stable: If d(V)d(t) d of V over d of t is negative definite on some domain D D containing (0,0) open paren 0 comma 0 close paren.

  • Stable (at non-linear level): If d(V)d(t) d of V over d of t is negative semi-definite.

Theorem: Given an autonomous system with critical point (0,0) open paren 0 comma 0 close paren, assume ∃V⁡(x,y) there exists V of open paren x comma y close paren continuous with continuous first partial derivatives, such that V⁡(0,0)=0 V of open paren 0 comma 0 close paren equals 0 and that in every neighbourhood of (0,0) open paren 0 comma 0 close paren there exists at least one point (x1,y1) open paren x sub 1 comma y sub 1 close paren where V⁡(x1,y1) V of open paren x sub 1 comma y sub 1 close paren is positive negative. If there exists some domain D D containing (0,0) open paren 0 comma 0 close paren where d(V)d(t) d of V over d of t is positive definite negative definite on D D, then (0,0) open paren 0 comma 0 close paren is an unstable critical point.

d(V)d(t)=x′⁢∂x(V)+y′⁢∂y(V) d of V over d of t equals x prime partial differential sub x V plus y prime partial differential sub y V

Fourier Analysis


Fourier Series of ƒ⁡(x) f with hook of x:

ƒ⁡(x)=a02+∑n=1∞(an⁢cos⁡(n⁢π⁢xL)-bn⁢sin⁡(n⁢π⁢xL)) f with hook of x equals the fraction with numerator a sub 0 and denominator 2 plus the sum from n equals 1 to infinity of open paren a sub n the cosine of open paren n pi x over L close paren minus b sub n the sine of open paren n pi x over L close paren close paren

where T=2⁢L T equals 2 L or ƒ⁡(x+2⁢L)=ƒ⁡(x) f with hook of open paren x plus 2 L close paren equals f with hook of x.

a0=1L⁢∫-LL(ƒ⁡(x)⁢d(x)) a sub 0 equals 1 over L the integral from negative L to L of open paren f with hook of x d of x close paren

an=1L⁢∫-LL(ƒ⁡(x)⁢cos⁡(n⁢π⁢xL)⁢d(x)) a sub n equals 1 over L the integral from negative L to L of open paren f with hook of x the cosine of open paren n pi x over L close paren d of x close paren

bn=1L⁢∫-LL(ƒ⁡(x)⁢sin⁡(n⁢π⁢xL)⁢d(x)) b sub n equals 1 over L the integral from negative L to L of open paren f with hook of x the sine of open paren n pi x over L close paren d of x close paren

Function is even if ƒ⁡(-x)=ƒ⁡(x) f with hook of negative x equals f with hook of x only cosine terms.

Function is odd if ƒ⁡(-x)=-ƒ⁡(x) f with hook of negative x equals negative f with hook of x only sine terms.

Partial Differential Equations


We have:

a2⁢∂x2(u)+β⁢u=∂t(u) a squared partial differential sub x squared u plus beta u equals partial differential sub t u

With non-homogeneous boundary conditions u⁡(0,t)=T1 u of open paren 0 comma t close paren equals T sub 1, u⁡(L,t)=T2 u of open paren L comma t close paren equals T sub 2, look for v⁡(x) v of x to solve:

a2⁢v″⁡(x)+β⁢v⁡(x)=0 a squared v double prime of x plus beta v of x equals 0, v⁡(0)=T1, v of 0 equals T sub 1 comma v⁡(L)=T2 v of L equals T sub 2

Then determine solution using:

u⁡(x,t)=v⁡(x)+w⁡(x,t) u of open paren x comma t close paren equals v of x plus w of open paren x comma t close paren

Identities


  • sinh⁡(x)=ex-e-x2 hyperbolic sine of x equals the fraction with numerator e to the x-th power minus e raised to the negative x power and denominator 2

  • cosh⁡(x)=ex+e-x2 hyperbolic cosine of x equals the fraction with numerator e to the x-th power plus e raised to the negative x power and denominator 2

  • cos⁡(x)=ei⁢x+e-i⁢x2 cosine x equals the fraction with numerator e raised to the i x power plus e raised to the negative i x power and denominator 2

  • sin⁡(x)=ei⁢x-e-i⁢x2⁢i sine x equals the fraction with numerator e raised to the i x power minus e raised to the negative i x power and denominator 2 i

  • cos⁡(μ)=0→μ=(2⁢n-1)⁢π2 cosine mu equals 0 right arrow mu equals open paren 2 n minus 1 close paren times pi over 2

  • 1-D Heat Eq: ∂t(u)=a2⁢∂x2(u) partial differential sub t u equals a squared partial differential sub x squared u

  • ∫(sin2⁡(a⁢x)⁢d(x))=x2-sin⁡(2⁢a⁢x)4⁢a the integral of open paren sine squared a x d of x close paren equals x over 2 minus sine 2 a x over 4 a

  • ∫(cos2⁡(a⁢x)⁢d(x))=x2+sin⁡(2⁢a⁢x)4⁢a the integral of open paren cosine squared a x d of x close paren equals x over 2 plus sine 2 a x over 4 a


More Cheat Sheets ->