数字信号处理实验六FIR滤波器设计

[n,wn,bta,ftype]=kaiserord([0.3 0.45 0.65 0.75],[0 1 0],[0.01 0.1087 0.01]);%用kaiserord函数估计出滤波器阶数n和beta参数

h1=fir1(n,wn,ftype,kaiser(n+1,bta),'noscale'); [hh1,w1]=freqz(h1,1,256); figure(1)

subplot(2,1,1)

plot(w1/pi,20*log10(abs(hh1))) grid

xlabel('归一化频率w：曾慧');ylabel('幅度/db'); subplot(2,1,2)

plot(w1/pi,angle(hh1)) grid

xlabel('归一化频率w：曾慧');ylabel('相位/rad');

图一

[n,wn,bta,ftype]=kaiserord([0.3 0.45 0.65 0.75],[0 1 0],[0.01 0.1087 0.01]);%用kaiserord函数估计出滤波器阶数n和beta参数

h1=fir1(n,wn,ftype,kaiser(n+1,bta),'noscale'); [hh1,w1]=freqz(

数字信号处理实验指导书（修订版1）

数字信号处理（本科）实验指导书

（修订版1）

戴 虹 编

2010年6月

1

数字信号处理实验指导书（修订版1）

目 录

实验一 信号、 系统及系统响应 实验二 FFT频谱分析 实验三 IIR数字滤波器设计 实验四 FIR数字滤波器设计

实验五 离散系统的时域及复频域分析

综合实验 双音多频（DTMF）通信设计的MATLAB仿真 附录 实验报告格式

2

数字信号处理实验指导书（修订版1）

实验一 信号、 系统及系统响应

一、 实验目的

1. 掌握典型序列的产生方法。

2. 掌握DFT的实现方法，利用DFT对信号进行频域分析。

3. 熟悉连续信号经采样前后频谱的变化，加深对时域采样定理的理解。

4. 分别利用卷积和DFT分析信号及系统的时域和频域特性，验证时域卷积定理。 二、实验环境

1． Windows2000操作系统 2． MATLAB6.0 三、实验原理 1. 信号采样

对连续信号xa(t)=Aesin(Ω0t)u(t)进行采样，采样周期为T,采样点0≤n<50，得采样序列xa(n)= 。 2. 离散傅里叶变换(DFT) 设序列为x(n),长度为N,则

X(e

2πMjωk

N?1-at

)=DFT[x(n)]=?n?0x(n) e-jωk

n

,

jωk

其中ωk=k(k=0,1,2,…

大连理工大学实验报告

学院（系）： 专业： 班级：

姓 名： 学号： 组： ___

实验时间： 实验室： 实验台： 指导教师签字： 成绩：

实验三 滤波器设计

一、实验结果与分析

IIR部分：

1.用buttord和butter函数，直接设计一个巴特沃兹高通滤波器，要求通带截止频率为0.6?，通带内衰减不大于1dB，阻带起始频率为0.4?，阻带内衰减不小于15dB，观察其频谱响应的特点： clc,clear

wp=0.6*pi/pi; ws=0.4*pi/pi; ap=1,as=15;

[N,wn]=buttord(wp,ws,ap,as); [bz,az]= butter(N,wn,'high'); [H,W]=freqz(bz,az);

plot(W,20*log10(abs(H))); grid on;

xlabel('频率

大连理工大学实验报告

学院（系）： 专业： 班级：

姓 名： 学号： 组： ___

实验时间： 实验室： 实验台： 指导教师签字： 成绩：

实验三 滤波器设计

一、实验结果与分析

IIR部分：

1.用buttord和butter函数，直接设计一个巴特沃兹高通滤波器，要求通带截止频率为0.6?，通带内衰减不大于1dB，阻带起始频率为0.4?，阻带内衰减不小于15dB，观察其频谱响应的特点： clc,clear

wp=0.6*pi/pi; ws=0.4*pi/pi; ap=1,as=15;

[N,wn]=buttord(wp,ws,ap,as); [bz,az]= butter(N,wn,'high'); [H,W]=freqz(bz,az);

plot(W,20*log10(abs(H))); grid on;

xlabel('频率

广东工业大学研究生课程考试试卷封面

学院： 专业： 姓名： 学号： 考试科目: 学生类别：

题号 1 2 3 4 5 6 7 8 合计 平时成绩 总评成绩 分数 考试时间: 第 周星期

( 年 月 日)

开课单位：

开课学期： 年秋季

任课教师：

说明：1.以上左栏学生填写,右栏老师填写; 2.学位课考试的答题均写在答题纸上,考查课若无课堂考试可不用答题纸，但仍应填写此页封面; 3.平时成绩根据任课教师所写的课程教学大纲规定评定； 4.任课教师评完分后从网上录入成绩，再打印成绩单一式两份连同试卷、答题纸交教务员保管（保管四年）。

1.设计内容

本设计是一个基于FPGA的一维数字信号处理算法的FIR的设计，设计使用Verilog语言编写FIR滤波器的模块，通过编译和综合，并通过MATLAB和modelsim仿真对比验证设计结果。

2.设计原理

有限冲击响应(FIR)滤波器和无限冲击响应(IIR)滤波器广泛应用于数字信号处理系统中。IIR数字滤波器方便简单,但它相位的线性,要采用全通网络进行相位校正。图象处理以及数据传输,都要求信道具有线性相位特性,而有限冲击响应(FIR)

本文分析了国内外数字滤波技术的应用现状与发展趋势，介绍了数字滤波器的基本结构，在分别讨论了IIR与FIR数字滤波器的设计方法的基础上，指出了传统的数字滤波器设计方法过程复杂、计算工作量大、滤波特性调整困难的不足，提出了一种利用MATLAB信号处理工具箱（Signal Processing Toolbox）快速有效的设计由软件组成的常规数字滤波器的设计方法。给出了使用MATLAB语言进行程序设计和利用信号处理工具箱的FDATool工具进行界面设计的详细步骤。利用MATLAB设计滤波器，可以随时对比设计要求和滤波器特性调整参数，直观简便，极大的减轻了工作量，有利于滤波器设计的最优化。本文还介绍了如何利用MATLAB环境下的仿真软件Simulink对所设计的滤波器进行模拟仿真。

1.1 数字滤波器的研究背景与意义

当今，数字信号处理[1] (DSP：Digtal Signal Processing)技术正飞速发展，它不但自成一门学科，更是以不同形式影响和渗透到其他学科：它与国民经济息息相关，与国防建设紧密相连；它影响或改变着我们的生产、生活方式，因此受到人们普遍的关注。

数字化、智能化和网络化是当代信息技术发展的大趋势，而数字化是智能化和网络化的

实验六、用窗函数法设计FIR 数字滤波器

一、实验目的：

（1）熟悉基本的窗函数，及其特点。

（2）掌握用窗函数法设计FIR 数字滤波器的原理和方法。 （3）熟悉线性相位FIR 数字滤波器特性。 二.实验原理

（一）FIR滤波器的设计

FIR滤波器具有严格的相位特性，这对于语音信号处理和数据传输是很重要的。目前FIR滤波器的设计方法主要有三种：窗函数法、频率取样法和切比雪夫等波纹逼近的最优化设计方法。本实验中的窗函数法比较简单，可应用现成的窗函数公式，在技术指标要求不高的时候是比较灵活方便的。它是从时域出发，用一个窗函数截取一个理想的hd(n)得到h(n)，以有限长序列h(n)近似理想的hd(n)：如果从频域出发，用理想的hd(ejw)在单位圆上等角度取样得到H（k），根据h(k)得到H(z)将逼近理想的Hd(z)这就是频率取样法。

（二）窗函数设计法

同其它的的数字滤波器设计的方法一样，用窗函数设计滤波器也是首先要对滤波器提出性能指标。一般是给定一个理想的频率响应hd(ejw)，使所设计的FIR滤波器的频率响应

h(ejw)去逼近所要求的理想的滤波器的响应hd(ejw)窗函数设计的任务在于寻找一个可实

现（有限长单位脉冲响应）的传递函数。

H(

数字信号实验报告

实验项目名称：用窗函数法设计FIR数字滤波器 所属课程名称： 数字信号处理 实 验 类 型 ： 综合型 指 导 教 师 ： 实 验 日 期 ： 2013.12.20 班 级 ： 学 号 ： 姓 名 ：

目录

一、实验目的································1

二、实验原理与方法······················1

三、实验内容与步骤·····················1

四、上级内容·······························2 1.程序 2.运行截图 3.注释

五、思考题·····························

计算机科学与工程学院

《数字信号处理》实验报告[4]

《数字信号处理》实验报告

2

计算机科学与工程学院 含载频成分。 (3)编程序调用 MATLAB 滤波器设计函数 ellipord 和 ellip 分别设计这三个椭圆滤波 器，并绘图显示其幅频响应特性曲线。 (4)调用滤波器实现函数 filter,用三个滤波器分别对信号产生函数 mstg 产生的信 号 st 进行滤波，分离出 st 中的三路不同载波频率的调幅信号 y1(n)、y2(n)和 y3(n), 并绘 图显示 y1(n)、y2(n)和 y3(n)的时域波形，观察分离效果。

三、信号产生函数 mstg 清单function st=mstg %产生信号序列向量 st,并显示 st 的时域波形和频谱 %st=mstg 返回三路调幅信号相加形成的混合信号，长度 N=1600 N=1600 %N 为信号 st 的长度。 Fs=10000;T=1/Fs;Tp=N*T; %采样频率 Fs=10kHz,Tp 为采样时间 t=0:T:(N-1)*T;k=0:N-1;f=k/Tp; fc1=Fs/10; %第 1 路调幅信号的载波频率 fc1=1000Hz, fm1=fc1/10; %第 1 路调幅信号的调制信号

FIR数字滤波器设计

摘要 数字滤波器的输入输出均为数字信号，信号通过数字滤波器后，可以改变

频率成分的相对比例或滤除某些频率成分。数字滤波器可以分为IIR数字滤波器

和FIR数字滤波器。文中只介绍FIR数字滤波器的设计，可以根据所给定的频率

特性直接设计FIR数字滤波器。FIR数字滤波器在保证幅度特性满足要求的同时，

能够做到严格的线性特性。本文采用了窗函数法、频率采样法以及基于firls函

数和remez函数的最优化方法设计FIR滤波器。在文中对FIR滤波器进行了详细

的理论分析，并且对应于每种方法都给出了设计实例。通过编写MATLAB语言程序，

运行程序，得到幅频和相频特性图。对于窗函数和firls函数设计的滤波器，还

通过建立Simulink系统模块进行仿真，观察滤波器滤波情况。

关键词 窗函数法；频率采样法；最优化方法；仿真

I

Design of FIR Digital Filter

Abstract Firstly, the paper gives a brief description of digital filter and MATLAB

language.The digital filter can be divided i