[X,fs,bits]=wavread('F:\情感语音库\录音\这不是我所预料的\女\高兴\he.wav');
N=fix(fs*0.015);
Xn=X(100+1:100+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);