数字信号处理实验

1、本科生实验报告实验课程 数字信号处理 学院名称 信息科学与技术学院 专业名称 信息工程 学生姓名 郑佳源 学生学号 201413010316 指导教师 杨斯涵 实验地点 5715 实验成绩 二 年 月 二 年 月填写说明1、 适用于本科生所有的实验报告（印制实验报告册除外）；2、 专业填写为专业全称，有专业方向的用小括号标明；3、 格式要求： 用A4纸双面打印（封面双面打印）或在A4大小纸上用蓝黑色水笔书写。 打印排版：正文用宋体小四号，1.5倍行距，页边距采取默认形式（上下2.54cm，左右2.54cm，页眉1.5cm，页脚1.75cm）。字符间距为默认值（缩放100%，间距：标准）；页码用

2、小五号字底端居中。 具体要求：题目（二号黑体居中）；摘要（“摘要”二字用小二号黑体居中，隔行书写摘要的文字部分，小4号宋体）；关键词（隔行顶格书写“关键词”三字，提炼3-5个关键词，用分号隔开，小4号黑体)； 正文部分采用三级标题；第1章 (小二号黑体居中，段前0.5行)1.1 小三号黑体（段前、段后0.5行）1.1.1小四号黑体（段前、段后0.5行）参考文献（黑体小二号居中，段前0.5行），参考文献用五号宋体，参照参考文献著录规则（GB/T 77142005）。目录实验一. 时域离散信号和系统的描述与分析4一.实验目的4二.实验原理4三.实验内容4四.实验结果4实验二. 时域离散信号和系统的

3、频域分析9一.实验目的9二.实验原理9三.实验内容9四.实验结果9实验三. 用FFT作谱分析12一.实验目的12二.实验原理12三.实验内容12四.实验结果12实验四. 用双线性变换法设计IIR数字滤波器14一.实验目的14二.实验原理14三.实验内容14四.实验结果14实验五.用窗函数法设计FIR数字滤波器17一实验目的17二实验原理17三.实验内容17四.实验结果17实验一. 时域离散信号和系统的描述与分析一. 实验目的掌握离散信号及其时间表示。二. 实验原理如果自变量t只在时间离散点上取值，那么称x(t)t=nTs为离散时间信号。这时的t=,-2Ts,Ts,0,Ts2Ts,=nTs。Ts

4、表示相邻两点之间的时间间隔，它的倒数fs-1/Ts，称为取样频率。一个时域离散线性非时变系统的输入/输出关系为y(n)=x(n)*h(n)=m=-xmh(n-m) （1-1）这里y(n)为系统的输出序列，x(n)为输入序列。h(n)、x(n)可以是无限长，也可以是有限长。为了计算机绘图观察方便，主要讨论有限长情况。如果h(n)和x(n)的长度分别为N和M，则y(n)的长度为L=N+M-1。这样，卷积运算就是序列移位、相乘和累加的过程，所以编程十分简单。三.实验内容1.MATLAB提供函数conv计算两序列x(n)和y(n)的卷积和。2.MATLAB提供函数randn产生随机数据，仿真“白噪声”

5、。键入程序并观察图形。3.用时间平均法检测某旋转构建的故障信号。4.给出系统：y(n)=0.75y(n-1)+0.125y(n-2)=x(n)-x(n-1)的单位冲激响应和阶跃响应。四.实验结果1.MATLAB提供函数conv计算两序列x(n)和y(n)的卷积和。（1）源程序:dalta=zeros(1,5);dalta(1)=1;x=1,2,1,3;conv(x,dalta) 结果：ans = 1 2 1 3 0 0 0 0（2）源程序：dalta=zeros(1,5);dalta(4)=1;x=1,2,1,3;conv(x,dalta)结果：ans = 0 0 0 1 2 1 3 02.M

6、ATLAB提供函数randn产生随机数据，仿真“白噪声”。键入程序并观察图形。图 1-1 白噪声3.用时间平均法检测某旋转构建的故障信号。图1-2 M=10图1-3 M=50图1-4 M=100图1-5 M=500图1-6 M=10004.给出系统：y(n)=0.75y(n-1)+0.125y(n-2)=x(n)-x(n-1)的单位冲激响应和阶跃响应。图 1-7 单位冲激响应和单位阶跃响应分析：从实验结果中可以发现，随着转动周期的提高，信噪比会逐渐增大，产生的图像中噪声逐渐变少，图像曲线趋于光滑，信号变好。实验二. 时域离散信号和系统的频域分析一.实验目的1.熟悉信号的频域表示；2.掌握有限长

