基础音频处理的Matlab代码整理

1、裁剪音频

sec1 = 0; %剪切开始的时间
sec2 = 4; %剪切结束的时间
[Y1,F]=audioread('summer.wav');

Y_new(1:4*F)=Y(((F*sec1+1):F*sec2),:);
audiowrite('Y_new.wav',Y_new,F);

2、产生一段白噪声

len = 10000000;
vwhiteMat = zeros(len, 1);  
for m= 1:1
    vwhiteMat(:,m) = randn(len, 1);
    sigma = sqrt(vwhiteMat(:,m).'*vwhiteMat(:,m)/len);
    vwhiteMat(:,m) = vwhiteMat(:,m)/sigma;
end
% save('white');
audiowrite('vwhite.wav',vwhiteMat,8000);

3、产生一段正弦信号

Fs = 24000; % 采样频率
T = 1;     % 时间长度
n = Fs*T;  % 采样点数
f = 500;   % 声音频率
y = sin(2*pi*f*T*linspace(0,1,n));
y = fft(y)';
loglog((1:24000),abs(y));

4、音乐频谱图与波形图

[y,Fs] = audioread('airplane.wav');

N = length(y);
y = y(:,1); % 这里假设声音是双声道，只取