ukmuiik

Question

I am wondering how to plot complex poles/zeros into Bode plot (signal magnitude vs. circular frequency (or $jomega$)). Complex poles/zeros differ from simple poles/zeros in such way that complex ones include imaginary part + real part, while simple ones only real part.

Simple poles/zeros can be directly plotted into Bode plot, just by knowing their real value. While complex poles cannot be so easily plotted (I guess), since they include imaginary part.

Although x-axis of Bode plot is represented with quantity that includes imaginary part (notice $jomega$) I am not able to "insert" those complex poles/zeros into Bode plot.

I will show what I mean on following example of some transfer function: (result for this example can be verified by Wolfram Alpha)

$$H(s)=frac10^5*s(s+100)(s+10)^2*(s^2+400s+10^6)$$

Zeros from numerator are $s=0$ and $s=-100$. Poles from denominator are $s=-10$ (2nd order pole) and $s=-200*(-1-2isqrt6)$ , $s=-200*(-1+2isqrt6)$.

Graphed transfer function can be seen as:

enter image description here

Effect of simple poles and zeros on transfer function can be seen at $omega=10$ and $omega=100$, but effect of complex poles can be seen at $omega=1000$. I am curious whether the term $s=-200*(-1pm2isqrt6)$ equals to $-1000?$ Since a double pole is located at $omega=1000$.

The problem here is that I cannot "convert" complex number into a real one. Or can be transfer function with complex poles/zeros graphed only with the help of a program and not possible to do it by hand? But even a program must be based on some "knowledge" so that it can perform tasks like this one.

I see from the terms you use that you are mixing up different properties of poles and zero. Complex and double are different concepts, as are real and simple. Stop a while to look into these important questions. — Aug 3 at 18:28

score 2 · Accepted Answer · 2018-08-04 19:49:34Z

$H(jomega)$
In these three diagrams, $H(s)$ has been computed at three imaginary points $s=jomega$ (represented by the red points on the imaginary axis) around $omega=200times 2sqrt6$, the imaginary part of the zeros of $s^2+400s+10^6$.

In each of these diagrams it is represented the complex plane with poles and zeros of $H(s)$: poles with an asterisk and zeros with an 'o'.
So you can see the two zeros, both real, at $0$ and $-100$ and the poles, one real at $-10$ (really it is a double pole, even though it is not evident from the diagram) and the other two, complex conjugate, at the positions you have already computed.

Then from each of these zeros and poles, a line (a vector) has been drawn to the red point on the imaginary axis, where we want to compute $H(s)$. (To make things a little clearer lines from zeros are dashed, while those from poles are solid).

Then the modulus of $H(s)$ is the product of the length of the dashed lines divided by the product of the length of the solid lines, multiplied by the gain of $H(s)$ that in your case amounts to $10^5$.

Please note that moving the red point from the position in the first diagram to that of the third, all the vectors are affected by a variation of their length. But the verctor that has the greatest relative variation of its length, as you can guess, is the smallest vector, that is, the vector starting from the complex pole at $-200(-1+2isqrt6)$. So to a first approximation the other vectors can be considered constant and the overall change in the modulus of $H(s)$ around $s=j200times 2sqrt6$ is due to the variation of the smallest vector only.
During the move of the red point, as you can see, the length of small vector decreases reaching a mininum on the second diagram and then increases again.

This length is at the denominator of the formula for computing the modulus of $H(s)$, so you can conclude that in that same range the modulus of $H(s)$ increases, then reaches a maximum, and then decreases again. The maximum of $vert H(s)vert$ is not exactly reached at the same point where the length of the small vector reaches its minimum, because we made some approximation. The smaller is the smallest vector w.r.t. all the other the better is the approximation.

That is to say, the farther poles and zeros are w.r.t. $s=-200(-1+j2sqrt6)$ and the nearer this pole is to the imaginary axis (that is, the smaller is its real part w.r.t. its imaginary part), the nearer a maximum of $vert H(jomega)vert$ is to $omega=200times 2sqrt6$.

EDIT:

The variable $s=sigma+jomega$ is a complex number. A complex number is a point on a plane (like a real number is a point on a line). The real part a complex number $sigma$ is the $x$ coordinate of the point in the plane. The imaginary part $jomega$ is its $y$ coordinate.
The modulus of $s$ is the length of the line from the origin of the plane to $s$.

Just to show some example of using the complex model, let's compare how curves are represented parametrically on a plane using the traditional real coordinate system and the complex one.

