Transfer function nyquist diagram problem

6 vues (au cours des 30 derniers jours)
Jacob
Jacob le 19 Juil 2014
Réponse apportée : Craig le 17 Sep 2014
Why does case (1) appear to scale/wrap infinity so that it is visible on the plot, while case (2) does not? Both systems are marginally stable. This is most likely some theory I don't understand but is there a way to wrap the nyquist plot so that I can see what is going on with the rest of the plot for case (2)?
case (1):
numg=[1 5];
deng=conv([1 6 100],[1 4 25]);
G=tf(numg,deng);
Gzpk=zpk(G)
nyquist(G)
case (2):
numg=poly([-3 -5]);
deng=poly([0 -2 -4]);
G=tf(numg,deng);
Gzpk=zpk(G)
nyquist(G)
  1 commentaire
Azzi Abdelmalek
Azzi Abdelmalek le 19 Juil 2014
appear to scale/wrap infinity: What does that mean?

Connectez-vous pour commenter.

Réponses (2)

Star Strider
Star Strider le 19 Juil 2014
Modifié(e) : Star Strider le 19 Juil 2014
I am not certain that I understand your question, but the second one has a pole at infinity while the first on does not. Extending the imaginary axis in the second one would not likely provide any additional useful information.
According to the documentation, you can extend the frequency axis if you want to:
  • nyquist(sys,w) explicitly specifies the frequency range or frequency points to be used for the plot. To focus on a particular frequency interval, set w = {wmin,wmax}. To use particular frequency points, set w to the vector of desired frequencies. Use logspace to generate logarithmically spaced frequency vectors. Frequencies must be in rad/TimeUnit, where TimeUnit is the time units of the input dynamic system, specified in the TimeUnit property of sys.

Craig
Craig le 17 Sep 2014
If I understand what you mean by scale/wrap. The first system has finite magnitude for all values of jw. The second system has an integrator, and therefore the magnitude goes to inf as jw approaches 0. This is why you do not see part of the nyquist plot for the second system at freqeuncies near 0.

Tags

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by