afni.nimh.nih.gov

One possibility is to create a regressor that specifies the square-wave (e.g., 1s at the time points of square-wave and 0s for the rest) and then feed it to the option -stim_file in 3dDeconvolve.

Gang