I have the following transfer function : G = 1/(2*(s*1E-2)*(1+s*1E-2));
I am looking to draw the phase and magnitude plot for this function. I used the following script :
s = tf('s');
G = 1/(2*(s*1E-2)*(1+s*1E-2));
bode(G)
[mag,phase,wout] = bode(G);
Is my aproach correct and if it is, can someone explain to me what are the poles and the zeros for the function ?
My guess is that I have 2 poles : one in 2*E2 and one in 1*E2. Is this correct ?

 Réponse acceptée

Star Strider
Star Strider le 16 Mar 2020

0 votes

The pzplot function:
s = tf('s');
G = 1/(2*(s*1E-2)*(1+s*1E-2));
bode(G)
[mag,phase,wout] = bode(G);
figure
pzplot(G)
shows one pole at -100 and another pole at 0.
Note that there is an error in the coding for ‘G’ that I corrected. There is either a misplaced parenthesis or a missing multiplication operator. The poles remain unchanged regardless how I permute those corrections.

8 commentaires

Grigorie Vlad
Grigorie Vlad le 16 Mar 2020
Ok, but can you help me understand how the 2 constant is influencing the magnitude plot ?
Star Strider
Star Strider le 16 Mar 2020
It is (obviously) dividing the magnitude by 2, however it is not affecting the pole locations because it multiplies the enitre denominator.
Grigorie Vlad
Grigorie Vlad le 16 Mar 2020
Modifié(e) : Grigorie Vlad le 16 Mar 2020
Sorry, but I still not undertand. I ploted the function with and without the 2. All I see that changed is that the magniture plot started from 40 dB(without the 2 constant) and from 34 dB(with the 2 constant). I am having problems understanding how that came to be.
Very simply:
Calc = 20*log10(1/2)
produces:
Calc =
-6.0206
Grigorie Vlad
Grigorie Vlad le 16 Mar 2020
Modifié(e) : Grigorie Vlad le 16 Mar 2020
Ah I think I understand the magnitude plot. One more thing : the phase plot does not change at all, is that correct ?
Star Strider
Star Strider le 16 Mar 2020
That is correct. The phase plot should not change.
Grigorie Vlad
Grigorie Vlad le 16 Mar 2020
Ok, now it is all clear. Thank you very much for the answer and for your time !
Star Strider
Star Strider le 16 Mar 2020
As always, my pleasure!

Connectez-vous pour commenter.

Plus de réponses (0)

Produits

Version

R2019b

Community Treasure Hunt

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

Start Hunting!

Translated by