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=ps′t−∫ps′dt=ps′t−∫ps′t′dx=ps′t−∫pt′ds=ps′t−pt′s+∫sd(pt′)=ps′t−pt′s+∫(pt′)′sdx.
So we need to show ps′t−pt′s=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.
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=ps′t−∫ps′dt=ps′t−∫ps′t′dx=ps′t−∫pt′ds=ps′t−pt′s+∫sd(pt′)=ps′t−pt′s+∫(pt′)′sdx.
So we need to show ps′t−pt′s=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