python filter design

最新推荐文章于 2024-06-30 17:30:52 发布

好清闲

最新推荐文章于 2024-06-30 17:30:52 发布

阅读量1.7k

点赞数 1

本文链接：https://blog.csdn.net/xingyu19871124/article/details/53787018

版权

该博客展示了如何使用Python中的scipy库设计和应用一个低通滤波器。通过butter函数计算滤波器系数，然后利用lfilter进行滤波，最后通过freqz绘制频率响应，并通过示例数据展示滤波前后的信号效果。

摘要由CSDN通过智能技术生成

import numpy as np;
import matplotlib.pyplot as plt
import cv2;
import math;
from cmath import sin
from numpy.core.numeric import dtype
from scipy.signal import butter, lfilter, freqz
from scipy.fftpack import fft;
'''    
A few comments:


The Nyquist frequency is half the sampling rate.
You are working with regularly sampled data, so you want a digital filter, not an analog filter. This means you should not use analog=True in the call to butter, and you should use scipy.signal.freqz (not freqs) to generate the frequency response.
One goal of those short utility functions is to allow you to leave all your frequencies expressed in Hz. You shouldn't have to convert to rad/sec. As long as you express your frequencies with consistent units, the scaling in the utility functions takes care of the normalization for you.
Here's my modified version