Matlab pitch period estimation post-processing

fid=fopen('zhouqi.txt', 'rt');
zhouqi=fscanf(fid, '%f');
fclose(fid);

zhouqi0=medfilt1(zhouqi, 5);
zhouqi1=medfilt1(zhouqi0, 3);
zhouqi2=linsmooth(zhouqi0, 5);

w=[];
w=zhouqi;
w1=w-zhouqi2;
w1=linsmooth(w1, 5);
zhouqi3=w1+zhouqi2;

v=[];
v(1)=0;v(2)=0;v(3)=0;v(4)=0;
for i=1:(length(zhouqi) - 4)
    v(i+4)=zhouqi(i);
end
%v(:,1) To get v The first column element in the matrix 
%v(1,:) To get v The first row of elements in the matrix 
%v(:) matrix v Convert to a row vector , It's a line 
% for example ：
%>> A = rand(2,3)
%
%A = 
%0.6925    0.3965    0.7802
%0.5567    0.0616    0.3376
%>> B = A(:)
%B =
%0.6925
%0.5567
%0.3965
%0.0616
%0.7802

v=v(:);
v1=v-zhouqi2;
v1=medfilt1(v1, 5);
v1=linsmooth(v1, 5);
zhouqi4=v1+zhouqi2;

figure(1);
subplot(511);
plot(zhouqi);
xlabel(' frames ');
ylabel(' Sample number ');
title(' Original pitch period trajectory ')

subplot(512);
plot(zhouqi1);
xlabel(' frames ');
ylabel(' Sample number ');
axis([0, 360, 0, 150]);
title(' Combination of five point median smoothing and three-point median smoothing ')

subplot(513);
plot(zhouqi2);
xlabel(' frames ');
ylabel(' Sample number ');
axis([0, 360, 0, 150]);
title(' A combination of five point median smoothing and five point linear smoothing ')

subplot(514);
plot(zhouqi3);
xlabel(' frames ');
ylabel(' Sample number ');
axis([0, 360, 0, 150]);
title(' Quadratic smoothing algorithm ')

subplot(515),plot(zhouqi4);
xlabel(' frames ')
ylabel(' Sample number ')
axis([0,360,0,150])
title(' Quadratic smoothing algorithm with delay ')

function [y]=linsmooth(x, n)
    %disp(size(x));% x  yes 360 That's ok ,1 Columns of the matrix 
    win=hann(n);
    win=win/sum(win); %  normalization 
    %disp(size(win));
    disp(win);
    [r,c]=size(x);%size(x) return  [ Row number , Number of columns ]
    if min(r,c)~=1 % ~= The representative is not equal to 
       error('sorry, no matrix here!:')
    end

    if r==1 %  Row vector 
        len=c;
    else    %  Column vector 
       len=r;
       x=x.';%  Transpose the column vector to the row vector 
    end
    %disp(len);
    % Come here  x  It must be a row vector 
    y=zeros(len,1);
    if mod(n,2)==0 % For return  m  Divide  2  Remainder after 
       m=n/2;
       x = [ones(1,m)*x(1) x ones(1,m)*x(len)]';
    else
       m=(n-1)/2; % Experimental here n=5, m=2
       x = [ones(1,m)*x(1) x ones(1,m+1)*x(len)]';% Generate new array , In the array  x  Add in front 2 Elements , In the array  x  Back plus 3 Elements 
    end

    for k=1:len
        y(k) = win'*x(k:k+n-1); %  One row and five columns  *  Five rows and one column 
    end
end