function xintp=intp_coef(m,n)
xintp=[floor(m/2)-1:-1:1,0:floor(m/2)].'*ones(1,n-1);
xintp=xintp+[ones(floor(m/2),1);-ones(floor(m/2),1)]*([1:n-1]/n);
tmp=-1*ones(1,floor(m/2));
tmp=[fliplr(tmp.^[0:floor(m/2)-1]),tmp.^[0:floor(m/2)-1]];
xintp=(tmp.'*sin(pi*[1:n-1]/n)/pi)./xintp;
x1=[zeros(floor(m/2)-1,1);1;zeros(floor(m/2),1)];
x2=flipud(x1);
xintp=[x1,xintp,x2];

% y=hamming(m*(n+1)).';
% z=zeros(m,n+1);
% for i=1:m
% 	z(i,:)=fliplr(y((i-1)*(n+1)+1:i*(n+1)));
% end
% xintp=z.*xintp;