首页 > 其他分享 >数字信号处理实验三:IIR数字滤波器设计及软件实现

数字信号处理实验三:IIR数字滤波器设计及软件实现

时间:2024-05-28 20:04:57浏览次数:33  
标签:滤波器 Fs 数字 IIR apos pi 阻带

一、实验目的

1. 掌握MATLAB中进行IIR模拟滤波器的设计的相关函数的应用;

2. 掌握MATLAB的工具箱中提供的常用IIR数字滤波器的设计函数的应用;

3.掌握MATLAB的工具箱中提供的模拟滤波器转数字滤波器的相关的设计函数的应用。

二、实验内容

本实验为综合性实验项目,要求通过利用MATLAB软件编写程序,设计IIR数字滤波器,并绘制波形。

三、实验器材(设备、元器件)

安装有MATLAB软件的计算机1台

四、实验步骤

1、设计一个工作于采样频率80 kHz的巴特沃斯低通数字滤波器,要求通带边界频率为4 kHz,通带最大衰减为0.5 dB,阻带边界频率为20 kHz,阻带最小衰减为45 dB。直接调用MATLAB工具箱函数buttord和butter设计数字滤波器,并显示数字滤波器系统函数H(z)的系数。

2、设计一个工作于采样频率80 kHz的切比雪夫Ⅰ型低通数字滤波器,要求通带边界频率为4 kHz,通带最大衰减为0.5 dB,阻带边界频率为20 kHz,阻带最小衰减为 45 dB。 直接调用MATLAB工具箱函数cheb1ord和cheby1设计数字滤波器,并显示数字滤波器系统函数H(z)的系数。

3、用脉冲响应不变法设计一个巴特沃斯低通数字滤波器,指标要求采样频率80 kHz,通带边界频率为4 kHz,通带最大衰减为0.5 dB,阻带边界频率为20 kHz,阻带最小衰减为45 dB。编写程序先调用MATLAB工具箱函数buttord和butter设计过渡模拟低通滤波器,再调用脉冲响应不变法数字化转换函数impinvar, 将过渡模拟低通滤波器转换成低通数字滤波器H(z)。

4、用双线性变换法(bilinear函数)设计符合以下指标的数字低通滤波器:通带内角频率低于0.2π rad时,容许幅度误差在1dB以内;在角频率0.3π到π之间的阻带衰减大于15dB,采样周期为0.5s。

五、实验程序及结果分析

1. 设计巴特沃斯低通数字滤波器代码如下:

close all;

clear all;

% 设置参数

Fs = 80000 ;% 采样频率

T = 1/Fs; %采样间隔

wp = 2 * pi * 4000/Fs; %%通带模拟角频率

ws = 2 * pi * 20000/Fs;%%阻带模拟角频率

rp = 0.5;

rs = 45;


% 进行设计

 [N,wc] =  buttord(wp/pi,ws/pi,rp,rs);

 [B,A] = butter(N,wc)

% [Bz,Az] = impinvar(B,A,Fs)

%  绘图

figure;

subplot(2,2,1);

[H,w]= freqz(B,A,1000);

plot(w/pi,20*log10(abs(H)));

title('损耗函数曲线');

xlabel('\omega/\pi');

ylabel('幅度(dB)');

subplot(2,2,2);

plot(w/pi,angle(H));

title('相频特性曲线');

xlabel('\omega/\pi');

ylabel('相位');

axis([0 1 -2*pi 2*pi]);

实验结果如下:

2.设计切比雪夫Ⅰ型低通数字滤波器代码如下:

close all;

clear all;

Fs = 80000; %Hz 采样频率

T = 1/Fs;

wp = 2 * pi * 4000/Fs; %%通带模拟角频率

ws = 2 * pi * 20000/Fs;%%阻带模拟角频率

rp = 0.5;

rs = 45;

%获取阶数和截止频率

[ N,wc ] = cheb1ord(wp/pi, ws/pi, rp , rs);

%获得转移函数系数

