[y,fs,nbits]=wavread('OriSound.wav'); %语音信号采集
sound(y,fs,nbits); %回放语音信号便于比较效果
n = length (y) ; %计算语音信号长度
Noise=0.2*randn(n,1); %产生随机噪声信号 Noise
s=y+Noise; %将 Noise 添加到原始信号,得到污
染信号 s
sound(s); %回放污染信号 s
figure;
subplot(2,1,1); %绘制加噪信号时域波形
plot(s);
title('加噪语音信号的时域波形','fontweight','bold');
axis([ 00000 80000 -1 1]);
grid;
S=fft(s); %对 s 进行快速傅里叶变换得到频谱
subplot(2,1,2); %绘制加噪信号频域频谱
plot(abs(S));
title('加噪语音信号的时域波形','fontweight','bold');
axis([ 0 150000 0 4000]);
grid;
结果如下:
评论0
最新资源