重庆大学光电学院数字信号处理实验参考模板.docx

1、重庆大学光电学院数字信号处理实验参考模板实验题目离散时间信号分析实验时间2017年10月24日实验地点主教1118实验成绩实验性质验证性设计性综合性教师评语：出勤率好原理正确方案合理实验结果正确回答问题正确报告规范一、实验目的1掌握各种常用的序列，理解其数学表达式和波形表示。2掌握在计算机中生成及绘制数字信号波形的方法。3掌握序列的相加、相乘、移位、反转等基本运算及计算机实现与作用。4掌握线性卷积软件实现的方法。5掌握计算机的使用方法和常用系统软件及应用软件的使用。6通过编程，上机调试程序，进一步增强使用计算机解决问题的能力。二、实验内容1、用matlab编制信号产生子程序，产生以下典型信号序

2、列。2、利用MATLAB编程完成序列(5)和(6)的移位、反转、加法、乘法等运算，并绘制运算后序列的波形。3、利用MATLAB编制一个计算两个序列线性卷积的通用程序，计算上述两序列，并绘制卷积后序列的波形。 三、实验原理1序列的基本概念离散时间信号在数学上可用时间序列x(n)来表示，其中x(n)代表序列的第n个数字，n代表时间的序列。注意：x(n)只在n为整数时才有意义, n不是整数时无定义, 但不能认为是0。离散时间信号可以是由模拟信号通过采样得到，例如对模拟信号xa(t)进行等间隔采样，采样间隔为T，得到一个有序的数字序列xa(nT)就是离散时间信号，简称序列。2常用序列单位脉冲序列（单位

3、抽样）、单位阶跃序列、矩形序列、实指数序列、复指数序列、正弦型序列等。3序列的基本运算3.1 序列的加法：z(n)=x (n)+y(n)加法用于序列的合成。它把两个序列中位置序号相同的样本相加，形成新的样本序列。3.2 序列的乘法：z(n)=x(n)y(n)序列的乘法是一种非线性运算，它用于信号的调制。它把两个序列中位置序号相同的样本相乘，形成新的样本序列。4序列的卷积上式的运算关系称为卷积运算，式中代表两个序列卷积运算。两个序列的卷积是一个序列与另一个序列反褶后逐次移位乘积之和，故称为离散卷积，也称两序列的线性卷积。四个步骤：（1）反转：先将x(n)和h(n)的变量n换成m，变成x(m)和h

4、(m) ，再将h(m)以纵轴为对称轴反转成h(-m) 。（2）移位：将h(-m)移位n，得h(n-m) 。当n为正数时，右移n位；当n为负数时，左移n位。（3）相乘：将h(n-m)和x(m)的对应点值相乘。（4）求和：将以上所有对应点的乘积累加起来，即得y(n) 。四、使用仪器、器材计算机，MATLAB软件五、实验过程原始记录（实验源程序、图标、计算等）1、用matlab编制信号产生子程序，产生以下典型信号序列。2、利用MATLAB编程完成序列(5)和(6)的移位、反转、加法、乘法等运算，并绘制运算后序列的波形。3、利用MATLAB编制一个计算两个序列线性卷积的通用程序，计算上述两序列，并绘制

5、卷积后序列的波形。六、实验结果分析1如何产生方波信号序列和锯齿波信号序列？答：方波信号由square函数可得，锯齿波信号由循环语句产生。%方波T=0:0.001:2*pi; %方波的时间向量y=square(6*T); %方波w=6,周期为T=2*pi/6subplot(2,1,1);plot(T,y); %画方波title(方波);axis(T(1)-1 T(end)+1 -2 2); %坐标轴区域hold on%锯齿波x=-10:10;y=;for i=1:21 %rem为取余函数，此循环使得y(i)交替取值1和-1if rem(i,2)=0 y(i)=-1;else y(i)=1;end

6、endsubplot(2,1,2);plot(x,y);title(锯齿波)2、实验中所产生的正弦序列的频率是多少？是否是周期序列？答：实验中正弦序列的频率是10HZ，是周期序列。实验题目用FFT进行谱分析实验时间2017年11月21日实验地点主教1118实验成绩实验性质验证性设计性综合性教师评语：出勤率好原理正确方案合理实验结果正确回答问题正确报告规范一、实验目的（1）进一步加深DFT算法原理和基本性质的理解（因为FFT只是DFT的一种快速算法，所以FFT的运算结果必然满足DFT的基本性质）。熟悉FFT程序结构及编程方法。（2）熟悉应用FFT对确定信号进行谱分析方法，熟悉FFT算法原理和FF

