극이 주파수 응답과 관련되는 방법


16

나는 주파수 1에서 무한 응답이 있기 때문에 극 s = 1을 고려 하여 최근 에 오류빠졌다 . 그러나 응답은 단지 1이었다. 이제, 극이 주어지면 주파수 응답을 도출 할 수 있는가?

둘째, 이론에 따르면 극이 왼쪽 s- 평면에있을 때 시스템이 안정적이므로 시간이 지남에 따라 감소합니다. 하지만 기다려. "극점 (pole)"은 무한한 반응, 즉 시간의 성장을 의미합니까?

마지막으로 DSP에서 올바른 질문입니까? IMO, D 는 디지털을 나타내고 s- 도메인은 아날로그입니다. 게시물에 레이블을 지정할 S- 평면 또는 Laplace 변형 태그를 찾지 못했습니다.

업데이트 답변 주셔서 감사합니다. 극소하지만 근본적인 것-극 (및 0)과 주파수의 관계를 제외하고는 그것을 얻은 것 같습니다. 기본적으로, 왜 고유 값 (또는, 어떻게 당신이 전화를 할 수 있습니다 s 연산자 / 변수) 주파수와 관련? 어떻게 든 지수 성장 및 Laplace 변환과 관련이 있어야합니다. 나는 극이 고유 값 (특히 이산 재발의 경우)이라는 것을 이해합니다. 그러나 이것은 주파수와 어떻게 관련이 있습니까?


1
"DSP 스택 교환"이 아니라 "신호 처리 스택 교환"입니다. :)
endolith 2016 년

4
endoith가 언급했듯이 아날로그 신호 처리가 주제입니다. DSP.SE은 초기 출시에 대한 편법 이름 이었지만, signals.stackexchange.com은 지금뿐만 아니라 여기에 연결합니다.
datageist

극점과 주파수 간의 관계를 요청할 때 정확히 무엇을 의미합니까?
Sudarsan

1
분명히, 극이 주파수 응답을 결정하는 방법과 이유입니다.
Val

대답은 이미 추측 한 것입니다. 주파수 응답은 jω 축을 따라 이동할 때 시스템 응답의 크기입니다 . 이 시스템의 전달 함수의 인수 분해 한 경우 H(s) 의 제품으로 1/(spi)(szi) , 당신이 할 필요가 상기 크기 찾을 수 있습니다 s=jω 전송을위한을 기능과 이것은 분명히 요인과 시스템의 반응에 나타날 것이기 때문에 극점과 영점의 위치에 의해 결정됩니다.
Sudarsan

답변:


12

귀하의 질문에 실제로 3 가지 질문이 있다고 생각합니다.

Q1 : (선형시 불변) 시스템의 극점에서 주파수 응답을 도출 할 수 있습니까?

예, 최대 한도까지 가능합니다. 경우 s,i , i=1,,N, 전달 함수의 극, 당신은 같은 전달 함수를 작성할 수 있습니다

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

