基于MATLAB采集语音信号分析与处理

语音信号的采集和频谱分析：
[y，fs，bits]=wavread('voice'); %读取音频信息（双声道，16位，频率44100Hz）
sound(y，fs，bits); %回放该音频
Y=fft(y，4096); %进行傅立叶变换
subplot(211);
plot(y);
title('声音信号的波形');
subplot(212)
plot(abs(Y));
title('声音信号的频谱');

窗函数设计低通滤波器：
fp=1000;
fc=1200;
as=100;
ap=1;
fs=22000;
wp=2*fp/fs;
wc=2*fc/fs;
N=ceil((as-7.95)/(14.36*(wc-wp)/2))+1;
beta=0.1102*(as-8.7);
window=Kaiser(N+1，beta);
b=fir1(N，wc，window);
freqz(b，1，512，fs);
结果：

滤波：
[y，fs，bits]=wavread('voice');
d=filter(b，a，y);
D=fft(d);
subplot(211)
plot(d);
title('滤波后的声音波形')
subplot(212)
plot(abs(D))
title('滤波后的声音频谱')

回放：
sound(d，fs，bits)
与滤波之前相比，噪音明显降低了许多。

采样，
fft频谱分析;
randn() 叠加噪声。

1. 采集建议使用cooledit，直观好用。
2. 时域波形直接plot，频域波形可以plot(abs(fft()))和plot(angle(fft()))，分别是幅度和相位，时频域分析可以用spectrogram命令，具体参数可以help一下。
3. awgn可以按信噪比添加白噪声，或者可以按需求生成需要的噪声，然后相加
4. 除非噪声是窄带的，否则无法用线性系统实现。因为语音信号是基带信号，如果不限定噪声的频率范围，则有可能与语音信号重叠，需要较复杂的算法进行分离。