fir-delay.m

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

%整数延时n个单位
function y=delay(x,n,d,t);%x为信号，n为整数延时，d是分数延时，t是FIR滤波器的阶数。
x=[1,2,4,6,5,4,7];
n=2;
d=0.3;
t=4;
h=[zeros(1:n-1),1];
a=conv(x,h);
c=a(1:length(x));
%信号a分数延时
s=1,m=1;
for i=0:t;
    for k=0:t;
        if i~=k;
            s=s*(d-k);
            m=m*(i-k);
            h(i)=s/m;
        end;
        l=[h(0):h(t)];
        y=filter(l,1,c);
        plot(y);