应用matlab对语音信号进行频谱分析及滤波

《应用matlab对语音信号进行频谱分析及滤波》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、数字信号处理课程设计报告——用MatLab对语音信号进行频谱分析及滤波器设计学院：物联网工程学院班系：姓名：赵一畅QQ号：1320846809指导老师：13第一章实验任务录制一段个人自己的语音信号，并对录制的信号进行采样；画出采样后语音信号的时域波形和频谱图；给定滤波器的性能指标，采用窗函数法和双线性变换设计滤波器，并画出滤波器的频率响应；然后用自己设计的IIR滤波器和FIR滤波器对采集的信号进行滤波，画出滤波后信号的时域波形和频谱，并对滤波前后的信号进行对比，分析信号的变化；回放语音信号。第二章实验原

2、理2.1采样频率、位数及采样定理采样频率[2]，也称为采样速度或者采样率，定义了每秒从连续信号中提取并组成离散信号的采样个数，它用赫兹（Hz）来表示。采样频率的倒数是采样周期或者叫作采样时间，它是采样之间的时间间隔。通俗的讲采样频率是指计算机每秒钟采集多少个声音样本，是描述声音文件的音质、音调，衡量声卡、声音文件的质量标准。采样频率越高，即采样的间隔时间越短，则在单位时间内计算机得到的声音样本数据就越多，对声音波形的表示也越精确。采样位数可以理解为声卡处理声音的解析度。这个数值越大，解析度就越高，录制和

3、回放的声音就越真实。我们首先要知道：电脑中的声音文件是用数字0和1来表示的。所以在电脑上录音的本质就是把模拟声音信号转换成数字信号。反之，在播放时则是把数字信号还原成模拟声音信号输出。采样定理又称奈奎斯特定理[2]，在进行模拟/数字信号的转换过程中，当采样频率fs不小于信号中最高频率fm的2倍时，采样之后的数字信号完整地保留了原始信号中的信息，一般实际应用中保证采样频率为信号最高频率的5～10倍。2.2时域信号的FFT分析信号的频谱分析就是计算机信号的傅里叶变换[5]。连续信号与系统的傅里叶分析显然不便

4、于用计算机进行计算，使其应用受到限制。而FFT是一种时域和频域均离散化的变换，适合数值运算，成为用计算机分析离散信号和系统的有力工具。对连续信号和系统，可以通过时域采样，应用DFT进行近似谱分析。2.3IIR数字滤波器设计原理利用双线性变换设计IIR滤波器（巴特沃斯数字低通滤波器的设计）[6]，首先要设计出满足指标要求的模拟滤波器的传递函数Ha(s),然后由Ha(s)通过双线性变换可得所要设计13的IIR滤波器的系统函数H(z)。如果给定的指标为数字滤波器的指标，则首先要转换成模拟滤波器的技术指标，这里

5、主要是边界频率Wp和Ws的转换，对ap和as指标不作变化。边界频率的转换关系为∩=2/Ttan(w/2)。接着，按照模拟低通滤波器的技术指标根据相应设计公式求出滤波器的阶数N和3dB截止频率∩c；根据阶数N查巴特沃斯归一化低通滤波器参数表，得到归一化传输函数Ha(p)；最后，将p=s/∩c代入Ha(p)去归一，得到实际的模拟滤波器传输函数Ha(s)。之后，通过双线性变换法转换公式s=2/T((1-1/z)/(1+1/z))得到所要设计的IIR滤波器的系统函数H(z)。2.4FIR数字滤波器设计原理基于窗

6、函数的FIR数字滤波器的设计方法通常也称之为傅立叶级数法，是用一定宽度窗函数截取无限脉冲响应序列，获得有限长的脉冲响应序列，从而得到FIR滤波器。它是在时域进行的，由理想滤波器的频率响应推导出其单位冲激响应hd(n)，再设计一个FIR数字滤波器的单位冲激响应h(n)去逼近hd(n)，表示=由此得到的离散滤波器的系统传递函数Hd(z)为＝，该hd(n)为无限长序列，因此Hd(z)是物理不可实现的。为了使系统变为物理可实现的，且使实际的FIR滤波器频率响应尽可能逼近理想滤波器的频率响应，采用窗函数将无限脉冲

7、响应hd(n)截取一段h(n)来近似表示hd(n)，可得：h(n)=hd(n)w(n),从而有：式中N表示窗口长度，这样H(z)就是物理可实现的系统。并且从线性相位FIR滤波器的充要条件可知，为了获得线性相位FIR数字滤波器的冲激响应h(n)，那么序列h(n)应有τ=(N−1)/2的延迟。由于窗函数的选择对结果起着重要的作用，针对不同的信号和不同的处理目的来确定窗函数的选择才能收到良好的效果。第三章实验过程一、原始语音信号采样后的时域、频域分析1、源程序如下：时域波形：z1=wavread('F:St

8、udy数字信号处理实验Qnylan.wav')sound(z1,22050);%播放语音信号figure(1)plot(z1);title('时域图');13xlabel('时间（n）');ylabel('幅值（n）');%图见3-1频域波形：z2=wavread('F:Study数字信号处理实验Qnylan.wav');y1=z2(1:8192);y1=fft(y1);n=0:8191;sound(z1,22050);%播放