DSP实验二语音信号分析与处理,南京理工大学紫金学院实验报告,信号与系统.docVIP

实验二 语音信号分析与处理 学号 姓名 注：1）此次实验作为《数字信号处理》课程实验成绩的重要依据，请同学们认真、独立完成，不得抄袭。 2）请在授课教师规定的时间内完成； 3）完成作业后，请以word格式保存，文件名为：学号+姓名 4）请通读全文，依据第2及第3 两部分内容，认真填写第4部分所需的实验数据，并给出程序内容。 1. 实验目的 (1) 学会MATLAB的使用，掌握MATLAB的程序设计方法 (2) 掌握在windows环境下语音信号采集的方法 (3) 掌握MATLAB设计FIR和IIR滤波器的方法及应用 (4) 学会用MATLAB对语音信号的分析与处理方法 2. 实验内容 录制一段自己的语音信号，对录制的语音信号进行采样，画出采样后语音信号的时域波形和频谱图，确定语音信号的频带范围；使用MATLAB产生白噪声信号模拟语音信号在处理过程中的加性噪声并与语音信号进行叠加，画出受污染语音信号的时域波形和频谱图；采用双线性法设计出IIR滤波器和窗函数法设计出FIR滤波器，画出滤波器的频响特性图；用自己设计的这两种滤波器分别对受污染的语音信号进行滤波，画出滤波后语音信号的时域波形和频谱图；对滤波前后的语音信号进行时域波形和频谱图的对比，分析信号的变化；回放语音信号，感觉与原始语音的不同。 3. 实验步骤 1）语音信号的采集与回放 利用windous下的录音机或其他软件录制一段自己的语音（规定：语音内容为自己的名字，以wav格式保存，如wql.wav），时间控制再2秒之内，利用MATLAB提供的函数wavread对语音信号进行采样，提供sound函数对语音信号进行回放。 [y,fs,nbits]=wavread(file), 采样值放在向量y中，fs表示采样频率nbits表示采样位数。Wavread的更多用法请使用help命令自行查询。 2）语音信号的频谱分析 利用fft函数对信号进行频谱分析 3）受白噪声干扰的语音信号的产生与频谱分析 4）据语音信号的频带情况，设计FIR和IIR两种滤波器 5）用滤波器对受污染语音信号进行滤波 FIR滤波器fftfilt函数对信号进行滤波，IIR滤波器用filter函数对信号进行滤波 6）比较滤波前后信号的波形与频谱 7）回放滤波后的语音信号 4. 实验数据及实验程序 实验数据 1）原始语音信号的时域波形和频谱图 带限白噪声信号的时域波形和幅频特性 受污染语音信号的时域波形和频谱图 4）滤波器的频响特性图 FIR滤波器的频响特性图 IIR滤波器的频响特性图 5）滤波后语音信号的时域波形和频谱图 6）滤波前后的语音信号时域波形对比图和频谱对比图 实验程序： 实验主程序 clc;clear;close %[x,fs,bits]=wavread(录音的名字); [x,fs,bits]=wavread(录音的名字.wav); %x=wavread(录音的名字); sound(x,fs,bits); plot(x);title(声音波形图) figure(2) %y=fft(x,4096); y=fft(x); plot(abs(y));title(声音频谱图) fh=0.9;fl=0.25;n=1;length(x); y1=fh*sinc(fh*(n-5000))-fl*sinc(fl*(n-5000)); figure(5) subplot(1,2,1);plot(y1);title(窄带噪声波形图) y2=fft(y1); subplot(1,2,2);plot(abs(y2));title(窄带噪声频谱图) y3=y1+x; figure(6) subplot(1,2,1);plot(y1+x);title(混合信号波形图) y4=fft(y3); subplot(1,2,2);plot(abs(y4));title(混合信号频谱图) wp=0.5;ws=0.55; detaw=ws-wp; n=ceil(8*pi/detaw); wc=(wp+ws)/2; b1=fir1(n-1,wc/pi,hanning(n)); freqz(b1,1,41856) f1=fftfilt(b1,y3); plot(f1) f2=fft(f1); plot(abs(f2)) f11=filter (bz,az, y3); figure(8) subplot(1,2,1);plot(f11);title(恢复信号波形图) f22=fft(f11); subplot(1,2,2);plot(abs(f22));title(恢复信号频谱图) sound(f11,fs,bits); figu