数字信号处理课设题目详细分解步骤.docVIP

数字信号处理课程详细步骤分解 语音（音乐）信号滤波去噪的选题 课题具体内容 语音信号的采集　　要求学生利用Windows下的录音机，录制，时间在 s左右。然后在Matlab软件平台下，利用函数wavread对语音信号进行采样，记住采样频率和采样点数。通过wavread函数的使用，学生很快理解了采样频率、采样位数等概念。e:\yuyin.wav); % 输入参数为文件的全路径和文件名，输出的第一个参数是每个样本的值，fs是生成该波形文件时的采样率，bits是波形文件每样本的编码位数。 sound(x,fs,bits); % 按指定的采样率和每样本编码位数回放 N=length(x); % 计算信号x的长度 fn=2100; % 单频噪声频率，此参数可改 t=0:1/fs:(N-1)/fs; % 计算时间范围，样本数除以采样频率 x=x; y=x+0.1*sin(fn*2*pi*t); sound(y,fs,bits); % 应该可以明显听出有尖锐的单频啸叫声 语音信号的频谱分析要求学生首先画出语音信号的时域波形；然后对语音号进行快速傅里叶变换，得到信号的频谱特性，从而加深学生对频谱特性的理解。设计数字滤波器和画出其频率响应给出各滤波器的性能指标：　，As＝dB，Ap＝1 dB要求学生用窗函数法双线性变换法设计上面要求的滤波器。在Matlab中，可以利用设计FIR滤波器，可以利用函数buttecheby1、cheby2和ellip设计IIR滤波器；利用Matlab中的函数freqz画出各滤波器的频率响应。wsd=fsd/fs*2*pi;wsu=fsu/fs*2*pi; M=ceil(6.2*pi/dw)+1; % 计算汉宁窗设计该滤波器时需要的阶数 n=0:M-1; % 定义时间范围 w_ham=hanning(M); % 产生M阶的汉宁窗 hd_bs=ideal_lp(wcd,M)+ideal_lp(pi,M)-ideal_lp(wcu,M); % 调用自编函数计算理想带阻滤波器的脉冲响应 h_bs=w_ham.*hd_bs; % 用窗口法计算实际滤波器脉冲响应 [db,mag,pha,grd,w]=freqz_m(h_bs,1); % 调用自编函数计算滤波器的频率特性 IIR滤波器的设计可参考教材 用滤波器对信号进行滤波要求学生用自己设计的各滤波器分别对采集的信号进行滤波，在Matlab中，FIR滤波器利用函数fftfilt对信号进行滤波，IIR滤波器利用函数filter对信号进行滤波。IIR滤波器比较滤波前后语音信号的波形及频谱回放语音信号在Matlab中，函数sound可以对声音进行回放。其调用格式：sound (fs,bits)；可以感觉滤波前后的声音有变化。