speech.m

来自「举了一些语音信号处理的例子」· M 代码 · 共 61 行

files=dir('E:\MATLAB6p5\work\Readme\speech');
voice_file1='E:\MATLAB6p5\work\Readme\speech\h.wav';
voice_file2='E:\MATLAB6p5\work\Readme\speech\v.wav';

[x fs]=wavread('h');
figure(1)
subplot(2,2,1);
plot(x);
title('h.wav信号的波形');

z1=x(1:8000);
figure(1)
subplot(2,2,3)
plot(z1);
title('截取h.wav前8000个点的波形');

z2=x(8001:16000);
subplot(2,2,4);
figure(1)
plot(z2);
title('截取h.wav后8000个点的波形');

R1=corrcoef(z1,z2);   %计算二者相关性
subplot(2,2,2);
figure(1)
plot(R1);
title('相关曲线');

[y fs]=wavread('v');
figure(2);
subplot(2,2,1);
plot(y);
title('v.wav信号的波形');

z3=y(1:8000);
figure(2);
subplot(2,2,3);
plot(z3);
title('截取v.wav前8000个点的波形');

z1=x(1:8000);
figure(2)
subplot(2,2,4)
plot(z1);
title('截取h.wav前8000个点的波形');

R2=corrcoef(z3,z1);   %计算二者相关性
subplot(2,2,2);
figure(2)
plot(R2);
title('相关曲线');

covx=cov(x);       %计算协方差
pc=pcacov(covx);   %进行pca变换
figure(3);
subplot(2,1,1);
plot(pc);
title('x经pca变换后的波形');