[ B,A ] = cheby1(N,rp,wc/pi,'low')

%滤波

figure;

subplot(2,2,1);

[H,w]= freqz(B,A,1000);

plot(w/pi,20*log10(abs(H)));

title('损耗函数曲线');

xlabel('\omega/\pi');

ylabel('幅度(dB)');

subplot(2,2,2);

plot(w/pi,angle(H));

title('相频特性曲线');

xlabel('\omega/\pi');

ylabel('相位');

axis([0 1 -2*pi 2*pi]);

实验结果如下:

 

3.设计巴特沃斯低通数字滤波器代码如下:

close all;

clear all;

%1.数字滤波器的技术指标

rp = 0.5;%通带最大衰减

rs = 45;%阻带最小衰减

fp = 4000;%通带截止频率

fs = 20000;%阻带截止频率

Fs = 80000;%抽样间隔

T = 1/Fs;

%2.将数字指标转化成模拟滤波器技术指标(采用冲激响应不变法)

wap = 2*pi*fp;

was = 2*pi*fs;

%3.设计滤波器

[N,wac] = buttord(wap,was,rp,rs,'s');%计算阶数N和3dB截止频率wac

[z,p,k] = buttap(N);%创建巴特沃斯低通滤波器 z零点p极点k增益

[Bap,Aap] = zp2tf(z,p,k);%由零极点和增益确定归一化Han(s)系数

[Bbs,Abs] = lp2lp(Bap,Aap,wac);%将s/wc 代替s,去归一化

[B,A] = impinvar(Bbs,Abs,Fs);%模拟域到数字域——冲激响应不变法

[H1,w] = freqz(B,A);%根据H(z)求频率响应

%绘制数字滤波器频响幅度谱

figure(1);

f = w*Fs/(2*pi);

subplot(2,1,1);

plot(f,20*log10(abs(H1)));%绘制幅度响应

title('冲激响应不变法——巴特沃斯BLPF(幅度)');

xlabel('频率/Hz');

ylabel('H1幅值/dB');

subplot(2,1,2);

plot(f,unwrap(angle(H1)));% 绘制相位响应

xlabel('频率/Hz');

ylabel('角度/Rad');

title('冲激响应不变法——巴特沃斯BLPF(相位)'); 

实验结果如下:  

4.用双线性变换法设计数字低通滤波器代码如下:

clear all;

wp = 0.2*pi;

ws = 0.3*pi;

rp = 1;

rs = 15;

Fs = 1;

wp1=2*Fs*tan(wp/2);

ws1=2*Fs*tan(ws/2);

[N,Wn] = buttord(wp1,ws1,rp,rs,'s');

[Z,P,K] = buttap(N);

[Bap,Aap] = zp2tf(Z,P,K);

[b,a] = lp2lp(Bap,Aap,Wn);

[bz,az] = bilinear(b,a,Fs);

[H,W] = freqz(bz,az);

disp(bz);

disp(az);

subplot(2,1,1);

plot(W*Fs/pi,abs(H));

grid on;

xlabel('频率/Hz');

ylabel('幅度');

title('(a)');

subplot(2,1,2);

plot(W/pi,20*log10(abs(H)));

grid on;

xlabel('\omega/\pi');

ylabel('幅度 (dB)');

title('(b)');

 实验结果如下:

标签:滤波器,Fs,数字,IIR,apos,pi,阻带
From: https://blog.csdn.net/m0_74460387/article/details/139275079

