《自-《数字信号处理》实验指导书-新》由会员分享，可在线阅读，更多相关《自-《数字信号处理》实验指导书-新（20页珍藏版）》请在金锄头文库上搜索。
1、数字信号处理实 验 指 导 书陈秋妹 编 写适用专业： 电子信息工程 通信工程 电子科学与计算机技术 浙江理工大学信息电子学院010年7月前 言数字信号处理是信息电子,通信工程等本科专业及其他相近专业的一门专业必修课。通过本课程的学习,学生应掌握以下基本概念、理论和方法：采样定理、离散序列的变换、离散信号的频谱分析；离散系统的传递函数、频率响应、离散系统的基本分析方法；数字滤波器的设计理论、滤波器的软件实现;离散傅立叶变换理论、快速傅立叶变换方法;有限字长效应。为了使学生更好地理解和深刻地把握这些知识,并在此基础上，训练和培养学生掌握离散系统的基本概念和分析方法,数字滤波器的设计和实现，以及如
2、何利用快速傅立叶变换等DSP技术对数字信号进行分析、滤波等处理,设置了以下三个实验:（1）离散时间序列卷积和MATLAB实现;内容：使用任意的编程语言编制一个程序，实现两个任意有限序列的卷积和。目的:理解线性非移变系统I/关系和实现要求：掌握使用计算机实现数字系统的方法(2）FT算法的ML实现;内容:使用MALA编程语言编制一个程序，实现任意有限序列的FFT。目的：理解FFT算法的意义和实现要求:掌握使用计算机实现FT算法的方法（3)数字滤波器的设计；内容：使用MATLA编程语言编制一个程序，实现FIR或IIR滤波器的设计目的:理解数字滤波器的设计技术要求：掌握使用计算机进行数字滤波器设计的方
3、法()窗函数设计FI滤波器；内容：使用MATLA编程语言编制一个程序，实现I或IR滤波器的设计目的:理解数字滤波器的设计技术要求：掌握使用计算机进行数字滤波器设计的方法该实验指导书是参照该课程的教学大纲而编制的,适合于信息电子工程、通信工程等本科专业及其他相近专业。数字信号知识预备:一典型的离散信号1单位抽样信号.脉冲串序列p(n)将在时间轴上延迟k个抽样周期,得若k从变到，则的所有移位可形成一个无限长的脉冲串序列p(n).单位阶跃序列4.正弦序列.复正弦序列6指数序列二.离散信号的运算1信号的延迟给定离散信号x(n),若信号,分别定义为=x(nk)=x(n+k)序列下（n)在某一时刻k时的值
4、可用的延迟来表示，即=2.两个信号的相加和相乘x()=1(n）2(n）y(n）x1(）x2(）(n）=cx(n)3.信号时间尺度的变化y(）(ta)，式中a若a1或a1呢？4.信号的分解5.信号的变换三.信号的分类1连续时间信号和离散时间信号2.周期信号和非周期信号3确定性信号和随机信号4能量信号和功率信号5.一维信号、二维信号及多通道信号四噪声x(n)=(n)+u(n)x()=s(n)（)五离散时间系统的基本概念y(n)=x(n)例：一个离散时间系统的输入、输出关系是y(n)y（-1)+x(n）求所给系统的单位抽样响应。例:系统式中,为常数。求所给系统的单位抽样响应六.有关离散系统的几个重要
5、定义线性对a1（n)bx2（n)的响应是ay1()y2(n)移不变性T(n）=y()Tx（k)=y(n-)3.因果性稳定性七SI系统的输入输出关系例：令,，求和的线性卷积。实验一：离散时间序列卷积和MATB实现实验学时：2实验类型：验证实验要求：必修（一）实验目的：学会用MATLAB对信号与系统分析的方法,理解离散序列卷积和的计算对进行离散信号与系统分析的重要性。（二)实验原理：、离散时间序列1（k)和()的卷积和定义:(k)f1（k)*f（k)=2、在离散信号与系统分析中有两个与卷积和相关的重要结论：a、f(k) f(）(k)即离散序列可分解为一系列幅度由（k）决定的单位序列（k）及其平移序
6、列之积。b、对线性时不变系统，设其输入序列为f(k），单位响应为h(k),其零状态响应为y(k），则有:(k)= （三）实验内容：onvm用来实现两个离散序列的线性卷积。其调用格式是:y=con(x，h)若的长度为N,h的长度为M,则y的长度L=+M1。题一:令x（n）=，h(n)=，(n）=x（n)*h(n)，求(n)。要求用subplo和st画出x（n),(n),y(n）与n的离散序列图形。题二:已知序列（k)=2(k)=调用con()函数求上述两序列的卷积和题三：编写计算两离散序列卷积和f(k)=f1(k)*2(k）的实用函数dcon().要求该程序在计算出卷积和f(k)的同时,还绘出序
7、列f1(),f()和(k）的时域波形图,并返回f(k)的非零样值点的对应向量。fuctonf,k=dnv(1,2，k,k2)%f1(k),f2(k)及f(k)的对应序号向量分别为k1，k和k。题四:试用MALAB计算如下所示序列f1(k)与f2（k)的卷积和（)，绘出它们的时域波形，并说明序列1(k)与f2(k）的时域宽度与序列f(k)的时域宽度的关系。提示：可用上述dcon()的函数来解决。f1(k）f(k)=题五：已知某LI离散系统,其单位响应h(k)=e(k)-e(k-4）,求该系统在激励为f(k）=e（k）-e(-)时的零状态响应，并绘出其时域波形图。提示：可用dcv(）的函数来解决。
8、（四)实验报告：1根据实验原理,编写代码，得出实验结果,并画出波形图。归纳、总结实验结果。.心得体会及其他。实验二：F算法的ATLAB实现实验学时：3实验类型:验证实验要求:必修（一)实验目的：理解离散傅立叶变换时信号分析与处理的一种重要变换，特别是FF在数字信号处理中的高效率应用。(二)实验原理:1、有限长序列x(n)的DT的概念和公式：、FFT算法调用格式是X= fft(x)或 =fft（,）对前者，若x的长度是的整数次幂，则按该长度实现x的快速变换，否则，实现的是慢速的非2的整数次幂的变换;对后者，N应为的整数次幂,若x的长度小于，则补零，若超过N,则舍弃N以后的数据。Ifft的调用格式
9、与之相同。例：N=;=0：N-1;x=4 32 7 8 9 0;X=ft(）X=390 -0.772 + 62929i 0 5000i 4.7782 - 7071i 5.00 .782 7.771i 0 + 000i-10782 6.2929iXk与n的维数相同，共有8个元素。例2：0.5*sin(2*pi15*t)+2*in（2*pi*40*t)。采样频率fs10Hz,分别绘制N=28、124点幅频图。clf；=100;N128；%采样频率和数据点数n=0:N-；t=n/; %时间序列=0.sn（2*i*15*)+2*sin(2*pi*t); 信号=t(x,); 对信号进行快速ourier变
10、换mag=abs(y）; %求得ouir变换后的振幅f=n*s/N; 频率序列supt（2，,1),plot(f,ma)； %绘出随频率变化的振幅xabl(频率Hz);ylbel(振幅);tle(N8);gid n;sbplot(,2,）,lot(f(:N),mag(1：N2);%绘出quis频率之前随频率变化的振幅xlbel(频率/Hz);ylbl(振幅)；tl(N=128)；rion;%对信号采样数据为24点的处理fs10;N=10；n=0:N-1;tn/fs；x0.5*sn（2*i*15*)+2i(2*pi*t); %信号yfft(x,N)； %对信号进行快速Forier变换ma=a()
11、；求取Fore变换的振幅f=n*fs;sublo(2,，3),po(f，ma）; %绘出随频率变化的振幅label(频率/Hz);ylabl（振幅)；title(N=1024);grid on;subplot(,2,4)plot(（1:N/2),ma(1：N2);绘出Nyuist频率之前随频率变化的振幅xbel(频率/z);ylabel(振幅);ttle（N1024）;gdo;运行结果:f10Hz，Nuist频率为fs=50Hz。整个频谱图是以Nyqist频率为对称轴的。并且可以明显识别出信号中含有两种频率成分：15H和40Hz。由此可以知道FFT变换数据的对称性。因此用FF对信号做谱分析,只
12、需考察0yquit频率范围内的福频特性。若没有给出采样频率和采样间隔，则分析通常对归一化频率1进行。另外,振幅的大小与所用采样点数有关，采用12点和02点的相同频率的振幅是有不同的表现值，但在同一幅图中,40Hz与15Hz振动幅值之比均为4：1，与真实振幅05:2是一致的。为了与真实振幅对应，需要将变换后结果乘以2除以N。（三)实验内容：题一：若x(n）=cos(npi/6)是一个=1的有限序列，利用MATLAB计算它的T并画出图形。题二:一被噪声污染的信号,很难看出它所包含的频率分量，如一个由50Hz和10Hz正弦信号构成的信号,受均值随机噪声的干扰，数据采样率为00H，通过FFT来分析其信
13、号频率成分,用ALAB实现。题三:调用原始语音信号ml,对其进行FT变换后去掉幅值小于1的FFT变换值,最后重构语音信号。（要求有四幅语音信号的频谱图在同一图形窗口以便比较:分别是1、原始语音信号;2、FT变换;去掉幅值小于的FF变换值;4、重构语音信号)（四)实验报告:1根据实验原理，编写代码,得出实验结果,并画出波形图。归纳、总结实验结果。心得体会及其他。实验三：数字滤波器的设计实验学时：3实验类型：设计实验要求：必修(一）实验目的：掌握IR数字低通滤波器的设计方法。（二）实验原理:1、滤波器的分类滤波器分两大类：经典滤波器和现代滤波器。经典滤波器是假定输入信号中的有用成分和希望取出的成分各自占有不同的频带。这样,当通过一个线性系统(即滤波器）后可讲欲去除的成分有
