jz.m

来自「我的数字信号处理课程设计」· M 代码 · 共 26 行

fs=22050; 
x1=wavread('b.wav');
f=fs*(0:511)/1024;
t=0:1/22050:(size(x1)-1)/22050;　　%将所加噪声信号的点数调整到与原始信号相同
Au=0.03;　　　　　　　　　　　
d=[Au*cos(2*pi*5000*t)]';　　　　%噪声为5kHz的余弦信号
x2=x1+d;
sound(x2,22050);　　　　　　　　%播放加噪声后的语音信号
y2=fft(x2,1024);
figure(1)
plot(t,x2)
title('加噪后的信号');
xlabel('time n');
ylabel('fuzhi n');
figure(2)
subplot(2,1,1);
plot(f,abs(y1(1:512)));
title('原始语音信号频谱');
xlabel('Hz');
ylabel('fuzhi');
subplot(2,1,2);
plot(f,abs(y2(1:512)));
title('加噪后的信号频谱');
xlabel('Hz');
ylabel('fuzhi');