Как полюсы связаны с частотной характеристикой


16

Недавно я впал в заблуждение , рассматривая полюс s = 1, поскольку на частоте 1 есть бесконечный отклик. Тем не менее, отклик был только 1. Теперь, вы можете получить частотный отклик, учитывая полюса?

Во-вторых, теория говорит, что система устойчива, когда полюса находятся в левой s-плоскости и, таким образом, затухают во времени. Но ждать. Разве «полюс» означает бесконечный отклик - рост во времени?

Наконец, это правильный вопрос в DSP? IMO, D обозначает цифровой, тогда как s-домен является аналоговым. Я не нахожу теги s-плоскости или преобразования Лапласа для маркировки моего поста.

Обновление Спасибо за ответы. Кажется, что у меня есть это, кроме одной незначительной, но фундаментальной вещи - отношения полюсов (и нулей) с частотой. По сути, почему собственные значения (или, как вы называете оператор / переменную s ) связаны с частотой? Это должно быть как-то связано с экспоненциальным ростом и преобразованием Лапласа. Я вполне понимаю, что полюса оказываются собственными значениями (особенно для дискретных повторений). Но как это связано с частотой?


1
Это «Обмен стека обработки сигналов», а не «Обмен стека DSP». :)
эндолит

4
Да, как уже упоминалось, обработка аналоговых сигналов идет по теме. DSP.SE было подходящим названием для первоначального запуска, но сигналов .stackexchange.com теперь ссылки здесь также.
обработке данных

Что именно вы имеете в виду, когда спрашиваете об отношениях между поляками и частотами?
Сударсан

1
Очевидно, именно так и почему полюса определяют частотную характеристику.
Вэл

Ответ уже дан, я думаю. Частотный отклик - это величина отклика системы при перемещении вдоль оси jω . Если вы преобразовали функцию переноса системы H(s) в произведение 1/(spi) и (szi) , все, что вам нужно сделать, - это найти величину при s=jω для переноса. Функция, и это, очевидно, определяется расположением полюсов и нулей, поскольку именно они будут отображаться в факторизованном ответе системы.
Сударсан

Ответы:


12

Я думаю, что на самом деле есть 3 вопроса в вашем вопросе:

Q1: Могу ли я получить частотную характеристику с учетом полюсов (линейной инвариантной по времени) системы?

Да, вы можете, до постоянной. Если s,i , i=1,,N, являются полюсами передаточной функции, вы можете записать передаточную функцию как

(1)H(s)=k(ss,1)(ss,2)(ss,N)

Отметим, что s является комплексной переменной s=σ+jω , а частотная переменная ω соответствует мнимой оси комплексной s плоскости. Теперь нам нужно получить частотную характеристику от передаточной функции. Для стабильных систем это можно сделать просто путем оценки передаточной функции H(s) для s=jω . Итак, вы заменяете s на jω в (1) и все готово. Заметим, однако, что это верно только для стабильных систем (т.е. если область сходимости H(s) включает ось jω ).

Q2: Как стабильная система может иметь полюса?

Как вы уже знаете, для причинных и устойчивых систем все полюса должны лежать в левой полуплоскости комплексной s плоскости. В самом деле, значение передаточной функции H(s) уйдет в бесконечность при полюсе s=s , но частотная характеристика будет в порядке, потому что, если все полюса находятся в левой полуплоскости, на jω Оу (или справа от него). Если вы посмотрите на него во временной области, то каждый (простой) полюс имеет вклад est в импульсную характеристику системы. Если полюс находится в левой полуплоскости, это означает, что s=σ+jω has a negative real part σ<0. So

est=eσejω

is an exponentially damped function and does not grow but decays, because σ<0.

Q3: Does this question belong here?

Other community members have to judge whether this question belongs here. I think that it does. It is obviously not directly related to pure DSP, but DSP engineers very often also have to deal with analog signals and systems before AD conversion, so they also know about continuous system theory. Second, almost all DSP people (at least the ones with traditional training) got quite some exposure to general signals and systems theory, including continuous-time and discrete-time systems.

By the way, for discrete-time systems you get the Z-transform instead of the Laplace-transform, and your complex variable is now called z instead of s. The variable D that you've mentioned is defined as D=z1 and is mainly used in the coding literature. By its definition, it denotes a delay element, so D stands for "delay" (not "digital").

If you know that the left half-plane of the complex s-plane maps to the region inside the unit circle of the complex z-plane (i.e. |z|<1), and the jω-axis maps to the unit circle |z|=1, then almost everything you know about one of the two domains will easily carry over to the other domain.


