Message Boards Message Boards


Open loop transfer function 50/(s^3+2s^2+s+4) ?

Posted 1 year ago
4 Replies
6 Total Likes

When I examine open loop transfer function L(s)= 50/(s^3+2s^2+s+4) using Routh-Hurowitz as well as the Roots[L(s)] function in Mathematica, both show two poles in right half plane. Yet, when I do a NyquistPlot[L(s) ], there is no encirclement of -1. Why is that?

4 Replies

You can explicitly see that the closed-loop system has 2 unstable poles

In[25]:= csys = SystemsModelFeedbackConnect@tf;
TransferFunctionPoles[csys][[1, 1]] // N

Out[26]= {-4.47411, 1.23705 - 3.2464 I, 1.23705 + 3.2464 I}

Or use the Routh–Hurwitz criterion which shows two sign changes in the first column.

In[27]:= Control`StabilityArray[csys][[1, 1]] // TableForm

enter image description here

Oops! Never mind, my mistake. There should be zero encirclement's since Z=N+P.=0+2=2.

Hi Craig, Thanks for the response. Attached is the code.


Perhaps your syntax is incorrect. The documentation is a reliable friend.

ls = 50/(s^3 + 2 s^2 + s + 4)
tf = TransferFunctionModel[{{ls}}, s]

Solve[0 == 1/ls, s]
Reply to this discussion
Community posts can be styled and formatted using the Markdown syntax.
Reply Preview
or Discard

Group Abstract Group Abstract