《数字信号处理》实验指导书(正文)解读
数字信号处理实验指导书
三. 实验主要仪器设备和材料 计算机，MATLAB6.5 或以上版本 四. 实验方法、步骤及结果测试 关于 MATLAB 它是由美国的 Math Works 公司推出的一个科技应用软件，它的名字 是由 MATRIX(矩阵）和 LABORATORY（实验室）的前三个字母组合而 成 MATLAB 是一种高性能的、用于工程计算的编程软件，它把科学计 算、结果的可视化和编程都集中在一个使用方便的环境中 优势在于能很容易求解复数 数值问题，速度快且容易扩展创建新的 命令和函数 主要组成部分： （1）编程语言：以矩阵和数组为基本单位的编程语言 （2）工作环境：包括一系列的应用工具，提供编程和调试程序的环 境
I
目录
目录
前言 ......................................................................................................................................I 目录 .................................................................................................................................... II 实验一 学习使用 MATLAB............................................................................................. 1 一. 实验目的 ................................................................................................. 1 二. 实验内容和要求 ..................................................................................... 1 三. 实验主要仪器设备和材料 ..................................................................... 1 四. 实验方法、步骤及结果测试 ................................................................. 1 五. 实验报告要求 ....................................................................................... 23 六. 思考题 ................................................................................................... 23 实验二 用 FFT 做谱分析 ............................................................................................... 24 一. 实验目的 ............................................................................................... 24 二. 实验内容和要求 ................................................................................... 24 三. 实验主要仪器设备和材料 ................................................................... 25 四. 实验方法、步骤及结果测试 ............................................................... 25 五. 实验报告要求 ....................................................................................... 37 六. 思考题 ................................................................................................... 37 实验三 用双线性变换法设计 IIR 数字滤波器 ............................................................. 38 一. 实验目的 ............................................................................................... 38 二. 实验内容和要求 ................................................................................... 38 三. 实验主要仪器设备和材料 ................................................................... 38 四. 实验方法、步骤及结果测试 ............................................................... 38 五. 实验报告要求 ....................................................................................... 43 六. 思考题 ................................................................................................... 43 实验四 用窗函数法设计 FIR 数字滤波器 .................................................................... 44 一. 实验目的 ............................................................................................... 44 二. 实验内容和要求 ................................................................................... 44 三. 实验主要仪器设备和材料 ................................................................... 45 四. 实验方法、步骤及结果测试 ............................................................... 45 五. 实验报告要求 ....................................................................................... 51 六. 思考题 ................................................................................................... 51
《数字信号处理》上机实验指导书
《数字信号处理》上机实验指导书一、引言“数字信号处理”是一门理论和实验密切结合的课程，为了深入地掌握课程内容，应当在学习理论的同时，做习题和上机实验。
上机实验不仅可以帮助学生深入地理解和消化基本理论，而且能锻炼初学者的独立解决问题的能力。
所以，根据本课程的重点要求编写了四个实验。
第一章、二章是全书的基础内容，抽样定理、时域离散系统的时域和频域分析以及系统对输入信号的响应是重要的基本内容。
由于第一、二章大部分内容已经在前期《信号与系统》课程中学习完，所以可通过实验一帮助学生温习以上重要内容，加深学生对“数字信号处理是通过对输入信号的一种运算达到处理目的”这一重要概念的理解。
这样便可以使学生从《信号与系统》课程顺利的过渡到本课程的学习上来。
第三章、四章DFT、FFT是数字信号处理的重要数学工具，它有广泛的使用内容。
限于实验课时，仅采用实验二“用FFT对信号进行谱分析”这一实验。
通过该实验加深理解DFT的基本概念、基本性质。
FFT是它的快速算法，必须学会使用。
所以，学习完第三、四章后，可安排进行实验二。
数字滤波器的基本理论和设计方法是数字信号处理技术的重要内容。
学习这一部分时，应重点掌握IIR和FIR两种不同的数字滤波器的基本设计方法。
IIR滤波器的单位冲激响应是无限长的，设计方法是先设计模拟滤波器，然后再通过S~Z平面转换，求出相应的数字滤波器的系统函数。
这里的平面转换有两种方法，即冲激响应不变法和双线性变换法，后者没有频率混叠的缺点，且转换简单，是一种普遍应用的方法。
学习完第六章以后可以进行实验三。
FIR滤波器的单位冲激响应是有限长的，设计滤波器的目的即是求出符合要求的单位冲激响应。
窗函数法是一种基本的，也是一种重要的设计方法。
学习完第七章后可以进行实验四。
以上所提到的四个实验，可根据实验课时的多少恰当安排。
例如：实验一可根据学生在学习《信号与系统》课程后，掌握的程度来确定是否做此实验。
若时间紧，可以在实验三、四之中任做一个实验。
《数字信号处理》实验指导书
《数字信号处理》实验指导书实验一离散傅里叶变换一、实验目的(1) 熟悉Matlab的主要操作命令；掌握Matlab的基本使用方法，能够运用Matlab软件分别产生常见的连续信号和离散信号，并对其进行一定的运算。
(2) 理解离散傅立叶变换是信号分析与处理的一种重要变换，特别是FFT在数字信号处理中的高效率应用。
掌握DFT的理论，通过DFT对典型信号进行的频谱分析，加深对DFT的理解。
(3) 通过对同一信号，作不同点数的FFT，比较其对应的频谱，比较两者的异同点，加深理解信号频谱概念和谱分析的原理与方法，了解快速傅立叶变换（FFT）可以提高运算量的特点，并运用Matlab软件分别对离散周期信号和非周期信号进行谱分析，同时绘出幅度频谱和相位频谱。
二、预习要求1、掌握连续信号和离散信号的特点及其运算方法；2、熟悉Matlab的基本编程语言及其变量、数组、向量与矩阵和部分运算符的使用；3、掌握部分Matlab基本数学函数和作图函数的使用。
三、和实验相关的一些功能函数正弦信号：A*sin(w0*t+phi)，A*cos(w0*t+phi)，A*sin(omega*n+phi)；方波信号：square(w0*t)，square(w0*t，DUTY)，A*square(omega*n)；注意DUTY的取值情况。
指数信号：A*exp(a*t)；矩形脉冲信号：rectpulse(t)，rectpulse(t,w)；单位脉冲信号和单位阶跃信号：ones（1,n）和zeros（1,n）；基本信号运算函数：abs（幅值）、 angle（相角）。
四、实验原理1、有限长序列x(n)的DFT的概念和公式：N?1?kn?x(k)??x(n)WN?n?0?N?1?kn?x(n)?1x(k)WN??Nk?0?0?k?N?10?n?N?1《数字信号处理》实验指导书WN?e?j(2?/N)2、FFT算法调用格式是X= fft(x)或 X=fft(x,N)对前者，若x的长度是2的整数次幂，则按该长度实现x的快速变换，否则，实现的是慢速的非2的整数次幂的变换；对后者，N应为2的整数次幂，若x的长度小于N，则补零，若超过N，则舍弃N以后的数据。
《数字信号处理》实验指导书(完整)
《数字信号处理》实验指导书通信教研室安阳工学院二零零九年三月第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 表示。
《数字信号处理》实验指导书
数字信号处理实验指导书电子与信息工程学院二○一二年前言数字信号处理(DSP)研究数字序列信号的表示方法，并对信号进行运算，以提取包含在其中的特殊信息。
数字信号处理是一门技术基础课程，实验是该课程教学的重要内容，是理论联系实际的重要手段。
学生通过实验，可以验证和巩固所学的理论知识，掌握数字信号处理实验的基本技能，提高分析和解决实际问题的能力，培养认真、严谨、实事求是的工作作风。
我们根据当前通信类新课程体系的流行趋势，充分考虑通信工程类专业的特殊要求，编写了这门实验课程指导书。
在内容安排上，我们在自身的教学基础上，吸收了兄弟院校的先进经验。
我们把重点放在对学生理论联系实际、分析和解决问题能力的训练上，力求丰富实验内容，简化实验方法与步骤，化抽象为具体，让学生通过实验能够举一反三，融会贯通，提高信息处理和信息加工的能力，为以后在信息领域的发明和创造打下牢固的基础。
在实验的具体编排上，我们按照循序渐进的原则，逐步加深实验内容，注意前后实验之间的连贯性，强化基本实验技能的培养，保证实验内容的丰富性、生动性，增强学生对数字信号处理实验课程的兴趣。
目录实验一信号的谱分析 (1)实验二基-2FFT算法的软件实现 (6)实验三 IIR数字滤波器的设计 (12)实验四 FIR数字滤波器的设计 (16)实验一 信号的谱分析一、实验目的1、熟练掌握快速离散傅里叶变换（FFT ）的原理及用FFT 进行频谱分析的基本方法；2、熟悉连续信号经理想采样前后的频谱变化关系，加深对时域采样定理的理解；3、进一步了解离散傅里叶变换的主要性质及FFT 在数字信号处理中的重要作用。
二、基本原理1、离散傅里叶变换（DFT ）及其主要性质DFT 表示离散信号的离散频谱，DFT 的主要性质中有奇偶对称特性，虚实特性等。
通过实验可以加深理解。
例如：实序列的DFT 具有偶对称的实部和奇对称的虚部，这可以证明如下： 由定义∑-==10)()(N n knNW n x k X∑∑-=-=-=1010)2sin()()2cos()(N n N n kn N n x j kn N n x ππ ∑-=-=-10)()()(N n nk N NW n x k N X∑-=-=1)(N n kn NNnW Wn x∑-=-=10)(N n knN W n x∑∑-=-=+=1010)2sin()()2cos()(N n N n kn N n x j kn N n x ππ)(*)(k N X k X -=∴对于单一频率的三角序列来说它的DFT 谱线也是单一的，这个物理意义我们可以从实验中得到验证，在理论上可以推导如下： 设：)()2sin()(n R n N n x N π=其DFT 为：∑-=-=102)()(N n kn Njen x k X πkn Nj N n e n N ππ210)2sin(--=∑=kn N j N n n Nj nN j e e e j πππ21022)(21--=-∑-=∑-=+----=10)1(2)1(2)(21N n k n Nj k n N j e e j ππ从而∑-=-=-=10220)(21)0(N n n Nj nN j e e j X ππ∑-=--==-=10422)1(21)1(N n n Nj N j j N e j X π0)2(=X0)2(=-N X22)(21)1(102)2(2N j j N e e j N X N n n j n N N j =-=-=-∑-=--ππ以上这串式中)0(X 反映了)(n x 的直流分量，)1(X 是)(n x 的一次谐波，又根据虚实特性)1()1(X N X -=-，而其它分量均为零。
数字信号处理实验指导讲解
数字信号处理实验指导编写马令坤郭文强电气与电子信息学院实验一离散时间的信号和系统一、实验目的1、复习离散时间的信号和系统，复习离散时间重要类型的信号和它们的运算的实现。
2、复习离散时间信号理论中一些重要的结果，它们在数字信号处理中很有用。
二、实验原理1、典型序列单位采样序列；单位阶跃序列；实数指数序列；复数指数序列；正余弦序列；随机序列：MATLAB可用rand(1,N)和randn(1,N)来生成；周期序列。
2、序列的运算信号加；信号乘；改变比例；移位；折叠：fliplr(x)；取样和：sum(x(n1:n2)) 取样积：prod(x(n1:n2))；信号能量：sum(abs(x)^2)；信号功率：sum(abs(x)^2)/length(x)3、一些有用的结果单位采样合成：奇偶合成：几何级数：序列相关：卷积运算：∑∞-∞=-=kknkxnx)()()(δ)()()(nxnxn xoe+=1||,11<-→∑∞=aaann对∑∞-∞=-=nyxllnynxlr称为移位),()()(,),(yxconv∑∑==---=MmNkkmknyamnxbny01) ()()(差分方程： 在Matlab 中： 三、实验内容 1.典型序列的实现单位阶跃序列；实数指数序列；复数指数序列；正余弦序列；随机序列用rand(1,N)和randn(1,N)来生成； 2.序列的运算给定序列x1=[1 2 3 4 5 6 7 8 9], ns1=-4; x2=[9 8 7 6 5 4 3 2 1], ns2=4 求：1) x1+x2; 2) y3=x1×x2；3) y1=0.5×x1+0.8×x2；4) y2=0.3×x1(n)×δ(n-6)+0.8×δ(n-5)×x2(n); 5) x1和x2的反折序列; 6) x1(n)和x2(n)的功率； 7) y3=x1*x2 (线性卷积)；四、本实验用到的一些MATLAB 函数Stem(x,y),Plot(x,y): x 轴和y 轴均为线性刻度（Linear scale ） Xlable (‘x’), Ylable(‘y’) Title(‘x’)axis([xmin,xmax,ymin,ymax]):调整图轴的范围 Subplot(x,y,z):同时画出数个小图形在同一视窗中。
数字信号处理实验指导书
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．掌握各种常用的序列，理解其数学表达式和波形表示。
2．掌握在计算机中生成及绘制数字信号波形的方法。
3．掌握序列的相加、相乘、移位、反褶等基本运算及计算机实现与作用。
4．掌握线性卷积软件实现的方法。
5．掌握计算机的使用方法和常用系统软件及应用软件的使用。
6．通过编程，上机调试程序，进一步增强使用计算机解决问题的能力。
二、实验原理1．序列的基本概念离散时间信号在数学上可用时间序列来表示，其中代表序列的第n 个数字，n 代表时间的序列，n 的取值范围为∞<<∞-n 的整数，n 取其它值)(n x 没有意义。
离散时间信号可以是由模拟信号通过采样得到，例如对)(t x a 模拟信号进行等间隔采样，采样间隔为T ，得到一个{})(nT x a 有序的数字序列就是离散时间信号，简称序列。
2．常用序列常用序列有：单位脉冲序列（单位采样）)(n δ、单位阶跃序列)(n u 、矩形序列)(n R N 、实指数序列、复指数序列、正弦型序列等。
3．序列的基本运算序列的运算包括移位、反褶、和、积、标乘、累加、差分运算等。
4．序列的卷积运算∑∞∞-*=-=)()()()()(n h n x m n h m x n y上式的运算关系称为卷积运算，式中代表两个序列卷积运算。
两个序列的卷积是一个序列与另一个序列反褶后逐次移位乘积之和，故称为离散卷积，也称两序列的线性卷积。
其计算的过程包括以下4个步骤。
(1)反褶：先将)(n x 和)(n h 的变量n 换成m ，变成)(m x 和)(m h ，再将)(m h 以纵轴为对称轴反褶成)(m h -。
(2)移位：将)(m h -移位n ，得)(m n h -。
当n 为正数时，右移n 位；当n 为负数时，左移n 位。
(3)相乘：将)(m n h -和)(m x 的对应点值相乘。
(4)求和：将以上所有对应点的乘积累加起来，即得)(n y 。
《数字信号处理实验》指导书
《数字信号处理实验》实验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 转换为级联型结构。
《数字信号处理》实验指导书-167
《数字信号处理》实验指导书实验1 数字滤波器的设计实验序号：1 实验名称：数字滤波器的设计 适用专业：通信工程、电子信息工程 学 时 数：4学时一、实验目的1．掌握双线性变换法设计IIR 数字滤波器的具体设计方法及其原理和窗函数设计FIR 滤波器的设计原理与基本方法。
2．观察双线性变换的频域特性。
熟悉Butterworth 滤波器的频率特性。
3．了解各种不同窗函数对滤波器性能的影响。
4．熟悉Matlab 计算机编程。
二、实验原理1．用双线性变换法设计IIR 数字滤波器方法（1）设计思想：将模拟滤波器转换成数字滤波器的实质是，用一种从s 平面到z 平面的映射函数将Ha(s)转换成H(z)。
对这种映射函数的要求是：(1) 因果稳定的模拟滤波器转换成数字滤波器，仍是因果稳定的。
(2)数字滤波器的频率响应模仿模拟滤波器的频响，s 平面的虚轴映射z 平面的单位圆，相应的频率之间成线性关系。
脉冲响应不变法和双线性变换法都满足如上要求。
s 平面与z 平面之间满足以下映射关系：1111--+-=zz ss 平面的虚轴单值地映射于z 平面的单位圆上，s 平面的左半平面完全映射到z 平面的单位圆内。
双线性变换不存在混叠问题。
双线性变换时一种非线性变换)/ω(tg 2=Ω，这种非线性引起的幅频特性畸变可通过预畸而得到校正。
（2）以低通数字滤波器为例，将设计步骤归纳如下：·确定数字滤波器的性能指标：通带临界频率f p 、阻带临界频率f s ；通带内的最大衰减A p ；阻带内的最小衰减A s ；· 确定相应的数字角频率，ωp=2πfp ；ωs=2πfs ；·计算经过预畸的相应模拟低通原型的频率，)2/(ωtg =Ω；·根据Ωp 和Ωs 计算模拟低通原型滤波器的阶数N ，并求得低通原型的传递函数Ha(s)； ·用上面的双线性变换公式代入Ha(s)，求出所设计的传递函数H(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)这时候即可看到一个绘有正弦曲线的窗口弹出 另：每条语句后面加“；”表示不要显示当前语句的执行结果 不加“；”表示要显示当前语句的执行结果。
数字信号处理实验指导书
% Program P1_2
% Generation of a sinusoidal sequence
n = 0:40;
பைடு நூலகம்f = 0.1;
phase = 0;
A = 1.5;
arg = 2*pi*f*n - phase;
x = A*cos(arg);
clf;
% Clear old graph
附录A MATLAB系统的常用概念 .........................................................................28
附录B
信号处理工具箱函数...........................................................................33
分析，从而进一步研究它们的性质。 2．熟悉离散时间序列的 3 种表示方法：离散时间傅立叶变换（DTFT）,离
散傅立叶变换（DFT）和 Z 变换。
二．实验相关知识准备 1．用到的 MATLAB 命令 运算符和特殊字符： < > .* ^ .^ 语言构造与调试： error function pause 基本函数： angle conj rem 数据分析和傅立叶变换函数： fft ifft max min 工具箱： freqz impz residuez zplane
数字信号处理 实验指导书
长沙理工大学电气与信息工程学院
目
录
实验一：离散时间信号的时域分析........................................................................3 实验二：离散时间系统的时域分析........................................................................6 实验三：离散时间信号的频域分析........................................................................9 实验四：线性时不变离散时间系统的频域分析..................................................13 实验五： IIR数字滤波器的设计...........................................................................17 实验六： FIR数字滤波器的设计..........................................................................24
《数字信号处理》实验指导书
实验一 采样率对信号频谱的影响一、实验目的1．理解采样定理； 2．掌握采样频率确定方法； 3．理解频谱的概念； 4．理解三种频率之间的关系。
二、实验原理理想采样过程是连续信号x a (t )与冲激函数串M (t )的乘积的过程∑∞-∞=-=k skT t t M )()(δ (1))()()(ˆt M t x t xa a = (2) 式中T s 为采样间隔。
因此，理想采样过程可以看作是脉冲调制过程，调制信号是连续信号x a (t )，载波信号是冲激函数串M (t )。
显然)()()()()(ˆs k s ak s aa kT t kT xkT t t xt x-=-=∑∑∞-∞=∞-∞=δδ (3)所以，)(ˆt xa 实际上是x a (t )在离散时间kT s 上的取值的集合，即)(ˆs a kT x 。
对信号采样我们最关心的问题是，信号经过采样后是否会丢失信息，或者说能否不失真地恢复原来的模拟信号。
下面从频域出发，根据理想采样信号的频谱)(ˆΩj X a和原来模拟信号的频谱)(Ωj X 之间的关系，来讨论采样不失真的条件∑∞-∞=Ω-Ω=Ωk ssakj j X T j X )(1)(ˆ (4)上式表明，一个连续信号经过理想采样后，其频谱将以采样频率Ωs ＝2π／T s 为间隔周期延拓，其频谱的幅度与原模拟信号频谱的幅度相差一个常数因子1／T s 。
只要各延拓分量与原频谱分量之间不发生频率上的交叠，则可以完全恢复原来的模拟信号。
根据式(4)可知，要保证各延拓分量与原频谱分量之间不发生频率上的交叠，则必须满足Ωs ≥2Ω。
这就是奈奎斯特采样定理：要想连续信号采样后能够不失真地还原原信号，采样频率必须大于或等于被采样信号最高频率的两倍h s Ω≥Ω2，或者h s f f 2≥，或者2hs T T ≤(5) 即对于最高频率的信号一个周期内至少要采样两点，式中Ωh 、f s 、T h 分别为被采样模拟信号的最高角频率、频率和最小周期。
数字信号处理实验指导书
《数字信号处理》实验指导书安阳工学院电子信息与电气工程学院目录实验一离散时间信号的表示及运算 (1)实验二离散时间LTI系统的时域分析 (10)实验三 z变换及离散时间LTI系统的z域分析 (17)实验四离散傅立叶变换及其快速算法 .................... ....... .. (29)实验五 IIR数字滤波器的MATLAB实现 (46)实验六 FIR数字滤波器的MATLAB实现.............................. .. (54)附录 MATLAB主要命令函数表....................... . (61)实验一 离散时间信号的表示及运算一、实验目的1、学会运用MATLAB 表示的常用离散时间信号；2、学会运用MATLAB 实现离散时间信号的基本运算。
二、实验条件装有MATLAB 软件的PC 机 三、 实验原理1、离散时间信号在MATLAB 中的表示离散时间信号是指在离散时刻才有定义的信号，简称离散信号，或者序列。
离散序列通常用)(n x 来表示，自变量必须是整数。
离散时间信号的波形绘制在MATLAB 中一般用stem 函数。
stem 函数的基本用法和plot 函数一样，它绘制的波形图的每个样本点上有一个小圆圈，默认是空心的。
如果要实心，需使用参数“fill”、“filled”，或者参数“.”。
由于MATLAB 中矩阵元素的个数有限，所以MATLAB 只能表示一定时间范围内有限长度的序列；而对于无限序列，也只能在一定时间范围内表示出来。
类似于连续时间信号，离散时间信号也有一些典型的离散时间信号。
①单位取样序列单位取样序列)(n δ，也称为单位冲激序列，定义为)0()0(01)(≠=⎩⎨⎧=n n n δ （1-1） 要注意，单位冲激序列不是单位冲激函数的简单离散抽样，它在n =0处是取确定的值1。
在MATLAB 中，冲激序列可以通过编写以下的impDT.m 文件来实现，即function y=impDT(n)y=(n==0); %当参数为0时冲激为1,否则为0 调用该函数时n 必须为整数或整数向量。
《数字信号处理》实验指导书-新
《数字信号处理》实验指导书陈秋妹编写适用专业：电子信息工程通信工程电子科学与计算机技术浙江理工大学信息电子学院2010年7月前言《数字信号处理》是信息电子，通信工程等本科专业及其他相近专业的一门专业必修课。
通过本课程的学习，学生应掌握以下基本概念、理论和方法：采样定理、离散序列的变换、离散信号的频谱分析；离散系统的传递函数、频率响应、离散系统的基本分析方法；数字滤波器的设计理论、滤波器的软件实现；离散傅立叶变换理论、快速傅立叶变换方法；有限字长效应。
为了使学生更好地理解和深刻地把握这些知识，并在此基础上，训练和培养学生掌握离散系统的基本概念和分析方法，数字滤波器的设计和实现，以及如何利用快速傅立叶变换等DSP技术对数字信号进行分析、滤波等处理，设置了以下三个实验：（1）离散时间序列卷积和MATLAB实现；内容：使用任意的编程语言编制一个程序，实现两个任意有限序列的卷积和。
目的：理解线性非移变系统I/O关系和实现要求：掌握使用计算机实现数字系统的方法（2）FFT算法的MATLAB实现；内容：使用MATLAB编程语言编制一个程序，实现任意有限序列的FFT。
目的：理解FFT算法的意义和实现要求：掌握使用计算机实现FFT算法的方法（3）数字滤波器的设计；内容：使用MATLAB编程语言编制一个程序，实现FIR或IIR滤波器的设计目的：理解数字滤波器的设计技术要求：掌握使用计算机进行数字滤波器设计的方法（4）窗函数设计FIR滤波器；内容：使用MATLAB编程语言编制一个程序，实现FIR或IIR滤波器的设计目的：理解数字滤波器的设计技术要求：掌握使用计算机进行数字滤波器设计的方法该实验指导书是参照该课程的教学大纲而编制的，适合于信息电子工程、通信工程等本科专业及其他相近专业。
数字信号知识预备：一．典型的离散信号 1．单位抽样信号⎩⎨⎧≠==0001)(n n n δ2．脉冲串序列p(n)将)(n δ在时间轴上延迟k 个抽样周期，得⎩⎨⎧≠==-kn kn k n 01)(δ若k 从∞-变到∞+，则的所有移位可形成一个无限长的脉冲串序列p(n)∑∞-∞=-=k k n n p )()(δ3．单位阶跃序列⎩⎨⎧<≥=0001)(n n n u 4．正弦序列 5．复正弦序列 6．指数序列二．离散信号的运算 1．信号的延迟给定离散信号x(n),若信号)(1n y ,)(2n y 分别定义为)(1n y ＝x(n-k) )(2n y =x(n+k)序列下x(n)在某一时刻k 时的值可用)(n δ的延迟来表示，即)(k x =)(n x )(k n -δ2.两个信号的相加和相乘x(n)=x1(n)+x2(n)y(n)=x1(n)x2(n)y(n)=cx(n)3.信号时间尺度的变化y(t)=x(t/a),式中a>0若a>1或a<1呢？4.信号的分解5．信号的变换三．信号的分类1．连续时间信号和离散时间信号2．周期信号和非周期信号3．确定性信号和随机信号4．能量信号和功率信号5．一维信号、二维信号及多通道信号四．噪声x(n)=s(n)+u(n)x(n)=s(n) u(n)五．离散时间系统的基本概念y(n)=T[x(n)]例：一个离散时间系统的输入、输出关系是y(n)=ay(n-1)+x(n)求所给系统的单位抽样响应。
《数字信号处理》实验指导书
的相角, Ai 就是极点 pi 到单位圆上的点 e jω 的矢量长度(距离),而θ i 就是该矢量 的相角,因此有:
M
∏ B e j(ψ1 +ψ 2 +⋅⋅⋅⋅+ψ M ) j
H (e jω ) =
j =1 N
= H (e jω ) e jϕ (ω )
∏ A e j(θ1+θ2 +⋅⋅⋅⋅+θ N ) i
(1) 设有直流信号 g(t)=1,现对它进行均匀取样，形成序列 g(n)=1。试讨 论若对该序列分别作加窗、补零，信号频谱结构有何变化。 四、实验过程及结果（含程序）
12
13
14
15
16
实验三 IIR 数字滤波器的设计
一、实验目的 (1)掌握双线性变换法及脉冲相应不变法设计 IIR 数字滤波器的具体设计 方法及其原理，熟悉用双线性变换法及脉冲响应不变法设计低通、高通和 带通 IIR 数字滤波器的计算机编程。 (2)观察双线性变换及脉冲响应不变法设计的滤波器的频域特性，了解双 线性变换法及脉冲响应不变法的特点。 (3)熟悉 Butterworth 滤波器、Chebyshev 滤波器和椭圆滤波器的频率特 性
《数字信号处理》
实验指导书
班级： 学号： 姓名： 苏州科技学院 电子教研室
实验一 信号、系统及系统响应
一、实验目的
(1) 熟悉 MATLAB 平台的使用，掌握离散信号、离散系统的 MATLAB 实现。 (2)掌握根据系统函数绘制系统零极点分布图的基本原理和方法。 (3)理解离散系统频率特性分析的基本原理，掌握根据系统函数零极点分布来分 析离散系统频率响应的几何矢量法。
17
变换类型 低通
Байду номын сангаас
数字信号处理实验指导书
1.5
1
0.5
0
-0.5
-1
-1.5
-2
0
2
4
6
8
10
12
14
16
18
20
（3）用 impz 函数 a1=[1,0.75,0.125]; b1=[1,-1]; impz(b1,a1,21);
Impulse Response 1.5
1
0.5
Amplitude
0
-0.5
-1
-1.5
-2
0
2
4
6
8
10 12 n (samples)
5
……
程序计算结果： I.
y[n] + 0.75 y[n − 1] + 0.125 y[n − 2] = x[n] − x[n − 1]
a. 单位冲激响应： (1) 用 filter 函数 a1=[1,0.75,0.125]; b1=[1,-1]; n=0:20; x1=[1 zeros(1,20)]; y1filter=filter(b1,a1,x1); stem(n,y1filter); title('y1filter'); xlabel('x'); ylabel('y');
3. 编制程序求解下列两个系统的单位冲激响应和阶跃响应，并绘出其图形。要求分
别用 filter、conv、impz 三种函数完成。
y[n] + 0.75 y[n − 1] + 0.125 y[n − 2] = x[n] − x[n − 1] y[n] = 0.25{x[n − 1] + x[n − 2] + x[n − 3] + x[n − 4]}
数字信号处理实验教材正文
实验一 离散时间信号与系统的实现与表示由于MATLAB 数值计算的特点，用它来分析离散时间信号与系统是很方便的。
一、离散时间信号在MA TLAB 中，可以用一个列向量来表示一个有限长度的序列，例如：(){}3,2,1,3,2,0,1,2-=n x在MA TLAB 中可以表示为： >>x = [2 1 0 2 3 -1 2 3];但是这种表示方法没有包含抽样的位置信息，要完全表示()n x ，要用x 和n 两个向量表示，例如：(){}3,2,1,3,2,0,1,2-=n x在MA TLAB 中可以表示为： >> n=[-4 -3 -2 -1 0 1 2 3]; >> x = [2 1 0 2 3 -1 2 3];由于内存有限，MA TLAB 无法表示无限序列。
下面介绍几种典型离散信号的表示方法。
1、单位脉冲序列()⎩⎨⎧≠==0001n n n δ（1）在MATLAB 中zeros(1,N)函数来实现有限区间的()n δ，例如：>>x= zeros(1,N); >>x(1)=1;（2）我们还可以通过编写impseq 函数来实现，该函数代码如下：（箭头表示第0个抽样点的位置）function [x,n]= impseq(n0,n1,n2) %产生单位脉冲序列x(n)=δ(n-n0) %[x,n]= impseq(n0,n1,n2) % 1,n=n0 % x = % 1,n ≠n0 %n=时间序列 n=n1:n2; x=[(n-n0)==0]该函数产生如下信号，取值区间为],[21n n ，201n n n ≤≤，得到序列()⎩⎨⎧≠==-00001n n n n n n δ例1：输入命令：>>[x,n]= impseq(0,-5,5)得到如下结果： x =0 0 0 0 0 1 0 0 0 0 0 n =-5 -4 -3 -2 -1 0 1 2 3 4 52、单位阶跃序列()⎩⎨⎧<≥=0001n n n u （1）在MATLAB 中ones(1,N)函数来实现有限区间的()n u ，例如：>>x= ones(1,N);（2）我们还可以通过编写impseq 函数来实现，该函数代码如下：function [x,n]= impseq(n0,n1,n2)%产生单位阶跃序列x(n)=u(n-n0);n1≤n0≤n2 %[x,n]= impseq(n0,n1,n2) % 1,n ≥n0 % x = % 1,n<n0 %n=时间序列 n=n1:n2; x=[(n-n0)>=0]该函数产生如下信号，取值区间为],[21n n ，201n n n ≤≤，得到序列()⎩⎨⎧<≥=-0001n n n n n n u例2：输入命令：>>[x,n]= impseq(0,-5,5)得到如下结果： x =0 0 0 0 0 1 1 1 1 1 1 n =-5 -4 -3 -2 -1 0 1 2 3 4 53、正弦序列sin 函数就可以产生正弦波，例如：>>n= 0:N-1; >>x=A*sin(n) 4、实指数序列()R a a n x n ∈=MATLAB 实现：>>n= 0:N-1; >>x=a..^n 5、复指数序列()()n j e n x ωσ+=MATLAB 实现：>>n= 0:N-1;>>x=exp((lu+j*w)*n); 6、随机序列rand(1,N)：产生[0,1]上均匀分布的随机序列；randn(1,N)：产生均值为0，方差为1的高斯随机序列，即白噪声序列。
数字信号处理实验指导书 正文
目录实验一熟悉MATLAB环境（2学时） (2)实验二求线性时不变系统的输出(2学时) (3)实验三时域及频域采样定理(2学时) (8)实验四零极点分布对系统频率响应的影响(2学时) ..... 错误！未定义书签。
实验五用DFT（FFT）对信号进行频谱分析(2学时) .. (13)实验六IIR滤波器的设计(2学时) (14)实验七FIR滤波器的设计(2学时) ................... 错误！未定义书签。
实验八数字音频信号的分析与处理(4学时) (22)附录： MATLAB基本操作及常用命令 (30)实验一 熟悉MATLAB 环境（2学时）一、 实验目的1．熟悉MATLAB 的主要操作命令。
2．学会用MATLAB 创建时域离散信号。
3．学会创建MATLAB 函数。
二、 实验原理参阅附录MATLAB 基本操作及常用命令。
三、 实验内容完成以下操作。
1．数组的加、减、乘、除运算。
输入A=[1 2 3 4]；B=[3 4 5 6]；计算：C=A+B ；D=A-B ；E=A.*B ；F=A./B ；G=A.^B ；并用stem 语句画出A 、B 、C 、D 、E 、F 、G 。
2．用MATLAB 实现以下序列(1)单位抽样序列(2)单位阶跃序列⎩⎨⎧<≥=000,0,1)n -(n n n n n u(3)矩形序列⎩⎨⎧≥<-≤≤=),0(,0)10(,1)(N n n N n n R N(4)正弦序列x (n )=5sin(0.5πn + π/4) (5)指数序列x (n )=exp(-0.5n )⎩⎨⎧≠==00,0,1)n -(n n n n n δ3．用MA TLAB 生成以下两个序列：)4(5)3(4)2(3)1(2)()(-+-+-+-+=n n n n n n x δδδδδ )3(2)2()1(2)()(-+-+-+=n n n n n h δδδδ并作以下运算，并绘制运算后序列的波形。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
矿产资源开发利用方案编写内容要求及审查大纲
矿产资源开发利用方案编写内容要求及《矿产资源开发利用方案》审查大纲一、概述
㈠矿区位置、隶属关系和企业性质。
如为改扩建矿山, 应说明矿山现状、
特点及存在的主要问题。
㈡编制依据
(1简述项目前期工作进展情况及与有关方面对项目的意向性协议情况。
(2 列出开发利用方案编制所依据的主要基础性资料的名称。
如经储量管理部门认定的矿区地质勘探报告、选矿试验报告、加工利用试验报告、工程地质初评资料、矿区水文资料和供水资料等。
对改、扩建矿山应有生产实际资料, 如矿山总平面现状图、矿床开拓系统图、采场现状图和主要采选设备清单等。
二、矿产品需求现状和预测
㈠该矿产在国内需求情况和市场供应情况
1、矿产品现状及加工利用趋向。
2、国内近、远期的需求量及主要销向预测。
㈡产品价格分析
1、国内矿产品价格现状。
2、矿产品价格稳定性及变化趋势。
三、矿产资源概况
㈠矿区总体概况
1、矿区总体规划情况。
2、矿区矿产资源概况。
3、该设计与矿区总体开发的关系。
㈡该设计项目的资源概况
1、矿床地质及构造特征。
2、矿床开采技术条件及水文地质条件。