René's URL Explorer Experiment

{"@context":"https://schema.org","@type":"DiscussionForumPosting","headline":"tf2ss() returns unstable system","articleBody":"I have the following transfer function:\r\n\r\n`_sys1.den=[[array([9.43513863e-11, 6.05312352e-08, 7.92752628e-07, 5.23764693e-06,\r\n         1.82502556e-05, 1.24355899e-05, 8.68206174e-06, 2.73818482e-06,\r\n         4.29133144e-07, 3.85554417e-08, 1.62631575e-09, 8.41098151e-12,\r\n         9.85278302e-15, 4.07646645e-18, 5.55496497e-22, 3.06560494e-26,\r\n         5.98908988e-31])]]`\r\n\r\n'_sys1.num=[[array([9.94004350e-13, 2.67602795e-11, 2.31058712e-10, 1.15119493e-09,\r\n         5.04635153e-09, 1.34066064e-08, 2.11938725e-08, 2.39940325e-08,\r\n         2.05897777e-08, 1.17092854e-08, 4.71236875e-09, 1.19497537e-09,\r\n         1.90815347e-10, 1.00655454e-11, 1.47388887e-13, 8.40314881e-16,\r\n         1.67195685e-18])]]'\r\n\r\nasking for poles:\r\n\r\n`_sys1.poles()=array([-6.28318525e+02+0.00000000e+00j, -6.28318531e+00+0.00000000e+00j,\r\n       -3.14159265e+00+5.44139809e+00j, -3.14159265e+00-5.44139809e+00j,\r\n       -1.19380521e-01+5.84842723e-01j, -1.19380521e-01-5.84842723e-01j,\r\n       -1.88495559e-01+0.00000000e+00j, -6.63225116e-02+9.92624461e-02j,\r\n       -6.63225116e-02-9.92624461e-02j, -1.00530965e-01+0.00000000e+00j,\r\n       -4.32233973e-03+0.00000000e+00j, -6.28318476e-04+0.00000000e+00j,\r\n       -6.28318592e-04+0.00000000e+00j, -6.29820167e-05+2.60763321e-07j,\r\n       -6.29820167e-05-2.60763321e-07j, -6.25315230e-05+0.00000000e+00j])`\r\n\r\nAll negative. System is stable. But then when converting _sys1 to state space using tf2ss() returns a system which have positive poles (unstable). I am using Slycot.\r\n\r\nI also repeated the same conversion but using scipy and it gives me a stable state space as a result.\r\n\r\n\r\npython-control version: '0.9.3.post2'\r\n","author":{"url":"https://github.com/pooyasf","@type":"Person","name":"pooyasf"},"datePublished":"2023-10-17T13:41:17.000Z","interactionStatistic":{"@type":"InteractionCounter","interactionType":"https://schema.org/CommentAction","userInteractionCount":3},"url":"https://github.com/935/python-control/issues/935"}