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]