preprocess1.m

该代码实现了am

clear all
filename='ring.wav'
[S00]=wavread(filename)
plot(S00)                                                    %display the waveform which loaded

n0=size(S00,1)
S01(1,1)=0
B=[1,-32735*2^(-15)]
A=[1,-1]
n=2:n0
S01(n,:)=filter(B,A,S00(n,:))                        %setover compensation

B=1
A=[1,-28180*2^(-15)]
S0(n,:)=filter(B,A,S01(n,:))                             %pre_aggravation

for k=1:n0/160
    S(k,:)=S0(1,1+(k-1)*160:k*160)
end                                                           %frame divide