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