Sturm-Liouville Orthogonality of Eigenfunctions

Amcote
Messages
16
Reaction score
0

Homework Statement


Consider the following Sturm-Liouville Problem:
\dfrac{d^2y(x)}{dx^2} + {\lambda}y(x)=0, \ (a{\geq}x{\leq}b)
with boundary conditions
a_1y(a)+a_2y{\prime}(a)=0, \ b_1y(b)+b_2y{\prime}(b)=0
and distinguish three cases:
a_1=b_1, a_2{\neq}0, b_2{\neq}0a_2=b_2=0, a_1{\neq}0, b_1{\neq}0a_21{\neq}0, a_2{\neq}0, b_1{\neq}0, b_2{\neq}0

a) Without determining the actual eigenvalues and eigenfunctions, prove, in detail, the orthogonality of non-degenerate (i.e. having different eigenvalues) eigenfunctions in all three cases.

Homework Equations


Suppose we have a solution with eigen functions y_j,y_m such that their eigen values are \lambda_j,\lambda_m, \lambda_m\neq{\lambda}_j
Then they are orthogonal if:
\int_{a}^{b}dxy_j(x)y_m(x)=0,

The Attempt at a Solution


So I am okay with doing this integral and proving the orthogonality of the two eigen functions. But what i do not understand is what a1, a2, b1 and b2 are and how I implement the special cases.

Basically I solve the integral and I am left with

[y_{j}p(x)y_m{\prime}(x)-y_m(x)p(x)y_j{\prime}(x)]^{b}_{a}-\int^{b}_{a}[y_{j}{\prime}(x)p(x)y_{m}{\prime}(x)-y_{m}{\prime}(x)p(x)y_{j}{\prime}(x)]
where clearly the second part is zero and the first part should be zero with boundary conditions. But like I said I have no idea what the a1, a2, b1, b2 mean.

Any help would be appreciated
 
Physics news on Phys.org
Amcote said:
I have no idea what the a1, a2, b1, b2 mean.
They just serve to identify three alternative sets of boundary conditions, which will give different solutions to the DE. Taking it just one step further, the alternative sets of boundary conds are:

(1) ##
a_1y(a)+a_2y{\prime}(a)=0, \ a_1y(b)+b_2y{\prime}(b)=0;
##
(2) ##
a_1y{\prime}(a)=0, \ b_1y{\prime}(b)=0
##
(3) ##
a_1y(a)+a_2y{\prime}(a)=0, \ b_1y(b)+b_2y{\prime}(b)=0
##
and none of the ##a_k,b_k## coefficients may be zero except ##a_1##.

You said the first term in your result 'should be zero with boundary conditions'. Have you done the calcs to test that?
 
A couple things I want to add
General Sturm Liouville equation:
\frac{d}{dx}p(x)\frac{dy}{dx}-s(x)y(x)+{\lambda}w(x)y(x)=0
from my understanding, in our problem p(x) = 1 and s(x)=0 and w(x)=1.
and so the last bit of that integral for us is:
[y_j(x)y{\prime}_m(x)-y_m(x)y{\prime}_j(x)]_a^b

And the reason I know this should be zero is because that entire integral should be 0 due to orthogonality. The problem I am having is, I guess, testing this using those boundary conditions. I don't know if I'm just missing something stupid but the constants (a1, b1, a2, b2) are completely throwing me off with this question.

Also I made a mistake writing the first case and third case, this is how it should be:
(1)a_1=b_1=0,a_2\neq0,b_2\neq0
(3)a_1\neq0,a_2\neq0,b_1\neq0,b_2\neq0
 
Amcote said:
I don't know if I'm just missing something stupid but the constants (a1, b1, a2, b2) are completely throwing me off with this question.
You're probably overthinking it. For example, if ##a_1 = 0## and ##a_2\ne 0##, you have ##a_2 y'(a) = 0##. Since ##a_2 \ne 0##, you can divide by it, so ##y'(a) = 0##.
 
vela said:
You're probably overthinking it. For example, if ##a_1 = 0## and ##a_2\ne 0##, you have ##a_2 y'(a) = 0##. Since ##a_2 \ne 0##, you can divide by it, so ##y'(a) = 0##.
Yeah that's what I ended up doing. And for case 3 I rearranged the boundary conditions so that y=\frac{a_2}{a_1}y\prime and similarly with the b2, b1 and was able to get it to zero.

Thanks for your help guys.
 
There are two things I don't understand about this problem. First, when finding the nth root of a number, there should in theory be n solutions. However, the formula produces n+1 roots. Here is how. The first root is simply ##\left(r\right)^{\left(\frac{1}{n}\right)}##. Then you multiply this first root by n additional expressions given by the formula, as you go through k=0,1,...n-1. So you end up with n+1 roots, which cannot be correct. Let me illustrate what I mean. For this...
Back
Top