《数字信号处理》实验指导书概要.docVIP

实验一 离散傅里叶变换 一、实验目的 (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的概念和公式： 2、FFT算法 调用格式是 X= fft(x) 或 X=fft(x,N) 对前者，若x的长度是2的整数次幂，则按该长度实现x的快速变换，否则，实现的是慢速的非2的整数次幂的变换；对后者，N应为2的整数次幂，若x的长度小于N，则补零，若超过N，则舍弃N以后的数据。Ifft的调用格式与之相同。 五、实验内容及实验步骤 １、使用MATLAB产生如下信号： （1）x(t)＝－2u(t－1)，－1t5 clear clear all;close all; t=-1:0.01:5;t0=1; u=[(t-t0)=0]; x=-2*u; plot(t,x);axis([-1,5,-2.5,0.5]); （2）x(n)＝2(（n+n0），n0＝2 clear all close all; n=-10:0.01:10; n0=2; t=n+n0; u=[t==0]; x=2*u; plot(n,x) 2、选择合适的变换区间长度N，用DFT对下列信号进行谱分析，画出幅频特性和相频特性曲线。 （1）x1(n)=2cos(0.2πn)R10(n); （2）x2(n)=sin(0.45πn)sin(0.55πn)R51(n); （3）x3(n)=2-|n|R21(n+10)。 六、问题与思考 对序列x(n)来说，做 N=64 点和 N=256 点的 DFT，得到的幅频特性相同吗？试比较异同点。 七、实验报告要求 1、独立完成以上实验内容； 2、采用标准实验报告纸书写实验报告； 3、实验报告内容包括：实验名称、实验目的，实验原理，所采用函数的基本说明、实验内容、程序清单（附主要语句说明或注释）、实验图表（打印或手绘）、归纳总结实验结果、心得体会及其他。 实验二 基于MATLAB的IIR数字滤波器的设计 一、实验目的 了解IIR数字滤波器的特点，了解IIR数字低通滤波器的设计方法和设计过程，掌握Matlab软件所提供的各种设计函数，运用Matlab软件分别实现基于脉冲响应不变法和双线性变换法的IIR滤波器设计。 二、预习要求 1、掌握IIR数字滤波器的特点及设计方法，重点掌握脉冲响应不变法和双线性变换法设计巴特沃斯的方法； 2、掌握部分Matlab基本数学函数、常用IIR数字滤波器设计和作图函数的使用； 3、掌握Matlab编程的基本方法。 三、实验所采用的功能函数 1、IIR滤波器阶数选择函数：buttord（巴特沃斯滤波器阶数选择函数） 调用格式：[N,wc]=buttord(wp,ws,Rp,As) (适用数字滤波器) [N,wc]=buttord(wp,ws,Rp,As, s) （适用模拟滤波器） 2、完整滤波器设计函数：butter（巴特沃斯滤波器设计） 调用格式：[B,A]=butter(N,wc,ftype) [B,A]=butter(N,wc,ftype