污染信号频谱分析模块程序

[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;