matlab设计模拟带通滤波器

matlab设计模拟带通滤波器简单记录下在matlab上如何设计出模拟的带通滤波器,包括:巴特沃斯滤波器、切比雪夫I型滤波器、切比雪夫II型滤波器、椭圆型滤波器。代码如下:%设计带通滤波器%巴特沃斯、切比雪夫I型、切比雪夫II型

大家好,又见面了,我是你们的朋友全栈君。

简单记录下在matlab上如何设计出模拟的带通滤波器,包括:巴特沃斯滤波器、切比雪夫I型滤波器、切比雪夫II型滤波器、椭圆型滤波器。
代码如下:

%设计带通滤波器
%巴特沃斯、切比雪夫I型、切比雪夫II型、椭圆型滤波器

clear all;

%wp和ws分别是通带和阻带的频率(截止频率)。当wp和ws为二元矢量时,为带通或带阻滤波器,这时求出的Wn也是二元矢量;当wp和ws为一元矢量时,为低通或高通滤波器:当wp<ws时为低通滤波器,当wp>ws时为高通滤波器。

%wp和ws为二元矢量
wp=[0.1*2*pi 0.15*2*pi]; %设置通带频率
ws=[0.05*2*pi 0.2*2*pi]; %设置阻带频率

Rp=1; %设置通带波纹系数
Rs=20; %设置阻带波纹系数

%巴特沃斯滤波器设计
[N,Wn]=buttord(wp,ws,Rp,Rs,’s’); %求巴特沃斯滤波器阶数,输出参数N代表满足设计要求的滤波器的最小阶数,Wn是等效低通滤波器的截止频率
%无论是高通、带通和带阻滤波器,在设计中最终都等效于一个截止频率为Wn的低通滤波器(我现在也不是很理解为啥是这样,毕竟我也是刚接触滤波器)
fprintf(‘巴特沃斯滤波器 N= %4d\n’,N); %显示滤波器阶数
[bb,ab]=butter(N,Wn,’s’); %求巴特沃斯滤波器系数,即求传输函数的分子和分母的系数向量
W=0:0.01:2; %设置模拟频率
[Hb,wb]=freqs(bb,ab,W); %求巴特沃斯滤波器频率响应
plot(wb/pi,20*log10(abs(Hb)),’b’); %作图
hold on

%切比雪夫I型滤波器设计
[N,Wn]=cheb1ord(wp,ws,Rp,Rs,’s’); %求切比雪夫I型滤波器阶数
fprintf(‘切比雪夫I型滤波器 N= %4d\n’,N); %显示滤波器阶数
[bc1,ac1]=cheby1(N,Rp,Wn,’s’); %求切比雪夫I型滤波器系数,即求传输函数的分子和分母的系数向量
[Hc1,wc1]=freqs(bc1,ac1,W); %求切比雪夫I型滤波器频率响应
plot(wc1/pi,20*log10(abs(Hc1)),’k’); %作图

%切比雪夫II型滤波器设计
[N,Wn]=cheb2ord(wp,ws,Rp,Rs,’s’); %求切比雪夫II型滤波器阶数
fprintf(‘切比雪夫II型滤波器 N= %4d\n’,N);%显示滤波器阶数
[bc2,ac2]=cheby2(N,Rs,Wn,’s’); %求切比雪夫II型滤波器系数,即求传输函数的分子和分母的系数向量
[Hc2,wc2]=freqs(bc2,ac2,W); %求切比雪夫II型滤波器频率响应
plot(wc2/pi,20*log10(abs(Hc2)),’r’); %作图

%椭圆型滤波器设计
[N,Wn]=ellipord(wp,ws,Rp,Rs,’s’); %求椭圆型滤波器阶数
fprintf(‘椭圆型滤波器 N= %4d\n’,N); %显示滤波器阶数
[be,ae]=ellip(N,Rp,Rs,Wn,’s’); %求椭圆型滤波器系数,即求传输函数的分子和分母的系数向量
[He,we]=freqs(be,ae,W); %求椭圆型滤波器频率响应
%作图
plot(we/pi,20*log10(abs(He)),’g’);
axis([0 max(we/pi) -30 2]);
legend(‘巴特沃斯滤波器’,’切比雪夫I型滤波器’,’切比雪夫II型滤波器’,’椭圆型滤波器’);
xlabel(‘角频率{\omega}/{\pi}’);
ylabel(‘幅值/dB’);
line([0 max(we/pi)],[-20 -20],’color’,’k’,’linestyle’,’–‘);%在画布上画线
line([0 max(we/pi)],[-1 -1],’color’,’k’,’linestyle’,’–‘);
line([0.2 0.2],[-30 2],’color’,’k’,’linestyle’,’–‘);
line([0.3 0.3],[-30 2],’color’,’k’,’linestyle’,’–‘);

