Graficas Señales Caudrada

Para la Señal original tendremos el código:
x = -0.5:0.01:0.5;
y = square(2*pi*.25*t,100);
plot(x,y)
tl = -1:0.01:-0.5;
yl = square(0*pi*.25*tl,0);
tr = 0.5:0.01:1;
yr =square(0*pi*.25*tl,0);
tg = [tl t tr];
yg = [yl y yr];
plot(tg,yg)
Usando este pulso de señal obtenemos el espectro de frecuencia complementado por el código anterior:
matlabFFT = figure;YfreqDomain = fft(yr);
grid
axis([0,100,0,120])
function [X,freq]=centeredFFT(x,Fs)
N=length(x);
if mod(N,2)==0
k=-N/2:N/2-1; % N even
else
k=-(N-1)/2:(N-1)/2; % N odd
end
T=N/Fs;freq=k/T;
y = square(2*pi*.25*t,100);
X=fft(x)/N;
X=fftshift(X);
[YfreqDomain,frequencyRange] = centeredFFT(y,Fs);
centeredFFT = figure;
function [X,freq]=positiveFFT(x,Fs)
N=length(x);k=0:N-1; %create a vector from 0 to N-1
T=N/Fs;
freq=k/T;
X=fft(x)/N;
cutOff = ceil(N/2);
X = X(1:cutOff);
freq = freq(1:cutOff);

Y obtenemos la siguiente grafica:

Empleado elmismo procedimiento pero variando el periodo de τ/T, variando τ y dejando a T constante se puede observar la variación de la señal en las siguientes graficas:
matlabFFT = figure;YfreqDomain = fft(yr);
grid
axis([0,100,0,120])
function [X,freq]=centeredFFT(x,Fs)
N=length(x);
if mod(N,2)==0
k=-N/2:N/2-1; % N even
else
k=-(N-1)/2:(N-1)/2; % N odd
end
T=N/Fs;freq=k/T;
y = square(2*pi*.25*t,125);
X=fft(x)/N;
X=fftshift(X);
[YfreqDomain,frequencyRange] = centeredFFT(y,Fs);
centeredFFT = figure;
function [X,freq]=positiveFFT(x,Fs)
N=length(x);k=0:N-1; %create a vector from 0 to N-1
T=N/Fs;
freq=k/T;
X=fft(x)/N;
cutOff = ceil(N/2);
X = X(1:cutOff);
freq = freq(1:cutOff);


matlabFFT =figure;
YfreqDomain = fft(yr);
grid
axis([0,100,0,120])
function [X,freq]=centeredFFT(x,Fs)
N=length(x);
if mod(N,2)==0
k=-N/2:N/2-1; % N even
else
k=-(N-1)/2:(N-1)/2; % N odd
end...