Pourquoi les filtres FIR sont-ils toujours stables même s'ils contiennent des pôles?

15
  • Comment se fait-il que les filtres FIR soient toujours stables?

  • Puisqu'ils contiennent des pôles, ne devraient-ils pas être plus affectés par les problèmes de stabilité que les autres?

user7277
la source
FIR est stable si tout son zéro est situé dans le cercle d'unité
dato datuashvili
2
Pas vrai: FIR est toujours stable et les zéros peuvent être où ils veulent, y compris en dehors du cercle unitaire. Exemple: le filtre [1 -6 11 -6] a des zéros à z = 1, 2 et 3
Hilmar
encore une fois, @ Hilmar, cela dépend de la façon dont le FIR est mis en œuvre. Les FIR implémentées en tant que RII tronqué (TIIR) peuvent ne pas être stables à l'intérieur. mis en œuvre comme un simple filtre FIR transversal, oui, qui est toujours stable. il est stable même s'il est implémenté en utilisant une "convolution rapide" (en utilisant une FFT et "overlap-add" ou "overlap-save"). et parfois lorsqu'il est implémenté en tant que filtre TIIR, il est stable (si l'IIF interne est stable). mais une FIR implémentée en tant que TIIR pourrait être instable en interne.
robert bristow-johnson

Réponses:

8

Les filtres FIR ne contiennent que des zéros et aucun pôle. Si un filtre contient des pôles, c'est IIR. Les filtres IIR souffrent en effet de problèmes de stabilité et doivent être manipulés avec précaution.

ÉDITER:

Après une réflexion plus approfondie et quelques gribouillages et recherches sur Google, je pense avoir une réponse à cette question des pôles FIR qui, espérons-le, sera satisfaisante pour les parties intéressées.

Commençant par la transformée en Z d'un filtre FIR apparemment sans polarité: Comme le montre la réponse de RBJ, les pôles FIR sont révélés en multipliant le numérateur et le dénominateur deH(z)parzN: H(z)=b0z N +b1z N - 1 +b2z N - 2 ++bN

H(z)=b0+b1z1+b2z2++bNzN1
H(z)zN Donnant ainsi nosNpôles à l'origine d'un filtre FIR général.
H(z)=b0zN+b1zN-1+b2zN-2++bNzN
N

Cependant, afin de le montrer, l'hypothèse de causalité est placée sur le filtre. En effet, si l'on considère un filtre FIR plus général où la causalité n'est pas supposée: Un nombre différent de pôles(N-k)apparaît à l'origine: G(z)=b0z N +b1z N - 1 +b2z N - 2 ++bN

g(z)=b0zk+b1zk-1+b2zk-2++bNzk-N1
(N-k)
g(z)=b0zN+b1zN-1+b2zN-2++bNzN-k