7、T子程序的应用。（3）学习用FFT对连续信号和时域离散信号进行谱分析的方法，了解可能出现的分析误差及其原因，以便在实际中正确应有FFT。二、实验内容(1）用matlab编程产生并画出信号x1(n)、x2(n)、x3(n)、x4(n)、x5(n)。（2）用matlab编制FFT函数对上述信号进行频谱分析，并画出上述信号谱图。三、实验原理一个序列x(n)的离散时间傅里叶变换就是它的频谱函数。四、使用仪器、器材计算机，MATLAB五、实验过程原始记录（实验源程序、图标、计算等）（1）用matlab编程产生并画出信号x1(n)、x2(n)、x3(n)、x4(n)、x5(n)。（2）用matlab编制F

8、FT函数对上述信号进行频谱分析，并画出上述信号谱图。六、实验结果分析1.说明FFT长度N改变对频谱的影响答：FFT随长度N增加，主瓣宽度减小，分辨率提高，主峰位置也较准确2.简要回答以下问题：在N=8时，x2(n)和x3(n)的幅频特性会相同吗？为什么？N=16呢？FFT在什么条件下也可以用来分析周期信号序列的频谱？如果正弦信号系统sin(2f0k),f0=0.1Hz，用16点FFT来做DFT运算，得到的频谱是信号本身的真实谱吗？为什么？（1）答：不相同。当N=8时，序列x1(n)和x2(n)中相同的元素值对应的 n 值是不同的，所乘的旋转因子的值也不同，因而得到的最终结果也是不同的。同理，N

9、=16时，所得的幅频特性也是不同的。（2）答：由于FFT算法对序列长度的要求是N=2M，M为正整数。所以，当周期信号序列一个周期的长度满足N=2M（M为正整数）的条件时，FFT可以用来分析周期信号的频谱。不是真实的频谱。因为序列的周期N=10不是2的整数次幂，所以不是真实的。实验题目用双线性变换法设计IIR数字滤波器实验时间2017年11月21日实验地点主教1118实验成绩实验性质验证性设计性综合性教师评语：出勤率好原理正确方案合理实验结果正确回答问题正确报告规范一、实验目的1. 掌握用双线性变换法设计IIR DF的原理及具体设计方法，熟悉用双线性变换法设计IIR DF的计算机编程。2. 观察

10、用双线性变换法设计的DF的频响特性，了解双线性变换法的特点。3. 熟悉用双线性变换法设计BW和CB型DF的全过程。二、实验内容1、读懂所给参考程序，熟悉程序的整体结构和功能。2、设计一个CB型低通DF，通带截频fp=3000Hz，衰耗满足Apmax=3dB，阻带截频fT=3400Hz,衰耗ATmin=31dB,取样频率fs=8000Hz。写出其模拟滤波器方程及数字滤波器方程。3、设计一个BW型低通DF，满足：通带截频fp=100Hz，衰耗满足Apmax=3dB，阻带截频fT=400Hz,衰耗ATmin=15dB,取样频率fs=2000Hz。写出其模拟滤波器方程及数字滤波器方程。4、设计一个BW

11、型高通DF，满足：通带截频fp=400Hz，衰耗满足Apmax=3dB，阻带截频fT=350Hz,衰耗ATmin=15dB,取样频率fs=1000Hz。写出其模拟滤波器方程及数字滤波器方程。5、设计一个CB型带通DF，满足：通带边界频率为100Hz500Hz，通带衰耗小于3dB，过渡带宽20Hz，阻带衰耗大于15dB，取样频率fs=2000Hz。写出其模拟滤波器方程及数字滤波器方程。三、实验原理 为了克服冲激响应不变法产生的频率混叠现象，这是从S平面到Z平面的标准变换zesT的多值对应关系导致的，为了克服这一缺点，产生了双线性变换法。双线性变换法的映射函数：固定其中一个变量，则另一个是线性的。

12、或者说，这种变换对于s和z是双线性的。 这说明模拟角频率和数字角频率的关系是非线性的，但是没出现混叠。在把变换为时产生了非线性畸变。为了克服它带来的问题，通常要使按上式预修正，以抵消畸变的影响。根据上述变换关系，可以写出：设计IIR数字滤波器的一般步骤：把给出的数字滤波器的性能指标转换为模拟滤波器的性能指标根据转换后的性能指标，通过滤波器阶数选择函数，来确定滤波器的最小阶数N和固有频率Wn；由最小阶数N得到低通滤波器原型；由固有频率Wn把模拟低通滤波器原型转换为低通、高通、带通、带阻滤波器运用脉冲响应不变法或双线性变换法把模拟滤波器转换成数字滤波器。四、使用仪器、器材计算机，matlab五、实

13、验过程原始记录（实验源程序、图标、计算等）1、读懂所给参考程序，熟悉程序的整体结构和功能。2、设计一个CB型低通DF，通带截频fp=3000Hz，衰耗满足Apmax=3dB，阻带截频fT=3400Hz,衰耗ATmin=31dB,取样频率fs=8000Hz。写出其模拟滤波器方程及数字滤波器方程。3、设计一个BW型低通DF，满足：通带截频fp=100Hz，衰耗满足Apmax=3dB，阻带截频fT=400Hz,衰耗ATmin=15dB,取样频率fs=2000Hz。写出其模拟滤波器方程及数字滤波器方程。4、设计一个BW型高通DF，满足：通带截频fp=400Hz，衰耗满足Apmax=3dB，阻带截频fT

14、=350Hz,衰耗ATmin=15dB,取样频率fs=1000Hz。写出其模拟滤波器方程及数字滤波器方程。5、设计一个CB型带通DF，满足：通带边界频率为100Hz500Hz，通带衰耗小于3dB，过渡带宽20Hz，阻带衰耗大于15dB，取样频率fs=2000Hz。写出其模拟滤波器方程及数字滤波器方程。六、实验结果分析1、与脉冲响应不变设计法相比较，简述双线性变换设计法的优缺点。双线性变换法优点：克服多值映射得关系，可以消除频率的混叠缺点：时域到频域的变换是非线性的，在高频处有较大的失真脉冲响应不变法的优点：1，模拟频率到数字频率的转换时线性的；2，数字滤波器单位脉冲响应的数字表示近似原型的模拟

15、滤波器单位脉冲响应，因此时域特性逼近好。缺点：会产生频谱混叠现象，只适合带限滤波器实验题目数字音效处理器的设计与实现实验时间2017年12月11日实验地点主教1118实验成绩实验性质验证性设计性综合性教师评语：出勤率好原理正确方案合理实验结果正确回答问题正确报告规范二、实验目的数字音效处理器，是用数字音频信号处理技术来提升和达到各种声音效果的电子装置或者音频变换系统，包括均衡器、数字混响器、立体声、特殊音效器等。其通过对音频进行增益处理、时间处理、频率处理，改变了原始音频的特征，从而达到某些需要特殊环境，如音乐厅、录音室、户外才能达到的特殊声音效果。通过项目实践，学生能够体会到数字信号处理技术

16、在数字音效处理方面的应用。在实用过程中比较处理前后的语音回放效果，学生能亲切地、身临其境地体会到实验的效果，能激发学生的学习兴趣，培养学生的创新精神、实践能力。二、实验内容1、 对于任意的语音信号进行采样，对采样后的语音信号进行频谱分析。2、 通过人机交互界面完成功能控制，如语音采集、播放、音效处理等按键、滑动条控制，要求界面简单易懂，功能丰富，尽量反应使用者的视觉和使用需求。3、 实现均衡算法，要求设计增益可调的滤波器，键入以分贝为单位的增益值或衰减值，改变对应中心频率处的幅值。4、 变声，将现场采集的音频信号，实现男女变声。5、 混音，要求将多路音频混合播出，可调单路音频的掩饰时间和衰减系

17、数。6、 其他（学生自由发挥）。三、实验原理1、频谱是频率谱密度的简称，是频率的分布曲线。复杂振荡分解为振幅不同和频率不同的谐振荡，这些谐振荡的幅值按频率排列的图形叫做频谱。频谱广泛应用于声学、光学和无线电技术等方面。频谱将对信号的研究从时域引入到频域，从而带来更直观的认识。把复杂的机械振动分解成的频谱称为机械振动谱，把声振动分解成的频谱称为声谱，把光振动分解成的频谱称为光谱，把电磁振动分解成的频谱称为电磁波谱，一般常把光谱包括在电磁波谱的范围之内。分析各种振动的频谱就能了解该复杂振动的许多基本性质，因此频谱分析已经成为分析各种复杂振动的一项基本方法。2、均衡器（Equalizer），是一种可

18、以分别调节各种频率成分电信号放大量的电子设备，通过对各种不同频率的电信号的调节来补偿扬声器和声场的缺陷，补偿和修饰各种声源及其它特殊作用，一般调音台上的均衡器仅能对高频、中频、低频三段频率电信号分别进行调节。在通信系统中，在系带系统中插入均衡器能够减小码间干扰的影响。3、通过自己发音，共振峰频率的改变是基本重采样实现的，从重采样原理知道，这也同时引发了基频的变化，为保证基频变化和共振频率变化的独立、互不相关，在基频移动是必须考虑抵消重采样带来的偏移，理论上只要基频检测足够精确，确保可以保证基频改变和共振峰频率改变间的互不相关，通过搬移和改变基频、语速，实现变声。4、混音器是处理音频的一种装置，分为软件类型和硬件类型。由于原理不同，软件类型的混音器和硬件类型的混音器的应用也大不相同，前者主要用于音频处理，后者主要用于音响设置。四、使用仪器、器材计算机，matlab五、实验过程原始记录（实验源程序、图标、计算等）1.语音载入和播放2.频谱3、混音器4、均衡器5、滤波器6、变声器7、效果图 -精心整理，希望对您有所帮助