Utilisez la fonction
filtfilt pour calculer la réponse à phase nulle d'un filtre, à l'aide d'un modèle de filtrage "ascendant-descendant".
1. Utilisez la fonction
butter pour obtenir les coefficients d'un filtre analogique Butterworth d'ordre 2.
2. Utilisez la fonction
iirlow pour obtenir les coefficients pour un filtre IIR lowpass avec fréquence de coupure de 0.25.
3. Utilisez les fonctions
gain et
arg pour définir la réponse d'intensité et de phase du filtre.
4. Redéfinissez la réponse d'amplitude et de phase du filtre pour toutes les valeurs f mais en termes de p.
5. Tracez la réponse d'amplitude, puis utilisez des marqueurs pour afficher son amplitude à la fréquence de coupure.
6. Tracez la réponse de phase, puis utilisez des marqueurs pour afficher sa valeur à la fréquence de coupure.
7. Définissez et tracez un signal d'impulsion d'échantillon.
L'entrée comprend une impulsion au point médian du signal.
8. Utilisez la fonction filtfilt pour calculer la sortie de phase nulle.
9. Utilisez la fonction
dft pour observer l'effet que le filtre a sur le signal en calculant les transformées de Fourier des deux signaux et tracez l'intensité et la phase de Y.
10. Utilisez la fonction mag pour montrer que la fonction de transfert, avec le taux des amplitudes de X et de Y, produit la forme de coupure Butterworth attendue.
11. Tracez la réponse d'amplitude, ainsi que la fonction de transfert.
12. Utilisez la fonction arg pour définir la différence de phase entre les arguments principaux de X et Y.
13. Tracez la réponse d'amplitude et la différence de phase entre X et Y.
• La fonction filtfilt a pour effet de doubler le filtrage de l'amplitude du signal, comme l'indique le tracé. Cependant, le décalage de phase est égal à zéro, contrairement à la phase non linéaire du filtre d'origine.
• Lors de l'interprétation du tracé de phase, n'oubliez pas que les phases de -2π et de 0 sont équivalentes. Lorsque l'amplitude de Y est très faible, par exemple pour les fréquences situées très au-dessus de la fréquence de coupure, la phase peut être insignifiante en raison du bruit numérique.