实验五 基于MATLAB的模拟滤波器设计

已于 2022-03-25 18:46:30 修改
阅读量2w 收藏 243
点赞数 42
分类专栏： 数字信号处理 文章标签： matlab 巴特沃斯模拟滤波器 切比雪夫型滤波器 椭圆滤波器 模拟域频率变换
于 2021-06-09 13:27:23 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_45696377/article/details/117739175
版权

目录

一、实验目的：

1．加深对模拟滤波器常用指标的理解；
2．学会模拟滤波器中的频率变换;
3．根据指标要求设计模拟滤波器，并进行信号的处理。

二、实验原理：

模拟滤波器的设计是其它滤波器设计的基础，其设计原理见课本。利用MATLAB设计模拟滤波器的调用函数见课本。常用的函数格式如下：

1.巴特沃斯模拟滤波器

[N,wc]=buttord(wp,ws,Ap,As,’s’)
[b,a]=butter(N,wc,’ftype’,’s’)

函数buttord的输入参数wp和ws（rad/s）分别表示滤波器的通带和阻带截频，Ap和As（dB）表示滤波器的通带和阻带衰减。‘s’表示所设计的是模拟滤波器。函数buttord返回的参数N是滤波器的阶数，wc为BW型滤波器3dB的截频。
利用函数butter获得BW型滤波器的分子多项式系数（b）和分母多项式系数（a）。
‘ftype’表示要设计滤波器的型号，当设计的为低通或带通滤波器是这个参数可以缺省。当ftype=high时，函数butter获得BW型高通滤波器的分子多项式系数（b）和分母多项式系数（a）；当当ftype=stop时，函数butter获得BW型带阻滤波器的分子多项式系数（b）和分母多项式系数（a）；
或者：

N=buttord(wp,ws,Ap,As,’s’); wc=wp/(10^(0.1*Ap)-1)^(1/2/N);  [b,a]=butter(N,wc, ’ftype’,’s’);

或者：

N=buttord(wp,ws,Ap,As,’s’); wc=ws/(10^(0.1*As)-1)^(1/2/N);  [b,a]=butter(N,wc, ’ftype’,’s’);

2.切比雪夫Ⅰ型滤波器

[N,wc]=cheb1ord(wp,ws,Ap,As,’s’);
[b,a]=cheby1(N,Ap,wc, ’ftype’,’s’);
或者，
N=cheb1ord(wp,ws,Ap,As,’s’);wc=wp; [b,a]=cheby1(N,Ap,wc, ’ftype’,’s’)
或者
N=cheb1ord(wp,ws,Ap,As,’s’); [b,a]=cheby1(N,Ap,wp, ’ftype’,’s’);

3. 切比雪夫Ⅱ型滤波器

[N,wc]=cheb2ord(wp,ws,Ap,As,’s’)
[b,a]=cheby2(N,As,wc, ’ftype’,’s’)

或者

N=cheb2ord(wp,ws,Ap,As,’s’);wc=ws; [b,a]=cheby2(N,As,wc, ’ftype’,’s’)

或者

N=cheb2ord(wp,ws,Ap,As,’s’); [b,a]=cheby2(N,As,ws, ’ftype’,’s’)

4.椭圆滤波器

[N,wc]=ellipord(wp,ws,Ap,As,’s’);
[b,a]=ellip(N,Ap,As,wc, ’ftype’,’s’);

或者

N= ellipord (wp,ws,Ap,As,’s’);wc=wp; [b,a]=ellip(N,Ap,As,wc, ’ftype’,’s’);

或者

N= ellipord (wp,ws,Ap,As,’s’); [b,a]=ellip(N,Ap,As,wp, ’ftype’,’s’);

5.模拟域频率变换

MATLAB提供了实现模拟域频率变换的函数，他们分别是：

（1）模拟低通到高通变换

[numt,dent] = lp2hp(num,den,W0)

高通中，一般W0=1。

（2）模拟低通到带通变换

 [numt,dent] = lp2bp(num,den,W0,B)

