%該段程序?qū)崿F(xiàn)將手機號碼語音x經(jīng)端點檢測出的第i個號的語音分幀提取LPC特征參數(shù)與碼本匹配，求出各幀所對應的索引號
%先進行端點檢測
function pipei(x)
%幅度歸一化到[-1,1]
x=double(x);
x=x/max(abs(x));
%常數(shù)設置
FrameLen=240;        %8k采樣，30ms對應240點
FrameInc=240;         %30ms對應楨移240點
amp1=2;
amp2=0.3;
zcr1=20;
zcr2=0.0005;
maxsilence=3; %3*30ms=90ms
minlen=6;	%6*30ms=180ms
status=0;
count=0;    %語音長度
silence=0;
%計算過零率
tmp1=enframe(x(1:length(x)-1),FrameLen,FrameInc);
tmp2=enframe(x(2:length(x)),FrameLen,FrameInc);
signs=(tmp1.*tmp2)<0;
diffs=(tmp1-tmp2)>0.02;
zcr=sum(signs.*diffs,2);
%計算短時能量
amp=sum(abs(enframe(filter([1-0.9375],1,x),FrameLen,FrameInc)),2);
%調(diào)整能量門限
amp1=min(amp1,max(amp)/4);
amp2=min(amp2,max(amp)/9);
d=zeros(22,1);
d(1)=1;
d(2)=1;
%開始端點檢測
for i=0:10
    status=0;
    count=0;
    for n=d(2*i+1):length(zcr)
        switch status
        case 0,					
            if amp(n)>amp2|zcr(n)>zcr2	%滿足此條件假定進入語音段
                count=count+1;
                status=1;
                silence=0;
            else				%靜音狀態(tài)
                status=0;
                count=0;
            end
        case 1,
            if amp(n)>amp2|zcr(n)>zcr2	%保持在語音段
                count=count+1;
                d(2*i+1)=max(n-count,1);
            else				%語音將結(jié)束
                silence=silence+1;
                if silence<maxsilence	%靜音還不夠長，尚未結(jié)束
                    count=count+1;
                elseif count<minlen	%語音長度太短，認為是噪聲
                    status=0;
                    silence=0;
                    count=0;
                else 			%語音結(jié)束
                    status=2;
                end
            end
        case 2,
            break;
        end
    end
    count=count-silence/2;
    d(2*i+2)=d(2*i+1)+count-0.5;
    d(2*i+3)=fix(d(2*i+2)+5);
end
%將手機號碼語音經(jīng)端點檢測出的第i個號的語音分幀提取LPC特征參數(shù)與碼本匹配，求出各幀所對應的索引號
load('index','i');
n=d(2*i+2)-d(2*i+1);
p=d(2*i+1)*240;
v1=x(p+1:p+240);
cx=lpc(v1,12)';
c3=cx(2:13);
p=p+240;
for i=1:n-1
    y=x(p+1:p+240);
    cx1=lpc(y,12)';
    c1=cx1(2:13);
    c3=horzcat(c3,c1);
    p=p+240;
end
load('codebook','cw');
d=disteu(c3,cw);
[m,index1]=min(d,[],2);
%對得到的索引號做一定的處理，減少匹配錯誤，方法是將連續(xù)相同的索引號中間出現(xiàn)跳變的索引號改為與前后索引號相同
if index1(1)~=index1(2);
    index1(2)=index1(1);
end
for i=2:length(index1)-1
    if index1(i-1)==index1(i+1)&&index1(i)~=index1(i+1)
        index1(i)=index1(i+1);
    end
end
l=length(index1);
if index1(l-1)~=index1(l)
    index1(l)=index1(l-1);
end
transindex=index1;
save('transmittedindex','transindex');