I think that frequency response involves complex conjugation in addition s in H(s) for s=jω.
Val

16

One thing that really helped me understand poles and zeros is to visualize them as amplitude surfaces. Several of these plots can be found in A Filter Primer. Some notes:

  • It's probably easier to learn the analog S plane first, and after you understand it, then learn how the digital Z plane works.
  • A zero is a point at which the gain of the transfer function is zero.
  • A pole is a point at which the gain of the transfer function is infinite.
  • Often there are zeros or poles at infinity, which aren't always included in descriptions of the transfer function, but are necessary to understand it.
  • The frequency response in the S plane happens along the jω axis only.
    • The origin is 0 Hz, or DC, and the cutoff frequency of filters increases radially away from the origin. Putting a pole at any point along a circle at a certain distance from the origin will produce the same cutoff frequency.
    • To increase the cutoff frequency of a filter, move the poles radially outward.
    • To increase the Q of a biquad filter, move the poles along the circle towards the jω axis, which keeps the cutoff frequency constant, but increases the effect that the pole has on the frequency response, making it more "peaky".
    • Moving poles along a circle keeps cutoff frequency constant but changes Q
  • If a zero appears on the jω axis, then the frequency response will drop to zero at that frequency; if you input a sine wave at that frequency, the output will be 0.
  • Если полюс появляется на оси jω, то импульсный отклик является осциллятором; любой импульс заставит его звонить вечно на этой частоте. Импульсы имеют конечную энергию, но отклик фильтра имеет бесконечную энергию, поэтому он имеет бесконечное усиление.

Простым примером является интегратор H (s) = 1 / s:

  • Эта функция равна 0, когда s бесконечно, поэтому она имеет ноль на бесконечности.
  • Эта функция равна бесконечности, когда s равен нулю, поэтому она имеет полюс в нуле.

Другими словами, он имеет бесконечный коэффициент усиления при постоянном токе (шаговая характеристика интегратора постоянно увеличивается), а коэффициент усиления уменьшается с увеличением частоты:

Bode plot of integrator

Перемещение полюса от начала координат вдоль мнимой оси в левую сторону плоскости S делает усиление при 0 Гц по оси jw снова конечным, и теперь у вас есть фильтр нижних частот:

enter image description here


+1, nice answer. But I don't understand what you mean by "Any point along a circle at a certain distance from the origin has the same frequency." Curves of constant frequency in the s-plane are lines parallel to the real axis. For circles with origin at s=0 you get σ2+ω2знак равносоNsT, where s=σ+jω.
Matt L.

He seems to confuse s-plane with z-plane
Val

@MattL.: Hmmm. I'm thinking of the poles of a Nth-order Butterworth filter being along a circle equidistant from the origin, for instance, or the poles of a biquad moving along a circle equidistant from the origin as you adjust the Q of the filter while keeping the frequency constant, or changing the cutoff of a filter by moving the poles closer to or away from the origin in a radial direction, or converting lowpass to highpass by inverting the poles about the unit circle. How should I reword this?
endolith

@Val: Cutoff frequency. I've already edited the post to correct it.
endolith

6
Val, No need for a douchy snarky comment to @endolith.
Spacey

1

I won't tell the full mapping from poles(1)/zeroes(0) to the frequency response but I think I can explain the connection between frequency and zero/infinite response, why do you have infinite/zero response at ejw=zzero/pole, i.e. what ejw has to do with z.

The general form of the linear system is

yn+a1yn1+a2yn2+=b0xn+b1xn1+b2xn2+,
which can be solved in z-from as
Y(z)=(b0+b1z+b2z2+)(1+a1z+a2z2+)X(z)=H(z)X(z)=(1z0z)(1z1z)(1p0z)(1p1z)X(z).

In the end, the series of binomial products (1z0z)11p0z can be considered as a series of systems, where first output, is the input for another.

I would like to analyze the effect of single pole and zero. Let's single out the first zero, considering it the transfer function so that the rest of H(z)X(z) is the input signal, Y(z)=(1z0z)Χ(z), which corresponds to some yn=b0xn+b1xn1. Let's take b0=b1=1 for simplicity. I mean that yn=xn+xn1.

What we want to determine the effect of the system H(z) upon harmonic signal. That is, the input is going to be test signal

