duanshiguolv.m

编写程序

close all
a1=zeros(1,620);K=450;R=zeros(1,250);
duanshnl=0;maxx=0;speech=zeros(1,1000);
sp1=wavread('h:\语音信号处理\a.wav');%11000采样频率
sp=wavread('h:\语音信号处理\f.wav');%11000采样频率
L=mean(sp);sp=sp-L;L2=mean(sp1);sp1=sp1-L2;
%以下程序计算短时过零率
m1=0;m2=0;
for i=1:2000
    if sp1(i)>0
        spp1(i)=1;
    else
        spp1(i)=-1;
    end
end

for i=1:2000
    if sp(i)>0
        spp(i)=1;
    else
        spp(i)=-1;   
    end
end

for i=2:2000
    m1=m1+0.5*abs(spp1(i)-spp1(i-1));
    m2=m2+0.5*abs(spp(i)-spp(i-1));   %  m2  >>  m1
end

plot(sp1(1:2000),'r');
hold on
plot(sp(1:2000));