Extract data from nyquist plot
6 visualizzazioni (ultimi 30 giorni)
Mostra commenti meno recenti
Andreas Bernatzky
il 15 Giu 2020
Commentato: Andreas Bernatzky
il 15 Giu 2020
Hello,
I am trying to extract the data of a nyquist plot and plot it into another figure (I have to do this for another API normally using nyquist would be ok for me).
If i try to plot the data again I get some strange behaviour compared to the original nyquist() plot. Anybody got an idea how to figure this out?
On the left is the nyquist plot created by matlab in
nyquist(G);
and on the right I replot the extracted data
plot(squeeze(re),squeeze(im));
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/316491/image.png)
Used lightweight example:
clear all;
close all;
G = tf([1],[1 0.5]);%observerd transfer funciton
G.OutputDelay = 8;
[MAG,PHASE,W] = bode(G);
figure %do "normal" nyquist plot and extract real and imaginary data
[re,im] = nyquist(G);%extract real and imaginary part
nyquist(G);
figure% do own nyquist plot with extracted data
plot(squeeze(re),squeeze(im));
Thanks for you help!
0 Commenti
Risposta accettata
Ameer Hamza
il 15 Giu 2020
Modificato: Ameer Hamza
il 15 Giu 2020
This code shows how you can extract the data from the graphic handles
G = tf(1,[1 0.5]);%observerd transfer funciton
G.OutputDelay = 8;
fig = figure; %do "normal" nyquist plot and extract real and imaginary data
nyquist(G);
hg = findall(fig, 'Type', 'hggroup');
l = hg.Children;
xdata1 = l(1).XData;
ydata1 = l(1).YData;
xdata2 = l(2).XData;
ydata2 = l(2).YData;
figure% do own nyquist plot with extracted data
plot(xdata1, ydata1, xdata2, ydata2);
3 Commenti
Ameer Hamza
il 15 Giu 2020
In your code, you can get a better resolution if you specify 'w' vector yourself
G = tf([1],[1 0.5]);%observerd transfer funciton
G.OutputDelay = 8;
[MAG,PHASE,W] = bode(G);
figure %do "normal" nyquist plot and extract real and imaginary data
w = 0:0.001:10;
[re,im] = nyquist(G, w);%extract real and imaginary part
nyquist(G);
figure% do own nyquist plot with extracted data
plot(squeeze(re),squeeze(im));
Più risposte (0)
Vedere anche
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!