t = linspace(0,20,20*60);

wn=2*pi/10;

a = 75*(sin(wn*t)+sin(4*wn*t))+228;
b = 0;
c = 75*(cos(wn*t)+cos(4*wn*t))+69;

xd = [a b c]';

plot(a,c)