function [xop] = x_operator(data)

sz = size(data);
ii = (sz(1)/2);

xop = [-ii:(ii-1)]' * ([1:sz]*0+1);
