See discussions, stats, and author profiles for this publication at: https://www.researchgate.

net/publication/269303397

FPGA-based match filter implementation in frequency domain using an

overlap-add method

Conference Paper · April 2014

DOI: 10.1109/SIU.2014.6830470

CITATIONS READS

5 781

5 authors, including:

Adnan ORDUYILMAZ Gokhan Kara

The Scientific & Technological Research Council of Turkey The Scientific & Technological Research Council of Turkey
41 PUBLICATIONS 211 CITATIONS 8 PUBLICATIONS 32 CITATIONS

SEE PROFILE SEE PROFILE

Mahmut Serin Alper Yildirim

The Scientific & Technological Research Council of Turkey Ankara Medipol Üniversitesi
27 PUBLICATIONS 147 CITATIONS 57 PUBLICATIONS 626 CITATIONS

SEE PROFILE SEE PROFILE

All content following this page was uploaded by Gokhan Kara on 05 December 2019.

The user has requested enhancement of the downloaded file.

 OVERLAP-ADD YÖNTEMİ KULLANILARAK FREKANS ALANINDA FPGA
TABANLI UYUMLU FİLTRE UYGULAMASI
FPGA-BASED MATCH FILTER IMPLEMENTATION IN FREQUENCY DOMAIN
USING AN OVERLAP-ADD METHOD

Adnan Orduyılmaz1, Gökhan Kara1, Mahmut Serin1, Alper Yıldırım1, Murat Efe2
1. TÜBİTAK-BİLGEM-İLTAREN 2. Elektrik- Elektronik Mühendisliği
{adnan.orduyilmaz,gokhan.kara,mahmut.serin,alper.yildirim} Bölümü
@tubitak.gov.tr Ankara Üniversitesi
[email protected]

ve dolayısıyla yeterli bellek kaynağının sağlanamamasıdır.

ÖZETÇE
Geleneksel yöntem işleyiş şeması Şekil 1’de sunulmuştur.
Bu bildiride, frekans alanında overlap-add yöntemi
kullanılarak gerçek zamanlı bir uyumlu filtre tasarımı
FGPA’de gerçeklenmiştir. Darbe sıkıştırma (DS) radarları ve
düşük olasılıklı algılanma (DOA) radarlarında alınan işaretin
gürültüye oranını (SNR) arttıran uyumlu filtre sayısal ortamda
tasarlanmıştır. Bu tasarım ara frekansta (2.5 GHz) örnekleme
yapabilen Xilinx Virtex 5 tabanlı işlemci kartlarında
uygulanmıştır. Overlap-add yönteminde paralel iki küçük FFT
bloğu kullanılarak uyumlu filtre tasarımı önerilmiştir. Farklı
darbe içi modülasyonları uygulanarak elde edilen uyumlu
filtre sonuçları sunulmuştur. Şekil 1 : Geleneksel yöntem

Bu çalışmada, overlap-add yöntemi kullanılarak radar almaç

ABSTRACT uyumlu filtresi gerçeklenmiştir. Bölüm 2’de overlap-add
yöntemi ve uyumlu filtre teorik olarak anlatılmıştır. Bölüm
In this research, a real time matched filter is implemented on
3’de FPGA üzerinde oluşturulan işaret oluşturma ve işleme
FPGA using an overlap-add method. The matched filter that
modülleri verilmiştir. Bölüm 4’de geliştirilen uyumlu filtre
increases the signal-to-noise ratio (SNR) for pulse
tasarımı farklı radar dalga formları için uygulanmış ve elde
compression and low probability intercept (LPI) radars is
sonuçlar Bölüm 5’de sunulmuştur.
implemented in digital domain. This design is implemented on
Xilinx Virtex-5 based processing board that samples in
2. TEORİ
intermediate frequency (2.5 GHz). In the overlap-add method,
we propose to design the matched filter by using two parallel Zaman alanında uyumlu filtre uygulaması radar darbe içi
FFT cores. Furthermore, the matched filter results are işareti örneklerinden oluşan filtre katsayıları kullanılarak
presented for different intra-pulse modulations. hedeften yansıyan işaret ile evrişimi ile gerçekleştirilir [3].
Eşitlik (1)’de zaman alanında uyumlu filtre uygulaması
sunulmuştur.
1. GİRİŞ N −1

