[x,fs]=audioread('000.wav');
sound(x,fs)
n=length(x);
x_p=fft(x,n);
f=fs*(0:n/2-1)/n;
figure(1)
subplot(2,1,1);
plot(x);
title('加噪语音信号采样后的时域波形'); xlabel('时间轴')
ylabel('幅值A')
subplot(2,1,2);
plot(f,abs(x_p(1:n/2)));
title('加噪语音信号采样后的频谱图'); xlabel('频率Hz');
ylabel('频率幅值');
L=length(x);
noise=0.1*randn(L,2);
x_z=x+noise;
sound(x_z,fs)
n=length(x);
x_zp=fft(x_z,n);
f=fs*(0:n/2-1)/n;
figure(2)
subplot(2,1,1);
plot(x_z);
title('加噪语音信号时域波形'); xlabel('时间轴')
ylabel('幅值A')
subplot(2,1,2);
plot(f,abs(x_zp(1:n/2))); title('加噪语音信号频谱图'); xlabel('频率Hz');
ylabel('频率幅值');
L=length(x) ;
noise=0.1*randn(L,2)
x_z=x+noise;
sound(x_z,fs)
t=[0:0.0001:10];
f1=3800;%噪声信号频率
y1=0.5*sin(2*pi*f1*t);%噪声信号
y2=y1(1,1:5000);%取噪声信号的5000点长度（1，1:5000）因为y1为列向量y=x+y2';%先对噪声信号转置再把噪声信号加入原始信号
Y=fft(y2,4096);%求噪声信号频谱
magY=abs(Y);%把傅里叶变换后的复数值取模
subplot(2,2,1);plot(y2);title('噪声信号波形');%绘制波形
subplot(2,2,2);plot(f,magY(1:2048));title('噪声信号频谱');%绘制波形
subplot(2,2,3);plot(y);title('加噪后信号波形');%绘制波形
Y1=fft(y,4096);%对加噪后的信号做4096点FFT变换
magY1=abs(Y1);%把傅里叶变换后的复数值取模
subplot(2,2,4);plot(f,magY1(1:2048));title('加噪后信号频谱');%绘制波形
(cos(6000/fs*pi*n)+cos(10000/fs*pi*n)+cos(15000/fs*pi*n))*0.5
