fe=podd(f); fe=podd(f,dim);
podd(f) returns the odd part of the periodic sequence f.
podd(f,dim) does the same along dimension dim.