Неадекватная частота спектра в fft

clear all% Очистка памяти
 
%% Параметры
fs=10% частота сигнала
%Tm=9;% Длина сигнала (с)
Fd=150;% Частота дискретизации (Гц)
FftL=256;% Количество линий Фурье спектра
ts=1/fs %их пер
td=1/Fd
 
%--Построение граффика одиночного импульса
NN=250;
 
time=0:td:NN*ts;% Массив отсчетов времени
n=1;
ch=0
x=0;
for pp=0:td:NN*ts;
x=x+td 
if (x < (2*ts/9))
y(n)=0; %0-2 0
end;
if (x >= (2*ts/9)) & (x<=(3*ts/9)) %2-3 1
y(n)=1;
end;
if (x > (3*ts/9)) & (x<(5*ts/9)) %3-5 0
y(n)=0;
end;
if (x >= (5*ts/9)) & (x<=(6*ts/9)) %5-6 1
y(n)=1;
end;
if (x > (6*ts/9)) & (x<=(9*ts/9)) %6-9 0
y(n)=0;
end;
n=n+1;
ch=ch+(td)
if (ch >= ts )
x=0+td; ch=0; end;
end;
 
subplot(4,1,1);
plot(time,y);
axis([0 9 -1 2]);
hold on;
xlabel('Время (с)');% Подпись оси х графика
ylabel('Амплитуда (Попугаи)');% Подпись оси у графика
title('Одиночный импульс');
 
pf=fft(y,FftL);
F1=0:Fd/FftL:Fd/2-1/FftL;% Массив частот вычисляемого спектра Фурье
pf=2*pf./FftL;% Нормировка спектра по амплитуде
pf(1)=pf(1)/2;% Нормировка постоянной составляющей в спектре
 
subplot(4,1,2);
stem(F1,abs(pf(1:length(F1))));% Построение спектра Фурье сигнала
title('Спектр последовательности MatLab');% Подпись графика