A set of two 2nd order linear ode

Click For Summary
SUMMARY

This discussion addresses the solution of a coupled set of second-order linear ordinary differential equations (ODEs) with constant coefficients, specifically the equations involving variables y1 and y2. The equations are given as: \(\frac{d^2 y_1}{d x^2} + a^2y_1 = -c y_2\) and \(\frac{d^2 y_2}{d x^2} - b^2y_2 = -d y_1\). The user seeks to find solutions under the constraint that the first derivatives at x=0 are zero. The discussion highlights the method of rewriting the system in vector form and finding eigenvalues, leading to the general solution involving exponential functions based on the eigenvalues derived from the matrix equation \((A-s^{2}I)\vec{v}=0\).

PREREQUISITES
  • Understanding of second-order linear ordinary differential equations (ODEs)
  • Familiarity with eigenvalues and eigenvectors in linear algebra
  • Knowledge of matrix operations and determinants
  • Basic concepts of boundary value problems in differential equations
NEXT STEPS
  • Study the method of solving second-order linear ODEs with constant coefficients
  • Learn about eigenvalue problems and their applications in differential equations
  • Explore the implications of complex eigenvalues on the solutions of ODEs
  • Research boundary conditions and their effects on the uniqueness of solutions
USEFUL FOR

Students and researchers in mathematics or engineering, particularly those focusing on differential equations, linear algebra, and boundary value problems.

omyojj
Messages
32
Reaction score
0
I'm doing an undergrad research job..
I have encountered the following coupled 2nd order linear ODE with constant coefficients a, b, c, d...
[tex]\begin{align} \frac{d^2 y_1}{d x^2} + a^2y_1 & = -c y_2 \\<br /> \frac{d^2 y_2}{d x^2} - b^2y_2 & = -d y_1 \end{align}[/tex]

In addition, I would like to impose constraints
[tex]\frac{d y_1}{d x}|_{x=0} = \frac{d y_2}{d x}_{x=0} = 0[/tex]
(y_1, y_2 are even in x)

if c were zero..the solution would be simple..
[tex]\begin{align} y_1 & = A \cos(ax) \\<br /> y_2 & = B \cosh(bx) + C \cos(ax) \end{align}[/tex]
but i cannot proceed further when c is not zero..
did I fall into the silly trap??
please help me find the solutions..
 
Physics news on Phys.org
Re-write this vector form:

[tex] \vec{y}''=A\vec{y}[/tex]

with y=[y1 y2]
and A=[-a^2 -c; -d b^2 ]

Like in an ordinary constant coefficient ODE, we suggest solutions of the form
[tex] \vec{y}=\vec{v}e^{sx}[/tex]

Substituting into the equation gives:

[tex] s^{2} \vec{v}e^{sx}=A\vec{v}e^{sx}[/tex]

or

[tex] (A-s^{2}I)\vec{v}=0[/tex]

For this algebric equation to have a non-trivial solution, the determinant of the matrix must equal to zero. You might recognize this as searching for eigenvalues, only here the eigenvalues are [tex]\lambda=s^{2}[/tex].

In the most simple case, you'll get two real positive eigenvalues with one eigenvector which,
for which you will have four basic solutions:

[tex]\vec{y}=K_{11}\vec{v_{1}}e^{\sqrt{\lambda_{1}}x}+K_{12}\vec{v_{1}}e^{-\sqrt{\lambda_{1}}x}+K_{21}\vec{v_{2}}e^{\sqrt{\lambda_{2}}x}+K_{22}\vec{v_{2}}e^{-\sqrt{\lambda_{2}}x}[/tex]

Off course other cases, will need more algebric manipulation, to pull out a real solution (f.e. if you get complex eigenvalues, your solutions will be periodic)
 

Similar threads

  • · Replies 5 ·
Replies
5
Views
2K
  • · Replies 7 ·
Replies
7
Views
2K
  • · Replies 2 ·
Replies
2
Views
4K
  • · Replies 3 ·
Replies
3
Views
2K
  • · Replies 2 ·
Replies
2
Views
3K
  • · Replies 5 ·
Replies
5
Views
2K
  • · Replies 1 ·
Replies
1
Views
3K
  • · Replies 2 ·
Replies
2
Views
2K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 5 ·
Replies
5
Views
2K