%整數(shù)延時n個單位
function y=delay(x,n,d,t);%x為信號，n為整數(shù)延時，d是分數(shù)延時，t是FIR濾波器的階數(shù)。
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分數(shù)延時
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);