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

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]

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.

Reply to this discussion
Community posts can be styled and formatted using the Markdown syntax.
Reply Preview
or Discard

Group Abstract Group Abstract