t=0:0.001:0.7;
tau=0.05;
num_realizz=10;
var_t0=0.03;
var_noise=0.015;
segnali=zeros(num_realizz,length(t));
for kk=1:num_realizz
    t0=normrnd(0.15,var_t0,1);
    if t0<0
        t0=0;
    end
noise=normrnd(0, var_noise,[1 length(t)]);
segnali(kk,:)=(t>t0).*(2*(t-t0).*exp(-(t-t0)/tau).*sin(2*pi*(t-t0)/0.22))+noise;

end

%sommo delle costanti in modo da rappresentare i segnali in una unica
%figura

costanti=0:0.1:0.1*(num_realizz-1);
costanti=repmat(costanti',1,length(t));
segnali_per_visualizz=segnali+costanti;
figure;
plot(t,segnali_per_visualizz)