Services & Resources / Wolfram Forums / MathGroup Archive
-----

MathGroup Archive 2012

[Date Index] [Thread Index] [Author Index]

Search the Archive

Re: Eigenvalue and eigenvectors of a 10x10 matrix

  • To: mathgroup at smc.vnet.net
  • Subject: [mg128086] Re: Eigenvalue and eigenvectors of a 10x10 matrix
  • From: Ray Koopman <koopman at sfu.ca>
  • Date: Fri, 14 Sep 2012 00:25:15 -0400 (EDT)
  • Delivered-to: l-mathgroup@mail-archive0.wolfram.com
  • Delivered-to: l-mathgroup@wolfram.com
  • Delivered-to: mathgroup-newout@smc.vnet.net
  • Delivered-to: mathgroup-newsend@smc.vnet.net
  • References: <k2s2o4$dpk$1@smc.vnet.net>

On Sep 13, 12:43 am, Redeemed <cakp... at gmail.com> wrote:
> I want to do the eigen analysis of the matrix below
> mat := {{0, 1, 0, 0, 0, 0, 0, 0, 0, 0}, {-(1 + K + K1), -0.1, K, 0, 0,
>      0, 0, 0, K1, 0},
>    {0, 0, 0, 1, 0, 0, 0, 0, 0, 0}, {K, 0, -(1 + 2 K + K1), -0.1, K, 0,
>      0, 0, K1, 0},
>    {0, 0, 0, 0, 0, 1, 0, 0, 0, 0}, {0, 0, K, 0, -(1 + 2 K + K1), -0.1,
>      K, 0, K1, 0},
>    {0, 0, 0, 0, 0, 0, 0, 1, 0, 0}, {0, 0, 0, 0, K,
>     0, -(1 + K + K1), -0.1, K1, 0},
>    {0, 0, 0, 0, 0, 0, 0, 0, 0, 1}, {K1, 0, K1, 0, K1, 0, K1,
>     0, -(1 + 4 K1), -0.1}};
>
> I kept getting a long solution with some Root [] and #1
> I do not know what I am doing wrong
> Any help,
> Its very urgent
>
> Thanks

Change .1 to 1/10

mat = {
{0, 1, 0, 0, 0, 0, 0, 0, 0, 0},
{-(1+K+K1), -1/10, K, 0, 0, 0, 0, 0, K1, 0},
{0, 0, 0, 1, 0, 0, 0, 0, 0, 0},
{K, 0, -(1+2K+K1), -1/10, K, 0, 0, 0, K1, 0},
{0, 0, 0, 0, 0, 1, 0, 0, 0, 0},
{0, 0, K, 0, -(1+2K+K1), -1/10, K, 0, K1, 0},
{0, 0, 0, 0, 0, 0, 0, 1, 0, 0},
{0, 0, 0, 0, K, 0, -(1+K+K1), -1/10, K1, 0},
{0, 0, 0, 0, 0, 0, 0, 0, 0, 1},
{K1, 0, K1, 0, K1, 0, K1, 0, -(1+4K1), -1/10}};

Eigenvalues[mat]

{(1/20)*(-1 + I*Sqrt[399]),
 (1/20)*(-1 - I*Sqrt[399]),
 (1/20)*(-1 - Sqrt[-399 - 2000*K1]),
 (1/20)*(-1 + Sqrt[-399 - 2000*K1]),
 (1/20)*(-1 - Sqrt[-399 - 800*K - 400*K1]),
 (1/20)*(-1 + Sqrt[-399 - 800*K - 400*K1]),
 (1/20)*(-1 - Sqrt[-399 - 800*K - 400*Sqrt[2]*K - 400*K1]),
 (1/20)*(-1 + Sqrt[-399 - 800*K - 400*Sqrt[2]*K - 400*K1]),
 (1/20)*(-1 - Sqrt[-399 - 800*K + 400*Sqrt[2]*K - 400*K1]),
 (1/20)*(-1 + Sqrt[-399 - 800*K + 400*Sqrt[2]*K - 400*K1])}

Eigensystem[mat]  returns similar output, but too much to give here.



  • Prev by Date: Re: Begginer's question: Saving a function for later use
  • Next by Date: Re: slide show questions
  • Previous by thread: Re: Eigenvalue and eigenvectors of a 10x10 matrix
  • Next by thread: FindMinimum function does not seem to work for simple problems