使用滤波器产生回声.doc

蒁%原始信号 蕿x1=wavread(ding.wav); %读取信号 螄figure;plot(x1);title(原始信号); 螀 sound(5*x1,40000); %对原始声音的回放 艿%采样信号 蚇[x,fs,bits]=wavread(ding.wav,[10000 40000]); 膄x=x(:,1); %只取单声道 薁sound(5*x,fs); %对声音的回放 莀n1=0:2000; 螅N=size(x,1); 薃figure;subplot(2,1,1);plot(x);title(采样后语音信号的时域波形); 芁Y=fft(x,2001); subplot(2,1,2); %对信号做2001点FFT变换 plot(n1(1:1000),Y(1:1000)); title(采样后语音信号的频谱图); %信号幅度 莁%延时的信号 肈[x,fs,bits]=wavread(ding.wav,[10000 40000]); 肃 x=x(:,1); %只取单声道 羁n1=0:2000; 腿z=[zeros(200,1);x]; figure(2); %信号的延时 芆subplot(2,1,1);plot(z); title(延时后的时域图); %画出延时后的信号时域图 蚆Z=fft(z,2001); 螂subplot(2,1,2);plot(n1(1:1000),Z(1:1000)); %延时后的信号频谱图 芀title(延时后的频谱图 ); 蕿sound(5*z,fs); 膅%混响后信号 蒂[x,fs,bits]=wavread(ding.wav,[10000 40000]); %读取语音信号 肇x=x(:,1); %只取单声道 蚇n1=0:2000; 薅z=[zeros(200,1);x]; %对语音信号进行延时 芃x=[x;zeros(200,1)]; %使语音信号与延时后信号同等长度 腿y1=x+z; %信号的混响 袅figure(3);subplot(2,1,1);plot(y1); title(混响的时域图); %混响时域图 羄Y1=fft(y1,2001); %对混响信号2001点的FFT变换 羃subplot(2,1,2);plot(n1(1:1000),Y1(1:1000)); title(混响的频谱图); %混响频谱图 sound(5*y1,fs); %回放混响后的信号 膀%单回声滤波器 膈[x,fs,bits]=wavread(ding.wav,[10000 40000]); %读取语音信号 蒃x=x(:,1); %只取单声道 螃n1=0:2000; 羈a=0.5; %a取小于等于1 芆z=[zeros(200,1);x]; %对语音信号进行延时 袃x=[x;zeros(200,1)]; %使语音信号与延时后信号同等长度 蒄y2=x+z*0.5; %信号经单回声滤波 聿figure;subplot(2,1,1);plot(y2);title(单回声滤波器时域图); 蚈Y2=fft(y2,2001); %对单回声信号做2001点的FFT变换 薆subplot(2,1,2);plot(n1(1:1000),Y2(1:1000)); %单回声信号频谱图 羀title(单回声滤波器频谱图); 肀sound(5*y2,fs,bits); %回放经单回声滤波器后的语音信号 螇%无限回声滤波器 羆[x,fs,bits]=wavread(ding.wav,[10000 40000]); %读取语音信号 蚀x=x(:,1);