MathGroup Archive 2012

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

Search the Archive

Re: trouble with obtaining eigenvalue of parametric matrix

  • To: mathgroup at smc.vnet.net
  • Subject: [mg128385] [mg128385] Re: trouble with obtaining eigenvalue of parametric matrix
  • From: "Dr. Wolfgang Hintze" <weh at snafu.de>
  • Date: Thu, 11 Oct 2012 23:58:26 -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: <k55o2o$e7h$1@smc.vnet.net>

On 11 Okt., 08:15, Mehrzad Sasanpoor <mz_sp... at yahoo.com> wrote:
> I have a parametric 13*13 matrix as following and I want to find their eigenvalues in term of these parameters.
>
> Eigenvalues[{{ec,
> tc, tc, ,0, 0, 0, 0, 0, 0, 0, 0, 0, 0}, {tc, ec, 0, tc, 0, 0, 0, 0, 0, 0, 0, 0,
> 0}, {tc, 0, ec, 0, tc, 0, 0, 0, 0, 0, 0, 0, 0}, {0, tc, 0, ec, 0, tc, zc, xc, 0,
> 0, 0, 0, 0}, {0, 0, tc, 0, ec, tc, 0, 0, 0, 0, 0, 0, 0}, {0, 0, 0, tc, tc, ec, 0,
> 0, 0, 0, 0, 0, 0} ,{0, 0, 0, zc, 0, 0, ez, 0, 0, 0, 0, 0, 0}, {0, 0, 0, xc, 0, 0,
> 0, ex, 0, xo, xy, xo, xy}, {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, ey, nyx, yy, yx, yy},
> {0, 0, 0, 0, 0, 0, 0, xo, nyx, ex, 0, 0, 0}, {0, 0, 0, 0, 0, 0, 0, xy, yy, 0, ey,
> 0, 0}, {0, 0, 0, 0, 0, 0, 0, xo, nyx, 0, 0, ex, 0}, {0, 0, 0, 0, 0, 0, 0, xy, yy,
> 0, 0, 0, ey}}]
>
> but I face with this error?
>
> Eigenvalues::matsq: Argument \
> {{ec,tc,tc,0,0,0,0,0,0,0,<<3>>},{tc,ec,0,tc,0,0,0,0,0,0,<<3>>},{tc,0,\
> ec,0,tc,0,0,0,0,0,<<3>>},<<5>>,{0,0,0,0,0,0,0,0,0,0,<<5>>},{0,0,0,0,0,\
> 0,0,xo,nyx,ex,<<3>>},<<3>>} at position 1 is not a nonempty square \
> matrix. >>
>
> what should I do?
>
> thanks a lot for your attention

There were simply some typos in your matrix.

I corrected (by brute force) the lines with a length different from 13
with the result:

In[73]:= m = {{ec, tc, tc, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}, {tc, ec, 0,
    tc, 0, 0, 0, 0, 0, 0, 0, 0, 0}, {tc, 0, ec, 0, tc, 0, 0, 0, 0, 0,
    0, 0, 0},
       {0, tc, 0, 0, ec, tc, zc, xc, 0, 0, 0, 0, 0}, {0, 0, tc, 0,
    ec, tc, 0, 0, 0, 0, 0, 0, 0}, {0, 0, 0, tc, tc, ec, 0, 0, 0, 0,
    0, 0, 0},
       {0, 0, 0, zc, 0, 0, ez, 0, 0, 0, 0, 0, 0}, {0, 0, 0, xc, 0, 0,
    0, ex, 0, xo, xy, xo, xy}, {0, 0, 0, 0, 0, 0, 0, 0, ey, nyx, yy,
    yx, yy},
       {0, 0, 0, 0, 0, 0, 0, xo, nyx, ex, 0, 0, 0}, {0, 0, 0, 0, 0,
    0, 0, xy, yy, 0, ey, 0, 0}, {0, 0, 0, 0, 0, 0, 0, xo, nyx, 0, 0,
    ex, 0},
       {0, 0, 0, 0, 0, 0, 0, xy, yy, 0, 0, 0, ey}};

Checkin the lengths using:
In[74]:= (Length[m[[#1]]] & ) /@ Range[13]
Out[74]= {13, 13, 13, 13, 13, 13, 13, 13, 13, 13, 13, 13, 13}

Looking at the appearance of the matrix using
In[75]:= MatrixForm[m]

An finally determining the Eigenvalues (I have supressed the very long
output)
In[76]:= Eigenvalues[m]

Hope this helps.

Regards,
Wolfgang



  • Prev by Date: Re: Integrating over 3D vector
  • Next by Date: Re: How to simplify hypergeometrics
  • Previous by thread: Re: trouble with obtaining eigenvalue of parametric
  • Next by thread: Re: trouble with obtaining eigenvalue of parametric matrix