Darbe sıkıştırma ve düşük olasılıklı algılanma radarlarının, y[n] = x[n]∗ h[n] = ∑ x[k ]h[n − k ] (1)
k =0
radar teknolojileri içinde en önemli özelliği, daha düşük bir
Burada y[n] uyumlu filtre çıkışını, x[n] filtrelenmesi istenen
tepe gücü ile nispeten daha uzun darbe süresi sayesinde verici
işaret dizisini, h[n] filtre katsayılarını ifade etmektedir.
ortalama gücünü yükseltebilmesi ve radar menzil
Uyumlu filtre frekans ekseninde iki işaretin frekans dönüşüm
çözünürlüğünü arttırabilmesidir [1]. Darbeli radarlarda, radar
katsayılarının çarpılıp, ters frekans dönüşümü yapılarak
almacında kullanılan uyumlu filtre uygulaması ortamda
gerçekleştirilir. Frekans alanında uyumlu filtre eşitlik (2)’de
bulunan gürültü ve kargaşayı bastırarak, maksimum kazancı
sunulmuştur.
sağlayıp hedeflerin daha rahat algılanabilmesi için geliştirilmiş
bir filtredir [1]. Uyumlu filtre uygulaması radar almacında y[n] = IFFT [ X [W ] × ( H [W ]) *]
analog ve sayısal tabanlı olarak gerçeklenebilmektedir [2].
= IFFT [FFT ( x[n]) × ( FFT (h[n])) *] (2)
Modern radar sistemlerinde uyumlu filtre sayısal tabanlı olarak
gerçeklenmektedir. Bu çalışmada sayısal tabanlı uyumlu filtre
uygulaması farklı radar dalga formları için uygulanmıştır. Burada X[W] ve H[W] sırasıyla x[n] ve h[n] işaret dizilerinin
Uyumlu filtrenin sayısal uygulamasında, radar ve hedef eko FFT’sini, y[n] uyumlu filtre çıkışını ifade etmektedir.
işaretlerinin FFT’si alınıp, birbiri ile çarpıldıktan sonra IFFT Darbe genişliği küçük olan radarlar için zaman alanında
dönüşümü gerçekleştirilir [1]. Bu yöntemde karşılaşılan en evrişim uygulanması daha kolay gerçeklenebilmektedir.
yaygın problem ise katsayılar arttıkça işlem sayısının artması, FPGA hazır IP çekirdekleri kullanılarak sınırlı sayıda filtre
katsayısı ile zaman alanında evrişim gerçekleştirilebilir [4].

978-1-4799-4874-1/14/$31.00 ©2014 IEEE

Fakat darbe genişliği arttıkça FPGA’de hem işlem yükü
artmakta, hem de filtre katsayı sınırı aşılmaktadır. Örneğin, N
örnek sayısına sahip bir sinyalin zaman alanında N2 işlem
yükü gerektirirken, frekans alanında (1/2)N[log2N] işlem
yükü gerekmektedir.

2.1. Overlap-Add Filtre

Overlap-add filtre metodunda, temel olarak sürekli gelen işaret
eşit boyuttaki parçalara ayırılır, her bir parçaya normal işaret
işleme algoritması uygulanır ve elde edilen sonuçların
birleştirilerek filtrelenmiş işaret elde edilir [5]. Şekil 2’de
overlap-add yönteminin zaman akış diyagramı sunulmuştur.

Şekil 3: FPGA’de uyumlu filtre tasarımı

3.1. Radar İşaret Simülatörü