xn=ejwnz1+ejwz+e2jwz2+=1/(1ejw)=X(z).
The response is going to be
yn=xn+xn1|xn=ejwn=ejwn+ejw(n1)=ejwn(1+ejw)
that is, 1+ejw is the transfer function or Y(z)=(1+z)(1ejwz)=(1+z)X(z).

Please note that 1+z basically says that output is sum of input signal plus shifted signal, since single z stands for single clock delay in time domain.

Now, as explained in, H(jw)=1+ejw=ejw/2(ejw/2+ejw/2)=ejw/22cos(w/2). Cosine makes it to behave like low-pass filter

{w=0H(j0)=12cos(0)=2w=πH(jπ)=ejπ/22cos(π/2)=0

It is also a good lesson that you get 2cosα=eiα+eiα because you will supply the real signals rather than complex imaginary ones in real life.

LTI with impulse response = {1,-1} is yn=xnxn|xn=ejwn=ejwn(1ejw) has transfer function of H(jw)=(1ejw)=ejw/2(ejw/2ejw/2)=ejw2sin(w/2), which has zero at w=0 since sin(0)=0 but it can be found from the frequency response

H(jw)=1ejw=0ejw=1=e0w=0.

After the textbooks, I can spot the surprising coincidence between transfer function H(z)=1±z and frequency response H(jw)=1±ejw. That is, z somehow corresponds to ejw, which is important for zero/pole analysis. I read it like

sine z-factor stands for a clock shift and yn=xn±xn1=0 means that next sample is ± previous one to get zero response, we need to have 1±z=0 in front of X(z). But, the frequency domain basis functions ejwn evolve by multiplying current value ejw(n1) with ejw every clock. Therefore, we have ejwn(1±ejw)=0 as condition for constant zero output. The latter 1±ejw matches perfectly with zero transfer function 1±z=0.

In general, single-zero LTI is given by yn=b0xn+b1xn1 or

Y(z)=(b0+b1z)X(z)=(b0+b1z)(1+x1z+x2z2+)=b0+(b0x1+b1x0)z+(b0x2+b1x1)z2+.
When b0+b1z=0, i.e. when z=b0/b1, whereas frequency response is,
yn(xn=ejwn)=b0ejwn+b1ejw(n1)=ejwn(b0+b1ejw)=ejwnb0(1z0ejw),

which goes to zero when 1z0ejw=0 or ejw=1/z0, which matches the computation for z if z=ejw. The only thing that bothers me is that fixed-amplitude complex exponential is not enough for the frequency (harmonic) basis. You cannot obtain arbitrary ratio 1/z0=ejw by choosing appropriate frequency w, a decaying harmonic signal is needed for that. That is weird because I have heard that any signal can be represented as sum of (constant amplitude) sines and cosines. But, anyway, we see that system zero stands for relationship between adjacent samples of input signal. When they are right, the output is identically 0 and we can choose such such frequency w so that zero z=1/z0=ejw.

Now, what about the poles? Let's single out a single pole a. The system has a from of yn=ayn1+(xn+xn1+), under assumption y0=0, has z-transform of Y(z)=X(z)/(1az).

The feedback a is equivalent to infinite impulse response 1,a,a2,z1+az+a2z2+=1/(1az). It says that response is infinite when z=1/a. What does it mean if we apply the test signal

xn=ejwnzX(z)=1+ejwz+e2jwz2+=1/(1ejwz)
to our system? We'll get Y(z)=11az11ejwz, or
yn=ejwn+aejw(n1)+a2ejw(n2)+=ejwn(1+aejw+a2e2jw+)=ejwn1aejw.
That is, frequency response is 1/(1aejw), which goes to infinity when ejw=1/a, the same as zpole above, ejw=zpole=1/a. But again, you can not always arrive at the pole 1/a adjusting the frequency w alone. The frequency basis functions must be decaying amplitude in general and look like (kejw)n.

That is, zeroes or poles of the transfer function H(z) happen to match the zeroes and poles of frequency response H(jw), which is really amazing. I noticed that this is related to the relation between adjacent samples, ejwn/ejw(n1)=ejw=1/zzero in case of zeroes. The fact that ejwn scales exponentially over time, along with the system with feedback a, also seems to be the key for matching between ejw and zpoles. It also seems important that you cannot simply look for the appropriate frequency of ejwn, the basis function must also have adjustable amplitude factor kn.

I would be happy if anybody could explain the same more condensely or more crisply.

Используя наш сайт, вы подтверждаете, что прочитали и поняли нашу Политику в отношении файлов cookie и Политику конфиденциальности.
Licensed under cc by-sa 3.0 with attribution required.