数字信号处理实验

[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