【含源代码】北邮dsp-MATLAB实验三梳状滤波器的应用

Dsp-matlab实验

实验三：梳状滤波器的应用

设

计

报

告

课题名称：梳状滤波器的应用

学生姓名：

班级：

班内序号：

学号：

日期：2015/06/15

目录

一、实验内容·········································

二、Matlab运行结果（含分析）·································

三、Matlab源代码····························

四、遇到的难题与解决方法····························

参考文献·························································

一、实验内容

录制一段自己的话音，时间长度及取样频率自定；对该段声音加入一次反射、三次反射和无穷多次反射。试验报告要求：

1、对试验原理的说明；

回声往往是原始声音衰减后的多个延迟叠加而组成的，因此回声可以用延迟单元来生成。X(n)表示原始声音信号，α为衰减系数，N为延迟周期，回声信号Y(n)=X(n)α*x(n-T）+α^2*x(n-2T)+……+α^N*x(n-NT). Z变换后的系统函数H（Z)可由梳状滤波器实现。MATLAB filter函数可用来仿真差分方程，本次实验用的就是这个函数。

2、在同一张图上，绘制原声音序列()

x n、加入一次反射后的声音序列

1()

x n、加入三次反射后的声音序列

3()

x n和加入无穷多次反射后的声音序列()

I

x n；

其中蓝色为原声音序列x(n)，粉红色为加入一次反射后的声音序列

x1(n)，绿色为加入三次反射后的声音序列x3(n)，红色为加入无穷多次反射后的声音序列x ∞(n)。

二、Matlab 运行结果（含分析）·

结合上述各序列，分析延时、衰减系数对回声效果的影响（提示：定量考察序列()x n 、1()x n 、3()x n 和()I x n 之间的区别）

延时不变时，衰减系数a 从零增大到1的过程中，回声效果由差变好再变差。a 很小时几乎听不到回声，a 在0.5±0.1时回声效果最明显，a 接近1时声音变得很不清晰，几乎不可识别。衰减系数不变时延时T 从零增大的过程中回声效果由差变好再变差。T 接近0时可以听到回声，但多次回声的层次感不清晰。0.1s1s 时回声感不强，感觉只是单纯地对输入声音的重复。

三、Matlab 源代码·

>> [x,fs]=audioread('a.wav');sound(x,fs);a=0.6;T=0.2;

y1=filter([1,zeros(1,T*fs-1),a],1,x);sound(y1,fs);wavwrite(y1,fs,'echo1.wav');

y2=filter([1,zeros(1,T*fs-1),a,zeros(1,T*fs-1),a^2,zeros(1,T*fs-1),a^3],1,x);

sound(y2,fs);wavwrite(y2,fs,'echo2.wav');y3=filter(1,[1,zeros(1,T*fs-1),a],x);sound(y3,fs);wavwri te(y3,fs,'echo3.wav');plot(y3,'m');

hold on;plot(y2,'r'); hold on;plot(y1,'g');hold on;plot(x,'b');

四、遇到的难题与解决办法

最开始遇到的问题是matlab 软件安装问题，因为电脑环境的特殊性尝试了多次才成功； 在建模过程中发现对实验原理因为学习时间过长有些不熟悉，于是翻书查阅复习，熟悉实验原理；

在实验过程中因为粗心，忘记保存，没有打符号等等之类问题使系统开始报错，细心调试之后成功建模

参考文献

《dsp-梳状滤波器的应用》来自百度文库《数字信号处理·第二版》科学出版社