修正自相关

修正自相关 fid=fopen('voice.txt','rt') b=fscanf(fid,'%f'); b1=b(1:320); N=160; A=[]; for k=1:160; sum=0; for m=1:N; sum=sum+b1(m)*b1(m+k-1); end A(k)=sum; end for k=1:160 A1(k)=A(k)/A(1); end figure(1) subplot(3,1,1) plot(A1); xlabel('延时k') ylabel('R(k)') length('N=160') axis([0,180,-1,1]); b2=b(1:240); N=120; B=[]; for k=1:120; sum=0; for m=1:N; sum=sum+b2(m)*b2(m+k-1); end B(k)=sum; end for k=1:120 B1(k)=B(k)/B(1); end figure(1) subplot(3,1,2) plot(B1); xlabel('延时k') ylabel('R(k)') length('N=120') axis([0,130,-1,1]); L=1; b3=b(1:160); N=60; C=[]; for k=1:60; sum=0; for m=1:N; sum=sum+b3(m)*b3(m+k-1); end C(k)=sum; end for k=1:60 C1(k)=C(k)/C(1); end figure(1) subplot(3,1,3) plot(C1); xlabel('延时k') ylabel('R(k)') length('N=60') axis([0,70,-1,1]);