Bu modül temel olarak darbeli bir radarın işareti evre uyumlu
Şekil 2: Overlap-add yöntemi. olarak ara frekansta üretebilmektedir [8]. Radar işaret
simülatörü farklı darbe içi modülasyona sahip dalga formları
Burada sürekli akan işaret N uzunluğunda parçalara (faz kodlu, doğrusal frekans modülasyonu (LFM), doğrusal
([x11,x12,...],[x21,x22,...]) ayrılmıştır. Daha sonra bu parçalara N olmayan frekans modülasyonu (NLFM) vb.) üretebilmektedir.
uzunluğunda sıfırlar eklenerek FFT’si alınır [6],[7]. FFT Örneklenmiş darbe içi işaret örnekleri öncelikle MATLAB
sonuçları filtre katsayıları ile çarpılarak IFFT alınıp, çıkışlar ortamında üretilir. Üretilen katsayılar FPGA’de bellek
([y11,y12,...],[y21,y22,...]) elde edilmiştir. İkinci çıkış işareti N bloklarına aktarılır. Belirlenen darbe tekrarlama aralığı (DTA)
uzunluğunda geciktirilerek, birinci ve ikinci işaret toplanır ve zamanında bellekten okunarak darbeli radar işareti üretilir.
sonuç işareti oluşturulur. Böylelikle overlap-add yöntemi
kullanılarak filtreleme sağlanır. 3.2. Hedef İşaret Simülatörü
Hedef işaret simülatöründe radar işaret simülatöründen alınan
3. UYGULAMA işaret bellek bloğuna kaydedilir. Daha sonra ara yüzden girilen
Uyumlu filtrenin kullanıldığı benzetim ortamında temel olarak menzil değerlerine göre bellek bloğundan darbe süresi
üç ana modül bulunmaktadır. Radar işaret simülatörü boyunca okunarak hedef işareti oluşturulur. Böylece belirli
modülünde ara frekans radar işaretleri üretilmektedir. Hedef menzillerde hedefler üretilir.
işaret simülatöründe hedef ekoları belirlenen sayıda ve menzil
3.3. Radar Sayısal Almaç
değerinde üretilmektedir. Radar sayısal almaç modülünde,
alınan işaret filtrelenmektedir. Şekil 3’de gerçeklenen uyumlu Radar sayısal almaç bloğunda öncelikle radar ve hedef
filtre tasarımının blok şeması verilmiştir. işaretleri 8 bit analog sayısal dönüştürücü (ASD) ile
örneklenir. Örneklenen işaretler frekans dönüşüm (FFT),
kompleks çarpım (X) ve ters frekans dönüşüm (IFFT)
bloklarından geçerek filtrelenir. Radar ve tehdit işaretlerinin
frekans dönüşümü overlap-add yöntemi uygulanarak Xilinx IP
çekirdekleri ile sağlanmaktadır. Alınan işaretlerin sürekli
şekilde işlenebilmesi için “Pipeline, Streaming I/O” FFT [4]
yöntemi uygulanmıştır.
Öncelikle radar işaretinin darbe genişliğinin iki katı boyutunda
FFT’si alınır. FFT sonucu bellek bloklarına kaydedilir. Radar
darbesinin FFT katsayıları, hedef işaretin FFT sonucuyla
senkronize şekilde bellek bloğundan her bir pencere boyunca
okunur. Yapılan çalışmada darbe genişliği 1000 ns (250 örnek)
seçilmiştir. Dolayısıyla frekans dönüşüm boyutu sıfırlar
eklendiği için minimum 512 olması gerekmektedir. Örnek bir
radar işareti, frekans dönüşüm sonucu ve bellek bloklarından 4. ANALİZLER
okunması Şekil 4‘de sunulmuştur.
Bu çalışmada uyumlu filtre uygulaması öncelikle bilgisayarda
MATLAB ve ModelSim ortamında benzetimi
gerçekleştirilmiş daha sonra Virtex-5 kartına aktarılmıştır.
Farklı darbe içi modülasyonlardan faz modülasyonu olarak
Barker kodları ve frekans modülasyonu olarak doğrusal
frekans modülasyonu seçilmiştir. Darbe içi modülasyon
olduğu durumlarda uyumlu filtrenin sağladığı işlem kazancı
benzetim sonuçlarında ve FPGA uygulamasında
hesaplanmıştır.
LFM modülasyonu için kullanılan radar parametre listesi
Tablo 1’de verilmiştir.
Şekil 4: Radar işareti ve frekans dönüşümü
Tablo 1: LFM radar parametre listesi
Hedef işaretinin frekans dönüşümü Şekil 2’de anlatıldığı üzere
Darbe Genişliği (PW) 2,5 µs
iki parçaya ayrılarak yapıldığı için, bellek bloklarından
Darbe Tekrarlama Aralığı (PRI) 100 µs
sonuçlar frekans dönüşüm boyutunun (NFFT) yarısı kadar
(NFFT/2) gecikmeli olarak oynatılmaktadır. Hedef işaretinin Bant Genişliği(BW) 40 MHz
overlap-add yöntemi uygulanarak yapılan frekans dönüşüm Ara Frekans (IF) 70 MHz
sonucu Zaman-Bant Genişliği Çarpımı (BT) 100
Şekil 5‘de sunulmuştur.
Bu değerler kullanılarak üretilen radar darbesi için uyumlu
filtre çıkışı benzetim ve gerçeklenme sonuçları sırasıyla Şekil
7 ve Şekil 8’de verilmiştir.

