重庆大学光电学院数字信号处理实验参考模板.docx
《重庆大学光电学院数字信号处理实验参考模板.docx》由会员分享，可在线阅读，更多相关《重庆大学光电学院数字信号处理实验参考模板.docx（33页珍藏版）》请在冰豆网上搜索。
重庆大学光电学院数字信号处理实验参考模板
实验题目
离散时间信号分析
实验时间
2017年10月24日
实验地点
主教1118
实验成绩
实验性质
□验证性□设计性√综合性
教师评语：
□出勤率好□原理正确□方案合理
□实验结果正确□回答问题正确□报告规范
一、实验目的
1．掌握各种常用的序列，理解其数学表达式和波形表示。
2．掌握在计算机中生成及绘制数字信号波形的方法。
3．掌握序列的相加、相乘、移位、反转等基本运算及计算机实现与作用。
4．掌握线性卷积软件实现的方法。
5．掌握计算机的使用方法和常用系统软件及应用软件的使用。
6．通过编程，上机调试程序，进一步增强使用计算机解决问题的能力。
二、实验内容
1、用matlab编制信号产生子程序，产生以下典型信号序列。
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.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（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编制一个计算两个序列线性卷积的通用程序，计算上述两序列，并绘制卷积后序列的波形。
六、实验结果分析
1．如何产生方波信号序列和锯齿波信号序列？
答：
方波信号由square函数可得，锯齿波信号由循环语句产生。
%方波
T=0:
0.001:
2*pi;%方波的时间向量
y=square（6*T）;%方波w=6,周期为T=2*pi/6
subplot（2,1,1）;
plot（T,y）;%画方波
title（'方波'）;
axis（[T
（1）-1T（end）+1-22]）;%坐标轴区域
holdon
%锯齿波
x=-10:
10;
y=[];
fori=1:
21%rem为取余函数，此循环使得y（i）交替取值1和-1
ifrem（i,2）==0
y（i）=-1;
else
y（i）=1;
end
end
subplot（2,1,2）;
plot（x,y）;
title（'锯齿波'）
2、实验中所产生的正弦序列的频率是多少？
是否是周期序列？
答：
实验中正弦序列的频率是10HZ，是周期序列。
实验题目
用FFT进行谱分析
实验时间
2017年11月21日
实验地点
主教1118
实验成绩
实验性质
□验证性□设计性√综合性
教师评语：
□出勤率好□原理正确□方案合理
□实验结果正确□回答问题正确□报告规范
一、实验目的
（1）进一步加深DFT算法原理和基本性质的理解（因为FFT只是DFT的一种快速算法，所以FFT的运算结果必然满足DFT的基本性质）。
熟悉FFT程序结构及编程方法。
（2）熟悉应用FFT对确定信号进行谱分析方法，熟悉FFT算法原理和FFT子程序的应用。
（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编制FFT函数对上述信号进行频谱分析，并画出上述信号谱图。
六、实验结果分析
1.说明FFT长度N改变对频谱的影响
答：
FFT随长度N增加，主瓣宽度减小，分辨率提高，主峰位置也较准确
2.简要回答以下问题：
①在N=8时，x2（n）和x3（n）的幅频特性会相同吗？
为什么？
N=16呢？
②FFT在什么条件下也可以用来分析周期信号序列的频谱？
如果正弦信号系统sin（2πf0k）,f0=0.1Hz，用16点FFT来做DFT运算，得到的频谱是信号本身的真实谱吗？
为什么？
（1）答：
不相同。
当N=8时，序列x1（n）和x2（n）中相同的元素值对应的n值是不同的，所乘的旋转因子的值也不同，因而得到的最终结果也是不同的。
同理，N=16时，所得的幅频特性也是不同的。
（2）
答：
由于FFT算法对序列长度的要求是 N=2^M，M为正整数。
所以，当周期信号序列一个周期的长度满足 N=2^M（M为正整数）的条件时，FFT可以用来分析周期信号的频谱。
不是真实的频谱。
因为序列的周期N=10不是 2 的整数次幂，所以不是真实的。
实验题目
用双线性变换法设计IIR数字滤波器
实验时间
2017年11月21日
实验地点
主教1118
实验成绩
实验性质
□验证性□设计性√综合性
教师评语：
□出勤率好□原理正确□方案合理
□实验结果正确□回答问题正确□报告规范
一、实验目的
1.掌握用双线性变换法设计IIRDF的原理及具体设计方法，熟悉用双线性变换法设计IIRDF的计算机编程。
2.观察用双线性变换法设计的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型高通DF，满足：
通带截频fp=400Hz，衰耗满足Apmax=3dB，阻带截频fT=350Hz,衰耗ATmin=15dB,取样频率fs=1000Hz。
写出其模拟滤波器方程及数字滤波器方程。
5、设计一个CB型带通DF，满足：
通带边界频率为100Hz～500Hz，通带衰耗小于3dB，过渡带宽20Hz，阻带衰耗大于15dB，取样频率fs=2000Hz。
写出其模拟滤波器方程及数字滤波器方程。
三、实验原理
为了克服冲激响应不变法产生的频率混叠现象，这是从S平面到Z平面的标准变换z＝esT的多值对应关系导致的，为了克服这一缺点，产生了双线性变换法。
双线性变换法的映射函数：
固定其中一个变量，则另一个是线性的。
或者说，这种变换对于s和z是双线性的。
这说明模拟角频率Ω和数字角频率ω的关系是非线性的，但是没出现混叠。
在把Ω变换为ω时产生了非线性畸变。
为了克服它带来的问题，通常要使Ω按上式预修正，以抵消畸变的影响。
根据上述变换关系，可以写出：
设计IIR数字滤波器的一般步骤：
把给出的数字滤波器的性能指标转换为模拟滤波器的性能指标
根据转换后的性能指标，通过滤波器阶数选择函数，来确定滤波器的最小阶数N和固有频率Wn；由最小阶数N得到低通滤波器原型；
由固有频率Wn把模拟低通滤波器原型转换为低通、高通、带通、带阻滤波器
运用脉冲响应不变法或双线性变换法把模拟滤波器转换成数字滤波器。
四、使用仪器、器材
计算机，matlab
五、实验过程原始记录（实验源程序、图标、计算等）
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=350Hz,衰耗ATmin=15dB,取样频率fs=1000Hz。
写出其模拟滤波器方程及数字滤波器方程。
5、设计一个CB型带通DF，满足：
通带边界频率为100Hz～500Hz，通带衰耗小于3dB，过渡带宽20Hz，阻带衰耗大于15dB，取样频率fs=2000Hz。
写出其模拟滤波器方程及数字滤波器方程。
六、实验结果分析
1、与脉冲响应不变设计法相比较，简述双线性变换设计法的优缺点。
双线性变换法
优点：
克服多值映射得关系，可以消除频率的混叠
缺点：
时域到频域的变换是非线性的，在高频处有较大的失真
脉冲响应不变法的优点：
1，模拟频率到数字频率的转换时线性的；
2，数字滤波器单位脉冲响应的数字表示近似原型的模拟滤波器单位脉冲响应，因此时域特性逼近好。
缺点：
会产生频谱混叠现象，只适合带限滤波器
实验题目
数字音效处理器的设计与实现
实验时间
2017年12月11日
实验地点
主教1118
实验成绩
实验性质
□验证性√设计性□综合性
教师评语：
□出勤率好□原理正确□方案合理
□实验结果正确□回答问题正确□报告规范
二、实验目的
数字音效处理器，是用数字音频信号处理技术来提升和达到各种声音效果的电子装置或者音频变换系统，包括均衡器、数字混响器、立体声、特殊音效器等。
其通过对音频进行增益处理、时间处理、频率处理，改变了原始音频的特征，从而达到某些需要特殊环境，如音乐厅、录音室、户外才能达到的特殊声音效果。
通过项目实践，学生能够体会到数字信号处理技术在数字音效处理方面的应用。
在实用过程中比较处理前后的语音回放效果，学生能亲切地、身临其境地体会到实验的效果，能激发学生的学习兴趣，培养学生的创新精神、实践能力。
二、实验内容
1、对于任意的语音信号进行采样，对采样后的语音信号进行频谱分析。
2、通过人机交互界面完成功能控制，如语音采集、播放、音效处理等按键、滑动条控制，要求界面简单易懂，功能丰富，尽量反应使用者的视觉和使用需求。
3、实现均衡算法，要求设计增益可调的滤波器，键入以分贝为单位的增益值或衰减值，改变对应中心频率处的幅值。
4、变声，将现场采集的音频信号，实现男女变声。
5、混音，要求将多路音频混合播出，可调单路音频的掩饰时间和衰减系数。
6、其他（学生自由发挥）。
三、实验原理
1、频谱是频率谱密度的简称，是频率的分布曲线。
复杂振荡分解为振幅不同和频率不同的谐振荡，这些谐振荡的幅值按频率排列的图形叫做频谱。
频谱广泛应用于声学、光学和无线电技术等方面。
频谱将对信号的研究从时域引入到频域，从而带来更直观的认识。
把复杂的机械振动分解成的频谱称为机械振动谱，把声振动分解成的频谱称为声谱，把光振动分解成的频谱称为光谱，把电磁振动分解成的频谱称为电磁波谱，一般常把光谱包括在电磁波谱的范围之内。
分析各种振动的频谱就能了解该复杂振动的许多基本性质，因此频谱分析已经成为分析各种复杂振动的一项基本方法。
2、均衡器（Equalizer），是一种可以分别调节各种频率成分电信号放大量的电子设备，通过对各种不同频率的电信号的调节来补偿扬声器和声场的缺陷，补偿和修饰各种声源及其它特殊作用，一般调音台上的均衡器仅能对高频、中频、低频三段频率电信号分别进行调节。
在通信系统中，在系带系统中插入均衡器能够减小码间干扰的影响。
3、通过自己发音，共振峰频率的改变是基本重采样实现的，从重采样原理知道，这也同时引发了基频的变化，为保证基频变化和共振频率变化的独立、互不相关，在基频移动是必须考虑抵消重采样带来的偏移，理论上只要基频检测足够精确，确保可以保证基频改变和共振峰频率改变间的互不相关，通过搬移和改变基频、语速，实现变声。
4、混音器是处理音频的一种装置，分为软件类型和硬件类型。
由于原理不同，软件类型的混音器和硬件类型的混音器的应用也大不相同，前者主要用于音频处理，后者主要用于音响设置。
四、使用仪器、器材
计算机，matlab
五、实验过程原始记录（实验源程序、图标、计算等）
1.语音载入和播放
2.频谱
3、混音器
4、均衡器
5、滤波器
6、变声器
7、效果图
---精心整理，希望对您有所帮助
升级会员 