（3）模拟低通到带阻变换

[numt,dent] = lp2bs(num,den,W0,B)

其中，num，den分别表示变换前模拟低通滤波器系统函数的分子多项式系数和分母多项式系数，W0和B为变换中的参数。numt，dent分别表示变换后模拟滤波器系统函数的分子多项式系数和分母多项式系数，

6.补充函数

（1）buttap函数用来返回设计的BW滤波器的零点、极点和增益。函数调用格式为：

[z,p,k]=buttap(N);

其中，z,p,k分别为滤波器系统函数H（S）的零点、极点和增益，N为滤波器的阶数。
（2）cheb1ap函数用来返回设计的CB Ⅰ型滤波器的零点、极点和增益，cheb2ap函数用来返回设计的CB Ⅱ滤波器的零点、极点和增益。其调用格式为：

[z,p,k]= cheb1ap (N，Rp);
[z,p,k]= cheb2ap (N，Rs);

其中，z,p,k分别为滤波器系统函数H（S）的零点、极点和增益，N为滤波器的阶数，Rp为滤波器在通带内的最大衰减值，Rs为滤波器在阻带内的最小衰减值。
（3）zp2tf函数，函数调用格式为：

[b,a]=zp2tf (z,p,k);

其中，z,p,k分别为滤波器系统函数H（S）的零点、极点和增益，b、a分别为滤波器系统函数H（s）的分子和分母多项式的系数，
（4）freqs函数用来求解模拟滤波器的频率响应。其函数调用格式为：

h=freqs（b,a,w）;
[h,w]=freqs（b,a,n）;

其中，b，a分别为滤波器系统函数H（s）的分子和分母多项式的系数，w表示频率点，n表示进行复频率响应的点数，其默认值为512。如果没有输出参数而直接调用freqs（b,a,w函数时，将直接画出复频率响应的幅频响应和相频响应。
freqs函数用来求解模拟滤波器的频率响应。

三、例题：

1．BW模拟低通滤波器的幅频特性曲线
程序如下：

clear all
for i=1:4  
switch i   
    case 1   
        N=2;    
    case 2     
        N=5;   
    case 3    
        N=10;   
    case 4    
        N=20;  
end  
[z,p,k]=buttap(N); 
 [b,a]=zp2tf(z,p,k); 
 [H,w]=freqs(b,a); 
 magH2=(abs(H)).^2;  
hold on 
plot(w,magH2);  
axis([0 2 0 1.1]);  
end  
xlabel('w/wc');  ylabel('/H(jw)/^2')

在这里插入图片描述

2、（1）设计一个巴特沃斯模拟低通滤波器，技术指标为：通带截频1000Hz，阻带截频1500Hz，通带波纹1dB，阻带衰减50dB。
（2）假设一个信号x（t）=sin(2pif1t)+sin(2pif2t)+sin(2pif3*t),其中f1=100Hz，f2=2000Hz，f3=2900Hz，信号的采样频率为10000Hz。将原信号与通过该滤波器的模拟信号进行比较。

clear all
wp=1000 *2*pi;
ws=1500*2*pi;
Ap=1;
As=50;
[N,wc]=buttord(wp,ws,Ap,As,'s'); %求的滤波器的最小阶数和截止频率
[b,a]=butter(N,wc,'s');           %设计模拟巴特沃斯滤波器
w=linspace(0,4000,1000)*2*pi; %设置绘制频率响应的频率点
H=freqs(b,a,w);      % 计算给定频率点的复数频率响应
figure(1)
plot(w/2/pi, 20*log10(abs(H)));
xlabel('频率/Hz');ylabel('振幅/dB');
grid on;
% 计算Ap,As

w1=[wp ws];
h=freqs(b,a,w1);
fprintf('Ap=%.4f\n',-20*log10(abs(h(1))));
fprintf('As=%.4f\n',-20*log10(abs(h(2))));

