Help implementing the subspace iteration- or Lanczos method in matlab
I am need to create an eigensolver using either the Subspace iteration method or the lanczoz method in matlab. The problem is related to dynamics in structures and the eigenfrequencies is wanted.
The structure is built with beams and bar elements, and the stiffness matrix "K" and the mass matrix "M" is made. I am having some big problems on getting started with the implementation on the eigensolvers...(Headbang) So help is deeply wanted!
I really appreciate any help!