clear

% *********************************************************************
T   = 1;
w0  = 2*pi/T;

d   = .001;
t   = 0:d:T;
phi = pi/2;
%phi = 0:d:2*pi;

% *********************************************************************
A1 = 10;
A2 = A1*1;
s1_t = A1*cos(w0*t);
s2_t = A2*cos(w0*t + phi);

r_t = s1_t + s2_t;
lrt = length(r_t);

for i=(1:lrt)
    deg = mod(i+floor(lrt/4),lrt)
    if r_t(i) > r_t(deg)
        x_t(i) = r_t(i);
    else
        x_t(i) = r_t(deg);
    end
end

plot(t,abs(x_t));