matlab语音信号采集与处理%加噪声fs=22050; x=wavread('E:\Windows XP.wav');f=fs*(0:511)/1024;Au=0.03;t=0:1/22050:(length(x)-1)/22050;d=[Au*cos(2*pi*5000*t)];x2=x+d;sound(x2,8000);y2=fft(x2,1024);figure(1);plot(t,x2)title('加噪后的信号'
来源:学生作业帮助网 编辑:作业帮 时间:2024/06/28 17:43:33
![matlab语音信号采集与处理%加噪声fs=22050; x=wavread('E:\Windows XP.wav');f=fs*(0:511)/1024;Au=0.03;t=0:1/22050:(length(x)-1)/22050;d=[Au*cos(2*pi*5000*t)];x2=x+d;sound(x2,8000);y2=fft(x2,1024);figure(1);plot(t,x2)title('加噪后的信号'](/uploads/image/z/6606174-30-4.jpg?t=matlab%E8%AF%AD%E9%9F%B3%E4%BF%A1%E5%8F%B7%E9%87%87%E9%9B%86%E4%B8%8E%E5%A4%84%E7%90%86%25%E5%8A%A0%E5%99%AA%E5%A3%B0fs%3D22050%3B+x%3Dwavread%28%27E%3A%5CWindows+XP.wav%27%29%3Bf%3Dfs%2A%280%3A511%29%2F1024%3BAu%3D0.03%3Bt%3D0%3A1%2F22050%3A%28length%28x%29-1%29%2F22050%3Bd%3D%5BAu%2Acos%282%2Api%2A5000%2At%29%5D%3Bx2%3Dx%2Bd%3Bsound%28x2%2C8000%29%3By2%3Dfft%28x2%2C1024%29%3Bfigure%281%29%3Bplot%28t%2Cx2%29title%28%27%E5%8A%A0%E5%99%AA%E5%90%8E%E7%9A%84%E4%BF%A1%E5%8F%B7%27)
matlab语音信号采集与处理%加噪声fs=22050; x=wavread('E:\Windows XP.wav');f=fs*(0:511)/1024;Au=0.03;t=0:1/22050:(length(x)-1)/22050;d=[Au*cos(2*pi*5000*t)];x2=x+d;sound(x2,8000);y2=fft(x2,1024);figure(1);plot(t,x2)title('加噪后的信号'
matlab语音信号采集与处理
%加噪声
fs=22050;
x=wavread('E:\Windows XP.wav');
f=fs*(0:511)/1024;
Au=0.03;
t=0:1/22050:(length(x)-1)/22050;
d=[Au*cos(2*pi*5000*t)];
x2=x+d;
sound(x2,8000);
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');
为什么总是出现错误,信号不能叠加.
matlab语音信号采集与处理%加噪声fs=22050; x=wavread('E:\Windows XP.wav');f=fs*(0:511)/1024;Au=0.03;t=0:1/22050:(length(x)-1)/22050;d=[Au*cos(2*pi*5000*t)];x2=x+d;sound(x2,8000);y2=fft(x2,1024);figure(1);plot(t,x2)title('加噪后的信号'
有两处错了:加%的部分
fs=22050;
x=wavread('E:\1.wav');
f=fs*(0:511)/1024;
Au=0.03;
t=0:1/22050:(length(x)-1)/22050;
d=[Au*cos(2*pi*5000*t)];
x2=x+[d' d']; %原来的命令为x2=x+d;x为m*2的矩阵,d为1*m的向量不能直接相加修改见%左边.
sound(x2,8000);
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(y2(1:512))); %plot(f,abs(y1(1:512))); y1在前面从未出现过,需要定义,这里为了测试把它改成了y2
title('原始语音信号频谱');
xlabel('Hz');
ylabel('fuzhi');
subplot(2,1,2);
plot(f,abs(y2(1:512)));
title('加噪后的信号频谱');
xlabel('Hz');
ylabel('fuzhi');