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 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?
Réponses:
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. Sis∞,i , i=1,…,N, sont les pôles de la fonction de transfert, vous pouvez écrire la fonction de transfert comme
Notez ques 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 plans 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 es∞t à 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 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 transformationZ 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=z−1 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 plans 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.
la source
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:
Un exemple simple est un intégrateur H (s) = 1 / s:
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:
É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:
la source
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 .e−jw=zzero/pole, e−jw z
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
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.(1−z0z)⋯11−p0z
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)=(1−z0z)X(z), yn=b0xn+b1xn−1. b0=b1=1 yn=xn+xn−1
Please note that1+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+e−jw=e−jw/2(ejw/2+e−jw/2)=e−jw/22cos(w/2) . Cosine makes it to behave like low-pass filter
It is also a good lesson that you get2cosα=eiα+e−iα because you will supply the real signals rather than complex imaginary ones in real life.
LTI with impulse response = {1,-1} isyn=xn−xn|xn=ejwn=ejwn(1−e−jw) has transfer function of H(jw)=(1−e−jw)=e−jw/2(ejw/2−e−jw/2)=e−jw2sin(w/2) , which has zero at w=0 since sin(0)=0 but it can be found from the frequency response
After the textbooks, I can spot the surprising coincidence between transfer functionH(z)=1±z and frequency response H(jw)=1±e−jw . That is, z somehow corresponds to e−jw , which is important for zero/pole analysis. I read it like
In general, single-zero LTI is given byyn=b0xn+b1xn−1 or
which goes to zero when1−z0e−jw=0 or e−jw=1/z0 , which matches the computation for z if z=e−jw . 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=e−jw 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=e−jw.
Now, what about the poles? Let's single out a single polea . The system has a from of yn=ayn−1+(xn+xn−1+⋯) , under assumption y0=0 , has z-transform of Y(z)=X(z)/(1−az) .
The feedbacka is equivalent to infinite impulse response 1,a,a2,…↔z1+az+a2z2+⋯=1/(1−az) . It says that response is infinite when z=1/a . What does it mean if we apply the test signal
That is, zeroes or poles of the transfer functionH(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(n−1)=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.
la source