bugGNU Octave - Bugs: bug #49801, [octave forge] (control) dlyap...

 
 

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

Submitted by:  None
Submitted on:  Thu 08 Dec 2016 03:00:58 AM UTC  
 
Category: Octave Forge PackageSeverity: 3 - Normal
Priority: 5 - NormalItem Group: None
Status: NoneAssigned to: None
Originator Name: Hernan J. GonzalezOriginator Email: -unavailable-
Open/Closed: OpenRelease: 4.2.0
Operating System: Any

Add a New Comment(Rich Markup)
   

You are not logged in

Please log in, so followups can be emailed to you.

 

Thu 08 Dec 2016 03:02:22 AM UTC, comment #1:

I forgot to add: tested with version 3.0.0 of control package

Anonymous
Thu 08 Dec 2016 03:00:58 AM UTC, original submission:

dlyap(C,D)

is supposed to solve

CXC' + D = X

(doc: https://octave.sourceforge.io/control/function/dlyap.html)

But it doesn't.

Example:

p=0.3, q=1-p;
C=[ q p 0 ; q 0 p ; q 0 0];
D = [0 0 0; 0 0 0 ; p p p];

X1=dlyap(C,D)
0.034394 0.029959 0.023144
0.029959 0.055090 0.021713
0.023144 0.021713 0.316853

CX1C'+D % this should coincide with X1

0.034394 0.029959 0.023144
0.029959 0.055090 0.021713
0.323144 0.321713 0.316853

See the huge discrepancy in the last row.

In Matlab 2016b it returns the correct value:
X1 =
0.5052 0.4708 0.3606
0.5382 0.5083 0.3833
0.6464 0.6232 0.5475

Anonymous

 

(Note: upload size limit is set to 16384 kB, after insertion of the required escape characters.)

Attach File(s):
   
   
Comment:
   

No files currently attached

 

Depends on the following items: None found

Items that depend on this one: None found

 

Carbon-Copy List
  • -unavailable- added by jwe (Updated the item)
  • -unavailable- added by nrjank (listed control package maintainer.)
  • -unavailable- added by None (Submitted the item)
  •  

    Do you think this task is very important?
    If so, you can click here to add your encouragement to it.
    This task has 0 encouragements so far.

    Only project members can vote.

     

    Please enter the title of George Orwell's famous dystopian book (it's a date):

     

     

    Follow 2 latest changes.

    Date Changed By Updated Field Previous Value => Replaced By
    Sun 13 Aug 2017 02:05:08 PM UTCjweSummarydlyap function (control package) does not work as expected=>[octave forge] (control) dlyap function does not work as expected
    Thu 08 Dec 2016 03:29:21 AM UTCnrjankCarbon-Copy-=>Added -unavailable-

    Back to the top


    Powered by Savane 3.1-cleanup1