matlab farrow,farrow结构滤波器matlab实现源码与仿真

clc;close all;clear

N = 29; % 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 (0,1)

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(1,L);      % impulse response vector

hvec=zeros(Nfil,L);  % impulse response coefficient matrix

magresp = zeros(Nfil,Npt);

phasdel = zeros(Nfil,Npt-1);

xvec=zeros(Nfil,1);     % fractional delay vector

P = 3; % polynomial order for FARROW structure (ca. 1-5)

C=zeros(P+1,N+1);      % polynomial coeff. matrix

wp = 0.85; % normalized bandwidth (0-1.0)

for i=1:Nfil

d=delay(i);

if d==0

d=d+0.0000001;   % a