Преобразование Фурье сигнала из файла *.BIN

clear all% Очистка памяти
 
%% Параметры
dat = fopen('2.BIN','r');
Signal=fread(dat); 
Tm=16777085/6.25e6;% Длина сигнала (с)
Fd=6.25e6;% Частота дискретизации (Гц)
FftL=1024;% Количество линий Фурье спектра
 
%% Генерация рабочих массивов
T=0:1/Fd:Tm;% Массив отсчетов времени
 
%% Спектральное представление сигнала
FftS=abs(fft(Signal,FftL));% Амплитуды преобразования Фурье сигнала
FftS=2*FftS./FftL;% Нормировка спектра по амплитуде
FftS(1)=FftS(1)/2;% Нормировка постоянной составляющей в спектре
 
%% Построение графиков
subplot(2,1,1);% Выбор области окна для построения
plot(T,Signal);% Построение сигнала
title('Сигнал');% Подпись графика
xlabel('Время (с)');% Подпись оси х графика
ylabel('Амплитуда');% Подпись оси у графика
 
F=0:Fd/FftL:Fd/2-Fd/FftL;% Массив частот вычисляемого спектра Фурье
figure% Создаем новое окно
subplot(2,1,1);% Выбор области окна для построения
plot(F,FftS(1:length(F)));% Построение спектра Фурье сигнала
title('Спектр сигнала');% Подпись графика
xlabel('Частота (Гц)');% Подпись оси х графика
ylabel('Амплитуда');% Подпись оси у графика