相关文章

  • 深度学习入门:3.手写数字分类
    这一章将基于Pytorch在只用全连接网络情况下实现手写数字识别,让大家基本了解怎么实现模型训练和预测使用。完整代码在最下面。一.MNIST数据集1.数据集介绍MNIST数据集包含60,000个训练样本和10,000个测试样本,每个样本都是28x28像素的灰度图像,表示一个0到9的手写......
  • WebGIS 智慧城市可视化合集 | 图扑数字孪生
    智慧城市可视化建设不仅提升了城市管理的科技含量和效率,还促进了城市可持续发展,提升了居民的生活质量。随着技术的不断发展和应用,智慧城市可视化建设将会更加丰富和完善,为城市发展带来更加广阔的前景。效果展示图扑应用自研HTforWeb产品搭建轻量化GIS智慧城市,一屏覆盖城市......
  • 2024提升数字思维能力加快企业数字化转型(74页PPT)
    方案介绍:本报告的价值在于为企业提供了一套系统的提升数字思维能力、加快数字化转型的理论框架和实践指南。通过本报告的学习和应用,企业可以更加清晰地认识到数字化转型的重要性和紧迫性,明确自身在数字化转型中的优势和不足,并找到适合自己的转型路径和策略。同时,本报告也为企......
  • CSP历年复赛题-P1179 [NOIP2010 普及组] 数字统计
    原题链接:https://www.luogu.com.cn/problem/P1179题意解读:统计l~r之间的整数包括多少个数字2。解题思路:枚举每一个数,对每一个数的每一位数字进行判断。100分代码:#include<bits/stdc++.h>usingnamespacestd;intl,r,ans;intmain(){cin>>l>>r;f......
  • 京准电子 | GPS北斗卫星时钟服务器助力数字化机场建设
    京准电子|GPS北斗卫星时钟服务器助力数字化机场建设京准电子|GPS北斗卫星时钟服务器助力数字化机场建设京准电子科技官微——ahjzsz摘要:某机场指挥调度、离港系统、航显广播等多个重要信息系统之间的时钟同步,对机场的正常运营和安全起到了重要作用,该文简单介绍该系统在某机......
  • 数字孪生——介绍
    数字孪生是一种先进的技术概念,其核心在于通过充分利用物理模型、传感器更新、运行历史等数据,在虚拟空间中完成对应实体装备的全生命周期过程的映射。以下是关于数字孪生的详细解释:定义:数字孪生是一种超越现实的概念,可以被视为一个或多个重要的、彼此依赖的装备系统的数字映射......
  • CDGA|像治理空气和水一样来治理数据:构建数字时代的新秩序
    在数字时代,数据已经渗透到我们生活的方方面面,成为推动社会发展的重要引擎。然而,随着数据的爆炸式增长,如何有效治理数据,确保其安全、可靠和合规,成为摆在我们面前的一大挑战。为了应对这一挑战,我们需要像治理空气和水一样来治理数据,构建数字时代的新秩序。治理数据如同治理空气......
  • 华贝甄选:通证经济生态引领者,数字化时代的引擎!
    在数字经济的浪潮中,华贝甄选科技有限公司凭借其卓越的技术实力和创新精神,成为通证经济领域的翘楚。华贝甄选致力于为客户提供全方位的通证经济解决方案,助力企业实现数字化转型。我们的团队由行业顶尖专家组成,拥有丰富的经验和深厚的技术功底。通证经济,作为数字经济的重要......
  • 突破带宽瓶颈,5SGSMD4E1H29C2G、5SGSMD4E1H29I2G、5SGSMD4E3H29I3G Stratix® V GS FPG
    StratixVGSFPGA——600-Mbps至12.5-Gbps收发器,适用于高性能数字信号处理(DSP)应用。中文参数:5SGSMD4E1H29C2GLAB/CLB数:135840逻辑元件/单元数:360000总RAM位数:19456000I/O数:360电压-供电:0.87V~0.93V安装类型:表面贴装型工作温度:0°C~85°C(TJ)封装/外壳:780-BBGA,FCBGA......
  • C# 在Excel中添加、应用或删除筛选器 (日期筛选、文本筛选、数字筛选)
    自动筛选器是Excel中的一个基本但极其有用的功能,它可以让你根据特定的条件来自动隐藏和显示你的数据。当有大量的数据需要处理时,这个功能可以帮你快速找到你需要的信息,从未更加有效地分析和处理相关数据。下面将介绍如何使用免费.NETExcel库在Excel中添加、应用和删除自动筛选......