基于MATLAB有噪声语音信号处理
一设计内容
选择一个语音信号作为分析的对象，或录制一段各人自己的语音信号，对其进行频谱分析；利用MATLAB中的随机函数产生噪声加入到语音信号中，模仿语音信号被污染，并对其进行频谱分析；设计FIR数字滤波器，并对被噪声污染的语音信号进行滤波，分析滤波后信号的时域和频域特征，回放语音信号。
二设计的具体实现
1语音信号的时频分析
利用MATLAB中的“wavread”命令来读入(采集)语音信号，将它赋值给某一向量。再对其进行采样，记住采样频率和采样点数。下面介绍wavread 函数几种调用格式。
(1)y=wavread(file)
功能说明：读取file所规定的wav文件，返回采样值放在向量y中。
(2)[y,fs,nbits]=wavread(file)
功能说明：采样值放在向量y中，fs表示采样频率(hz)，nbits表示采样位数。
接下来，对语音信号kalong.wav进行采样。其程序如下：
[y,fs,nbits]=wavered (‘kalong’);把语音信号加载入Matlab 仿真软件平台中然后，画出语音信号的时域波形，再对语音信号进行频谱分析。MATLAB提供了快速傅里叶变换算法FFT计算DFT的函数fft,其调用格式如下：Xk=fft(xn,N)
参数xn为被变换的时域序列向量，N是DFT变换区间长度，当N大于xn 的长度时，fft函数自动在xn后面补零。当N小于xn的长度时，fft函数计算xn 的前N个元素，忽略其后面的元素。
在本次设计中，我们利用fft对语音信号进行快速傅里叶变换，就可以得到信号的频谱特性。其程序如下：
[x,fc,nbits]=wavread ('kalong.wav');
sound(x,fc,nbits); %回放语音信号
N=length(x); %求出语音信号的长度
t=0:1/fc:(size(x)-1)/fc;
X=abs(fft(x,N)); %傅里叶变换
X=X(1:N/2);
f=fc/N*(0:1:(N/2)-1);
figure(1)
subplot(2,1,1);
plot(t,x);
title('原始信号波形');xlabel('s');
subplot(2,1,2);
plot(f,abs(X));%可见声音信号频谱在0～4e3范围之内
