[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;
文章来源地址https://www.toymoban.com/news/detail-494243.html
文章来源：https://www.toymoban.com/news/detail-494243.html
到了这里，关于污染信号频谱分析模块程序的文章就介绍完了。如果您还想了解更多内容，请在右上角搜索TOY模板网以前的文章或继续浏览下面的相关文章，希望大家以后多多支持TOY模板网！
