From what I understand, you want the falling edge of your unit step signal to coincide with that of the square wave.
If that is the case then the following fix should work:
unitstep(t>=3)=-1 % Instead of "unitstep(t>=1)=-1;".
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!