참고 s 복소 변수 s=σ+jω 및 상기 주파수 가변 ω 복소수의 허수 축에 대응하는 s -plane. 이제 전달 함수로부터 주파수 응답을 얻어야합니다. 안정적인 시스템 경우 s = j ω에 대해 전달 함수 H(s) 를 평가하여 간단히 수행 할 수 있습니다 . 따라서 (1)에서 sj ω로 바꾸면 끝났습니다. 참고 그러나,이 (안정적 시스템 만 사실이라고, 즉 융합의 지역의 경우 H는 ( Ss=jωsjωH(s) 에는jω 축이포함됩니다.

Q2 : 안정적인 시스템은 어떻게 극을 가질 수 있습니까?

아시다시피, 인과적이고 안정적인 시스템의 경우 모든 극은 복합 s 평면의 왼쪽 절반 평면에 있어야합니다 . 실제로 전달 함수 H(s) 의 값은 극점 s=s 에서 무한대로 진행 되지만 주파수 응답은 정상입니다. 모든 극이 왼쪽 절반 평면에 있으면 극점이 없습니다. jω 축 (또는 오른쪽). 시간 영역에서 살펴보면 각 (단순) 극은 시스템의 임펄스 응답 에 est 기여를 갖습니다 . 극이 왼쪽 절반 평면에있는 경우 s=σ+jω 는 음의 실수 부σ<0 갖습니다. 그래서

est=eσejω

σ<0 이므로 지수 적으로 감쇠 된 함수이며 커지지 않지만 쇠퇴 합니다.

Q3 : 이 질문이 여기에 속합니까?

다른 커뮤니티 회원은이 질문이 여기에 속하는지 판단해야합니다. 나는 그렇게 생각합니다. 순수한 DSP와 직접 관련이있는 것은 아니지만 DSP 엔지니어는 AD 변환 전에 아날로그 신호 및 시스템을 다루어야하는 경우가 많으므로 지속적인 시스템 이론에 대해서도 알고 있습니다. 둘째, 거의 모든 DSP 사람들 (적어도 전통적인 훈련을받은 사람들)은 연속 시간 및 이산 시간 시스템을 포함한 일반적인 신호 및 시스템 이론에 상당히 노출되었습니다.

그건 그렇고, 이산 시간 시스템 의 경우 Laplace-transform 대신 Z -transform 을 얻습니다. 복잡한 변수는 이제 s 대신 z 입니다 . 언급 한 변수 DD = z - 1 로 정의되며 주로 코딩 문헌에서 사용됩니다. 정의에 따라 지연 요소를 나타내므로 D 는 "지연"( "디지털"이 아님)을 나타냅니다.sDD=z1D

complex s 평면의 왼쪽 절반 평면이 complex z 평면 의 단위 원 내부 영역에 매핑 됨을 알고있는 경우 (예 : |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.
  • S 평면의 주파수 응답은 jω 축을 따라서 만 발생합니다.
    • 원점은 0Hz 또는 DC이며 필터의 차단 주파수는 원점에서 반경 방향으로 멀어집니다. 원점에서 특정 거리에있는 원을 따라 임의의 지점에 극을 놓으면 동일한 차단 주파수가 생성됩니다.
    • 필터의 차단 주파수를 높이려면 극을 방사형 바깥쪽으로 움직입니다.
    • 바이 쿼드 필터의 Q를 높이려면 극을 jω 축쪽으로 원을 따라 이동하여 차단 주파수를 일정하게 유지하지만 극이 주파수 응답에 미치는 영향을 증가시켜 더 "피크"하게 만듭니다.
    • Moving poles along a circle keeps cutoff frequency constant but changes Q
  • jω 축에 0이 표시되면 해당 주파수에서 주파수 응답이 0으로 떨어집니다. 해당 주파수에서 사인파를 입력하면 출력은 0이됩니다.
  • jω 축에 극이 나타나면 임펄스 응답은 오실레이터입니다. 임펄스가 있으면 해당 주파수에서 영원히 울립니다. 임펄스에는 한정된 에너지가 있지만 필터의 응답에는 무한한 에너지가 있으므로 무한 이득을 갖습니다.

간단한 예는 적분기 H (s) = 1 / s입니다.

  • s가 무한한 경우이 함수는 0과 같으므로 무한대에는 0이 있습니다.
  • 이 함수는 s가 0 일 때 무한대와 같으므로 극점이 0입니다.

즉, DC에서 무한 이득을 가지며 (적분기의 스텝 응답이 계속 증가 함) 주파수가 증가함에 따라 이득이 감소합니다.

Bode plot of integrator

가상 축을 따라 원점에서 멀리 떨어진 S 평면의 왼쪽으로 이동하면 jw 축에서 0Hz의 게인이 다시 유한하게되고 이제 저역 통과 필터가 있습니다.

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=const, 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.

당사 사이트를 사용함과 동시에 당사의 쿠키 정책개인정보 보호정책을 읽고 이해하였음을 인정하는 것으로 간주합니다.
Licensed under cc by-sa 3.0 with attribution required.