[bug #68196] [octave forge] (control) Wrong phase in bode plot

Please use the bug tracker to post updates to a bug report.  The mailing list is intended as a read-only notification stream.  Info posted to this mailing list address won't appear in the tracker database where it is most useful.

Follow-up Comment #8, bug #68196 (group octave):

Too many files, too many variations...
 With a Torsten's "bode.m" modified for `sqrt(eps)` tolerance, attached as
(file bode_1.m)
and with test script (file t1.m) 

$ cat t1.m 
w = 2*pi*[0.001:0.001:50];
Ts = 1/100;
z = tf('z',Ts);
G = Ts/2 * (z+1)/(z-1);
alpha = 2*pi;
Gf = alpha/(z-(1-alpha));
Gnew = Gf*G;
[M,P,W] = bode(Gnew, w);
P(1)*pi/180
H = freqresp(Gnew,w);
angle(H)(1)


I see:

octave:1> pkg load control
octave:2> t1
ans = -1.5708
ans = -1.5708
octave:3> bode(G, w)
octave:4> legend("off")

with plot (file `bode_1.png`) attached.

I think results and the plot are as expected.

Dmitri.
-- 



(file #58417, file #58418, file #58419)

    _______________________________________________________

Additional Item Attachment:

Name: bode_1.m                       Size: 8.0KiB
    <https://file.savannah.gnu.org/file/bode_1.m?file_id=58417>

Name: bode_1.png                     Size: 65KiB
    <https://file.savannah.gnu.org/file/bode_1.png?file_id=58418>

Name: t1.m                           Size: 204B
    <https://file.savannah.gnu.org/file/t1.m?file_id=58419>


    AGPL NOTICE

These attachments are served by Savane. You can download the corresponding
source code of Savane at
https://savannah.gnu.org/source/savane-f290f6b25beb8cb99bbe243a6cd2c5fef79ffcde.tar.gz


    _______________________________________________________

Reply to this item at:

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

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

-----BEGIN PGP SIGNATURE-----

iHUEABYIAB0WIQQk97aszIMMAvLLwm6qLAuaBUf3TgUCacrzJwAKCRCqLAuaBUf3
TvEaAP9V9+oZNyd/cJJHGS15sexY9H8CMxy3h1tvEoceE7xGQwEAnWqkMtXSFrau
iIkMy7P3yK+A4Fu+LNCfWl0jWJbONg4=
=3z+0
-----END PGP SIGNATURE-----