function [H,h] = circul(hin,symmetric) if(symmetric == 1) h = [hin, 0, hin([length(hin):-1:2])]; else h = hin; end hrev = [h(1), h(length(h):-1:2)]; H = toeplitz(h,hrev);