Ana Sayfa / Yazılar / Yapı Sağlığı İzleme SHM-01
Örnekleme, Nyquist ve Takma: Kayıttan Sonra Düzeltilemeyen Hata
Ölçüm zincirinin ilk halkası ve tek geri dönüşü olmayan kararı. Yanlış seçilmiş bir örnekleme frekansı, verinin içine sonradan ayrılamayacak bir yanlış gömer — ve bunun hiçbir yazılımsal çaresi yoktur.
1. Ölçüm zinciri
Yapı sağlığı izlemede bir moda ulaşana kadar veri şu halkalardan geçer:
Bu ayrım yazının bütün konusudur. Pencereyi yanlış seçtiyseniz spektrumu yeniden hesaplarsınız. Örnekleme frekansını yanlış seçtiyseniz ölçümü tekrarlamaktan başka çareniz yoktur.
2. Örnekleme ne yapar?
Sürekli bir ivme sinyalini eşit aralıklarla okuyup sayıya çeviririz. Saniyedeki okuma sayısı örnekleme frekansı fs'dir. Okumalar arasında olan biteni bilmeyiz — ve bilmediğimizi de bilmeyiz. Sayısal veri, aradaki boşluk hakkında hiçbir iz taşımaz.
Sezgisel olarak: bir titreşimin çevrimini tanıyabilmek için çevrim başına en az iki nokta gerekir. Bir tepe, bir çukur. Bundan azıyla, hızlı bir salınımı yavaş bir salınımdan ayırt etmenin matematiksel bir yolu yoktur.
3. Nyquist sınırı ve katlama
Bu sezgi Nyquist–Shannon örnekleme teoremiyle kesinleşir. fs ile örneklenen bir sinyalde doğru temsil edilebilen en yüksek frekans Nyquist sınırıdır:
Görünen frekans katlama bağıntısıyla bulunur:
4. Uygulama: katlanan bileşenler
Aşağıdaki tablo fs = 50 Hz için hesaplanmıştır; sayılar sayfa açılırken motordan üretilir. Nyquist sınırı —'dir.
| Gerçek frekans | Spektrumda görünen | Durum |
|---|---|---|
| hesaplanıyor… | ||
Dikkat edilmesi gereken satır 45 Hz'lik olandır. Yüksek frekanslı bir gürültü kaynağı — bir makine, bir pompa, elektrik şebekesi harmoniği — spektrumda düşük frekanslı bir yapısal mod gibi görünür. Spektruma bakan mühendisin elinde onu ayırt edecek hiçbir bilgi yoktur.
5. Neden sonradan düzeltilemez
Bu yüzden çözüm sayısal değil analogdur: örneklemeden önce, sinyal hâlâ süreklidir, Nyquist üstü içeriği bir alçak geçiren süzgeçle bastırırsınız. Buna takma önleyici (anti-aliasing) süzgeç denir ve ivmeölçer ile ADC arasında donanımda bulunur.
Pratikte süzgeç ideal değildir; kesme frekansından sonra bastırma tedricidir. Bu yüzden örnekleme frekansı genellikle ilgilenilen en yüksek frekansın iki katından belirgin biçimde yüksek seçilir — süzgecin geçiş bandına yer bırakmak için.
6. Çakışma: en kötü durum
Katlamanın en sinsi hâli, birden fazla bileşenin aynı görünen frekansa düşmesidir.
Bu durumda spektrumda tek bir tepe görürsünüz ve o tepe, birbiriyle hiçbir ilgisi olmayan üç ayrı fiziksel olayın toplamıdır. Genliği yorumlamak anlamsızdır; hatta tepenin varlığını bir moda atfetmek de yanlıştır.
7. Kayıt süresi ayrı bir karardır
Örnekleme frekansı ne kadar yükseğe bakabileceğinizi belirler. Kayıt süresi ise ne kadar ince ayırt edebileceğinizi belirler. Bunlar bağımsız iki karardır ve sık karıştırılır.
Sayısal örnek: 0,02 Hz aralıkla duran iki modu ayırmak için en az 50 saniyelik kayıt gerekir; örnekleme frekansını 100'den 1000 Hz'e çıkarmak bu konuda hiçbir şey değiştirmez. Bu ayrım SHM-02'de ayrıntılı işlenir.
8. Ölçüm düzeni kurarken
- İlgilendiğiniz en yüksek modu belirleyin. Yapı sağlığı izlemede genellikle ilk üç–beş mod yeterlidir; bina periyotları için bu birkaç Hz'lik bir bölgedir.
- Örnekleme frekansını bunun en az iki katı, pratikte belirgin biçimde üstünde seçin. Süzgecin geçiş bandına yer bırakın.
- Takma önleyici süzgecin var olduğunu doğrulayın. Bu bir yazılım ayarı değildir; cihazın teknik belgesinde aranır. Yoksa hiçbir örnekleme frekansı sizi kurtarmaz.
- Kayıt süresini çözünürlük ihtiyacından belirleyin, örnekleme frekansından değil.
- Ortamdaki yüksek frekanslı kaynakları not edin. Makine, jeneratör, trafik. Katlandıklarında nereye düşeceklerini önceden hesaplayın — bu yazının tablosu tam olarak bunun içindir.
9. Kendini sınama
- fs = 200 Hz ile örnekleme yapıyorsunuz. 130 Hz'lik bir bileşen spektrumda hangi frekansta görünür?
- Yukarıdaki durumda örnekleme frekansını 400 Hz'e çıkarmak sorunu çözer mi? Peki kaydı iki katına uzatmak?
- Bir spektrumda 5 Hz'de beklenmedik bir tepe görüyorsunuz. Bunun katlanmış bir gürültü olup olmadığını nasıl sınarsınız? (İpucu: örnekleme frekansını değiştirip tekrar ölçün — gerçek bir mod yerinde kalır, katlanan bileşen yer değiştirir.)
- Takma önleyici süzgeç neden sayısal değil analog olmak zorundadır?
- 0,01 Hz aralıklı iki modu ayırmak için en az ne kadar kayıt gerekir?
Kaynaklar
- Nyquist–Shannon örnekleme teoremi. Türetme evrenseldir ve hiçbir ülke yönetmeliğine ait değildir.
- Katlama bağıntısı fgörünen = |f − fs·round(f/fs)|, spektrumun periyodikliğinden doğrudan çıkar.
Bu yazıdaki bütün sayılar engine/dsp.js tarafından üretilir. Doğrulama, üçüncü taraf bir çıktıyla karşılaştırma yoluyla değil, kapalı çevrimle yapılır: bilinen bir sinyali kendimiz üretir, zincirden geçirir ve motorun bilinen değerleri geri bulup bulmadığını ölçeriz.
Takma, kayıttan sonra hiçbir sayısal işlemle düzeltilemez. Önlemi donanımdadır.