Реализация фильтра Опции

[_Anatoliy]

Всем доброго времени суток.Есть задача - сделать фильтр согласованный с входным сигналом.Входной сигнал представляет собой радиоимпульс ЛЧМ с полосой 300МГц и длительностью 25 мкс. При оцифровке сигнала,например,с частотой 1ГГц получается окно для фильтрации 25000 отсчётов.Ясно что такой фильтр ни в какую FPGA не поместится,это если делать в лоб,во времеННой области.Как делать в частотной области - вообще не понятно,ведь время прихода сигнала не известно.Может есть какие хитрые приёмы реализации подобных фильтров?

[petrov]

Как делать в частотной области - вообще не понятно,ведь время прихода сигнала не известно.

Чем не устраивает быстрая свёртка? Коэффициенты фильтра можно считать на лету.

Может есть какие хитрые приёмы реализации подобных фильтров?

Хитрые приёмы связаны с хитрыми CAZAC последовтаельностями, посмотрите последовательности Chu если так ЛЧМ хочется.
Для некоторых последовательностей согласованные фильтры очень простые, например Golay Complementary Sequences.

[_Anatoliy]

Чем не устраивает быстрая свёртка? Коэффициенты фильтра можно считать на лету.
Хитрые приёмы связаны с хитрыми CAZAC последовтаельностями, посмотрите последовательности Chu если так ЛЧМ хочется.
Для некоторых последовательностей согласованные фильтры очень простые, например Golay Complementary Sequences.

Спасибо за ответ,начал поиск литературы.

Как я понял - быстрая свёртка это и есть то что я назвал в топовом сообщении "в частотной области". А здесь мне не понятен один момент:для FFT я набираю массив данных длиной 2^N,но может статься что в этом массиве окажется например 30% импульса ЛЧМ,следующие 70% попадут уже в следующий набор данных.Это как то скажется на качестве корреляции?Мне сдаётся что вместо одного пика на выходе согласованного фильтра я получу два пика(для первого и второго набора данных).

[alex_os]

Спасибо за ответ,начал поиск литературы.

Как я понял - быстрая свёртка это и есть то что я назвал в топовом сообщении "в частотной области". А здесь мне не понятен один момент:для FFT я набираю массив данных длиной 2^N,но может статься что в этом массиве окажется например 30% импульса ЛЧМ,следующие 70% попадут уже в следующий набор данных.Это как то скажется на качестве корреляции?Мне сдаётся что вместо одного пика на выходе согласованного фильтра я получу два пика(для первого и второго набора данных).

Нужно использовать результаты предидидущих блоков. Как вариант см. код ниже. Используется тот факт что свертка суммы последовательностей есть сумма сверток. Хорошо написано было про это кажется у Рабинера с Голдом.

Код

N = 13;
x = randn(1, 7*N);  % signal
h = randn(1,N);     % impulse response of a filter

yref = conv(x, h);

y = zeros(1, length(x)+length(h)-1);

% size of FFT is 2*N !!!
H = fft([h, zeros(1,N)]);
prev_hx = zeros(1, N-1);

for i=1:7
    
    X = fft([x((i-1)*N+1 : i*N), zeros(1,N)]);
    hx = ifft(H.*X);
    % прибавляем предидущий хвост
    hx(1:N-1) = hx(1:N-1) + prev_hx;
    prev_hx = hx(N+1:end-1);
    % запоминаем "хвост"  
    y((i-1)*N+1 : i*N) = hx(1:N);
end

figure(1);
% Результаты вычисленные двумя способами идентичны
plot(1:7*N, yref(1:7*N), 1:7*N, y(1:7*N))