使用脉冲响应不变法设计的切比雪夫I型IIR滤波器.ppt

1、使用脉冲响应不变方法设计并行Chebyshev I型滤波器，语音信号滤波器去噪，描述点，点1:记录语音信号和频谱分析点2:记录噪声声音信号的频谱分析点3:绘制设计滤波器并行结构点4:绘制滤波器并行结构图点5:噪声信号过滤和频谱分析，设计为脉冲响应不变的并行Chebyshev I 使用MATLAB的wavread(x，fs，bits)函数读取，绘制语音信号的波形，分析频谱特征，返回，录制的语音信号：点1，使用MATLAB的源程序见下图。 右：x，fs，bits=wav read(ldd . wav)；%输入参数是文件的完整路径和文件名，输出的第一个参数是“每个采样的值”，fs是生成波形文件时的采

2、样率，bits是波形文件中每个采样的编码位数。Sound(x、fs、bits)；% N=以指定的采样率和每个采样的编码位播放length(x)。%计算信号x的长度t=0333691/fs :(n-1)/fs；%计算时间范围，抽样数抽样频率x=xx=ABS(FFT(x)；%源和去噪信号的FFT转换，振幅谱X=X(13360n/2)；%裁剪前半部分deltaf=fs/N；%计算的光谱的光谱线间距f=0: deltaf : fs/2-del taf；%计算频谱频率范围，点2，内容：向原始语音信号添加单频率3000Hz的噪声信号，绘制频谱图，观察噪声去除前后的频谱变化，使用声音(y，fs，bits)函

3、数在添加噪声后播放语音对比前后的变化，返回，噪声后语音：通过Sound(y、fs、bits)播放此录音时，听到的是尖锐的单频率噪声。要点3，内容：数字滤波器的设计指标包括：通过带阻频率fn=3000，fp=fn-200=2800Hz，带阻频率fc=fn-50=2950HZ，带通波纹Rp=1，带阻波纹Rp=1fs=8000Fn=3000FP=fn-200；fc=fn-50；%定义通过频带和关闭频率Rp=1；As=16%定义通过频带波形和频带衰减WP=(FP/fs)* 2 * pi；ws=(fc/fs)* 2 * pi；%计算相应的数字频率T=1。%定义采样间隔ome gap=WP/T；ome g

4、as=ws/T；%截止频率线性变换Ripple=10(-Rp/20)；Attn=10(-As/20)；将%转换为绝对格式cs，ds=AFD _ chb1 (omegap，omegas，RP，as)；%计算过滤器顺序和截止频率b，a=imp_invr(cs，ds，T)；%脉冲响应不变方法是滤波系统函数分子分母系数db，mag，pha，grd，w=freqz_m(b，a)；%过滤器达到指定的性能delta=1，zeros(1，99)。Ha=filter(b，a，delta)；%计算脉冲响应，点3，程序运行结果过滤器响应也右：命令行上wp/pi，ws/pi验证过滤器是否满足设计标准：wp/pi 0.7000 ws/pi=0.7375，点4，内容：完成过滤器后，在表单命令行上C；返回，结果并行结构系数图：点4，可以根据并行结构系数和过滤器顺序绘制过滤器的结构图右侧：点5，内容：如果使用设计的过滤器过滤噪声信号，则过滤器调用方法如下：返回，y_fil=filter(b，a，y)；% IIR过滤器过滤信号y。其中b，a是上面设计的过滤器参数。y _ fil=ABS(FFT(y _ fil)；对于% FFT变换，振幅谱Y _ fil=Y _ fil(13360n/2)；%裁剪的前半部分，使用sound(y_fil，fs，bits)语句过滤的语音播放：要点5，右图表示通过运行程