function TZSB(y,a)
%ph1=zeros(1,5000-a);
%ph2=zeros(1,5000-a);
z=hilbert(y);
    for n=1:5000-a;
     %   temp0=0;temp1=0;
       % for k=n:n+a/2-1;
        %    temp0=temp0+z(k);
          %  temp1=temp1+z(k+a/2);
        %end
        %DHWT(n)=(abs(temp0-temp1))/sqrt(a); %#ok<AGROW>
         ph1(n+1)=abs(angle(z(n)));
         ph2(n+1)=ph1(n+1)+ph2(n);
    end
         ph0=ph2(5000-a+1)/(5000-a+1);
     for n=1:5000-a;
         ph4(n)=ph1(n)-ph0;%NL
         ph3(n)=ph4(n)*ph4(n);
     end
          for n=1:5000-a-1;
         ph(n+1)=ph4(n+1)+ph4(n);
         ph3(n+1)=ph3(n+1)+ph3(n);
          end
     ph=(ph4(5000-a)*ph4(5000-a))/((5000-a)*(5000-a));
canshu3=sqrt((ph3(5000-a)/(5000-a))-ph)

  %    h=5000-a;
  %  k=1:h-1;
  %  subplot(312);plot(k,ph1(k));
  %  l=1:h-1;
  %  subplot(313);plot(l,ph2(l));
 %grid on