function amostrador(t, r, Ts)

% sinal r, frequência máxima do sinal wr, período de amostragem desejado Ts
disc=0.001;
Ts=0.5; ws=2*pi/Ts;
t=0:disc:10; %20s
r=sin(2*pi/10*t);%+sin(3*t)+sin(1/2*t+pi/2); sinal com frequência wr=pi rad/ms --> Tr=2s;
cont=(t(end)/disc)/(Ts/disc);
Impulso=[]; 
for i=1:cont 
    Impulso=[Impulso 1 zeros(1,Ts/disc-1)]; 
end
Impulso=[Impulso 1];
rs=r.*Impulso; % sinal amostrado
plot(t,r,t,rs);
%   
Filtro=[ones(1,Ts/disc) zeros(1,length(t)-Ts/disc)];
y=conv(rs,Filtro);
plot(t,rs,t,y(1:length(t)));
