语音信号的滤波与频谱分析.docVIP

语音信号的滤波与频谱分析 生物医学信号处理大作业 题目：语音信号的滤波与频谱分析 学生姓名 学院名称 精密仪器与光电子工程 专 业 学 号 一、实验目的 语音信号的滤波与频谱分析 录制自己的一段语音：“天津大学精密仪器与光电子工程学院生物医学工程X班XXX, College of precision instrument and opto-electronics engineering, biomedical engineering”，时间控制在15秒到30秒左右；利用wavread函数读入语言信号，记住采样频率。 二、实验过程 （1）求原始语音信号的特征频带：可以分别对一定时间间隔内，求功率谱（傅里叶变换结果取模的平方）并画出功率谱。 （2）根据语音信号频谱特点，设计FIR或IIR滤波器，分别画出滤波器幅频和相频特性曲线。说明滤波器特性参数。用设计的滤波器对信号滤波，画出滤波后时域波形。用sound函数回放语音信号。 （3）求出特征频段语音信号随时间变化的曲线（每隔0.05秒求一次功率谱，连接成曲线）。 （4）选做：语谱图：横轴为时间，纵轴为频率，灰度值大小表示功率谱值的大小。(提示，可以采用spectrogram函数) （1）读入语音文件并画出其时域波形和频域波形，实现加窗fft并求出其功率谱。 clc clear all; close all; ); %因为语音信号主要集中在低频段，所以这里只需要显示低频段即可，取（1:3000） end 加窗后的时域波形 0.80.60.40.20-0.2-0.4-0.6-0.8012345678x 1095 第1~2秒语音功率谱5000450040003500300025002000150010005000050010001500200025003000 5~6秒语音功率谱18x 104第16141210864200500100015002000250030007~8秒语音功率谱10x 105第9876543210050010001500200025003000第11~12秒语音功率谱12x 10510864200500100015002000250030004第15~16秒语音功率谱7x 106543210050010001500200025003000第3~4秒语音功率谱3.5x 10432.521.510.50050010001500200025003000 5第9~10秒语音功率谱7x 106543210050010001500200025003000 13~14秒语音功率谱6x 105第5432100500100015002000250030004第17~18秒语音功率谱2.5x 1021.510.50050010001500200025003000 频谱分析：人说话的频率基本集中在1200Hz以内的低频段，但我这里可以在高频段（2000Hz左右）可以观察到部分能量，这应该是在录音的时候电脑本身的噪音以及录音设备的误差造成 （2）根据功率谱线，大致可观察到语音信号在150~400和500~650之间有两个波峰，因此取这两个频宽作为我语音信号的特征频带，根据要求选择汉宁窗作为滤波器滤波。 1）汉宁窗滤波器的幅频特性显示 %汉宁窗 f1=150; %通带上下两个频率 f2=400; %带通滤波器的通带范围 w1=2*pi*f1/Fs; %归一化 w2=2*pi*f2/Fs; w=[w1,w2]; N=ceil(4*pi/(2*pi*200/Fs));%求出所需滤波器的阶数 Windows=HANNING(N); b1=FIR1(N-1,w, Windows); %带通滤波器 figure;freqz(b1,1,512);title(‘汉宁窗带通（150~400）滤波器的频率响应’); %数字滤波器频率响应 f1=500; f2=650; %带通滤波器的通带范围 w1=2*pi*f1/Fs; w2=2*p