7、序列离散傅立叶变换DFT。二.实验原理有限长离散序列傅立叶变换DFT的定义为：时域N点长序列X(n)的离散傅立叶变换是频域的N点的长序列X（k）。它们的关系是：X（k）=DFTx(n)=n=0N-1xnWnkn,k=0,1,2,N-1 (2-1)X（n）=IDFTx(k)=1Nn=0N-1XkWn-kn,n=0,1,2,N-1 (2-2)式中WN=e-j2N三.实验内容1.求周期序列：xn=,0,1,2,3,0,1,2,3,0,1,2,3,的DFS表示。2.求周期序列x（n）=R5n，周期分别为N=20和N=60时的|x(k)|。3.已知序列x(n)=1，2，2，1，y1，-1，-1，1，分别

8、用MATLAB中求解线性卷积的函数，循环卷积和DFT求解其线性卷积。4.已知序列x(n)=R4（n），求其DTFT及4点和8点的DFT。5.设x（n）=(-0.9)n,-5n5,求其DFT和IDFT，探讨其周期性。四.实验结果1.求周期序列：x (n)=,0,1,2,3,0,1,2,3,0,1,2,3,的DFS表示。图 2-1 DFS表示2.求周期序列x （n）=(R_5 ) （n），周期分别为N=20和N=60时的|x (k)|。 图 2-2 离散傅立叶级数3.已知序列x(n)=1，2，2，1，y1，-1，-1，1，分别用MATLAB中求解线性卷积的函数，循环卷积和DFT求解其线性卷积。 图

9、 2-3 线性卷积4.已知序列x(n)=R_4 （n），求其DTFT及4点和8点的DFT。图2-4 DTFT与DFT图2-5 DFT与IDFT分析：在时域和频域之一中进行卷积操作，对应于另一个域中的乘法。在一个域中进行“采样”的操作相当于把原信号与一个无限长的冲激序列相乘。那么，在另一个域中，相应的效果就是与无限长冲激序列的傅里叶变换进行卷积，结果就是周期延拓，并且还会产生混叠现象。实验三. 用FFT作谱分析一.实验目的1.进一步加深DFT算法原理和基本性质的理解；2.熟悉FFT算法原理和FFT子程序的应用；3.了解利用FFT进行信号消噪；4.学习用FFT对连续信号和时域离散信号进行谱分析的方

