信号与系统simulink声音的录制播放滤波

综合训练②

实验内容：利用matlab中simulink工具，从计算机声卡中录取音频文件。在matlab中利用该音频文件，产生带回音的音频文件，要求回音产生的次数、间隔、大小可以调整，并通过simulink从声卡播放混有回音的音频文件。再设计去回音的系统，消除音频中的回音，在通过simulink从声卡播放消去回音的音频文件。

实验步骤：

一、利用simulink从声卡录取音频文件，并产生回声。

分析：录取音频文件，则需要利用一个接收信号的装置和一个储存信号的装置。若要产生回声则需要一个信号的衰减和一个信号的延迟共同作用。

如图所示，其中话筒型装置用来收集外部声音信号，通过一个放大器和一个延迟器，再与原信号共同叠加变成回声信号，其中右下角的装置用来收集声音信号。喇叭装的装置用来输出声音信号。（理论上也可以采用负反馈的方式来使原信号产生回声）

经过回声叠加以后产生的音频信号如图所示。

二、消除音频中的回声

分析：若要消除音频中的回声，则需要做一个滤波器，则先应该把其延迟和回声大小大概估计，然后将延迟的信号过滤掉。[y,fs,nbits]=wavread('output');

figure(1)

r=xcorr(y);

plot(r);

grid on;

title('y的自相关函数');

[u,v]=max(r);

r1=r;

r1(v-100:v+100,1)=0;

[u1,v1]=max(r1);

N=v-v1;%利用自相关函数求出回声延迟

a=[1,zeros(1,384),0.5];

b=1;

z=filter(1,a,y);

t=[0:441343];

figure(2),plot(t,z,'r',t,y,'b'),grid on;

title('带回声的音频信号，回声衰减a=0.5');

xlabel('t');

sound(z,44100);

以上程序利用wavread读取声音信号，然后再对回声延迟和回声衰减进行估计。

此为过滤前的输出波形：

然后此为过滤后的输出波形：

以下为滤波程序：

load lineup.mat sound(y,8192);

a=[1,zeros(1,999),0.5]; b=1;

z=filter(1,a,y);

t=linspace(0,0.855,7000);

figure(1),plot(t,0.5*z,'r',t,0.5*y,'b'),grid on;

title('带回声的信号与回声消除之后的信号进行对比'); xlabel('t');

sound(z,8192);

三、通过simulink从声卡播放消去回音的音频文件

分析：将上述已经过滤了的波形进行保存.mat，然后通过与simulink中自带的喇叭相连，就可以播放出声音来。

Loui荣誉巨献