Python移植MATLAB的巴特沃斯滤波器，椭圆滤波器

瀛海的鲲

已于 2022-11-18 17:22:52 修改

阅读量1.2k

点赞数 2

分类专栏： Python 文章标签： matlab python 开发语言

于 2022-11-18 17:10:52 首次发布

本文链接：https://blog.csdn.net/X888885304/article/details/127925907

版权

本文介绍了如何将MATLAB的巴特沃斯滤波器和椭圆滤波器移植到Python环境中，分别展示了MATLAB的原始示例和代码，以及Python的实现代码，但未涉及滤波器的详细原理。

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

1、MATLAB示例

2、MATLAB代码历程

MATLAB 巴特沃斯滤波器

clear all;
fp=3500; fs=5500; Rp=2.5; Rs=25;
[n,fn] =buttord(fp,fs,Rp,Rs,'s');  
Wn=2*pi*fn;
[b, a] =butter (n,Wn,'s');  
f=0:100:10000;
s=j*2*pi*f;
Hs=polyval(b,s)./polyval(a,s); 
figure(1);
subplot (2,1,1);
plot(f, 20*log10(abs(Hs)));
axis([0 10000 -40 1]);
xlabel('频率Hz') ; ylabel('幅度dB');
grid on;
subplot (2,1,2); plot(f,angle(Hs));
xlabel('频率Hz'); ylabel('相角rad');
disp('滤波器阶数和截止频率:')
grid on;
n, fn, Wn, b, a