数字信号处理实验指导书
数字信号处理实验指导书(M)
数字信号处理实验电子信息科学与技术实验室2007年7月目录实验一离散时间信号的时域表示 (3)实验二离散信号的卷积和 (6)实验三离散傅立叶变换及其特性验证 (8)实验四信号处理中FFT的应用 (11)实验五离散系统的Z域分析 (15)实验六无限冲激响应(IIR)数字滤波器的三种结构 (19)实验七冲激响应不变法IIR数字滤波器设计 (23)实验八双线性变换法IIR数字滤波器设计 (26)实验一 离散时间信号的时域表示一、实验目的1、熟悉Matlab 命令，掌握离散时间信号－序列的时域表示方法。
2、掌握用Matlab 描绘二维图像的方法。
3、掌握用Matlab 对序列进行基本的运算和时域变换的方法。
二、实验原理与计算方法（一）序列的表示方法 序列的表示方法有列举法、解析法和图形法，相应的用Matlab 也可以有这样几种表示方法，分别介绍如下：1、列举法 在Matlab 中，用一个列向量来表示一个有限长序列，由于一个列向量并不包含位置信息，因此需要用表示位置的n 和表示量值的x 两个向量来表示任意一个序列，如：例1.1：>>n=[-3,-2,-1,0,1,2,3,4]; >>x=[2,1,-1,0,1,4,3,7];如果不对向量的位置进行定义，则Matlab 默认该序列的起始位置为n=0。
由于内存有限，Matlab 不能表示一个无限序列。
2、解析法对于有解析表达式的确定信号，首先定义序列的范围即n 的值，然后直接写出该序列的表达式，如：例1.2：实现实指数序列nn x )9.0()(=，100≤≤n 的Matlab 程序为：>>n=[0:10]; >>x=(0.9).^n;例 1.3：实现正余弦序列)5.0sin(2)31.0cos(3)(n n n x πππ++=，155≤≤n 的Matlab 程序为：>>n=[5:15];>>x=3*cos(0.1*pi*n+pi/3)+2*sin(0.5*pi*n); 3、图形法在Matlab 中用图形法表示一个序列，是在前两种表示方法的基础上将序列的各个量值描绘出来，即首先对序列进行定义，然后用相应的画图语句画图，如：例1.4：绘制在例1.1中用列举法表示的序列的图形，则在向量定义之后加如下相应的绘图语句：>>stem(n,x);此时得到的图形的横坐标范围由向量n 的值决定，为-3到4，纵坐标的范围由向量x 的值决定，为-1到7。
《数字信号处理》实验指导书(完整)
《数字信号处理》实验指导书通信教研室安阳工学院二零零九年三月第1章 系统响应及系统稳定性1.1 实验目的● 学会运用MATLAB 求解离散时间系统的零状态响应；● 学会运用MATLAB 求解离散时间系统的单位取样响应；● 学会运用MATLAB 求解离散时间系统的卷积和。
1.2 实验原理及实例分析1.2.1 离散时间系统的响应离散时间LTI 系统可用线性常系数差分方程来描述，即∑∑==-=-Mj jN i i j n x b i n y a 00)()( （1-1） 其中，i a （0=i ，1，…，N ）和j b （0=j ，1，…，M ）为实常数。
MATLAB 中函数filter 可对式（13-1）的差分方程在指定时间范围内的输入序列所产生的响应进行求解。
函数filter 的语句格式为y=filter(b,a,x)其中，x 为输入的离散序列；y 为输出的离散序列；y 的长度与x 的长度一样；b 与a 分别为差分方程右端与左端的系数向量。
【实例1-1】 已知某LTI 系统的差分方程为)1(2)()2(2)1(4)(3-+=-+--n x n x n y n y n y试用MATLAB 命令绘出当激励信号为)()2/1()(n u n x n=时，该系统的零状态响应。
解：MATLAB 源程序为>>a=[3 -4 2];>>b=[1 2];>>n=0:30;>>x=(1/2).^n;>>y=filter(b,a,x);>>stem(n,y,'fill'),grid on>>xlabel('n'),title('系统响应y(n)')程序运行结果如图1-1所示。
1.2.2 离散时间系统的单位取样响应系统的单位取样响应定义为系统在)(n 激励下系统的零状态响应，用)(n h 表示。
数字信号处理实验指导书
注意此书用的时候N要先付值数字信号处理实验指导书目录前言 (1)第一章MATLAB基础知识 (1)第二章MATLAB基本数值运算 (4)第三章MATLAB的图形处理功能 (8)第四章MATLAB的程序设计 (11)第五章常用数字信号处理函数 (16)第六章MATLAB在数字信号处理中的应用 (23)实验一常见离散信号的MATLAB产生和图形显示 (33)实验二离散系统的频率响应分析和零、极点分布 (37)实验三序列线性卷积、圆周卷积的计算及其关系的研究 (39)实验四利用DFT分析信号的频谱 (41)实验五信号时间尺度变换的研究 (43)实验六快速傅里叶变换及其应用 (47)实验七IIR滤波器的实现与应用 (56)实验八FIR滤波器的实现与应用 (61)第一章MATLAB基础知识§1-1 MA TLAB软件简介MATLAB，Matrix Laboratory的缩写，是由Mathworks公司开发的一套用于科学工程计算的可视化高性能语言，具有强大的矩阵运算能力。
它集数值分析、矩阵运算、信号处理和图形显示于一体，构成了一个界面友好的用户环境，在这个环境中，问题与求解都能方便地以数学的语言（主要是矩阵形式）或图形方式表示出来。
与大家常用的Fortran 和C等高级语言相比，MA TLAB的语法规则更简单，更贴近人的思维方式，被称为“草稿纸式的语言”。
§1-2 MA TLAB应用入门1．MATLAB的安装与卸载MATLAB软件在用户接口时具有较强的亲和力，其安装过程比较典型，直接运行光盘中的安装向导支撑程序SETUP.exe，按其提示一步步选择即可。
MATLAB自身带有卸载程序，在其安装目录下有uninstall子目录，运行该目录下uninstall.exe的即可；也可以通过Windows系统的安装卸载程序进行卸载。
2．MATLAB的启动与退出MATLAB安装完成后，会自动在Windows桌面上生成一个MA TLAB图标，它是指向安装目录下\bin\win32\matlab.exe的链接，双击这个图标即可来到MATLAB集成环境的基本窗口；也可以在开始菜单的程序选项中选择MATLAB 快捷方式；还可以在MA TLAB的安装路径的bin子目录中双击可执行文件matlab.exe。
数字信号处理指导书
实验一卷积运算1. 实验目的(1) MATLAB中序列的表示；(2) 序列的图形显示；(3) 序列的卷积计算。
2. 实验原理与方法(1) 信号在MATLAB中的表示方法MATLAB中用两个参数向量来表示有限长序列x(n)，一个是x(n)中各点的样值向量，一个是各点的位置向量。
两个向量长度相等，假设位置向量的第m 个元素的值为k，则样值向量的第m个元素的值即为x(k)。
(2) 序列的图形显示MATLAB中可调用stem函数来显示序列，其具体形式为：stem(X,Y)stem(...,'fill')stem(...,LineSpec)(3) 序列的卷积运算卷积和是离散信号与系统分析的有效方法和工具，两个序列x(n)和h(n)的卷积和定义为：∑∞-∞=-= =mmnhmxnhnxny)()()(*)()(利用MATLAB求离散序列卷积和的专用函数conv可以实现离散信号卷积和的计算。
其具体形式为：w = conv(u,v)3. 实验内容及步骤(1) 熟悉MATLAB造作环境，复习时域离散信号和系统的相关知识。
(2) 编写实验程序，产生以下序列并显示其图形：14234()()403()3470()c o s 01543()s in774x n R n n n x n n n x n n n x n nn ππ=-≤≤⎧⎪=-≤≤⎨⎪⎩=≤≤=-≤≤其它(3) 编制程序，计算x 2(n)*x 1(n)、x 3(n)*x 1(n)、x 4(n)*x 1(n)，并显示其计算结果。
(4) 手动计算上述卷积和，并与程序运行结果进行比较。
4．实验方式及要求每人一台安装有Matlab7.0的计算机，在计算机上编程仿真。
一人一组，独立完成。
5. 思考题脚本文件与函数文件编写上有什么区别？如何利用函数文件完成任意两序列的卷积运算？6. 实验报告要求(1) 简述实验目的及实验原理。
(2) 按实验步骤附上实验过程中离散序列的时域波形，并对所得结果进行分析和解释。
《数字信号处理》实验指导书
R6(n)
0.5
0 0.2
0
5
10
15
20
25 n
30
35
40
45
ቤተ መጻሕፍቲ ባይዱ
50
Magnitude
0.1
0 4
-3
-2
-1
0 Frequency(rad)
1
2
3
Phase
2 0 -2 -4 -3 -2 -1 0 Frequency(rad) 1 2 3
k＝0 时的直流分量及其合成的波形：
0.5
the Kth harmonic
the Kth harmonic
0
the Kth harmonic
0 5 10 15 20 25 n 30 35 40 45 50
0
-0.5
-0.5
0
5
10
15
20
sum of the first K+1 harmonics
1
sum of the first K+1 harmonics
25 n
30
35
0
-0.5
-0.5
0
5
10
15
20
sum of the first K+1 harmonics
1
sum of the first K+1 harmonics
25 n
30
35
40
45
50
1
0.5
0.5
0 0 5 10 15 20 25 n 30 35 40 45 50
0 0 5 10 15 20 25 n 30 35 40 45 50
数字信号处理实验指导书
j j
括幅频特性和相频特性）曲线。并将其和第 4 步中得到的结果进行比较。
七. 实验报告内容与要求
1． 简述实验目的、实验原理及实验方法和步骤。 2． 对各实验所得结果进行分析和解释。 3． 打印程序清单和要求的各信号波形。 4． 总结实验中的主要结论。 5． 简要回答思考题。
八. 思考
1． 信号的频域特性即信号的傅立叶变换利用 MATLAB 程序如何实现？ 2． 信号的频域特性即频率响应函数 H (e ) 利用 MATLAB 程序如何求取？
4
X (e j ) FT [ x(n)]
n
x ( n) e
j n
（2.1）
序列和信号的傅立叶变换是ω的连续函数， 而计算机只能计算出有限个离散频率点的 函数值。因此在取得频谱函数后，应该在 0～2π之间取许多点，计算这些点的频谱函数 的值，并取它们的包络，该包络才是需要的频率特性。当然，点数取得多一些，该包络才
y(n) 0.05 x(n) 0.05 x(n 1) 0.9 y(n 1) 的响应 y2 (n) ，并绘出 y2 (n) 的时域特性曲
线。
( n) ，并绘出 y1 ( n) 的 5． 利用卷积函数 conv （） 求信号 x1 ( n) 通过系统 h1 (n) 的响应 y1
j 能接近真正得频率特性。通常对 X (e ) 在［0，2π］上取模 X (e ) ，绘出幅频特性曲
j
线进行观察分析。系统的频域特性，通常是指求系统频率响应函数 H (e ) ，即系统单位 脉冲响应 h(n)的傅里叶变换。 对于线性时不变时域离散系统，当系统的输入序列为 x(n) ，系统的单位脉冲响应为 为 h(n) ，则线性时不变系统的输出序列为
《数字信号处理实验》指导书
《数字信号处理实验》实验1 常用信号产生实验目的：学习用MATLAB编程产生各种常见信号。
实验内容：1、矩阵操作：输入矩阵：x=[1 2 3 4;5 4 3 2;3 4 5 6;7 6 5 4]引用 x的第二、三行；引用 x的第三、四列；求矩阵的转置；求矩阵的逆；2、单位脉冲序列：产生δ（n）函数；产生δ（n-3）函数；3、产生阶跃序列：产生U（n）序列；产生U（n-n0）序列；4、产生指数序列：x(n)=0.5n⎪⎭⎫⎝⎛4 35、产生正弦序列：x=2sin(2π*50/12+π/6)6、产生取样函数：7、产生白噪声：产生[0，1]上均匀分布的随机信号：产生均值为0，方差为1的高斯随机信号：8、生成一个幅度按指数衰减的正弦信号：x(t)=Asin(w0t+phi).*exp(-a*t)9、产生三角波：实验要求：打印出程序、图形及运行结果，并分析实验结果。
实验2 利用MATLAB 进行信号分析实验目的：学习用MATLAB 编程进行信号分析实验内容：1数字滤波器的频率响应：数字滤波器的系统函数为：H(z)=21214.013.02.0----++++z z z z ， 求其幅频特性和相频特性：2、离散系统零极点图：b =[0.2 0.1 0.3 0.1 0.2];a=[1.0 -1.1 1.5 -0.7 0.3];画出其零极点图3、数字滤波器的冲激响应：b=[0.2 0.1 0.3 0.1 0.2];a=[1.0 -1.1 1.5 -0.7 0.3];求滤波器的冲激响应。
4、 计算离散卷积：x=[1 1 1 1 0 0];y=[2 2 3 4];求x(n)*y(n)。
5、 系统函数转换：（1）将H(z)=)5)(2)(3.0()1)(5.0)(1.0(------z z z z z z 转换为直接型结构。
（2）将H （z ）=3213210.31.123.7105.065.06.11-------+--+-zz z z z z 转换为级联型结构。
《数字信号处理》实验指导书(全)
数字信号处理实验指导书电子信息工程学院2012年6月目录实验一离散信号产生和基本运算 (3)实验二基于MATLAB的离散系统时域分析 (7)实验三基于ICETEK-F2812-A 教学系统软件的离散系统时域分析 (9)实验四基于MATLAB 的FFT 算法的应用 (16)实验五基于ICETEK-F2812-A 的FFT 算法分析 (18)实验六基于ICETEK-F2812-A 的数字滤波器设计 (20)实验七基于ICETEK-F2812-A的交通灯综合控制 (24)实验八基于BWDSP100的步进电机控制 (26)实验一离散信号产生和基本运算一、实验目的（1）掌握MATLAB最基本的矩阵运算语句。
（2）掌握对常用离散信号的理解与运算实现。
二、实验原理1.向量的生成a.利用冒号“：”运算生成向量，其语句格式有两种：A=m:nB=m:p:n第一种格式用于生成不长为1的均匀等分向量，m和n分别代表向量的起始值和终止值，n>m 。
第二种格式用于生成步长为p的均匀等分的向量。
b.利用函数linspace()生成向量，linspace()的调用格式为：A=linspace(m,n)B=linspace(m,n,s)第一种格式生成从起始值m开始到终止值n之间的线性等分的100元素的行向量。
第二种格式生成从起始值m开始到终止值n之间的s个线性等分点的行向量。
2.矩阵的算术运算a．加法和减法对于同维矩阵指令的A+BA-B对于矩阵和标量（一个数）的加减运算，指令为：A+3A-9b．乘法和除法运算A*B 是数学中的矩阵乘法，遵循矩阵乘法规则A．*B 是同维矩阵对应位置元素做乘法B=inv（A）是求矩阵的逆A/B 是数学中的矩阵除法，遵循矩阵除法规则A．/B 是同维矩阵对应位置元素相除另'A表示矩阵的转置运算3.数组函数下面列举一些基本函数，他们的用法和格式都相同。
sin(A),cos(A),exp(A),log(A)(相当于ln)sqrt(A)开平方 abs(A)求模 real(A)求实部 imag(A)求虚部 式中A 可以是标量也可以是矩阵 例： 利用等差向量产生一个正弦值向量 t=0:0.1:10 A=sin(t) plot(A)这时候即可看到一个绘有正弦曲线的窗口弹出 另：每条语句后面加“；”表示不要显示当前语句的执行结果 不加“；”表示要显示当前语句的执行结果。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
本科数字信号处理实验指导书基于EXP-III（TMS320F2812）实验教学平台编写：姚晓通兰州交通大学电工电子实验中心第一章实验系统介绍一、系统概述EL-DSP-EXP III 教学系统是一种综合的教学实验系统，采用模块化分离式结构，使用灵活，方便用户二次开发。
通过“E_LAB”和“TECH_V”扩展总线，可以扩展声、光、机、电等不同领域的控制对象。
客户可根据自己的需求选用不同类型的CPU适配板，我公司所有CPU适配板是完全兼容的，用户在不需要改变任何配置情况下，更换CPU 适配板即可做TI公司的不同类型的DSP的相关试验。
现已开发的CPU板类型有’C5000系列的：5402、5409、5410、5416，’C2000系列的2407、2812。
注：我公司将陆续推出VC5509、ARM系列CPU板，最新产品信息请向总公司或各地分公司咨询。
系统组成框图实验箱仿真器第二章算法实验指导实验一快速傅立叶变换（FFT）算法实验一．实验目的1．加深对DFT算法原理和基本性质的理解；2．熟悉FFT算法原理和FFT子程序的应用；3．学习用FFT对连续信号和时域信号进行谱分析的方法，了解可能出现的分析误差及其原因，以便在实际中正确应用FFT。
二．实验设备计算机，CCS 2.0 版软件，EXP3实验箱，DSP仿真器，导线三．基本原理1．离散傅立叶变换DFT的定义：将时域的采样变换成频域的周期性离散函数，频域的采样也可以变换成时域的周期性离散函数，这样的变换称为离散傅立叶变换，简称DFT。
2．FFT是DFT的一种快速算法，将DFT的N2步运算减少为（N/2）log2N步，极大的提高了运算的速度。
3．旋转因子的变化规律。
4．蝶形运算规律。
5．基2FFT算法。
四．实验步骤1．复习DFT的定义、性质和用DFT作谱分析的有关内容；2．复习FFT算法原理与编程思想，并对照DIT-FFT运算流程图和程序框图，了解本实验提供的FFT子程序；3．阅读本实验所提供的样例子程序；4．运行CCS软件，对样例程序进行跟踪，分析结果；记录必要的参数。
5．填写实验报告。
6．提供样例程序实验操作说明A．实验前准备用导线连接“模拟信号源单元”中2号孔接口“信号源1”和“A/D单元”的2号孔接口“ADIN2”；信号源1选择正弦波，并适当调节波形的频率和幅值，“A/D单元”JP3中的4打到“ON”。
S23全部置OFF；B．实验正确完成计算机、DSP仿真器和实验箱的连接后，系统上电。
启动CCS 2.0，Project/Open打开“exp3_01_fft”子目录下“Example_281x_fft.pjt”工程文件；双击“Example_281x_fft.pjt”及“Source”可查看各源程序；加载“Example_281x_fft.out”；在中断程序中，“k++；”处设置断点；单击“Run”运行程序，或按F5运行程序；程序将运行至断点处停止；2) 用View / Graph / Time/Frequency打开一个图形观察窗口；设置该观察图形窗口变量及参数；采用双踪观察在启始地址分别为px和pz，长度为128，数值类型为16位整型，这两段存储单元中分别存放的是经A/D转换后的输入信号和对该信号进行FFT变换的结果；单击“Animate”运行程序，或按F10运行；调整模拟信号源的频率、幅值、种类观察窗口并观察输入信号波形及其FFT变换结果；单击“Halt”暂停程序运行，关闭窗口，本实验结束。
实验结果：在CCS2.0环境，同步观察输入信号波形及其FFT变换结果；五．思考题1．对于不同的N，幅频特性会相同吗？为什么？2．FFT进行谱分析，可以应用的什么方面？六．实验报告要求1．简述实验原理及目的；2．结合实验中所给定典型序列幅频特性曲线，与理论结果比较，并分析说明误差产生的原因以及用FFT作谱分析时有关参数的选择方法；3．总结实验所得主要结论。
七．程序参数说明void kfft(pr,pi,n,k,fr,fi,l,il)：基2快速傅立叶变换子程序，n为变换点数，应满足2的整数次幂，k为幂次（正整数）；八．子程序流程图：实验二有限冲击响应滤波器（FIR）算法实验一、实验目的1．掌握用窗函数法设计FIR数字滤波器的原理和方法；2．熟悉线性相位FIR数字滤波器特性；3．了解各种窗函数对滤波特性的影响。
二、实验设备计算机，CCS 2.0 版软件，EXP3实验箱，DSP仿真器，导线三、实验原理1．有限冲击响应数字滤波器的基础理论；2．模拟滤波器原理（巴特沃斯滤波器、切比雪夫滤波器、贝塞尔滤波器）；3．数字滤波器系数的确定方法。
四、实验步骤1．复习如何设计FIR数字滤波；阅读本实验原理，掌握设计步骤；2．阅读本实验所提供的样例子程序；3．运行CCS软件，对样例程序进行跟踪，分析结果；4．填写实验报告。
5．样例程序实验操作说明A．实验前准备用导线连接“模拟信号源单元”中2号孔接口“信号源1”和“A/D单元”的2号孔接口“ADIN2”；“A/D单元”JP3中的4打到“ON”。
S23其中一个置ON，选择混频，调节信号源的频率和幅值，使混频得到满意的波形。
B．实验启动CCS 2.0，Project/Open打开“EXP3_02_fir”子目录下“Example_281x_fir.pjt”工程文件；双击“Example_281x_fir.pjt”及“Source”可查看各源程序；加载“Example_281x_fir.out”；在中断程序中，在“k++；”处设置断点；单击“Run”运行程序，程序将运行至断点处停止；用View / Graph / Time/Frequency打开一个图形观察窗口；设置观察图形窗口变量及参数为：采用双踪观察启始地址分别为px和py，长度为256，这两个数组中分别存放的是经A/D转换后的输入混叠信号（输入信号）和对该信号进行FIR滤波的结果；单击“Animate”运行程序，或按F10运行程序；调整观察窗口并观察滤波结果；单击“Halt”暂停程序运行，激活“Example_281x_fir.c”的编辑窗口；实验程序说明：该程序为51阶FIR低通滤波器算法程序，采用矩形窗函数实现，数组h 和xmid长度均为51，fs为采样频率，fstop为滤波器截止频率，可以修改以上参数来改变滤波器性能。
重新“Rebuild All ”后，重新加载，单击“Animate ”，可得到不同的实验结果。
实验结果：在CCS2.0环境，同步观察输入信号波形及其FIR 低通滤波结果。
五、思考题 1． 如果给定通带截止频率和阻带截止频率以及阻带最小衰减，如何用窗函数法设计线性相位低通滤波器？写出设计步骤。
2． 定性说明本实验中，3dB 截止频率的理论值在什么位置？是否等于理想低通的截止频率？3． 如果要求用窗函数法设计带通滤波器，且给定上下边带截止频率，试求理论带通的单位脉冲响应。
六、实验报告要求1． 简述实验目的及理论。
2． 自己设计一串数据应用样例子程序，进行滤波。
3． 总结设计FIR 滤波器的主要步骤。
4． 描绘出输入、输出数组的曲线七、FIR 程序参数说明系统函数 ∑=−=Mk k k Z b z H 0)(对应的常系数线性差分方程：()∑=−=Mk k k n x b n y 0)( 程序参数说明： 输入信号：输入信号经A/D 转换后，写入数组px ，长度256，16位无符号整型； 输出信号：FIR 低通滤波器输出，写入数组py ，长度256，16位无符号整型。
八、程序流程图：实验三无限冲击响应滤波器（IIR）算法实验一、实验目的1．熟悉设计IIR数字滤波器的原理与方法；2．掌握数字滤波器的计算机仿真方法；3．通过观察对实际信号的滤波作用，获得对数字滤波的感性认识。
二、实验设备计算机，CCS 2.0 版软件，EXP3实验箱，DSP仿真器，导线三、实验原理1．无限冲击响应数字滤波器的基础理论；2．模拟滤波器原理（巴特沃斯滤波器、切比雪夫滤波器、贝塞尔滤波器）；3．双线性变换的设计原理。
四、实验步骤1．复习有关巴特沃斯滤波器设计和用双线性变换法设计IIR数字滤波器的知识；2．阅读本实验所提供的样例子程序；3．运行CCS软件，对样例程序进行跟踪，分析结果；4．填写实验报告。
5．样例程序实验操作说明A．实验前准备用导线连接“模拟信号源单元”中2号孔接口“信号源1”和“A/D单元”的2号孔接口“ADIN2”；“A/D单元”JP3中的4打到“ON”。
S23其中一个置ON，选择混频，调节信号源的频率和幅值，使混频得到满意的波形。
B．实验启动CCS 2.0，Project/Open打开“Algorithm”目录中“exp3_03_iir”子目录下“Example_281x_iir.pjt”工程文件；双击“Example_281x_iir.pjt”及“Source”可查看各源程序；加载“Example_281x_iir. out”；在中断程序中，在“k++；”处设置断点；单击“Run”运行程序，程序将运行至断点处停止；用View / Graph / Time/Frequency打开一个图形观察窗口；采用双踪观察在启始地址分别为px和py，长度为256的单元中数值的变化，这两个数组分别存放的是经A/D转换的混叠信号和对该信号进行IIR低通滤波后的输出信号；单击“Animate”运行程序，或按F10运行程序；调整观察窗口，并观察滤波结果；单击“Halt”暂停程序运行，激活“Example_281x_iir.c”的编辑窗口；该iir低通滤波器滤波性能参数为：采样频率为250KHz，通带内最大允许衰减3dB，阻带内最小衰减大于30dB，过渡带宽度约为50KHz；通带上限频率：20KHz；阻带下限截止频率：70KHz。
可以修改以上参数的归一化参数“nlpass”和“nlstop”来改变滤波器性能。
修改“Example_281x_iir.c”程序中“nlpass”和“nlstop”参数可改变IIR低通滤波器的滤波性能。
重新“Rebuild All”后，加载，单击“Animate”，可得到不同的结果。
实验结果：在CCS2.0环境下，同步观察输入信号及其IIR低通滤波结果。
五、思考题1．试述用双线性变换法设计数字滤波器的过程？2．实验中，计算每个二阶滤波器的输出序列时，如何确定计算点数？3．对滤波前后的信号波形，说明数字滤波器的滤波过程与滤波作用。
六、实验报告要求1． 简述IIR 滤波器的基本原理；2． 对比FIR 滤波器与IIR 滤波器的异同；3． 描绘出输入、输出数组的曲线。
七、IIR 程序参数说明系统函数：对应的常系数线性差分方程：()∑=−+=Nk k k n y a n x n y 0)()( 程序参数说明： void biir2lpdes(double fs, double nlpass, double nlstop, double a[], double b[])：IIR 低通滤波器参数设计子程序参数说明：fs ：采样频率； nlpass ：通带上限频率归一化参数； nlstop ：阻带下限截止频率归一化参数；设置时，采样频率对应为1，应使“nlpass ”和“nlstop ”两参数均要小于0.5，且“nlpass ”要比“nlstop ”小0.2，否则，将不能满足阻带的最大衰减大于30dB 。