clf;

n=0:10;
t=0:.01:10;

f1 = 1/4;
F1 = 1/4;
A1 = 1;
Th1 = pi/4;

x1 = A1*cos(2*pi*f1*n + Th1);
xa1 = A1*cos(2*pi*F1*t + Th1);

subplot(211)
stem(n,x1);
hold on
plot(t,xa1);
plot([0 10],[0 0]);   % Plot axis
grid
ti=title('CT and DT sinusoids');
set(ti,'FontSize',18);
yl=ylabel('x_1');
set(yl,'FontSize',14);
set(yl,'FontWeight','bold');

f2 = 5/4;
F2 = 5/4;
A2 = 1;
Th2 = pi/4;

x2 = A2*cos(2*pi*f2*n + Th2);
xa2 = A2*cos(2*pi*F2*t + Th2);

subplot(212)
stem(n,x2);
hold on
plot(t,xa2);
plot([0 10],[0 0]);   % Plot axis
grid
xl=xlabel('n (sample)');
set(xl,'FontSize',14);
set(xl,'FontWeight','bold');
yl=ylabel('x_2');
set(yl,'FontSize',14);
set(yl,'FontWeight','bold');

tt=text(9.8,-1.5,'8/25/00');
set(tt,'FontSize',8);
set(tt,'FontName','Courier');

print -dps two_freq_ct_dt.ps
