Matlab分析音频

来源：互联网发布：mac 207 编辑：程序博客网时间：2024/05/17 05:02

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');

0 0