A parametric curve on a plane is represented by a system of two equations (because in a plain a point has two coordinates, $x$ and $y$) in one parameter $t$ (because one is the geometrical dimension of a curve)
beginequation
begincases
x=x(t)\
y=y(t)
endcases
endequation

For example a straight line passing through the point $(x_0,y_0)$ with a slope of $theta$ radians w.r.t. the $x$-axis is given by these equations:
beginequation
begincases
x=x_0 + tcos theta \
y = y_0 +tsin theta
endcases
endequation
where $t$ runs on $mathbbR$.

A circle of radius $R$ centered at $(x_0,y_0)$ has equations:
beginequation
begincases
x=x_0 + Rcos t \
y = y_0 +Rsin t
endcases
endequation
where $t$ runs in $[0,2pi[$.

Using the complex model, you can think that every point in the plane have a single complex coordinate instead of two reals.
Then a parametric curve is given by only one equation, like this:
beginequation
s=s(t)
endequation

And again the same straight line as before has equation:
$$s = s_0 + te^jtheta$$
where $t$ runs on $mathbbR$ and $s_0 = x_0+j y_0$.

While the same circle as before has equation:
$$s=s_0 + Re^jt$$
where $t$ runs in $[0,2pi[$ and $s_0 = x_0+j y_0$.

Now a word on the funtions.
A real function of a real variable can be plotted as a curve in a 2D diagram: one axis is used for the variable, the other for the values of the function. A complex function of a complex variable like $H(s)$ cannot be plotted because it would require a 4D diagram: two axes for the variable $s$, and two for the complex values of the function!
On the contrary, a real function of a complex variable like $vert H(s)vert$ can be plotted as a surface on a 3D diagram: two axes for the variable $s$, one for the value of the function.

If you have a function as simple as this $$H(s) = s$$
you cannot plot it. But if you are interested only on the modulus $vert H(s)vert$ then you can. Think of what is the shape of the surface. In the origin of the plane it is$vert H(0)vert=0$. On every points of the unit circle centered at the origin of the plane it is $vert H(1e^jomega)vert = 1$ (do you see the unit circle here? I used its parametric equation, that you must know as of now)

On all points of the circle of radius $2$ centered at the origin of the plane it is $vert H(2e^jomega)vert = 2$, ...

So the shape of the surface of this simple $vert H(s)vert$ is an inverted cone with the vertex in the origin of the plane.

The shape of the modulus of $H(s)=s-a$ is an inverted cone with vertex in the point $a$ of the plane.

You can now guess how to construct the surface described by the modulus of $H(s)=frac1s$. It is $frac1vert svert$, that is, in every point $s$ it is the reciprocal of the length of the line from $0$ (the origin of the plane) to the point of complex coordinate $s$. In $0$ in particular $H(s)$ is not defined, because it goes to $infty$ (the length of the line is zero). It is again a surface with a cylindrical symmetry centered on the axis dedicated to the function values, the $z$-axis. If you make an intersection between this surface and a plane that contains the $z$-axis and the $x$-axis you will have a broken diagram at $x=0$ where the left and right branches both go to $infty$ and the left curve decreases to zero for $x$ that goes from $0$ to $-infty$, while the right curve decreases to zero for $x$ that goes from $0$ to $+infty$. This same diagram is obtained if the intersection is made by the surface and the plane containing the $z$-axis and the $y$-axis. And in general the diagram is the same for any plane containing the $z$-axis whatever is its orientation (this is the meaning of cylindrical symmetry).

If you have $H(s)=frac1s-a$, its modulus is analogous to the previous one, only that it is translated parallel to the plane to bring the axis of symmetry at the point $a$ because it is $vert H(s)vert=frac1vert s-avert$ that is the reciprocal of the lenght of the line segment joining $a$ and $s$.

If you have a complex function of a complex veriable $s$ that is a fraction of polynomial of the first degree $s-a$ like for instance:
$$H(s)=Kfrac(s-z_1)(s-z_2)s(s-p_1)(s-p_2)$$
its modulus is given by:
$$vert H(s)vert=vert Kvertfracvert s-z_1vert vert s-z_2vertvert svert vert s-p_1vert vert s-p_2vert$$
that is, the product of the length of the line segments joining $s$ to each zero $z_i$ divided by the product of the length of the line segments joining $s$ to each pole $p_i$ multiplied by the constant $vert Kvert$.
Please note that if you have (at numerator or at denominator) a polynomial raised to some integer positive power $(s-a)^2$ you need to consider that line segment as many times as is expressed by the exponent (in this case 2 times). Note also that in this case $a$ is said a multiple zero or multiple pole, to contrast this case with that of the simple zero or simple pole (that is when $s-a$ is not raise to power greater than the first).

The plane that I represented in the diagrams attached to the answer is the plane where $s$ lives. On each of this points $vert H(s)vert$ takes a value (ploes excluded). In particular its zeros are the points in the plane where $vert Hvert$ is zero. In your case there are only real zeros, that is, two points on the $x$-axis where $vert H(s)vert$ is zero. Poles are points where $vert H(s)vert$ goes to $infty$. In your case there are 3 poles: one real, that is, on the $x$-axis and 2 complex conjugate, that is, they are neither on the $x$-axis nor on the $y$-axis and they are placed simmetrically w.r.t. the $x$-axis.

Now you are interested to $vert H(jomega)vert$, that is, only on the values that $vert H(s)vert$ takes on the imaginary axis (= the $y$-axis). You can obtain it by intersecting the surface represented by $vert H(s)vert$ with the plane containing the $z$-axis and the $y$-axis. So $vert H(jomega)vert$ is a function that can be plotted as a curve on a 2D diagram, because you need one axis for $omega$ and the other for the values of $vert H(jomega)vert$.

In my original answer I tried to let you see why you have the evidence of the presence of a pole at $omegaapprox 1000$. Now if have followed me you should understand that the pole in not at $$omegaapprox 1000$$ but you have only a perturbation of the plot of $vert H(jomega)vert$ at $omegaapprox 1000$ because of the presence of a pole in a point whose imaginary part (that is, $y$-coordinate) is $omegaapprox 1000$.

Hmm... I am really struggling to understand at least few things you wrote, but am unable to do so. I feel like these two answers are at least ten levels above my current understanding of this topic we are discussing. I feel like don't even understand the foundations regrading complex numbers, complex plane, Bode plot, etc. — Aug 4 at 16:07
I am so confused right now I don't even know where to start learning this stuff. I feel like there is a ton of literature behind your writing. — Aug 4 at 16:08
I just "bumped" into this topic, without any lessons on complex numbers, etc. before doing so, so there might be few reasons I am not understanding even a bit of your answer. — Aug 4 at 16:10
I try to give you some notion, adding them on the answer. Hoping that it helps you. — 2 days ago

score 2 · Answer 2 · 2018-08-03 03:04:41Z

I believe you have trouble with understanding what the Bode plot actually is. I suggest returning to the definition. Let's look at the Bode magnitude plot. It is the graph of the function:
$$mid H(s = jomega) mid $$
The abscissa is logarithmic and the ordinate is the magnitude given in decibels. In order to understand what is happening let's look at a simple example of a transfer function:
$$ H(s) = frac1s^2 + 1$$
We obviously have two poles at: $s = j$ and $s = -j$. Let us now think how we could graph this in the described way. We could say:
$$A(omega) = 20 log mid H(j omega) mid $$
This makes it simple for us to draw the plot by hand in order to get a better understanding of what is happening. So we have:
$$A(omega) = 20 log fracmid 1 midmid -omega^2 + 1 mid = -20 log sqrt(-omega^2 + 1)^2 $$
For $omega leq 1$ we have(obviously we are only considering positive frequencies):
$$A(omega) = -20 log (-omega^2 + 1)$$
and for $omega > 1$ we have:
$$A(omega) = -20 log (omega^2 - 1)$$
It doesn't really matter where we include the pole $s = j omega$ for $omega = 1$ we get that as $omega$ is getting close to $1$ our Bode magnitude plot is going towards infinity, as the graph shows:
Bode plot
Finally concerning your question about your particular example whether(where $s_c$ is your complex pol) "$s = s_c$ equals to $-1000$" or not, I suggest doing the same thing I did for the previous simple example. If we look at your transfer function, we can isolate the part with the complex poles:
$$H_c(s) = frac1s^2 + 400 s + 10^6$$
If I haven't messed up when taking the complex modulus we end up with:
$$ frac1sqrt160000 omega^2 + (10^6 - omega^2)^2 $$

so for $omega = 200 sqrt23 approx 1000$ we do have an extreme and that's why we can at that point see influence.

score 2 · Accepted Answer · 2018-08-04 19:49:34Z