help-octave
[Top][All Lists]
Advanced

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: eig from Matlab to Octave


From: Riccardo Corradini
Subject: Re: eig from Matlab to Octave
Date: Mon, 17 Sep 2007 16:32:18 +0200 (CEST)

Am I right if  I use
[D, RCOND] =  det (A)
istead for Matlab rcond(A)?
Finally
I used this function I found in the mailing list archive
## 2001-03-07 Ross Lippert
##   * adapted from rlab, using a slightly different similarity transform
## 2001-03-08 Paul Kienzle
##   * cleanup; add Rolf Fabians comment about direct csf computation

function [U, T] = rsf2csf (U, T)

  if (nargin != 2)
    usage ("[U,T] = rsf2csf (U, T)");
  endif

  ## Find complex unitary similarities to zero the subdiagonal elements.
  n = columns (T);
  if (n > 1)
    idx = find (diag (T, -1))';
    for m = idx
      k = m:m+1;
      d = eig (T (k, k));
      cs = [ T(m+1,m+1)-d(1), -T(m+1,m) ];
      cs = cs / norm (cs);
      G = [ conj(cs); cs(2), -cs(1) ];
      T (k, m:n) = G' * T (k, m:n);
      T (1:m+1, k) = T (1:m+1, k) * G;
      U (:, k) = U (:, k) * G;
      T (m+1, m) = 0;
    endfor
  endif

endfunction 

Why is not included in Octave-forge?
Thanks
Riccardo

David Bateman <address@hidden> ha scritto:
Riccardo Corradini wrote:
> Hi all,
> I am trying to use the e4 toolbox for matlab in octave 2.9.13
> but unfortunately
> I found the following instruction
> [W,d] = eig(T,Y);
> not compatible with eig of octave
> Are there any hints concerning this issue?
> Thanks a lot to all
> Respectfully
> Riccardo
>
> ------------------------------------------------------------------------
> ------------------------------------------------------------------------
>
Use the "qz" function

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




L'email della prossima generazione? Puoi averla con la nuova Yahoo! Mail
reply via email to

[Prev in Thread] Current Thread [Next in Thread]