Farrow结构设计分数延时滤波器

clc;close all;clear

N = 23; % filter order odd better
L = N+1;             % filter length;
Npt = 256;           % no. of frequency points for plots
w = （0:1:Npt-1）/Npt; % frequenc scan （01）

delay = [0 0.1 0.2 0.3 0.4 0.5];   % delay range x=0..0.5
Nfil = length（delay）; % number of filters

h = zeros（1L）;      % impulse response vector
hvec=zeros（NfilL）;  % impulse response coefficient matrix
magresp = zeros（NfilNpt）; 
phasdel = zeros（NfilNpt-1）;
xvec=zeros（Nfil1）;     % fractional delay vector

P = 2; % polynomial order for FARROW structure （ca. 1-5）
C=zeros（P+1N+1）;      % polynomial coeff. matrix

wp = 0.8; % normalized bandwidth （0-1.0）

for i=1:Nfil
    d=delay（i）;
    if d==0
        d=d+0.0000001;   % avoid sin（0）/0;
    end
    xvec（i）=d;
    if