Matlab分析音频

clc;clear all;

%%%%% 读入文件 %%%%%

x = wavread('F:\4.wav',[100001,200000]);

y = wavread('F:\4.wav',[200001,300000]);

z = wavread('F:\4.wav',[300001,400000]);

N=100000;

signal1=x(:,2);%获取右声道

signal2=y(:,2);

signal3=z(:,2);

figure(1);

subplot(311);plot(signal1);title('音频波形');xlabel('Time/s');

subplot(312);plot(signal2);xlabel('Time/s');

subplot(313);plot(signal3);xlabel('Time/s');

%%%%%% 均值 %%%%%%%

 figure(2);

signal4=[signal1,signal2,signal3];

signal5=mean(signal4,2);

plot(signal5);title('均值');xlabel('Time/s');

%%%%% 方差 %%%%%%

figure(3);

signal16=std(signal4,1,2);

signal13=(signal16).^2;

plot(signal13);title('方差');xlabel('Time/s');

%%%%% 自相关 %%%%

figure(4);

[signal7,N]=xcorr(signal1);

[signal8,N]=xcorr(signal2);

[signal9,N]=xcorr(signal3);

subplot(311);plot(N,signal7);title('自相关函数');xlabel('时间间隔/s');

subplot(312);plot(N,signal8);xlabel('时间间隔/s');

subplot(313);plot(N,signal9);xlabel('时间间隔/s');

%%%%% 互相关 %%%%

figure(5);

[signal10,N]=xcorr(signal1,signal2);

[signal11,N]=xcorr(signal2,signal3);

[signal12,N]=xcorr(signal3,signal1);

subplot(311);plot(N,signal10);title('互相关函数');xlabel('时间间隔/s');

subplot(312);plot(N,signal11);xlabel('时间间隔/s');

subplot(313);plot(N,signal12);xlabel('时间间隔/s');