[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: Good eigenvalue/eigenvector test problems for eigs?
From: |
David Bateman |
Subject: |
Re: Good eigenvalue/eigenvector test problems for eigs? |
Date: |
Fri, 22 Sep 2006 17:42:13 +0200 |
User-agent: |
Thunderbird 1.5.0.5 (X11/20060817) |
Sebastien Loisel wrote:
>
> ## Symmetric real PD
> A = spdiags([ones(n,1),4*ones(n,1),ones(n,1)],[-2,0,2],n,n);
>
>
> This should not be ill conditioned. A-2I is diagonally dominant and
> symmetric so all its eigenvalues are positive. 6I-A is also diagonally
> dominant and symmetric so its eigenvalues are all positive. As a
> result, the eigenvalues of A are between 2 and 6, this is a pretty
> tight bound. The condition number is uniformly bounded by 3.
>
> Sebastien Loisel
This was my thinking as well, the the problem should be easy. By try the
code I supplied under matlab and see what you get for "det(A - lamba *
eye(n))".
D.
--
David Bateman address@hidden
Motorola Labs - Paris +33 1 69 35 48 04 (Ph)
Parc Les Algorithmes, Commune de St Aubin +33 6 72 01 06 33 (Mob)
91193 Gif-Sur-Yvette FRANCE +33 1 69 35 77 01 (Fax)
The information contained in this communication has been classified as:
[x] General Business Information
[ ] Motorola Internal Use Only
[ ] Motorola Confidential Proprietary
Re: Good eigenvalue/eigenvector test problems for eigs?, Sebastien Loisel, 2006/09/22