FFT, Pencereleme ve Sönüm: "Her Zaman Hann Kullan" Tavsiyesi Nerede Yanlış

1. Kapalı çevrim: doğrulamayı nasıl kuruyoruz

Bir hesap aracının doğruluğunu göstermenin alışılmış yolu, üçüncü taraf bir yazılımın çıktısıyla karşılaştırmaktır. Sinyal işlemede buna ihtiyaç yoktur, çünkü girdiyi kendimiz üretebiliriz.

Buna ek olarak dönüşümün kendi matematiksel özellikleri sınanır. En temizi Parseval eşitliğidir — zaman ve frekans alanındaki enerjinin eşit olması gerekir:

2. Çözünürlük 1/T'dir, başka bir şey değil

SHM-01'de örnekleme frekansının üst sınırı belirlediğini gördük. Frekans çözünürlüğünü belirleyen ise kayıt süresidir:

Bunun pratik sonucu keskindir: birbirine 1/T'den yakın iki mod ayrılamaz. Çakık modlar, burulmalı-ötelenmeli eşleşmeler ve simetrik yapılardaki yakın frekans çiftleri bu sınıra takılır. Çözüm tek yerdedir — daha uzun kayıt.

3. Sıfır dolgu çözünürlük değildir

FFT'nin verimli çalışması için uzunluğun 2'nin kuvveti olması istenir; eksik kalan kısım sıfırla doldurulur. Bu işlem spektrumdaki nokta sayısını artırdığı için sık sık "çözünürlük arttı" diye yorumlanır. Yanlıştır.

Sıfır dolgunun gerçek faydası, tek bir tepenin konumunu ve genişliğini daha iyi örneklemektir. Bu, bölüm 9'daki yarı-güç ölçümünde işe yarar.

4. Sızıntı ve pencerelemenin gerekçesi

FFT, kaydın sonsuz tekrarlandığını varsayar. Kaydın başı ve sonu birbirine uymuyorsa bu varsayım bir süreksizlik üretir ve süreksizlik spektruma sızıntı olarak yayılır: bileşenin enerjisi kendi bininde kalmaz, komşu binlere dağılır.

Sızıntının en görünür sonucu genliğin düşük okunmasıdır. Aşağıdaki tabloda, bin merkezine oturmayan saf bir sinüzoidin genliği ölçülmüştür — gerçek genlik 1,0000'dır.

Buraya kadar tablo, ders kitaplarındaki "sızıntıyı azaltmak için pencereleme yapın" tavsiyesini doğrular. Sıradaki bölüm bu tavsiyenin nerede çöktüğünü gösteriyor.

5. Bulgu: geçici sinyalde Hann sönümü yarıya indirir

Şimdi aynı zinciri serbest sönümlü bir kayda uygulayalım. Sinyali kendimiz ürettiğimiz için gerçek sönüm oranını biliyoruz: ζ = %2,00. Yarı-güç yöntemiyle geri okumaya çalışıyoruz.

Bu, ihmal edilebilir bir yanlılık değil. Sönümü %2 yerine %1 olarak raporlamak, yapının enerji yutma kapasitesini yarı yarıya küçük göstermek demektir — ve sönüm, deprem yanıtı hesabında doğrudan taban kesme kuvvetine giren bir büyüklüktür.

6. Neden? Enerjinin durduğu yer

Sebep, pencerenin biçiminde ve sinyalin biçiminde birlikte yatıyor.

Serbest sönümlü bir yanıtın bütün enerjisi kaydın başındadır. Sinyal t = 0'da en büyüktür ve üstel olarak söner. Hann penceresi ise uçlarda tam olarak sıfıra iner ve ortada bire çıkar. Yani bu iki eğri çarpıldığında pencere, sinyalin en bilgi taşıyan bölümünü söndürür.

Sonuçta elde edilen dalga, artık üstel olarak sönen bir sinyal değildir: önce yükselen sonra alçalan, yapay biçimde simetrikleştirilmiş bir şeydir. Spektral tepesi de gerçek Lorentz biçiminden farklıdır ve bant genişliği daralır — dolayısıyla ζ küçük çıkar.

