How to extract phase information from freqs?

2 views (last 30 days)
I'm plotting using this:
freqs(butter_tf_num,butter_tf_den, scale)
But I want to adjust the X axis to another unit and set limits. So I do:
[H,W] = freqs(butter_tf_num,butter_tf_den, scale);
Freq_Hz = W(:)/2/pi;
subplot(3,2,1)
% Magnitude
plot(Freq_Hz, abs(H))
subplot(3,2,2)
% Phase
plot(Freq_Hz, atand(imag(H)./real(H)))
axis([600 1000 -200 200])
But the result is different... Why?

Accepted Answer

Rick Rosson
Rick Rosson on 3 Mar 2016
Edited: Rick Rosson on 4 Mar 2016
Instead of
atand(imag(H)./real(H))
please try using
angle(H)*180/pi

More Answers (0)

Tags

Community Treasure Hunt

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

Start Hunting!

Translated by