Article View: sci.math.num-analysis
Article #288Re: Eigenvalue,Eigenvector (part II)
From: Ferhat F. Hatay
Date: Tue, 14 Apr 1992 14:42
Date: Tue, 14 Apr 1992 14:42
73 lines
2802 bytes
2802 bytes
I am having difficulties to reach Eurepan dudes (dudesses!) by mail,
somewhat they keep rejecting my mail. So, please forgive my overuse
of network to communicate the following -deep thoughts- to them.
I don't claim any liability for my suggestions, anyone is welcomed
to find any holes in my reasoning. feel free to poke my failing eyes.
In article <17229@sophia.inria.fr> falzon@sophia.inria.fr writes:
>
> Does anybody have the general forms of eigenvalues and eigenvectors of an NxN
> pseudo-tridiagonal matrix of the form :
>
> For example with N=4;
>
>
> [(1+1/a) -1/a 0 0 ]
> A = [-1/a (1+2/a) -1/a 0 ]
> [ 0 -1/a (1+2/a) -1/a ]
> [ 0 0 -1/a (1+1/a)]
>
> Thanks in advance.
>
>
> Fred.
>-----------------------------------------------------------------------------------------
Hmm, The problem adds to its lure when we consider it with
some boundary conditions. It stops being circulant and
turns to a headache. Assuming that this is a solution
to a second order pde like Helmholtz or Poisson, let me
suggest the following:
For dirichlet boundary conditions we can fold them
onto (or into?) a system of order of N-2. The
resulting -nonsingular- system (matrix) will again
have a circulant structure and a different right
hand side.
For neumann boundary conditions: I use my eigenvalue solver!
I still think there is an analytical solution for the problem
mentioned here. If we solve the eigenproblem for a matrix A:
A = Q D Qt where D is diagonal matrix of eigenvalues.
Q A Qt = D
what would be the solution for
B = A + S ? where S is a diagonal matrix.
let's try:
Q B Qt = Q A Qt + Q S Qt
Q B Qt = D + Q S Qt
then what? eigenvectors of matrix B are Q same as the original
matrix A and eigenvalues are given by D + S?
Do you see any holes in the reasoning above? It was rather
spontaneous so I should have probably checked myself.
cheers,
Ferhat
--
Ferhat F. Hatay || Disclamer: |||||||||||||||||
Aerospace Engineering Sciences || I speak for myself unless ||
University of Colorado, Boulder || otherwise indicated. ||
e-mail: hatay@spot.colorado.edu |||||||||||||||||||||||||||||||
Technical Details
Message-ID:
<1992Apr14.144201.2845@ucsu.Colorado.EDU>Path:
arc1-ca.pugleaf.net!sparky!uunet!boulder!ucsu!spot.Colorado.EDU!hatayReferences:
<17229@sophia.inria.fr>Show Full Headers
Complete Article Headers
Newsgroups: sci.math.num-analysis Path: arc1-ca.pugleaf.net!sparky!uunet!boulder!ucsu!spot.Colorado.EDU!hatay From: hatay@spot.Colorado.EDU (Ferhat F. Hatay) Subject: Re: Eigenvalue,Eigenvector (part II) Message-ID: <1992Apr14.144201.2845@ucsu.Colorado.EDU> Sender: news@ucsu.Colorado.EDU (USENET News System) Nntp-Posting-Host: spot.colorado.edu Organization: University of Colorado, Boulder References: <17229@sophia.inria.fr> Date: Tue, 14 Apr 1992 14:42:01 GMT Lines: 73