function [t, x] = ondaquadra(A,Tp,d,T)
% input:
%  A  = ampiezza dell'onda quadra
%  Tp = periodo dell'onda quadra
%  d  = duty cicle dell'onda quadra d = 2D/Tp
%  T  = passo di campionamento 
% output:
%  t = tempi del periodo [0, Tp) campionati a passo T
%  x = onda quadra in [0, Tp)
%

% asse dei tempi
t = [0:T:Tp-1e-10]; 
% onda quadra
D = d*Tp/2;
x = A*(t<D)+A*(t>=Tp-D);