def resampleMatrix(Tso, Tsf, o, f):
    retval = numpy.zeros((f, o))
for i in range(f):
    retval[i] = numpy.arange(Tsf * i / Tso, Tsf * i / Tso - o, -1.0)
return numpy.sinc(retval)
resampling, interpolating matrix