Eigendecomposition problem?

• To: mathgroup at smc.vnet.net
• Subject: [mg126921] Eigendecomposition problem?
• From: Guido <gbellomo at gmail.com>
• Date: Mon, 18 Jun 2012 05:42:02 -0400 (EDT)
• Delivered-to: l-mathgroup@mail-archive0.wolfram.com

```Hi! I'm having some trouble with this. In the code I define two
matrices, then calculate its eigenvalues/eigenvectors, and finally I
plot the absolute values of the components of the eigenvectors. The
results are consistent with (I'm quite sure) the analytical ones but
only below "t approx 30". Above this value, (only) the first plot
corresponding to one of the vectors shows some strange behaviour.

Any idea? Thanks!

This is the code:

Clear[a, b, w, g, t]
\[Rho]t = {{a, b*Exp[-I*w*t - g*t]}, {b\[Conjugate]*Exp[I*w*t - g*t],
1 - a}};
{e1, e2} = Eigenvalues[\[Rho]t];
{v1, v2} = Eigenvectors[\[Rho]t];

a = 0.6`1000; b = 1`1000 + 2`1000*I; w = 1`1000; g = 1`1000;
v1n = Normalize[N[v1, 1000]]; v2n = Normalize[N[v2, 1000]];
Abs1 = {Abs[v1n[[1]]], Abs[v1n[[2]]]};
Abs2 = {Abs[v2n[[1]]], Abs[v2n[[2]]]};

Plot[Abs1, {t, 0, 50}, Filling -> Bottom]
Plot[Abs2, {t, 0, 50}, Filling -> Bottom]

```

• Prev by Date: Re: Export points 3D plots
• Next by Date: Re: Add scale bar to Image
• Previous by thread: Re: NIntegrate through discontinuities
• Next by thread: Bug in Sort !?!?