声音信号的采集、加噪,再滤波
[y,Fs,bits]=wavread('D:\q.wav');%读出信号,采样率和采样位数。
y=y(:,1);%取单声道作分析
yl=length(y)%求语音信号长度
yy=fft(y,yl);%傅里叶变换
t=[0:1/8000:4 zeros(1,yl-32001)]';
m=0.07*sin(10000*pi*t);%产生噪声
n=y+m;%加入噪声
nl=length(n)%求语音信号长度
nn=fft(n,nl);%傅里叶变换
figure(1);
subplot(2,1,1);plot(n);title('噪声信号波形')
subplot(2,1,2);plot(y);title('原信号波形')
figure(2);
subplot(2,1,1);plot(abs(nn));title('噪声信号频谱');
subplot(2,1,2);plot(abs(yy));title('原信号频谱');
sound(n,fs)
1