delay.m

本次的算法是相位周期展开聚类

%整数延时n个单位
function y=delay(x,n,d);%x为信号，n为整数延时，d是分数延时
%n=10;
%d=0.5;
%snr=10;
%x=randn(1,30);
%[sx,fs,nibits]=wavread('d:\aaa\ttt.wav');
%x=x';
%x=awgn(x,snr);
if d==0
    h=[zeros(1,n),1];
    y=conv(x,h);
    %c=c1(1:(length(x)-n));
    else 
    h=[zeros(1,n),1];
    c=conv(x,h);
    %c=c1(1:(length(x)-n));
   
    %c=a(1:length(x));
%信号a分数延时
%s=1,m=1,t=4;
%for i=1:t;
    %for k=1:t;
        %if i~=k;
            %s=s*(d-k);
           % m=m*(i-k);
            %h(i)=s/m;
       % end
    %end
%end
h1=lagrange(10,d);
%h1=[(d-1)*(d-2)*(d-3)*(d-4)/24,-d*(d-2)*(d-3)*(d-4)/6,d*(d-1)*(d-3)*(d-4)/4,d*(d-1)*(d-2)*(d-4)/6,d*(d-1)*(d-2)*(d-3)/24];
y=filter(h1,1,c);
end