运行结果如下:

巴特沃斯滤波器 N= 4
切比雪夫I型滤波器 N= 3
切比雪夫II型滤波器 N= 3
椭圆型滤波器 N= 2

<span role="heading" aria-level="2">matlab设计模拟带通滤波器

 

 注:
1、求各种滤波器的传输函数的分子和分母向量系数的函数,例如:

[bc1,ac1]=cheby1(N,Rp,Wn,’s’); %求切比雪夫I型滤波器系数,即求传输函数的分子和分母的系数向量

这里的函数cheby1()中的参数为

cheby1(N,Rp,Wn,’type’,’s’);

其中注意’type’可为高通:’high’或带阻:‘stop’,如果不填的话,默认是低通或者带通,至于是低通还是带通,主要看你前面的wp和ws是二元矢量还是一元矢量:二元矢量对应带通,一元矢量对应低通。

2、matlab中的line()函数简单介绍:
https://blog.csdn.net/qq_36554582/article/details/83352155
————————————————
版权声明:本文为CSDN博主「逸凌123」的原创文章,遵循 CC 4.0 BY-SA 版权协议,转载请附上原文出处链接及本声明。
原文链接:https://blog.csdn.net/qq_36554582/article/details/83350837

版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 举报,一经查实,本站将立刻删除。

发布者:全栈程序员-用户IM,转载请注明出处:https://javaforall.cn/155297.html原文链接:https://javaforall.cn

【正版授权,激活自己账号】: Jetbrains全家桶Ide使用,1年售后保障,每天仅需1毛

【官方授权 正版激活】: 官方授权 正版激活 支持Jetbrains家族下所有IDE 使用个人JB账号...

(0)
blank

相关推荐

  • wireshark的安装教程(暖气片安装方法图解)

    百度Wireshark后,下载,进行安装。开始安装,许可,安装选项,附加任务,安装位置,必须安装WinPcap,才能捕获数据包,USBPcap,开始安装,安装WinPcap,协议,驱动,安装NPF驱动时出错,Wireshark必须有NPF驱动才能工

  • Socket 非阻塞模式下connect 返回EINPROGRESS(115)错误[通俗易懂]

    Socket 非阻塞模式下connect 返回EINPROGRESS(115)错误[通俗易懂]今天再测试socket的时候,发现一个很奇怪的问题,就是客户端再connect的时候第一次connect总是会返回-1,errno是115,往往第二次连接就可以成功了。但是对于服务端来说,第一次连接已经成功返回了。后来想想可能跟自己的设置socket是非阻塞的有关系,后来吧socket设置成阻塞的,问题确实就没有了。后来有反复尝试了非阻塞的。我先把服务器关闭,让客户端连接,可以发现从打出来的e…

  • matlab函数plot函数_动态变量

    matlab函数plot函数_动态变量引子对于真实系统或者仿真平台,数据是增量式的产生的。Matlab除了强大的矩阵运算外,还具有强大的数据可视化库。由于静态画图的方法较多,本文只针对增量式数据流的动态显示。本文主要介绍几种Matlab实现数据的动态显示方法。方法主要有两种:holdonset函数与drawnow函数组合holdon方法1.方法介绍此种方法比较原始,适合于即时数据,原理是先画上一帧,接着保留原始图像,

    2022年10月10日
  • 服务自启动设置_服务器设置启动项

    服务自启动设置_服务器设置启动项一、linux下的服务分类RPM包默认安装的服务这些服务是通过RPM包安装的,可以被服务管理命令识别。又分为两种子分类:1)独立的服务:就是独立启动的意思,这类型的服务可以自行启动,而不

  • voliate关键字的作用[通俗易懂]

    voliate关键字的作用[通俗易懂]一、内存可见性基于缓存一致性协议,当用voliate关键字修饰的变量改动时,cpu会通知其他线程,缓存已被修改,需要更新缓存。这样每个线程都能获取到最新的变量值。二、基于内存屏障的防止指令重排用voliate修饰的变量,可以防止cpu指令重排序。底层的实现方式是基于4种内存屏障:读读、读写、写读、读读屏障。…

  • netty通信框架_王国风云2控制台代码

    netty通信框架_王国风云2控制台代码写在前面所属章节链接第二章2.1BIO通信aaaTimeServer代码importjava.io.IOException;importjava.net.ServerSocket;importjava.net.Socket;publicclassTimeServer{ publicstaticvoidmain(String[]args)throwsIOException{ intport=8081; if(args

发表回复

您的电子邮箱地址不会被公开。

关注全栈程序员社区公众号