How to find the max Phase lead provided by a Compensator?
12 views (last 30 days)
Show older comments
If I have a compensator such as s+1/s+3, is there a command to find the max phase lead angle provided by it?
0 Comments
Answers (2)
David Goodmanson
on 1 Dec 2023
Edited: David Goodmanson
on 1 Dec 2023
Hi Gidel,
It's not a single command and it's numeric, but is
w = 0:.001:10; maxangle = (180/pi)*max(angle((i*w+1)./(i*w+3))) % degrees
maxangle = 30.0000
sufficient?
0 Comments
Sam Chak
on 1 Dec 2023
Hi @Gidel
sys = tf([1 1], [1 3])
h = bodeplot(sys);
setoptions(h, 'FreqUnits', 'Hz', 'MagVisible', 'off');
grid on
[mag, phase, wout] = bode(sys);
maxAngle = max(phase)
2 Comments
Sam Chak
on 2 Dec 2023
Hi @Gidel
Since your original question has been answered, perhaps you would also like to determine the compensator's magnitude at the frequency where the maximum phase angle occurs.
% Transfer function of the Lead Compensator
sys = tf([1 1], [1 3])
h = bodeplot(sys); grid on
% Extract data points from Bode plot
[mag, phase, wout] = bode(sys);
% Max phase angle
maxAngle = max(phase)
% Frequency at which the max phase angle occurs, ω_max
idx = find(phase == maxAngle);
wmax = wout(idx)
% Compensator's magnitude at ω_max
magwmax = mag(idx)
% Check mag dB
magdB = 20*log10(magwmax)
See Also
Categories
Find more on Get Started with Control System Toolbox in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!