Comment les pôles sont liés à la réponse en fréquence


16

Je suis récemment tombé dans l'erreur , considérant le pôle s = 1 car il y a une réponse infinie à la fréquence 1. Pourtant, la réponse n'était que de 1. Maintenant, pouvez-vous dériver la réponse en fréquence, étant donné les pôles?

Deuxièmement, la théorie dit qu'un système est stable lorsque les pôles sont dans le plan s gauche et, par conséquent, se désintègrent dans le temps. Mais attendez. Est-ce que «pôle» signifie la réponse infinie - la croissance dans le temps?

Enfin, est-ce une bonne question dans DSP? IMO, D signifie numérique tandis que le domaine s est analogique. Je ne trouve pas de balises de transformation s-plane ou Laplace pour étiqueter ma publication.

mise à jour Merci pour les réponses. Il semble que je l'ai, sauf la seule chose mineure mais fondamentale - la relation des pôles (et des zéros) avec la fréquence. Au fond, pourquoi (ou valeurs propres, comment appelez-vous le s opérateur / variable) liée à la fréquence? Cela devrait être en quelque sorte lié à la croissance exponentielle et à la transformation de Laplace. Je comprends très bien que les pôles se trouvent être des valeurs propres (en particulier pour les récurrences discrètes). Mais, comment est-ce lié à la fréquence?


1
C'est "échange de pile de traitement du signal", pas "échange de pile DSP". :)
endolith

4
Oui, comme mentionné plus haut, le traitement du signal analogique est un sujet. DSP.SE était un nom approprié pour le lancement initial, mais signaux.stackexchange.com propose désormais des liens ici également.
datageist

Que voulez-vous dire exactement lorsque vous demandez la relation entre les pôles et les fréquences?
Sudarsan

1
Évidemment, c'est comment et pourquoi les pôles déterminent la réponse en fréquence.
Val

Je suppose que la réponse a déjà été donnée. La réponse en fréquence est la magnitude de la réponse du système lorsque vous vous déplacez le long de l' axe jω . Si vous avez factorisé la fonction de transfert du système H(s) en produit de 1/(spi) et (szi) , tout ce que vous devez faire est de trouver l'amplitude à s=jω pour le transfert La fonction et celle-ci est évidemment déterminée par l'emplacement des pôles et des zéros, car ce sont ceux qui apparaissent dans la réponse du système factorisé.
Sudarsan

Réponses:


12

Je pense qu'il y a en fait 3 questions dans votre question:

Q1: Puis-je dériver la réponse en fréquence compte tenu des pôles d'un système (invariant dans le temps linéaire)?

Oui, vous pouvez, jusqu'à une constante. Si s,i , i=1,,N, sont les pôles de la fonction de transfert, vous pouvez écrire la fonction de transfert comme

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

Notez que s est une variable complexe s=σ+jω , et la variable de fréquence ω correspond à l'axe imaginaire du plan s complexe . Maintenant, nous devons obtenir la réponse en fréquence de la fonction de transfert. Pour les systèmes stables, cela peut simplement être fait en évaluant la fonction de transfert H(s) pour s=jω . Vous remplacez donc s par jω dans (1) et vous avez terminé. Notez cependant que cela n'est vrai que pour les systèmes stables (c'est-à-dire si la région de convergence de H(s) inclut l’axejω ).

Q2: Comment un système stable peut-il avoir des pôles?

Comme vous le savez déjà, pour les systèmes causaux et stables, tous les pôles doivent se trouver dans le demi-plan gauche du plan s complexe . En effet, la valeur de la fonction de transfert H(s) ira à l'infini à un pôle s=s , mais la réponse en fréquence sera OK, car si tous les pôles sont dans le demi-plan gauche, il n'y a pas de pôles sur le jω -axis (ou à sa droite). Si vous le regardez dans le domaine temporel, chaque pôle (simple) a une contribution de est à la réponse impulsionnelle du système. Si le pôle est situé dans le demi-plan gauche, cela signifie que s=σ+jω a une partie réelle négativeσ<0 . Donc

est=eσejω

est une fonction exponentiellement amortie et ne croît pas mais décroît, car σ<0 .

Q3: Cette question appartient-elle ici?

Les autres membres de la communauté doivent juger si cette question appartient ici. Je pense que oui. Il n'est évidemment pas directement lié au DSP pur, mais les ingénieurs DSP doivent très souvent également traiter les signaux analogiques et les systèmes avant la conversion AD, de sorte qu'ils connaissent également la théorie des systèmes continus. Deuxièmement, presque toutes les personnes DSP (au moins celles ayant une formation traditionnelle) ont été assez exposées aux signaux généraux et à la théorie des systèmes, y compris les systèmes à temps continu et à temps discret.

Soit dit en passant, pour les systèmes à temps discret, vous obtenez la transformation Z au lieu de la transformation Laplace, et votre variable complexe est maintenant appelée z au lieu de s . La variable D que vous avez mentionnée est définie comme D=z1 et est principalement utilisée dans la littérature de codage. Par sa définition, il désigne un élément de retard, donc D signifie "retard" (pas "numérique").

Si vous savez que le demi-plan gauche du plan s complexe correspond à la région à l'intérieur du cercle unité du plan z complexe (c.-à-d. |z|<1 ), et que l' axe jω correspond au cercle unité |z|=1 , alors presque tout ce que vous savez sur l'un des deux domaines sera facilement transféré sur l'autre domaine.


Je pense que la réponse en fréquence implique une conjugaison complexe en plus s dans H (s) pour s = jω.
Val

16

Une chose qui m'a vraiment aidé à comprendre les pôles et les zéros est de les visualiser comme des surfaces d'amplitude. Plusieurs de ces tracés se trouvent dans A Filter Primer . Quelques notes:

  • Il est probablement plus facile d'apprendre d'abord le plan S analogique, puis après l'avoir compris, découvrez comment fonctionne le plan Z numérique.
  • Un zéro est un point auquel le gain de la fonction de transfert est nul.
  • Un pôle est un point où le gain de la fonction de transfert est infini.
  • Il y a souvent des zéros ou des pôles à l'infini, qui ne sont pas toujours inclus dans les descriptions de la fonction de transfert, mais sont nécessaires pour la comprendre.
  • La réponse en fréquence dans le plan S se produit uniquement le long de l'axe jω.
    • L'origine est 0 Hz, ou DC, et la fréquence de coupure des filtres augmente radialement loin de l'origine. Mettre un pôle à n'importe quel point le long d'un cercle à une certaine distance de l'origine produira la même fréquence de coupure.
    • Pour augmenter la fréquence de coupure d'un filtre, déplacez les pôles radialement vers l'extérieur.
    • Pour augmenter le Q d'un filtre biquad, déplacez les pôles le long du cercle vers l'axe jω, ce qui maintient la fréquence de coupure constante, mais augmente l'effet que le pôle a sur la réponse en fréquence, le rendant plus "pointu".
    • Moving poles along a circle keeps cutoff frequency constant but changes Q
  • Si un zéro apparaît sur l'axe jω, alors la réponse en fréquence tombera à zéro à cette fréquence; si vous entrez une onde sinusoïdale à cette fréquence, la sortie sera 0.
  • Si un pôle apparaît sur l'axe jω, alors la réponse impulsionnelle est un oscillateur; toute impulsion le fera sonner pour toujours à cette fréquence. Les impulsions ont une énergie finie, mais la réponse du filtre a une énergie infinie, donc elle a un gain infini.

Un exemple simple est un intégrateur H (s) = 1 / s:

  • Cette fonction est égale à 0 lorsque s est infini, elle a donc un zéro à l'infini.
  • Cette fonction est égale à l'infini lorsque s est nul, elle a donc un pôle à zéro.

En d'autres termes, il a un gain infini à DC (la réponse en échelon d'un intégrateur est en constante augmentation), et le gain diminue à mesure que la fréquence augmente:

