amp = 0.5*pi;
bias = pi;
T = [0:0.05:40]';
n = size(T,1);
TETA = amp*prbs(n,80,120,0.8,1)+bias;
plot(T,TETA);
