interpolate.m

滤波器的设计以及零极点对应频率响应分析

%function is=interpolate(s)
% s is the signal

%*****************************
%******* this is for interpolate******
s=x0;
f=find(s==0);
l=length(f);
lt=l;
is=s;
i=2;
while i<l-10
    t=f(i+1)-f(i);
    if t==1
        j=i+1;
        while j<l-10
           if (f(j+1)-f(j))==1
                    j=j+1;
                    continue;
           end
           break;
           end
    if j<(length(s)-4)
       is(f(i):f(j))=linspace(is(f(i)-1),is(f(j)+1),(f(j)-f(i)+1));
        i=j;
    else
         i=j;
        break;
    end
    continue;
    end
    is(f(i))=(s(f(i)-1)+s(f(i))+s(f(i)+1))/2;
    i=i+1;  
end
 
s=is;
%********* this is for smooth*********
n=9;
 %n is the smooth parameters
ix(1:(n-1))=s(1:(n-1));
l=length(s);
for i=n:(l-n)
ix(i)=(s(i-4)+s(i+3)+s(i-2)+s(i-1)+s(i)+s(i+1)+s(i+2)+s(i-3)+s(i+4))/9;
end