Bode plot of integrator

Éloigner le pôle de l'origine, le long de l'axe imaginaire dans la main gauche du plan S, rend à nouveau le gain à 0 Hz sur l'axe jw fini, et maintenant vous avez un filtre passe-bas:

enter image description here


+1, belle réponse. Mais je ne comprends pas ce que vous entendez par «tout point le long d'un cercle à une certaine distance de l'origine a la même fréquence». Les courbes de fréquence constante dans le plan sont des lignes parallèles à l'axe réel. Pour les cercles d'origine à s = 0, vous obtenez σ 2 + ω 2 = c o n s t , où s = σ + j ω . ss=0σ2+ω2=consts=σ+jω
Matt L.

Il semble confondre avion S et avion Z
Val

@MattL .: Hmmm. Je pense aux pôles d'un filtre Butterworth d'ordre N se trouvant le long d'un cercle équidistant de l'origine, par exemple, ou aux pôles d'un biquad se déplaçant le long d'un cercle équidistant de l'origine lorsque vous ajustez le Q du filtre tout en gardant la constante de fréquence, ou changer la coupure d'un filtre en rapprochant ou éloignant les pôles de l'origine dans une direction radiale, ou en convertissant le passe-bas en passe-haut en inversant les pôles autour du cercle unitaire. Comment devrais-je reformuler cela?
endolith

@Val: fréquence de coupure . J'ai déjà édité le post pour le corriger.
endolith

6
Val, pas besoin d'un commentaire snarky douchy à @endolith.
Spacey

1

Je ne dirai pas le mappage complet des pôles (1) / zéros (0) à la réponse en fréquence mais je pense que je peux expliquer le lien entre la fréquence et la réponse zéro / infinie, pourquoi avez-vous une réponse infinie / zéro à c'est -à- dire ce que e - j w a à voir avec z .ejw=zzero/pole,ejwz

La forme générale du système linéaire est ce qui peut être résolu en z-from comme Y ( z ) = ( b 0 + b 1 z + 2

yn+a1yn1+a2yn2+=b0xn+b1xn1+b2xn2+,
Y(z)=(b0+b1z+b2z2+)(1+a1z+a2z2+)X(z)=H(z)X(z)=(1z0z)(1z1z)(1p0z)(1p1z)X(z).

Au final, la série de produits binomiaux peut être considéré comme une série de systèmes, où la première sortie est l'entrée d'un autre.(1z0z)11p0z

Je voudrais analyser l'effet du monopôle et du zéro. Distinguons le premier zéro, en le considérant comme la fonction de transfert de sorte que le reste de soit le signal d'entrée, Y ( z ) = ( 1 - z 0 z ) Χ ( z ) , qui correspond à certains y n = b 0 x n + b 1 x n - 1 . Prenons b 0H(z)X(z)Y(z)=(1z0z)Χ(z),yn=b0xn+b1xn1.b0=b1=1yn=xn+xn1

xn=ejwnz1+ejwz+e2jwz2+=1/(1ejw)=X(z).
yn=xn+xn1|xn=ejwn=ejwn+ejw(n1)=ejwn(1+ejw)
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.

En utilisant notre site, vous reconnaissez avoir lu et compris notre politique liée aux cookies et notre politique de confidentialité.
Licensed under cc by-sa 3.0 with attribution required.