数字信号处理课程设计1要点解析.pptVIP

数字信号处理课程设计 ——语音信号的去噪 董芳 朱浩然 张传祥 Contents 噪声信号 1 巴特沃兹滤波器 2 kaiser窗滤波器 3 滤波结果 4 Noise 高频余弦噪声 加性高斯白噪声 噪声信号 Contents 噪声信号 1 巴特沃兹滤波器 2 kaiser窗滤波器 3 滤波结果 4 Procedure 1 2 频率预畸 3 设置参数 Ωp=2*4500; Ωs=2*5500; Rp=3dB; As=50dB; 计算N值 N=ceil(log10((10^(As/10)-1)/(10^(Rp/10)-1))/(2*log10(OmegaS/OmegaP))) Procedure 4 5 6 归一化原型 双线性变换 去归一化 Butterworth Contents 噪声信号 1 巴特沃兹滤波器 2 kaiser窗滤波器 3 滤波结果 4 Procedure 1 设置参数 fp=4500;fc=5500;As=80;Ap=2; 2 模拟域转换成数字频率 wc=2π*fc/fs; wp=2π*fp/fs; 3 计算β 0.1102(α－0.87) α50 β= 0.5842(α－21)0.4+0.07886(α－21) 0 α21 Procedure 4 窗的长度点数N N=ceil((As-8)/2.285/wdel) 5 求kaiser函数 wn= kaiser(N+1,beta) 计算得β=7.9856相当于布莱克曼窗 Contents 噪声信号 1 巴特沃兹滤波器 2 kaiser窗滤波器 3 滤波结果 4 Results Appendix pause； 暂停后再执行后面的语句，使语音信号播放不会重叠 y=(:,1); 因为读取的语音为双声道与噪声的矩阵维度不一致，所以取y第一列其维度一致 Appendix OmegaP1=2*fs*tan(wp/2); OmegaS1=2*fs*tan(ws/2); 频率预畸：为了克服临界频率点的频率非线性畸变，必须采用频率预畸 [bz,az]=bilinear(b,a,fs) 双线性变换法求数字滤波器系统函数的系数 z=filter(bz,az,y_z); 一维数字滤波器 Appendix b=fir1(N,ws,wn); fir1是用窗函数法设计线性相位FIR DF，N为阶数，wc为截止频率，wn为指定的窗 x=fftfilt(b,y_z); 频域滤波，只对FIR滤波器有效 Thank you 数字信号处理课程设计