帮我看看这程序是何种算法？解释下 谢谢

gumuxiao

[X,fs,bits]=wavread('F:\情感语音库\录音\这不是我所预料的\女\高兴\he.wav');
N=fix(fs*0.015);
Xn=X(761+1:761+N);%取韵母段某一帧
Xn1=Xn;
Xmax=max(Xn);       %求一帧信号最大值
for i=1:N
   if((Xn(i)>-0.38*Xmax)&(Xn(i)<0.38*Xmax))
      Xn(i)=0;
   end
   if(Xn(i)> 0.38*Xmax)
      Xn(i)=Xn(i)-0.38*Xmax;
   end
   if(Xn(i)<-0.38*Xmax)
      Xn(i)=Xn(i)+0.38*Xmax;
  end
end  

%求削波前的AMDF
for k=1:N
    F1(k)=0;
    for m=1:N-k
        F1(k)=F1(k)+abs(Xn1(m+k)-Xn1(m));
    end
end
%求削波后的AMDF
for k=1:N
    F2(k)=0;
    for m=1:N-k
        F2(k)=F2(k)+abs(Xn(m+k)-Xn(m));
    end
end
subplot(221);
plot(Xn1);
subplot(223);
plot(F1);
subplot(222);
plot(Xn);
subplot(224);
plot(F2);