10、法，了解可能出现的分析误差及其原因，以便在实际中正确应用FFT。二.实验原理快速傅立叶变换：取N为正整数，且N=2E，E为正整数，将N点序列x(n)按n奇偶分成两组，分别计算它们的DFT，则复数乘法次数为2（N2）2=N2/2,省去一半，既然将序列分为两半可使乘法次数减半，何乐而不为，一直分下去，最后得到不能再分的序列组，每个序列组只有两点，然后依次计算两点的DFT，并合而得到整个序列的DFT，这就是所谓“基-2按时间抽取的FFT算法”的基本思想。三.实验内容试用FFT方法消噪。四.实验结果1.试用FFT方法消噪。源程序：n=1024;t=1: n;x=randn(size(t);Y=fft(

11、x, n);Y(200:824)=zeros(size(1:625);z=ifft(Y,n);for i=1:nt=i/n;t1=(t-1/8)2; t3=(t-3/8)2;t4=(t-4/8)2; t6=(t-6/8)2;t7=(t-7/8)2; m=640*pi;s(i)=cos(2*pi*280*t)*(exp(-m*t1)+exp(-m*t3)+exp(-m*t4)+exp(-m*t6)+exp(-m*t7);endy=s+z;y1=fft(y,n);s1=fft(s,n);k=1:1024;subplot(5,2,1);plot(k,y);axis(0,n,-3,3);xlabel(

12、t);ylabel(t);title(频域原信号);Y=fft(y,n);subplot(5,2,2);plot(k,Y);axis(0,n,-70,70);xlabel(t);ylabel(t);title(时域原信号);Y(924:1024)=zeros(size(1:101);y=ifft(Y,n);subplot(5,2,3);plot(k,y);axis(0,n,-3,3);xlabel(t);ylabel(t);title(频域第一次消噪);Y=fft(y,n);Y=fft(y,n);subplot(5,2,4);plot(k,Y);axis(0,n,-70,70);xlabel(

13、t);ylabel(t);title(时域第一次消噪);Y(1:100)=zeros(size(1:100);y=ifft(Y,n);subplot(5,2,5);plot(k,y);axis(0,n,-3,3);xlabel(t);ylabel(t);title(频域第二次消噪);Y=fft(y,n);subplot(5,2,6);plot(k,Y);axis(0,n,-70,70);xlabel(t);ylabel(t);title(时域第二次消噪);Y(824:924)=zeros(size(1:101);y=ifft(Y,n);subplot(5,2,7);plot(k,y);axis

14、(0,n,-3,3);xlabel(t);ylabel(t);title(频域第三次消噪);Y=fft(y,n);subplot(5,2,8);plot(k,Y);axis(0,n,-70,70);xlabel(t);ylabel(t);title(时域第三次消噪);Y(100:200)=zeros(size(1:101);y=ifft(Y,n);subplot(5,2,9);plot(k,y);axis(0,n,-3,3);xlabel(t);ylabel(t);title(频域第四次消噪);Y=fft(y,n);subplot(5,2,10);plot(k,Y);axis(0,n,-70,

15、70);xlabel(t);ylabel(t);title(时域第四次消噪);图3-1 FFT方法消噪分析：现实生活中，有用信号的频率通常比较低，而噪声的频率通常比较高，所以可以将高频的部分变为零来去除噪声。而通过上面的降噪处理，第四次左右噪声基本消失。实验四. 用双线性变换法设计IIR数字滤波器一.实验目的1.熟悉用双线性变换法设计IIR数字滤波器的原理和方法。2.掌握数字滤波器的计算机仿真方法。3.通过观察对实际心电图信号的滤波作用，获得数字滤波的感性知识。二.实验原理IIR数字滤波器与FIR数字滤波器的区别：1、单位响应IIR滤波器的单位脉冲响应为无限长，网络中有反馈回路。FIR滤波器的

16、单位脉冲响应是有限长的，一般网络中没有反馈回路。FIR滤波器的系统函数一般是一个有理分式，分母多项式决定滤波器的反馈网络。2、幅频特性IIR数字滤波器幅频特性精度很高，不是线性相位的，可以应用于对相位信息不敏感的音频信号上；FIR数字滤波器的幅频特性精度较之于IIR数字滤波器低，但是线性相位，就是不同频率分量的信号经过fir滤波器后他们的时间差不变，这是很好的性质。3、实时信号处理FIR数字滤波器是有限的单位响应也有利于对数字信号的处理，便于编程，用于计算的时延也小，这对实时的信号处理很重要。三.实验内容1.设计满足指标的模拟巴特沃斯滤波器。2.设计切比雪夫I型高通数字滤波器。3.用双线性变换

17、法设计一个巴特沃斯低通IIR数字滤波器。四.实验结果四、实验结果1.设计满足指标的模拟巴特沃斯滤波器。图4-1 巴特沃斯滤波器图4-2 巴特沃思滤波器幅频特性曲线图 4-3 产生结果2.设计切比雪夫I型高通数字滤波器。图4-4 产生结果图4-5 切比雪夫I型高通数字滤波器幅频曲线3.用双线性变换法设计一个巴特沃斯低通IIR数字滤波器。图4-6 心电图分析：用matlab实现数字滤波器设计，可通过修改滤波器的参数十分方便地改变滤波器的特性，实用性很高。实验五.用窗函数法设计FIR数字滤波器一 实验目的1.掌握用窗函数法设计FIR数字滤波器的原理和方法。2.熟悉线性相位FIR数字滤波器特性。3.了

18、解各种窗函数对滤波器特性的影响。二实验原理不同的窗函数对信号频谱的影响是不一样的，这主要是因为不同的窗函数，产生泄漏的大小不一样，频率分辨能力也不一样。信号的截断产生了能量泄漏，而用FFT算法计算频谱又产生了栅栏效应，从原理上讲这两种误差都是不能消除的，但是我们可以通过选择不同的窗函数对它们的影响进行抑制。(矩形窗主瓣窄，旁瓣大，频率识别精度最高，幅值识别精度最低；布莱克曼窗主瓣宽，旁瓣小，频率识别精度最低，但幅值识别精度最高)。设计滤波器时，通常给定的幅度特性是分段恒定的、且在频带的边界有不连续点。逼近这样的特性，需用无限冲激响应(IIR)数字滤波器。为采用有限冲激响应(FIR)实现，可以选合适的、有限时宽窗序列对无限冲激响应序列加权，以构成有限冲激响应数字滤波器实现给定特性。因为时间域内两个序列相乘，频率域内等效为两个序列离散傅里叶变换的卷积。因此，窗函数与给定滤波特性的卷积可达到以有限冲激响应逼近给定特性的目的。三.实验内容1.应用频率采样法设计相位低通FI