function s = sirene(f,f1,beta,T,Fe)
%        s = sirene(f,f1,beta,T,Fe)
% s(t) = sin(2*pi*f*t+beta*sin(2*pi*f1*t)/f1);
%---------------------------------------------

dt = 1/Fe;
t = 0:dt:T;
n = length(t);

s = sin(2*pi*f*t+beta*sin(2*pi*f1*t)/f1);
s = .9*s/max(abs(s));

