function [t, x]=ondaquadra(A0,T0,d0,T)
% A0= ampiezza dell'onda quadra;
% T0= periodo dell'onda quadra;
% d0= duty cicle dell'onda quadra= 2*T1/T0 
%T=passo di campionamento T
% x=onda quadra in [0, T0)
passo=T;
t=[0:passo:T0]'; %asse dei tempi
%segnale x(t)
x=A0*ones(length(t),1);
x(d0*T0/2/passo:T0/passo-d0*T0/2/passo)=zeros(size(x(d0*T0/2/passo:T0/passo-d0*T0/2/passo)));