《使用滤波器产生回声》由会员分享，可在线阅读，更多相关《使用滤波器产生回声（6页珍藏版）》请在金锄头文库上搜索。
1、蒁%原始信号蕿 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);%
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
3、:1000);% 延时后的信号频谱图芀 title( 延时后的频谱图 );蕿 soun d(5*z,fs);膅%混响后信号蒂x,fs,bits=wavread(di ng.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(混响的时域图);%混响时域图羄丫仁fft(y1,2001);%对混响信号20
4、01点的FFT变换羃 subplot(2,1,2);plot(n1(1:1000),Y1(1:1000);title(混响的频谱图);% 混响频谱图soun d(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;sub
5、plot(2,1,1);plot(y2);title(单回声滤波器时域图蚈 Y2=fft(y2,2001);%对单回声信号做2001点的FFT变换%单回声信号频谱图薆 subplot(2,1,2);plot (n 1(1:1000),Y2(1:1000);袇 Bz1=a,0,0,0,0,0,0,0,0,0,1;%分子的系数羀 title( 单回声滤波器频谱图 );肀 sound(5*y2,fs,bits);%回放经单回声滤波器后的语音信号螇%无限回声滤波器羆x,fs,bits=wavread(di ng.wav,10000 40000); % 读取语音信号蚀 x=x(:,1);%只取单声道袈
6、n1=0:2000;袅 a=0.5;%a 取小于等于 1莅 Bz=0,0,0,0,0,0,0,0,0,0,1;%分子的系数蒁 Az=1,0,0,0,0,0,0,0,0,0,-a;%分母的系数罿 yy1=filter(Bz,Az,x);%FIR 滤波器进行滤波芈 figure(5);subplot(2,1,1);plot(yy1);%无限回声滤波器时域波形袄title(无限个回声滤波器时域波形);賺 YY仁 fft(yy1,2001);%经无限回声滤波器后的信号做2001 点的 FFT 变换羀 subplot(2,1,2);plot (n 1(1:1000),YY1(1:1000);%无限回声滤
7、波器频谱图莆title(无限个回声滤波器频谱图);芄 sound(5*yy1,fs,bits);%回放经无限回声滤波器后的语音信号羂 pause(2);%延时 2 秒螈%全通结构的混响器 螈x,fs,bits=wavread(di ng.wav,10000 40000); % 读取语音信号蚃 x=x(:,1);%只取单声道蚂 n1=0:2000;莆 Az1=1,0,0,0,0,0,0,0,0,0,a;%分母的系数莂 yy2=filter(Bz1,Az1,x); figure(6);%FIR 滤波器进行滤波羁 subplot(2,1,1);plot(yy2); title( 全 通 应 的 时
8、域 图 );% 全 通 结 构 的 混 响 器 时 域 波 形YY2=fft(yy2,2001);%经全通结构的混响器后的信号做2001 点的 FFT 变换subplot(2,1,2);plot(n1(1:1000),YY2(1:1000);% 全通结构的混响器频谱图羅 title( 全通应的频谱图 ); sound(5*yy2,fs,bits);%回放经全通结构的混响器后的语音信号以下无正文仅供个人用于学习、研究；不得用于商业用途。For personal use only in study and research; not for commercial use.仅供个人用于学习、研究；不
9、得用于商业用途Nur f u r den pers?nlichen f u r Studien, Forschung, zu kommerziellen Zwecken verwendet werden.Pour l e tude et la recherche uniquementa des fins personnelles; pasa des fins commerciales.仅供个人用于学习、研究；不得用于商业用途to刀 bKO gHrogefiko Topieno 刀 b3ymm 日ca6yqeHu耳 c 刀 egoB u Huetgo 刀冰HbiUCnO/B3OBaTBC 刃b KOMMepqeckuxue/HX.
