怎么往信号上叠加噪音?

guofengji

高手帮看下:
fs=22050;
x1=wavread('d:\ding.wav');
f=fs*(0:511)/1024;
%t=0:1/22050:(size(x1)-1)/22050;　　%将所加噪声信号的点数调整到与原始信号相同
d=0.04*randn(size(x1));
%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');
错误显示:
??? sound(x2,22050);　　　　　　　　%播放加噪声后的语音信号
                    |
Error: Missing variable or function.
加噪前语音是每个电脑自带的ding.wav,噪音是余弦5kHz,或用 附件-娱乐-录音机 录为.wav再叠加也可以,我的邮箱是guofengji@gmail.com
谢谢了先.