timeshift.m

用b-splines变换的方法实现分数延迟滤波器matlab仿真程序

function [y]=timeshift(x,tau) 
% the time-shift tau must be [0,1]
h1=[1 -3 3 -1]*(tau)^3/6;
h2=[0 1 -2 1]*(tau)^2/2;
h3=[0 1 0 -1]*(tau)/2;
h4=[0 1 4 1]/6;
hs=h1+h2+h3+h4;
x=conv(hs,x);
h=[1 4 1 ];
r=roots(h);
y1=rfilter(x,r(1));
y2=filter(1,[1 -r(2)],y1);
y=6*y2;
y=y(2:end-2);
% Copyright 2008 J.T. Olkkonen and H. Olkkonen