数字信号处理告.docVIP

课 程 设 计 报 告 课程名称 数字信号处理课程设计 课题名称 数字滤波器设计及在语音信号分析中的应用 专 业 信息与计算科学 班 级 0801班 学 号 200810010112 姓 名 吴丽萍 指导教师 彭 祯，颜国风，曹炎炎 2010年 12 月 湖南工程学院 课 程 设 计 任 务 书 课程名称 数字 课 题 专业班级 0801 学生姓名 学 号 指导老师 彭祯 审 批 任务书下达日期 20 年 月 日 任务完成日期 20 年 月 日 y=wavread(‘wuliping.wav’); [y,fs,n]=wavread(‘wuliping.wav’);对语音信号采样。fs是采样频率，n是采样点数。 sound(y,fs);听采样后的语音 yk=fft(y); yk_abs=abs(yk);angle(yk); 8.2、含噪语音信号的合成 （1） 加白噪音代码如下： n1=randn(592704,2); yn=y+0.01*n1; plot(yn);title(‘加噪后的信号’) 加白噪后信号 （2）加单频噪音代码如下： t=0:1:592703; n2=[sin(2*pi*592704*t)]; ym=y(:,2)+50*n2; 加噪后的语音信号 8.3、数字滤波器的设计及滤波 （1）双线性变换法设计带阻滤波器 Fs=100; fp1=2; fp2=45; fc1=5; fc2=44; wpz=[2*pi*fp1/Fs,2*pi*fp2/Fs] wsz=[2*pi*fc1/Fs,2*pi*fc2/Fs] wp=2*Fs*tan(wpz/2); ws=2*Fs*tan(wsz/2); rp=10; rs=35; [N,wc]=buttord(wp,ws,rp,rs,s); [B,A]=butter(N,wc,s); % freqs(B,A); [Bz,Az]=bilinear(B,A,Fs); figure(1); freqz(Bz,Az,512,Fs); x=filter(Bz,Az,yn); X=fft(x); Y=fft(yn); subplot(2,2,1),plot(yn); title(滤波前信号波形); subplot(2,2,2),plot(x); title(滤波后信号波形); subplot(2,2,3),plot(abs(Y)); title(滤波前信号频谱); subplot(2,2,4),plot(abs(X)); title(滤波后信号频谱); （2）设计巴特沃斯低通数字滤波器，通带截至频率Wp=HZ，阻带截至频率为Ws=HZ，通带衰减不大于0.5DB，阻带衰减不小于40DB，抽样频Fs=HZ。 代码为： fp = ; fs = 58000; Fs = 44000; wp = 2*pi*fp/Fs; ws = 2*pi*fs/Fs; wp = tan(wp/2); ws = tan(ws/2); % 通带最大衰减为0.5dB，阻带最小衰减为40dB [N, wn] = buttord(wp, ws, 0.5, 40, s); % 模拟低通滤波器极零点 [z, p, k] = buttap(N); % 由极零点获得转移函数参数 [b, a] = zp2tf(z, p, k); % 由原型滤波器获得实际低通滤波器 [B, A] = lp2lp(b, a, wp); [bz, az] = bilinear(B, A, .5); [h, w] = freqz(bz, az, 256, Fs); figure plot(w, abs(h)) grid on ym_butter=filter(w,abs(h),ym); plot(ym_butter); sound(10*ym_butter); Butterworth低通