Şekil 5: Hedef işareti ve frekans dönüşümü

Bu durumda çıkan radar ve hedef işaretlerinin dönüşüm Şekil 7: LFM işareti uyumlu filtre benzetim sonucu
sonuçları çarpılacaktır. Zaman alanındaki evrişim frekans
alanındaki kompleks çarpmaya denk gelmektedir. Dolayısıyla
frekans dönüşüm sonuçlarının gerçek ve sanal (real and
imaginary) parçaları Eşitlik (3)’e göre çarpılacaktır.

X (W ) = a + bi ⇒ a = re{ X (W )}, b = im{ X (W )}

Y (W ) = c + di ⇒ c = re{Y (W )}, d = im{Y (W )} (3)
X (W )* Y (W ) = (a + bi ) × (c − di )

Denklem (3)’de X(W) radar işareti FFT katsayıları ve Y(W)

hedef işareti FFT katsayıları olarak kullanılmıştır. Böylece Şekil 8: LFM işareti uyumlu filtre gerçeklenme sonucu
hedef işareti frekans alanında filtrelenmiş olmaktadır.
Kompleks çarpım sonuçları ise ters frekans dönüşümü (IFFT) Şekil 7’de sunulan uyumlu filtre benzetim sonucunda yan-lob
yapılarak zaman alanında uyumlu filtre sonuçları elde seviyesi -13.4 dB olarak görülmektedir. Şekil 8’de ana-lob ve
edilmiştir. Darbe içi modülasyon bulunmayan radar darbe yan-lob voltaj seviyeleri ölçülerek ve -20log(Vana-lob/ Vyan-lob)
işareti için uyumlu filtre FPGA benzetim sonuçları Şekil 6‘de hesabıyla -13,38 dB yan-lob seviyesi hesaplanmıştır.
gösterilmiştir. Barker-13 modülasyonu için kullanılan radar parametre listesi
Tablo 2’de verilmiştir.
Tablo 2: Barker kod radar parametre listesi

Darbe Genişliği 2,6 µs

Darbe Tekrarlama Aralığı 100 µs
Ara Frekans 70 MHz
Şekil 6: Frekans alanındaki uyumlu filtre uygulaması Barker Kod Uzunluğu(N) 13
 Bu değerler kullanılarak üretilen radar darbesi için uyumlu Tablo 4: Kaynak Kullanımı
filtre çıkışı benzetim ve gerçeklenme sonuçları sırasıyla Şekil
9 ve Şekil 10’da verilmiştir.
Device Utilization Summary (estimated values) [-]
Avail Utiliz
Logic Utilization Used
able ation
Number of Slice Registers 29270 44800 65%
Number of Slice LUTs 27099 44800 60%
Number of fully used LUT-FF pairs 20904 35465 58%
Number of Block RAM/FIFO 20 148 13%
Number of BUFG/BUFGCTRLs 2 32 6%
Number of DSP48Es 102 128 79%

