function[]=AddNoise();
[x,fs,bits]=wavread('bluetooth_input.wav');%播放原始信號
N=length(x);%返回采樣點數(shù)
t=(1:N)/fs;
df=fs/N;%采樣間隔
n1=1:N/2;
f=(n1-1)*df;%頻帶寬度
y1=awgn(x,45);%加入高斯白噪聲，信噪比為45
wavplay(y1,fs);%播放加噪信號
figure(4);
plot(y1,'LineWidth',2),grid on;%信號的時域波形
title('加噪后信號的時域波形');
xlabel('時間/t');
ylabel('幅值/A');
y2=fft(y1);
figure(5);
plot(f,20*log10(abs(y2(n1)))),grid on;
title('加噪后信號的頻譜圖');
xlabel('頻率/Hz');
ylabel('幅值/db');
grid;