fs=10000;
dt=1/fs;   %模拟信号采样间隔
f1=100;f2=2000;f3=2900;
t=0:dt:0.1;
x=sin(2*pi*f1*t)+sin(2*pi*f2*t)+sin(2*pi*f3*t);
H=[tf(b,a)];  %滤波器在MATLAB中的表示
y=lsim(H,x,t); % 模拟输出
figure(2)
subplot(2,1,1)
plot(t,x);
xlabel ('(a)输入信号');
subplot(2,1,2)
plot(t,y);
xlabel('(b)输出信号');

在这里插入图片描述

		图1  模拟低通滤波器的幅度响应函数

在这里插入图片描述

			  图2 滤波器的输入输出函数

Ap=0.9053； As=50.0000
2.设计一个巴特沃斯模拟带通滤波器，技术指标为：通带频率：1000-2000Hz，两侧过渡带宽500Hz，通带波纹1dB，阻带衰减100dB。
方法一：

clear all
wp1=1000*2*pi;
wp2=2000*2*pi;
ws1=500*2*pi;
ws2=2500*2*pi;
B=wp2-wp1;
W0=sqrt(wp1*wp2);
Ap=1;
As=100;
ws11=(ws1^2-W0^2)/(B*ws1);     %带通到低通的频率变换
ws22=(ws2^2-W0^2)/(B*ws2);
wp=1;
ws=min(abs(ws11),abs(ws22));    %低通滤波器的阻带截频
[N,wc]=buttord(wp,ws,Ap,As,'s'); %求的滤波器的最小阶数和截止频率
[b,a]=butter(N,wc,'s');           %设计模拟巴特沃斯滤波器
[c,d]=lp2bp(b,a,W0,B);    %将低通变换为带通
w=linspace(1,3000,1000)*2*pi; %设置绘制频率响应的频率点
H=freqs(c,d,w);      % 计算给定频率点的复数频率响应
plot(w/2/pi,abs(H));
xlabel('频率/Hz');ylabel('振幅/dB');
grid on;
% 计算Ap,As
w1=[wp1 wp2 ws1 ws2];
h=freqs(c,d,w1);
fprintf('Ap1=%.4f\n',-20*log10(abs(h(1))));
fprintf('Ap2=%.4f\n',-20*log10(abs(h(2))));
fprintf('As1=%.4f\n',-20*log10(abs(h(3))));
fprintf('As1=%.4f\n',-20*log10(abs(h(4))));

运行结果：
在这里插入图片描述

			图3 带通滤波器的幅度响应函数

Ap1=0.9713；Ap2=0.9716；As1=244.2668；As1=100.0000

方法二：

程序如下：
clear all
wp=[1000 2000]*2*pi;
ws=[500 2500]*2*pi;
Ap=1;
As=100;
[N,wc]=buttord(wp,ws,Ap,As,'s'); %求的滤波器的最小阶数和截止频率
[b,a]=butter(N,wc,'s');           %设计模拟巴特沃斯滤波器
w=linspace(1,3000,1000)*2*pi; %设置绘制频率响应的频率点
H=freqs(b,a,w);      % 计算给定频率点的复数频率响应
plot(w/2/pi,abs(H));
xlabel('频率/Hz');ylabel('振幅/dB');
grid on;
% 计算Ap,As
w1=[wp ws];
h=freqs(b,a,w1);
fprintf('Ap1=%.4f\n',-20*log10(abs(h(1))));
fprintf('Ap2=%.4f\n',-20*log10(abs(h(2))));
fprintf('As1=%.4f\n',-20*log10(abs(h(3))));
fprintf('As1=%.4f\n',-20*log10(abs(h(4))));

在这里插入图片描述

			图3 带通滤波器的幅度响应函数

Ap1=0.9724；Ap2=0.9721；As1=244.2646；As1=99.9998

四、作业：

1.假设一个信号x（t）= sin(2pif1t)+0.5cos(2pif2t)，其中f1=20Hz，f2=100Hz。
（1）请设计一个模拟滤波器能把f2滤除掉，请写出程序，并画出原信号与原信号通过滤波器的输出信号的图形。

