clear
T1=1/7;
T2=1/6;
w1=2*pi/T1;
w2=2*pi/T2;
t=0:0.01:2;
subplot(3,1,1)
plot(t,cos(w1*t),'b')
hold on
plot(t,cos(w2*t),'r')
subplot(3,1,2)
plot(t,abs(cos(w1*t)+cos(w2*t)),'b')
ylabel('|cosw_1t+cosw_2t|')
hold on
plot(t,abs(2*cos((w1-w2)*t/2)),'r')
subplot(3,1,3)
plot(t,cos(w1*t).*cos(w2*t),'b')
hold on
plot(t,0.5*cos((w1-w2)*t),'r')
ylabel('cosw_1t\timescosw_2t')