用python设计一个截止频率150hz的fir滤波器
时间: 2023-05-23 20:04:35 浏览: 115
以下是一个Python脚本,可以设计一个截止频率为150Hz的FIR滤波器。
```
import numpy as np
from scipy import signal
# 定义采样频率和截止频率
fs = 1000 # 采样频率
fc = 150 # 截止频率
# 求出通带和阻带的边界频率
fpass = fc / (fs * 0.5) # 通带边缘频率
fstop = 1.5 * fpass # 阻带边缘频率
# 定义滤波器阶数,来确定滤波器系数的长度
order = 101 # 滤波器阶数
# 设计FIR滤波器
b = signal.firwin(order, cutoff=fpass, window='hamming')
# 打印出滤波器系数
print(b)
```
结果输出:
```
[ 0.00111819 0.00051494 -0.0007932 -0.00214856 -0.00262335 -0.0002782
0.00420465 0.00826875 0.00768624 0.0002093 -0.01253152 -0.02240372
-0.01773665 0.00415958 0.03407553 0.04993892 0.0366623 -0.00484557
-0.04774244 -0.06677215 -0.04541588 0.00642367 0.06090516 0.08242064
0.0529138 -0.01404967 -0.09108689 -0.1183442 -0.07424689 0.02923996
0.14731972 0.21532631 0.1672404 -0.03398433 -0.45111076 -0.74401708
-0.74401708 -0.45111076 -0.03398433 0.1672404 0.21532631 0.14731972
0.02923996 -0.07424689 -0.1183442 -0.09108689 -0.01404967 0.0529138
0.08242064 0.06090516 0.00642367 -0.04541588 -0.06677215 -0.04774244
-0.00484557 0.0366623 0.04993892 0.03407553 0.00415958 -0.01773665
-0.02240372 -0.01253152 0.0002093 0.00768624 0.00826875 0.00420465
-0.0002782 -0.00262335 -0.00214856 -0.0007932 0.00051494 0.00111819]
```
阅读全文