梳妆滤波器实现

1、梳状滤波器实现摘要本文中根据实际使用产生回声的“梳状滤波器”进行MATLAB实现。分别从时 域和频域两个角度仿真模拟梳状滤波器的效果。时域采用“延时拼凑”的方法，直 观地产生回声效果。频域模拟典型梳状传输函数来实现。关键字梳状滤波器时域频域传输函数 MATLAB提出问题梳状滤波器它是由许多按一定频率间隔相同排列的通带和阻带，只让某些特定频率范围的信号通过。梳状滤波器其特性曲线象梳子一样，故称为梳状滤波器。那么如何在MATLAB环境下实现梳状滤波器呢？解决方案.差分方程仿真梳状滤波器回声可以认为是由原始声音衰减后的多个延迟叠加组成的,因此回声可以用延迟单元来生成。x(n)表示原始声音信号，x(n

2、-N)是原始信号在时间轴上 延时N后的信号，a为衰减系数，T为延迟周期，回声信 号yn =xn +axn -T + . +akxnkT。图解如右：因此，构建差分方程能从时域的角度实现梳状滤波的效果。.构建递归数字滤波器的传输函数按照书中门9。提供的传输函数H(C) = i 一。:进行试验取一段音频文件(hello.wav )作为输入xn,观察检测输出效果yn具体实现1.差分方程仿真MATLAB代码如下：x, fs = wavread(hello.wav);% sound(x, fs);a = 0.6; T = 0.2;y1 = filter(1, zeros(1,T*fs-1), a, 1,

3、x);% sound(y1, fs); wavwrite(y1, fs, hello.wav );y2 = filter(1, zeros(1,T*fs-1), a, zeros(1,T*fs-1), aA2, zeros(1,T*fs-1), aA3,1, x);y3 = filter(1, 1, zeros(1,T*fs-1), a, x);% sound(y3, fs); wavwrite(y3, fs, hello.wav );plot(y3, m);hold on; plot(y2, r);hold on; plot(y1, g);hold on; plot(x, b);*代码引入函

4、数解释：.y=filter(p,d,x)用来实现差分方程 d表示差分方程输出 y的系数，p表示输入x的系数，x表示输入序列。实现差分方程，先从简单的说起：filter(1,2,1,1,2,3,4,5),实现 yk=xk+2*xk-1y1=x1+2*0=1(x1之前状态都用 0)y2=x2+2*x1=2+2*1=4. zeros(M,N)生成M行N列的0矩阵. y,fs,nbits= wavread(wavFilel);y就是被载入的音频信号；fs是采样频率，比如说16000就是每秒16000次；nbit是采样精度，比如说16就是指16位精度的采样；.hello.wav为事先用Cool Edit

5、录制好的发音为“hello”的音频文件，并导入进 document/MATLAB 中*绘图结果：1, zeros(1,T*fs-1)i,-；-右的枳系数糠车飞行数为；一r i科3口小加来定义 x n，0 x n L 0 xn - 2 . ),-1cax n - m在 xn处有输入，在延时m 后又有增益为a的输入axn-m,就从时域角度组合起来了梳状滤波器的回声效果。其他同理。2.递归数字滤波器B=1 zeros(1,7);A =1 zeros(1,6)-0.5;H,w=freqz(B,A,400,whole);Hf=abs(H);Hx=angle(H);clffigure(1)plot(w,H

6、f)title(离散系统幅频特性曲线)岛髓质统HI期带打曲理x,fs=audioread(hello.wav)y=filter(B,A,x)MWB=1 zeros(1,6498);岛髓质统HI期带打曲理x,fs=audioread(hello.wav)y=filter(B,A,x)MWB=1 zeros(1,6498);A =1 zeros(1,6497)-0.5figure(2)plot(w,Hx)title(离散系统相频特性曲线sound(y,fs) sound(y,fs)解释：根据书中提供的梳状滤波器的传输函数，定义其在MATLAB中的表达方式，并用filter函数来进行，信号输入滤波器进行处理后输出的过程。总结上次教训是对函数