基于MATLAB的音频处理

>> [y,fs,bits]=wavread('D:\Matlab 2012b\bin\冰雨.wav');
sound(y,fs);
n=length(y)



>>  y_p=fft(y,n);
f=fs*(0:1456128/2-1)/1456128;
figure(1)
subplot(2,1,1);
plot(y);
title('原始语音信号采样后时域波形')
xlabel('t')
ylabel('A')
subplot(2,1,2);
plot(f,abs(y_p(1:1456128/2)));
title('原始语音信号采样后频谱图');
xlabel('频率Hz');
ylabel('频率幅值');
L=length(y)



>> noise=0.2*randn(L,2);
y_z=y+noise;
sound(y_z,fs)
n=length(y);
n=length(y)