对语音信号进行及处理精编.docVIP

 PAGE \* MERGEFORMAT 17 设计目的 1.进一步巩固数字信号处理的基本概念、理论、分析方法和实现方法；使自身对信号的采集、处理、传输、显示和存储等有一个系统的掌握和理解； 2.增强应用Matlab语言编写数字信号处理的应用程序及分析、解决实际问题的能力； 3．培养自我学习的能力和对相关课程的兴趣； 二、设计过程 1、语音信号的采集 采样频率，也称为采样速度或者采样率，定义了每秒从连续信号中提取并组成离散信号的采样个数，它用赫兹（Hz）来表示。 采样位数可以理解为声卡处理声音的 HYPERLINK /view/194157.htm \t _blank 解析度。这个数值越大，解析度就越高，录制和回放的声音就越真实 采样定理又称奈奎斯特定理，在进行模拟/数字信号的转换过程中，当 HYPERLINK /view/82683.htm \t _blank 采样频率fs不小于信号中最高频率fm的2倍时，采样之后的数字信号完整地保留了原始信号中的信息，一般实际应用中保证采样频率为信号最高频率的5～10倍。 利用Windows下的录音机，录制了一段发出的声音，内容是“数字信号”，时间在3 s内。接着在D盘保存为WAV格式，然后在Matlab软件平台下．利用函数wavread对语音信号进行采样，并记录下了采样频率和采样点数，在这里我们还通过函数sound引入听到采样后自己所录的一段声音。 [x1,fs,bits]=wavread(E:\数字信号.wav); %读取语音信号的数据，赋给变量x1，返回频率fs 44100Hz，比特率为16 。 2 、语音信号的频谱分析 首先画出语音信号的时域波形； 程序段： x=x1(60001:1:120000); %截取原始信号60000个采样点 plot(x) %做截取原始信号的时域图形 title(原始语音采样后时域信号); xlabel(时间轴 n); ylabel(幅值 A); （2）然后用函数fft对语音号进行快速傅里叶变换，得到信号的频谱特性； y1=fft(x,6000); %对信号做N=6000点FFT变换 figure(2) subplot(2,1,1),plot(k,abs(y1)); title(|X(k)|); ylabel(幅度谱); subplot(2,1,2),plot(k,angle(y1)); title(arg|X(k)|); ylabel(相位谱); （3）产生高斯白噪声，并且对噪声进行一定的衰减，然后把噪声加到信号中，再次对信号进行频谱特性分析，从而加深对频谱特性的理解； d=randn(1,60000); %产生高斯白噪声 d=d/100; %对噪声进行衰减 x2=x+d; %加入高斯白噪声 3、设计数字滤波器 (1)IIR低通滤波器性能指标通带截止频，阻带截止频率，通带最大衰减，阻带最小衰减。 (2)FIR低通滤波器性能指标通带截止频率，阻带截止频率， 通带衰减1dB，阻带衰减100dB。 IIR高通滤波器的设计指标，，，阻带最小衰减，通带最大衰减。 （4）FIR高通滤波器的设计指标，，，阻带最小衰减，通带最大衰减。 用自己设计的各滤波器分别对采集的信号进行滤波，在Matlab中，FIR滤波器利用函数fftfilt对信号进行滤波，IIR滤波器利用函数filter对信号进行滤波。比较滤波前后语音信号的波形及频谱，在一个窗口同时画出滤波前后的波形及频谱。 在Matlab中，函数sound可以对声音进行回放。其调用格式：sound(x，fs，bits)；可以感觉滤波前后的声音有变化。 结果分析 1、原始语音信号采样后的时域波形及FFT变换后频谱 图 1.原始信号时域图 图 2.原始信号的频谱图 2、加入噪声后信号的时域波形及FFT变换后频谱 图 3.加入噪声后的时域信号 分析：由图可以看出加入噪声后有明显的不一样，运行sound(x2 ,fs);播放加入高斯噪声后的语音，可以明显听出兹兹的噪声 图 4.加入噪声后的频谱图 3、IIR滤波器及原始信号通过IIR后的时域波形频谱变化 图 5. IIR低通滤波器 图 6. IIR低通滤波器滤波前后时域波形 分析：经过滤波器后的信号和原始信号很近似 图 7. IIR低通滤波器滤波前后的频谱 分析：从图7可以看出，经过IIR低通滤波器滤波后将高频部分滤除了。 4、FIR滤波器及原始信号通过FIR后的时域波形频谱变化 图 8. FIR低通滤波器 图 9.FIR低通滤波器滤波前后时域波形 图 1