数字信号处理梳妆滤波器应用.doc

信息与通信工程学院
数字信号处理实验报告

实验题目:梳状滤波器的应用

班 级： 2009211114
姓 名： 廖宁
学 号： 09210411
日 期： 2011.6

实验要求：

录制一段自己的话音，时间长度及取样频率自定；对该段声音加入一次反射、三次反射和无穷多次反射。试验报告要求：
1、 对试验原理的说明；
2、 在同一张图上，绘制原声音序列、加入一次反射后的声音序列、加入三次反射后的声音序列和加入无穷多次反射后的声音序列；
3、 结合上述各序列，分析延时、衰减系数对回声效果的影响（提示：定量考察序列、、和之间的区别）。

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

2.
蓝色为原始声音信号，红色为一次反射信号，黄色为三次反射信号，绿色为无穷次反射信号。

3. 从原始声音信号到无穷次反射信号，信号的幅值按照衰减系数的幂在衰减，反射次数越高，衰减越厉害；延时造成信号的水平位移，是回声产生的必要条件。延时只要大于人的听觉0.1s,人能明显感觉到回声。存在一定的延时，减小衰减系数，回声的效果越好。

4. 附源程序代码：

[x,fs]=wavread('1.wav');
sound(x,22050)
figure(1);
plot(x,'b');
hold on;%取衰减系数=0.2，延时=0.2s，0.2*22050=4410；
y=filter([1,zeros(1,4410),0.2],1,x);
sound(y,22050)
plot(y,'r')
hold on;
y1=filter([1,zeros(1,13230),0.008],1,x);
sound(y1,22050);
plot(y1,'y');
hold on;
y2=filter([1,zeros(1,44100),0.2^10],1,x);
sound(y2,22050);
plot(y2,'g');