Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

Orthonormalization using Matlab

  1. Oct 29, 2007 #1
    I'm using matlab's [V,D]=eig(A,B) function to find the eigenvectors and eigenvalues given two full matrices of A and B.

    I know the eigenvectors that I get are not orthonormalized, so how do I do this?

    Let's say I'm solving a simple Sturm-Liouville problem like [tex]\phi''(x)}+\lambda\sigma(x)\phi(x) = 0[/tex] where [tex]\sigma(x) = 1 - x^{2}[/tex].

    The general solution that I have by formulae is

    [tex]\phi_{n}(x)\cong\frac{1}{\sigma^{1/4}}sin[\lambda_{n}^{1/2}\int\sigma(s)^{1/2}ds], \lambda_{n}\cong\frac{(n\pi)^{2}}{(\int\sigma(s)^{1/2}ds)^{2})}[/tex]

    When I compare the graph of the eigenfunction from my formula to the numerical eigenfunction I got, they are quite similar except it looks like it is missing some weighting function.
  2. jcsd
Know someone interested in this topic? Share this thread via Reddit, Google+, Twitter, or Facebook

Can you offer guidance or do you also need help?
Draft saved Draft deleted

Similar Threads - Orthonormalization using Matlab Date
How to program this in Matlab Feb 28, 2018
Using Maxima to plot error in Fourier series Nov 5, 2017
Matlab Solving equation with integration using MATLAB Aug 4, 2017
What Software do you use to do your math manipulations? May 21, 2017