function P = fir_polyfit(M, R, hp, Nb, offset)

Np = length(hp);
order = M - 1;

for k=1:Nb,
    P(k,:) = polyfit((0:R)/R, hp(offset+(k-1)*R+1:offset+k*R+1), order);
end;