clear all
wp=52*2*pi;
ws=80*2*pi;
Ap=1;
As=50;
[N,wc]=buttord(wp,ws,Ap,As,'s'); %求的滤波器的最小阶数和截止频率
[b,a]=butter(N,wc,'s');           %设计模拟巴特沃斯滤波器
w=linspace(0,500,1000)*2*pi; %设置绘制频率响应的频率点
H=freqs(b,a,w);      % 计算给定频率点的复数频率响应
figure(1)
plot(w/2/pi, 20*log10(abs(H)));
xlabel('频率/Hz');ylabel('振幅/dB');
grid on;
% 计算Ap,As

w1=[wp ws];
h=freqs(b,a,w1);
fprintf('Ap=%.4f\n',-20*log10(abs(h(1))));
fprintf('As=%.4f\n',-20*log10(abs(h(2))));

fs=10000;
dt=1/fs;   %模拟信号采样间隔
f1=20;f2=100;
t=0:dt:0.1;
x=sin(2*pi*f1*t)+0.5*cos(2*pi*f2*t);
H=[tf(b,a)];  %滤波器在MATLAB中的表示
y=lsim(H,x,t); % 模拟输出
figure(2)
subplot(2,1,1)
plot(t,x);
xlabel ('(a)输入信号');
subplot(2,1,2)
plot(t,y);
xlabel('(b)输出信号');

在这里插入图片描述

（2）请设计一个模拟滤波器能把f1滤除掉，请写出程序，并画出原信号与原信号通过滤波器的输出信号的图形。
在这里插入图片描述

在这里插入图片描述
在这里插入图片描述

2.假设一个信号x（t）=sin(2pif1t)+sin(2pif2t)+sin(2pif3*t),其中f1=200Hz，f2=1500Hz，f3=2900Hz，信号的采样频率为10000Hz。
（1）请设计一个模拟滤波器将f2保留，f1和f3滤除。

clear all
wp1=1000*2*pi;
wp2=2000*2*pi;
ws1=500*2*pi;
ws2=2500*2*pi;
B=wp2-wp1;
W0=sqrt(wp1*wp2);
Ap=1;
As=100;
ws11=(ws1^2-W0^2)/(B*ws1);     %带通到低通的频率变换
ws22=(ws2^2-W0^2)/(B*ws2);
wp=1;
ws=min(abs(ws11),abs(ws22));    %低通滤波器的阻带截频
[N,wc]=buttord(wp,ws,Ap,As,'s'); %求的滤波器的最小阶数和截止频率
[b,a]=butter(N,wc,'s');           %设计模拟巴特沃斯滤波器
[c,d]=lp2bp(b,a,W0,B);    %将低通变换为带通
w=linspace(1,3000,1000)*2*pi; %设置绘制频率响应的频率点
H=freqs(c,d,w);      % 计算给定频率点的复数频率响应
plot(w/2/pi,abs(H));
xlabel('频率/Hz');ylabel('振幅/dB');
grid on;
% 计算Ap,As
w1=[wp1 wp2 ws1 ws2];
h=freqs(c,d,w1);
fprintf('Ap1=%.4f\n',-20*log10(abs(h(1))));
fprintf('Ap2=%.4f\n',-20*log10(abs(h(2))));
fprintf('As1=%.4f\n',-20*log10(abs(h(3))));
fprintf('As1=%.4f\n',-20*log10(abs(h(4))));

在这里插入图片描述

（2）请设计一个模拟滤波器将f2滤除，f1和f3保留。

更多相关文章点这里哦

数字信号处理----全套Matlab实验报告

通信工程(信息类，电子类，电气工程，自动化，计算机，软件工程，机电，等相关专业)全套学习指导

在这里插入图片描述
答疑资料qq群:1007576722

Gopher-毛
关注
  • 42
    点赞
  • 踩
  • 243
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 4
    评论