Ainsi, je conclus ce qui suit:

  • (Répondre à la question d'origine) En général, un filtre FIR a des pôles, bien que toujours à l'origine du plan Z. Parce qu'ils ne sont jamais au-delà du cercle unitaire, ils ne menacent pas la stabilité d'un système FIR.
  • NkNth(k=0)N
  • H(z)=z-1=1z
Kenneide
la source
2
Les filtres IIR ne sont pas très dangereux en fait.
user7358
19

z=0 .

parce que tous les pôles sont situés à l'intérieur du cercle unitaire, le filtre FIR est ostensiblement stable.

ce n'est probablement pas le filtre FIR auquel l'OP pense, mais il existe une classe de filtres FIR appelés filtres IIR tronqués (TIIR) qui peuvent avoir un pôle sur ou à l'extérieur du cercle unitaire qui est annulé par un zéro au même endroit. l'exemple le plus simple est la somme mobile ou la moyenne mobile. mais, du point de vue des E / S, ces filtres TIIR sont FIR.

mais je ne garantirais pas naïvement la "stabilité". en utilisant le langage du système de contrôle, le filtre TIIR n'est pas "complètement observable" et peut sembler stable parce que sa réponse impulsionnelle semble de longueur finie, mais à l'intérieur des états du filtre pourrait aller en enfer, et avec une précision numérique finie, cette instabilité interne finira par apparaître à la sortie.

nous devons nous désabuser de l'idée que "les filtres FIR n'ont pas de pôles" . ce n'est pas vrai.

robert bristow-johnson
la source
Pouvez-vous montrer mathématiquement que les filtres FIR ont des pôles, parce que je ne le vois pas.
Jim Clay
Le meilleur exemple de FIR à pôles est le filtre en cascade à peigne intégré (CIC). Il commence par un simple filtre de moyenne mobile (coefficients comme 1, 1, 1, 1) et le réécrit récursivement - introduisant ainsi un pôle. Voir lien . Celles-ci sont souvent implémentées sur des FPGA comme première étape de la conversion descendante car sous leur forme récursive, elles sont assez bon marché à implémenter par calcul. Voir la documentation Graychip à titre d'exemple. Ils sont généralement mis en œuvre en virgule fixe pour maintenir la stabilité.
David
1
Je pense que nous devrons accepter d'être en désaccord - le résumé de l'article original de Hogenauer se lit comme suit: "Une classe de filtres numériques à réponse impulsionnelle finie en phase linéaire (FIR) pour la décimation (diminution du taux d'échantillonnage) et l'interpolation (augmentation du taux d'échantillonnage) sont présentés."
David
4
NthN
2
@JimClay, un filtre CIC à somme mobile ou moyenne mobile est très certainement un filtre FIR. son IR est F. il n'est normalement pas implémenté comme un filtre FIR transversal, mais il pourrait certainement l'être si vous vouliez le payer avec MIPS.
robert bristow-johnson
14

"Pouvez-vous montrer mathématiquement que les filtres FIR ont des pôles, parce que je ne le vois pas." - Jim Clay

peut-on supposer que cette FIR est causale?

NN+1

la réponse impulsionnelle finie: h[n]=0n>N, n<0

fonction de transfert de la FIR:

H(z)=n=+h[n]zn=n=0Nh[n]zn=n=0NzNh[n]zNn=zNn=0Nh[Nn]zn=n=0Nh[Nn]znzN=h[N]+h[N1]z+h[N2]z2++h[1]zN1+h[0]zN(z0)N

il vous suffit de factoriser le numérateur et vous saurez où sont les zéros. mais il est assez évident où tous les pôles sont pour un filtre FIR. et il y a autant de pôles que l'ordre du filtre FIR. notez que ces pôles n'affectent pas la réponse en fréquence. sauf pour la phase.

robert bristow-johnson
la source
6
Je me suis trompé. Merci pour l'explication.
Jim Clay
1

Un peu par définition, en fait. Étant donné que vous entrez de l'énergie finie et que le filtre ne fournira au maximum qu'un multiple de l'énergie fournie (sa réponse impulsionnelle a une énergie finie), le signal résultant aura au maximum un multiple de l'énergie fournie. Il ne peut pas résonner et donc dégénérer, comme le peuvent les filtres IIR. C'est aussi derrière la réponse de Kenneides.

user7358
la source
oui, et c'est aussi faux que la réponse de Kenneide.
robert bristow-johnson
2
Cher Robert Bristow-Johnson, veuillez nous éclairer, mortels. D'où vient le filtre FIRH(z)=1avoir un poteau?
user7358
2
ok, bon point là-bas. le filtre FIR ou IIR d'ordre 0 (parfois appelé "fil") n'a pas de pôles ni de zéros (sauf si vous voulez penser àH(z)=1=zz où le pôle et le zéro s'annulent). Je me suis trompé.
robert bristow-johnson
L'unité tarde-t-elle H(z)=zavoir un poteau?
user7358
1
eh bien, c'est une avance d' unité . mais le délai unitaire,H(z)=z-1 n'ont un seul pôle àz=0.
robert bristow-johnson
1

Personne n'a vraiment expliqué pourquoi les pôles d'un filtre FIR sont amovibles, j'ai donc essayé de répondre à cette question ci-dessous.

Les filtres FIR auront des pôles amovibles à l'origine, car la limitation de leur réponse impulsionnelle l'exige. C'est autour du pôle, il est possible de définir la fonction de façon à ce qu'elle soit toujours holomorphe (différentiable en tout point de son domaine).

C'est un théorème de Riemann que si un signal est différenciable en chaque point de son domaine (sauf pour un nombre fini de points), alors il existe un voisinage autour de ces points spéciaux où la fonction est bornée. Les implications sont à double sens dans ce théorème, donc parce que les filtres FIR doivent avoir une réponse impulsionnelle bornée, alors la réponse impulsionnelle doit être différenciable en tout point du cercle unitaire. Ainsi, le signal peut être étendu de manière cohérente afin qu'il n'y ait pas de singularités (c'est-à-dire que les pôles sont amovibles).

C'est pourquoi la transformée en z d'un filtre FIR ne contient aucun pouvoir négatif de z.

Tom Kealy
la source
1
Tom, je pense que la transformée en Z d'une FIR réalisable ne contient que des pouvoirs négatifs dez. eh bien, seuls les pouvoirs non positifs dez.
robert bristow-johnson
@ robertbristow-johnson désolé, oui. Je pensais à une fonction génératrice. Cependant, je ne pense pas que la réponse ci-dessus change sous l'action dez -> z-1.
Tom Kealy