online advertising
Loading [MathJax]/jax/output/HTML-CSS/jax.js

Sunday, January 3, 2016

Sturm–Liouville theory - show that the linear operator is self-adjoint.

The goal of this post is to document my derivation that the linear operator for the Sturm Liouville theory is actually self-adjoint

The differential equation is

ddx(p(x)ddxy)+(λw(x)+q(x))y=0

Rearranging, we can write this as an eigenvalue problem as follow:

ddx(pddxy)+qy=λwy

1w(ddx(pddxy)+qy)=λy

Define the operator L=1w(ddx(pddx)+q), it is trivial to verify the operator is linear.

Consider the function space with inner product defined as <a,b>=abwdx, we wanted to verify the operator L is self-adjoint in this space, in other words, <As,t>=<s,At>.

<As,t>=1w(ddx(pddxs)+qs)twdx=(ddx(pddxs)+qs)tdx=ddx(pddxs)tdx+qstdx.

Similarly we can write the other side as

<s,At>=1w(ddx(pddxt)+qt)swdx=(ddx(pddxt)+qt)sdx=ddx(pddxt)sdx+qtsdx.

So all we needed to prove is (ps)tdx=(pt)sdx, we do that using integration by part:

(ps)tdx=td(ps)dx=pstpsdt=pstpstdx=pstptds=pstpts+sd(pt)=pstpts+(pt)sdx.

So we need to show pstpts=0. This should come from the boundary conditions, which is not specified here. The boundary of a Sturm Liouville problem should satisfy this. For example, p(x)=0 for the Legendre's problem at the boundary.

No comments:

Post a Comment