5. SONUÇ
Şekil 9 : Barker kodlu işareti uyumlu filtre benzetim Bu çalışmada, frekans alanında uyumlu filtre tasarımı overlap
sonucu add metodu kullanılarak FPGA’de gerçeklenmiştir, teorik
sonuçlar ile pratikte gerçekleştirilen sonuçların birbirine çok
yakın olduğu gözlemlenmiştir. Gerçeklenen tasarım farklı faz
ve frekans modülasyonları içeren radar darbeleri için
uygulanmıştır. Overlap-add metodu kullanılarak tasarlanan
uyumlu filtre sonucunda elde edilen ana-lob ve yan-lob voltaj
seviyelerinin teorik olarak hesaplanan değerler ile aynı olduğu
görülmüştür. Overlap-add metodu ile kullanılarak gelecek
çalışmalarda, daha fazla kaynak içeren FPGA kartlarında
tasarımın gerçeklenmesi ve tasarımın en iyilenmesi
hedeflenecektir. Overlap-add metodu kullanılarak bu
çalışmada büyük kaynak tüketen FFT bloklarının yerine
paralel küçük FFT blokları kullanılmış, kaynak tüketiminden
kazanç sağlanmıştır. Ayrıca kaynak tüketimi çok olan hazır
FFT IP çekirdekleri yerine daha az kaynak tüketecek FFT
bloğu tasarlanacaktır.

Şekil 10: Barker kodlu işaretin uyumlu filtre 6. KAYNAKÇA

gerçeklenme sonucu [1] Serin, M., Yılmaz, A.E., Navruz, I., “A comparison on
the new electronic attack techniques against pulse
Şekil 9’de sunulan uyumlu filtre benzetim sonucunda yan-lob compression radars”, Proc. of Signal Processing and
seviyesi -22,28 dB olarak görülmektedir. Şekil 10’de ana-lob Communications Applications Conference (SIU), 2010
ve yan-lob voltaj seviyeleri ölçülerek ve -20log(Vana-lob/ Vyan- IEEE 18th, Diyarbakir, 22-24 April 2010.
lob) hesabıyla -19,13 dB yan-lob seviyesi hesaplanmıştır. [1] Yu,J., Huang, D., Pei, N., Zhao, S., Guo, J., Xu, Y.,
Farklı zaman-bant genişliği çarpımı sahip LFM ve Barker “CORDIC-based design of matched filter weighted
kodları için benzetim ve gerçeklemeler uygulanmış olup elde algorithm for pulse compression system”, ICSP, 2012.
edilen sonuçlar Tablo 3’de sunulmuştur. [2] M. I. Skolnik, “Radar Handbook”, 3rd edition, McGraw
Tablo 3: Farklı darbe içi dalga formları için uyumlu Hill, 2008.
filtre sonuçları [3] Agrawal, N., “Multichannel Pulse Compression System
for SAR Data Based on Xilinx Virtex-7 FPGA”.
[4] Xilinx Inc. “LogiCORE Fast Fourier Transform v7.1”,
Benzetim Gerçekleme Matematiksel
2010.
dB dB hesaplama
dB [5] Smith, W., S., “The Science and Engineer’s Guide to
LFM(BT = 20) 14.47 12.61 13.4 Digital Sinal Processing”, 2. Edition, California
LFM(BT = 50) 13.55 15.24 13.4 Technical Publishing, 1999.
LFM(BT = 100) 13.43 13.56 13.4 [6] Özdil, Ö, İspir, M., Onat, E.,Yıldırım, A., “FPGA Tabanlı
Barker(N = 3) 9.542 8.91 9.54 Overlap-add Filtre Geliştirilmesi”, 21. SIU , Kıbrıs 2013,.
Barker(N = 7) 16.9 14.7 16.9 [7] Özdil, Ö, İspir, M., Onat, E.,Yıldırım, A.,
Barker(N = 13) 22.28 19.13 22.27 “Implementation of FPGA-Based FFT Convolution”, IET
Radar ,2013.
Tablo 3’de elde edilen sonuçlar [2]’de verilen değerler ile [8] Orduyilmaz, A., Kara, G., Ispir, M., Yildirim, A.,
uyumlu olduğu görülmektedir. “Electronic Attack Techniques Validation Environment”,
Tablo 4’de Virtex-5 Xilinx işlemci kartında tüketilen kaynak SIU, 2013 21st, Haspolat, 24-26 April 2013.
kullanımı sunulmuştur.

View publication stats