function B=flattop1(A,Tx)
u=ones(1,Tx); % define a single flat top pulse
B=kron(A,u); % perform kronecker product