专栏目录
使用Matlab设计数字滤波器，从原理到代码
乐观的lishan的博客
12-05 8万+
本文是在Matlab下设计数字滤波器的一些总结，减少了理论的说明，重在代码的实现方法，给出了普遍使用的数字滤波器代码，修改参数即可移植使用。主要内容包括：滤波器的特性介绍、FIR、IIR数字滤波器的生成方法，包括冲激响应不变法和双线性映射法的实现。最后介绍了Matlab自带的滤波器生成工具的使用方法和扫频信号的生成，可以避开复杂的理论计算，直接设计需要的滤波器，提高工作效率。
模拟滤波器设计简介
02-23
模拟滤波器设计是模拟电子电路设计中很重要的一个部分。本文首先综述了模拟滤波器的一些基本内容，然后着重介绍了有源滤波器的设计方法、电路结构、性能评价和器件选择。
matlab 模拟滤波器,基于MATLAB的模拟滤波器设计
weixin_39897070的博客
03-21 540
第12卷 第3期2004年9月 山东交通学院学报JO UR NAL OF SHANDON G JIAOTONG UNIVERS ITY Vol.12No.3 Sep.2004基于M A TLAB 的模拟滤波器设计姜春玲1,史 玲2(1.泰山学院物理系,山东泰安 271021;2.山东交通学院信息工程系,山东济南 250023)摘要:模拟滤波器是现代测控系统中的重要部件,传统设计方法的设计过程繁琐。...
模拟滤波器的设计
weixin_43389179的博客
05-18 3889
参考书籍：模拟滤波器与电路设计手册 【美】Arthur B. Williams 常用的滤波器频率响应是可以用两个含s的负数多项式比来决定，其实s=jw（w为角频率） 极点：分母多项式的根 零点：分子多项式的根 频率响应的幅值为频率点到零点的矢量长度之积除以到极点的矢量长度之积；相位为频率点到零点的矢量相位和减去到极点的矢量相位和。 当频率远小于某零点/极点的模时，该零点/极点对频率响应的影响可以忽...
模拟滤波器工程设计
pp_0604的博客
12-08 422
使用matlab验证滤波器零、极点对幅频响应的影响
数字滤波器(三)--模拟滤波器的设计
魔法师の屋子★↓☆
06-01 1万+
模拟滤波器的设计引前言1. 模拟滤波器的系统函数2 巴特沃斯滤波器2.1 设计指标2.2 设计流程3. 切比雪夫Ⅰ型滤波器 引 数字滤波器(一)–IIR与FIR的基本结构与MATLAB实现 数字滤波器(二)–最小相位延时系统和全通系统 前言 IIR滤波也就是无限长单位冲击响应滤波器，基本结构和MATLAB实现可以参考本专栏第一篇博客：数字滤波器(一)–IIR与FIR的基本结构与MATLAB实现。 IIR滤波器的设计步骤主要由下图所示的三个大步骤组成。首先我们要设计出模拟原型滤波器，然后将滤波器数字化为相应的
基于MATLAB的滤波器的设计ppt文件
04-20
**基于MATLAB的滤波器设计** 在信号处理领域，滤波器扮演着至关重要的角色。它们被用于去除噪声，提取有用信号，以及改善信号质量。MATLAB（矩阵实验室）作为一个强大的数学计算软件，提供了丰富的工具箱来实现各种...
维纳滤波器的设计及Matlab仿真实现.doc
06-25
4. 仿真与验证：在Matlab等环境中，可以通过生成随机信号，模拟信道传输和噪声添加，然后应用设计的滤波器，最后比较输出信号y(n)与期望响应d(n)的相似度，验证滤波器的效果。 在给定的Matlab代码中，信号s(n)首先...
基于matlab数字滤波器的设计.pdf
12-12
MATLAB（矩阵实验室）是一种交互式的程序设计环境，专为数值计算和符号计算而设计，它支持多种科学和工程计算任务，包括数据分析、算法开发、用户界面设计等。 3.2 MATLAB产生的历史背景： MATLAB起源于20世纪70...
基于matlab的卡尔曼滤波器设计与实现
06-23
卡尔曼滤波器是一种在信号处理和控制理论领域广泛应用的估计方法，特别是在处理随机变化且存在噪声的数据时。...通过模拟实验和对比不同滤波效果，可以不断优化滤波器性能，以获得更精确的估计结果。
低通滤波器 matlab 实现
08-29
采用自带的buttord 和butter函数实现对信号的低通滤波
Analog Filters Using MATLAB（使用 MATLAB 进行模拟滤波）MATLAB 代码
06-05
使用MATLAB的模拟滤波器全面概述了频率选择模拟滤波器的理论和现代设计方法，并介绍了如何使用MATLAB选择和设计模拟滤波器。该材料涵盖具有集总电路和分布式电路元件的经典无源滤波器的合成和设计，以及基于运算放大器，跨导体和电流传输器的有源滤波器。
基于Matlab的FIR滤波器设计实验报告
06-08
基于Matlab的FIP滤波器的设计实验报告
滤波器设计实例及Matlab代码
12-10
数字滤波器设计原理公式及其实例，包括各种窗函数、带通滤波、高通、低通滤波；有限脉冲（FIR）滤波器
不同模拟滤波器设计matlab代码
11-04
模拟滤波器设计matlab代码，低通、高通、带通、带阻，巴特沃斯滤波器
matlab卡尔曼滤波器激光雷达SLAM算法设计
04-17
《基于Matlab的卡尔曼滤波器与激光雷达SLAM算法设计详解》 在自动驾驶领域，定位和建图（Simultaneous Localization and Mapping，简称SLAM）是核心技术之一，它使得车辆能够在未知环境中自主导航。其中，卡尔曼...
用matlab设计滤波器实验报告,数字信号出来实验报告--matlab滤波器设计
weixin_33941859的博客
03-21 1757
数字信号出来实验报告--matlab滤波器设计 广 西 工 学 院 实 验 报 告 用 纸 实验名称 IIR数字滤波器的设计 实验成绩 指导老师 陈艳 系(院) 计算机工程系 班级 学号 学生姓名 一、实验目的 1、掌握双线性变换法和脉冲响应不变法设计IIR数字滤波器的具体设计方法及原理，熟悉用双线性变换法及脉冲响应不变法设计低通、高通和带通IIR数字滤波器的计算机编程。 2、观察双线性变换及脉冲...
IIR数字滤波器设计（模拟原型低通滤波器设计）
最新发布
weixin_43388156的博客
04-18 2149
前面介绍了IIR滤波器的设计流程，设计IIR数字滤波器需要先设计相应的模拟滤波器，这里将对巴特沃斯滤波器、切比雪夫滤波器及椭圆滤波器三种模拟滤波器的设计方法进行介绍。ωp\omega_pωp​：通带截止频率；ωs\omega_sωs​：阻带介质频率；δp\delta_pδp​：通带波动；δs\delta_sδs​：阻带波动。
常用模拟低通滤波器的设计~经典 IIR 滤波器之巴特沃斯滤波器
热门推荐
yishuihanq的博客
05-29 1万+
巴特沃斯 (Butterworth) 滤波器又称为“最平”的幅频响应滤波器，该滤波器的特点是在通带内具有最大平坦的幅度且随频率升高呈单调减小。。。2.1、buttap 函数。。。2.2、buttord 函数。。。2.3、butter 函数。。。
iir数字滤波器设计matlab实验
04-06
B'IIR数字滤波器设计MATLAB实验是一种通过MATLAB软件进行IIR数字滤波器设计与模拟的实验。IIR数字滤波器是一种基于时域信号的滤波器，具有滤波器阶数低、滤波器品质高等优点。在实验过程中，可以通过MATLAB自带的滤波器设计工具箱，选择不同类型的滤波器、设置滤波器参数、进行滤波器性能分析等。通过这种实验可以更加深入地理解数字滤波器的设计与应用，为研究数字信号处理打下基础。
评论 4
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

Gopher-毛

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付：¥1
获取中
扫码支付

您的余额不足，请更换扫码支付或充值

打赏作者

实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值