7. Doğru kural: sinyalin türüne bakın

Ortam titreşimi kayıtları durağan olduğu için, operasyonel modal analizde Hann doğru tercihtir. Bu yüzden "her zaman Hann" tavsiyesi çoğu SHM bağlamında zararsız görünür — ta ki serbest sönüm kaydından sönüm çıkarmaya çalışılana kadar.

8. Tepe seçme ve parabolik ara değer

Modal frekanslar spektrumdaki yerel maksimumlardır. Ancak gerçek frekansın tam bir bin merkezine düşmesi için hiçbir sebep yoktur; genellikle iki bin arasına düşer. Tepe komşularıyla birlikte bir parabole oturtularak konumu bin çözünürlüğünün altına indirilir.

9. Yarı-güç bant genişliğiyle sönüm

Sönümlü bir tepenin genişliği, sönüm oranıyla doğrudan ilişkilidir. Genliğin tepe değerinin 1/√2'sine düştüğü iki frekans arasındaki bant:

Yöntemin iki koşulu vardır ve ikisi de sonuçla birlikte taşınmalıdır:

• Hafif sönüm. ζ büyüdükçe yaklaşım bozulur; %10 üzerinde sonucu bir üst sınır olarak okumak gerekir.
• Yeterli çözünürlük. Bant genişliği yalnızca birkaç bin kadarsa ölçüm anlamsızdır. Motorumuz bant içindeki bin sayısını raporlar ve dörtten azsa uyarır.

İkinci koşul beklenmedik bir sonuç doğurur: sönümü ölçmek için gereken kayıt süresi, sönümün kendisiyle ters ilişkilidir. Küçük ζ dar bir tepe demektir, dar tepe ince çözünürlük gerektirir, ince çözünürlük uzun kayıt gerektirir. Az sönümlü yapıda sönümü ölçmek, çok sönümlü yapıdakinden daha uzun kayıt ister.

10. Uçtan uca doğrulama

Zincirin tamamı, iki modlu bir sentetik sinyalle sınanır: bilinen frekanslar 1,25 ve 3,75 Hz.

11. Yaygın hatalar

1. Pencereyi sinyalin türüne bakmadan seçmek. Bu yazının bütün konusu. Geçici kayıtta Hann, sönümü yarıya indirir.
2. Sıfır dolguyu çözünürlük artışı sanmak. Nokta sayısı artar, ayırt etme yeteneği artmaz.
3. Örnekleme frekansını artırarak çözünürlük beklemek. f_s üst sınırı, T çözünürlüğü belirler.
4. Yarı-güç sonucunu bin sayısını kontrol etmeden raporlamak. Bant genişliği 2–3 bin kaplıyorsa sonuç gürültüdür.
5. Parabolik ara değeri mod ayırmakla karıştırmak. Görünmeyen iki tepe, ara değerle görünür hâle gelmez.
6. Genlik düzeltmesini atlamak. Tarayıcı pencere genliği düşürür; tutarlı kazanca bölünmeden okunan genlik sistematik olarak küçüktür.

12. Kendini sınama

1. 60 saniyelik bir kaydın frekans çözünürlüğü nedir? Bu kayıtla 0,01 Hz aralıklı iki mod ayrılabilir mi?
2. Çekiç darbesiyle uyarılmış bir kirişin serbest sönüm kaydından sönüm çıkaracaksınız. Hangi pencereyi seçersiniz ve neden?
3. Aynı kirişi rüzgârda sürekli titreşirken kaydettiniz. Pencere seçiminiz değişir mi?
4. ζ = %0,5 olan bir modun bant genişliği f₀ = 2 Hz'de kaç Hz'dir? Bunu 8 bin ile örneklemek için en az ne kadar kayıt gerekir?
5. Parseval eşitliğinin sağlanması, spektrumun doğru yorumlandığını da gösterir mi? Cevabınızı gerekçelendirin.