[Octave-bug-tracker] [bug #49801] [octave forge] (control) dlyap function does not work as expected

[Torsten Lilge]

Follow-up Comment #2, bug #49801 (project octave):

I think, the matrix D in your discrete-time Lyapunov equation must be a
Hermitian matrix, i.e., a complex square matrix that is equal to its own
conjugate transpose. This is not the case and explains why the entries (3,1)
and (3,2) exactly differ by the related entries in D which would be 0 for a
Hermitian matrix.

When changing D to

D = [0 0 p;
     0 0 p;
     p p p];

the result is as expected.

    _______________________________________________________

Reply to this item at:

  <https://savannah.gnu.org/bugs/?49801>

_______________________________________________
  Message sent via Savannah
  https://savannah.gnu.org/