n=0:0.002:4;
y=sin(100*pi*n);
subplot(4,1,1);%做出圖像 
axis([0 length(y) -1 1]);
plot(y);
x=awgn(y,0.05);
subplot(4,1,2);%做出圖像 
plot(x);
axis([0 length(x) -1 1]);
%-------------------------------------------------------------------------- 
Framelen = 200;%幀長 
Framemove = 80;%幀移 
%--------------------------------------------------------------------------
A1 = enframe(x(1:end-1), Framelen, Framemove);%采樣信號中第一值到第N-1個信號值 
A2 = enframe(x(2:end) , Framelen, Framemove);%采樣信號中第二值到第N個信號值 
signs = (A1.*A2)<0;%signs為符號數組,用于存儲相鄰兩個采樣點的符號是否相同,即是否 
%穿過零電平 
W = abs(A1 -A2)>0.02;%W用來度量相鄰兩個采樣點之間的距離,如果距離大于門限0.02(經驗值), 
%則為1,否則為0. 
ZC = sum(signs.*W, 2);%單幀過零次數 
ZCR=ZC/(20*0.02);%短時過零率
subplot(4,1,3);%做出圖像 
plot(ZCR);%單體坐標 
axis([0 length(ZCR) 0 max(ZCR)+100])%坐標橫縱軸設置 
xlabel('ZCR');%坐標名