Akademik araştırmalar, anlaşılır dil

Verianla | Akademik Araştırmalardan Türkçe Ekonomi ve Bilim İçerikleri

27 Eylül 2026, Pazar
VERİANLABağımsız bilim yayıncılığı
Menüyü aç veya kapat
...
Home / Uygulamalı Bilimler / Mühendislik / Yüksek İrtifa İnsansız Hava Araçlarında RINS/CNS Bütünleşik Navigasyon Sistemi İçin Tam Parametreli Bir Kalibrasyon Yöntemi
Mühendislik

Yüksek İrtifa İnsansız Hava Araçlarında RINS/CNS Bütünleşik Navigasyon Sistemi İçin Tam Parametreli Bir Kalibrasyon Yöntemi

Bu araştırma, yüksek irtifada uzun süre görev yapan insansız hava araçları için dönel modülasyonlu ataletsel navigasyon sistemi ile göksel navigasyon sistemini aynı platformda birleştiren çift eksenli bir RINS/CNS mimarisinin bütün temel hata parametrelerini tek bir kalibrasyon sürecinde belirlemeyi amaçlıyor.

10/08/2026  Veri Anla 26 görüntüleme
Yüksek İrtifa İnsansız Hava Araçlarında RINS/CNS Bütünleşik Navigasyon Sistemi İçin Tam Parametreli Bir Kalibrasyon Yöntemi

Bu araştırma, yüksek irtifada uzun süre görev yapan insansız hava araçları için dönel modülasyonlu ataletsel navigasyon sistemi ile göksel navigasyon sistemini aynı platformda birleştiren çift eksenli bir RINS/CNS mimarisinin bütün temel hata parametrelerini tek bir kalibrasyon sürecinde belirlemeyi amaçlıyor. Sistemde ataletsel ölçüm birimi (IMU) iki eksenli bir indeksleme/gimbal mekanizmasına yerleştiriliyor; keşif amacıyla kullanılan kamera ise ayrıca yıldız sensörü olarak değerlendiriliyor. Böylece ayrı bir yıldız izleyici ve ek kalibrasyon donanımı gereksiniminin azaltılması hedefleniyor.

Çalışmanın temel yeniliği yalnız jiroskop ve ivmeölçer bias değerlerini kalibre etmek değildir. Araştırmacılar aynı hata durum vektöründe IMU bias, ölçek faktörü ve çapraz eksen/kurulum hatalarını; gimbal eksenlerinin ortogonallik bozukluklarını; iç ve dış eksen encoder açı hatalarını; kameranın dış yönelim parametrelerini (EOP); odak uzaklığı, ana nokta ve skew gibi iç parametrelerini (IOP) ve Brown–Conrady radyal/teğetsel lens distorsiyon parametrelerini birlikte tahmin etmektedir.

Bu amaçla 47 boyutlu bir hata-durum modeli ve yıldız görüntüsündeki piksel artıklarına dayanan Extended Kalman Filter yapısı geliştirilmiştir. IMU verileri 100 Hz'de yayılırken yıldız-piksel güncellemeleri 2 Hz'de kullanılmaktadır. Kalibrasyon yolu, farklı hata gruplarını aynı anda gelişi güzel uyarmak yerine belirli segmentlerde belirli parametrelerin gözlemlenebilirliğini artıracak biçimde tasarlanmıştır.

Üç temel ayrıştırma mekanizması kullanılmaktadır: iç eksenin simetrik taramalarında tek/çift (odd/even) simetriyle iç encoder hatasının ayrılması; dış eksenin farklı referans konumlarıyla ölçüm Jacobian'ının bazının değiştirilmesi; dış encoder hatasının ise üçgen dalga biçimli dış eksen hareketine frekans etiketi uygulanarak diğer quasi-DC hata bileşenlerinden ayrılması.

Saf simülasyonda IMU, gimbal ve kamera parametrelerinin büyük bölümü kendileri için bilgi taşıyan ilgili hareket segmentinde yakınsamış ve çalışma sonunda hata seviyeleri genel olarak yaklaşık %5'in altında kalmıştır. HIL/semi-physical doğrulamada gerçek IMU ve gimbal encoder verileri kullanıldığında sonuçlar daha gürültülü hale gelmiş; araştırmacılar son hata seviyelerini genel olarak yaklaşık %6–16 aralığında raporlamıştır.

Ancak HIL doğrulaması gerçek bir yüksek irtifa İHA uçuşu değildir. IMU ve gimbal verilerinin fiziksel donanımdan gelmesine rağmen yıldız ölçümleri SAO yıldız kataloğundan, gerçek yönelim ve görüş alanı kısıtları kullanılarak sentetik biçimde oluşturulmuştur. Dolayısıyla çalışma, yöntemin laboratuvar/yarı-fiziksel uygulanabilirliğini desteklemekte; gerçek gökyüzü görüntüsü, uçuş titreşimi, termal gradyanlar ve uzun süreli navigasyon hatası altında tam uçuş doğrulamasını henüz göstermemektedir.

Türkiye açısından anlamı: Türkiye'deki yüksek irtifa İHA, uzun havada kalış platformu, GNSS'siz navigasyon, elektro-optik keşif sistemi ve ataletsel navigasyon Ar-Ge çalışmalarında benzer bir birleşik kalibrasyon yaklaşımı yöntemsel açıdan değerlendirilebilir. Özellikle görev kamerasının aynı zamanda göksel referans sensörü olarak kullanılması, platformdaki sensör ve kütle tekrarını azaltma potansiyeli taşır. Ancak bu araştırmadaki kalibrasyon rotası, IMU kalitesi, yıldız görünürlüğü ve hata büyüklükleri belirli bir simülasyon/HIL konfigürasyonuna aittir; gerçek Türkiye hava sahası görevleri için atmosferik koşullar, kamera optiği, gimbal yapısı, uçuş titreşimi, zaman senkronizasyonu ve gerçek yıldız görünürlüğüyle yeniden doğrulama gerekir.

Neden RINS ile CNS birlikte kullanılıyor?

Strapdown ataletsel navigasyon sistemi dışarıdan sinyal almadan konum, hız ve yönelim hesabı yapabildiği için GNSS'nin karıştırıldığı veya erişilemediği görevlerde önemlidir. Ancak jiroskop ve ivmeölçerlerdeki çok küçük bias, ölçek veya eksen hizalama hataları zamanla entegre edilerek navigasyon hatasının büyümesine neden olur.

Rotational-Modulation Inertial Navigation System (RINS), IMU'yu kontrollü biçimde döndürerek sabit veya yavaş değişen sensör hatalarını sinüzoidal biçimlere taşımayı ve zaman içinde ortalama etkilerini azaltmayı amaçlar.

Celestial Navigation System (CNS) ise yıldız doğrultularını mutlak yönelim referansı olarak kullanabilir. Böylece ataletsel sistemin zamanla büyüyen yönelim hataları göksel ölçümlerle sınırlandırılabilir.

Çalışmanın yaklaşımı bu iki mekanizmayı birleştirerek yalnız navigasyon yapmak değil, aynı hareketleri sensör sisteminin kendi kalibrasyonu için de kullanmaktır.

Sistemin fiziksel mimarisi nasıl?

Kaynak çalışmanın sayfa 5'teki Şekil 1'i sistemin temel yapısını gösterir. Ortak bir mekanik platform üzerinde:

  • iç gimbal ekseni,
  • dış gimbal ekseni,
  • IMU,
  • optik lens/kamera

birlikte bulunmaktadır.

İç gimbal nominal olarak gövde \(+X_b\) ekseni etrafında \(\alpha\) açısıyla; dış gimbal ise \(+Z_b\) ekseni etrafında \(\beta\) açısıyla hareket etmektedir.

Aynı sayfadaki Şekil 2, Dünya/nav koordinat sistemi, gövde sistemi ve sensör/gimbal sistemleri arasındaki dönüşümleri göstermektedir. Bu koordinat zinciri kritik önemdedir; çünkü encoder, gimbal ortogonallik, IMU kurulum ve kamera boresight hatalarının tamamı yıldız doğrultusunun tahmin edilen piksel konumuna ulaşmasına kadar aynı dönüşüm zincirinde birbirine karışabilmektedir.

Gimbal hataları nasıl modelleniyor?

Gerçek mekanik sistemde iki gimbal ekseninin kusursuz biçimde birbirine dik olduğu veya IMU eksenlerinin mekanik dönme eksenleriyle tam çakıştığı varsayılamaz.

Çalışma gimbal geometrisini dört küçük ortogonallik parametresiyle ifade etmektedir:

\[ \theta_x,\;\theta_y,\;\theta_z,\;\phi \]

Ayrıca iç ve dış encoder ölçümlerindeki açı hataları:

\[ \delta\alpha,\qquad\delta\beta \]

olarak tanımlanmıştır.

Gerçek eksen açıları:

\[ \tilde{\alpha}=\alpha+\delta\alpha \] \[ \tilde{\beta}=\beta+\delta\beta \]

biçiminde yazılmaktadır.

Birinci dereceden lineerleştirme sonucunda araştırmacılar üç önemli gözlem elde ediyor:

  • \(\delta\alpha\) en güçlü biçimde iç eksen taramalarıyla uyarılır.
  • \(\delta\beta\) dış eksen hareketi olmadan güçlü biçimde gözlemlenemez.
  • Gimbal non-orthogonality parametreleri ölçümde quasi-DC bir bozucu gibi davranır.

IMU hata modeli neleri içeriyor?

Jiroskop ölçüm modeli kaynakta:

\[ \tilde{\omega}^{s} = (I+\Delta K_g)\omega^s + \Delta A_g\omega^s + \varepsilon^s + n_g \]

olarak verilmektedir.

İvmeölçer için:

\[ \tilde{f}^{s} = (I+\Delta K_a)f^s + \Delta A_a f^s + \nabla^s + n_a \]

kullanılır.

Burada:

  • \(\varepsilon^s\): jiroskop sabit bias'ları,
  • \(\nabla^s\): ivmeölçer sabit bias'ları,
  • \(\Delta K_g,\Delta K_a\): ölçek faktörü hataları,
  • \(\Delta A_g,\Delta A_a\): çapraz eksen/kurulum hataları,
  • \(n_g,n_a\): ölçüm gürültüleridir.

Kamera neden yalnız dış yönelimle kalibre edilmiyor?

Çalışmada kamera basit bir yönelim sensörü olarak değerlendirilmemektedir. Kamera iç geometrisi de kalibrasyon durum vektörüne dahil edilmiştir.

İç parametre grubu:

\[ \eta_{IOP} = [f_x,f_y,c_x,c_y,s,k_1,k_2,p_1,p_2]^T \]

olarak tanımlanmıştır.

Burada \(f_x,f_y\) piksel cinsinden odak uzaklıkları; \(c_x,c_y\) ana nokta; \(s\) pixel skew; \(k_1,k_2\) radyal ve \(p_1,p_2\) teğetsel distorsiyon parametreleridir.

Brown–Conrady distorsiyon modeli

Normalize edilmiş ideal görüntü koordinatları \(x,y\) için:

\[ r^2=x^2+y^2 \]

ve distorsiyonlu koordinatlar:

\[ x_d = x+x(k_1r^2+k_2r^4) +2p_1xy +p_2(r^2+2x^2) \] \[ y_d = y+y(k_1r^2+k_2r^4) +p_1(r^2+2y^2) +2p_2xy \]

olarak modellenmektedir.

Bu denklem S5 segmentinin neden kamerayı yalnız görüntü merkezinde değil kenar ve köşe bölgelerinde de yıldızlarla uyarması gerektiğini açıklar. \(r\) küçük kaldığında özellikle radyal distorsiyon parametrelerine ilişkin bilgi zayıf kalmaktadır.

Bir yıldız gözlemi bütün hata gruplarına nasıl bağlanıyor?

Araştırmanın temel birleşik hata denklemlerinden biri yıldız birim vektöründeki küçük hatayı dört bileşene ayırmaktadır:

\[ \delta s^n \approx -[s^n\times]\phi^n - [s^n\times]C^n_b\delta\theta_{EOP} + C^n_b(\delta C^b_s)R_{BC}s^c + C^n_bR_{BC}J_i\delta\eta_{IOP} \]

Bu ifade aynı yıldız piksel artığının:

  • navigasyon/gövde yönelim hatasını,
  • kamera EOP hatasını,
  • gimbal encoder ve non-orthogonality hatalarını,
  • kamera IOP ve lens distorsiyon hatalarını

aynı anda taşıyabildiğini göstermektedir.

Dolayısıyla temel problem yalnız çok sayıda parametre tahmin etmek değildir; birbirine çok benzer ölçüm imzaları oluşturan parametrelerin birbirinden ayrılabilmesidir.

Neden 47 boyutlu hata durumu kullanılıyor?

Çalışmadaki birleşik hata-durum vektörü toplam 47 parametreden oluşmaktadır.

Parametre grubuDurum sayısı
Navigasyon hız + tutum hataları6
IMU bias'ları6
IMU ölçek faktörleri6
IMU çapraz eksen/kurulum terimleri12
Gimbal non-orthogonality + encoder hataları6
Kamera IOP + distorsiyon8 + skew
Kamera EOP3

Kaynak denkleminde IOP bloğu \(f_x,f_y,c_x,c_y,k_1,k_2,p_1,p_2,s\) olmak üzere dokuz parametre içerir ve bütün bloklarla birlikte toplam durum boyutu 47'dir.

EKF hangi hızlarda çalışıyor?

Önerilen hata-durum EKF yapısında:

  • IMU propagation: 100 Hz,
  • yıldız-piksel güncellemesi: 2 Hz,
  • HIL gimbal encoder akışı: 1 kHz

olarak verilmiştir.

Her görülebilir yıldız için tahmin edilen piksel ile ölçülen piksel arasındaki artık filtreye ölçüm olarak girer.

Görüntü artıklarında aykırı değerlerin filtreyi bozmasını azaltmak için \(\chi^2\) gating, residual whitening ve geometrik görünürlük kontrolleri uygulanmaktadır.

Statik segmentlerde ZUPT neden var?

S0 ve S6 statik bekleme bölümlerinde opsiyonel zero-velocity update uygulanabilmektedir:

\[ z_{zupt} = \hat v^n-0 \]

Bu ek ölçüm özellikle ivmeölçer bias ve ilgili hata durumlarının gözlemlenebilirliğini desteklemektedir.

Üç hata ayrıştırma mekanizması

1. İç eksende tek/çift simetri

Dış eksen \(\beta\) sabitken iç eksende \(+\alpha\) ve \(-\alpha\) simetrik geçişler uygulanmaktadır.

\[ z^{odd} = \frac{1}{2} [z(+\alpha,\beta)-z(-\alpha,\beta)] \] \[ z^{even} = \frac{1}{2} [z(+\alpha,\beta)+z(-\alpha,\beta)] \]

İç encoder hatasının etkisi \(\alpha\)'ya göre tek simetri gösterirken EOP, dış encoder ve non-orthogonality terimleri birinci mertebede daha çok çift/quasi-DC bileşende kalmaktadır.

Böylece odd kanalın:

\[ z^{odd}\rightarrow\delta\alpha \]

parametresini kuvvetli biçimde ayırması amaçlanır.

2. Basis switching

İç eksen taraması yalnız \(\beta=0^\circ\) konumunda yapılırsa bazı Jacobian kolonları birbirine çok benzer hale gelebilir.

Bu nedenle dış eksen:

\[ \beta=0^\circ,\quad+45^\circ,\quad-45^\circ \]

referanslarına taşınmaktadır.

Yeni dış eksen konumu ölçüm hassasiyet matrisinin geometrik bazını değiştirir. Böylece özellikle gimbal parametreleri ile sabit kamera EOP hataları arasındaki korelasyonun azaltılması hedeflenir.

3. Dış encoder için frekans etiketleme

Dış eksende rate-limited üçgen hareket uygulanarak \(\delta\beta\) parametresinin dinamik imzası quasi-DC terimlerden ayrılmaktadır.

Kaynakta kullanılan demodülasyon mantığı:

\[ D_{\beta}[z] = \frac{2}{T} \int_{t_0}^{t_0+T} z(t)\, \mathrm{sgn}(\dot{\beta}(t))dt \]

şeklindedir.

Dış encoder bileşeni bu etiketlemeyle korunurken sabit EOP ve non-orthogonality bileşenlerinin bastırılması amaçlanmaktadır.

Gözlemlenebilirlik nasıl ölçülüyor?

Her hareket segmenti için Fisher Information Matrix oluşturulmakta ve Singular Value Decomposition ile bilgi içeriği değerlendirilmektedir.

Segment \(S\) için:

\[ W_S = \sum_{k\in S} \Phi^T_{k,0} H^T_k R^{-1}_k H_k \Phi_{k,0} \]

tanımlanmıştır.

Belirli bir parametre grubu için normalize gözlemlenebilirlik indeksi:

\[ NOI_S(\Theta) = \frac{ \sigma_{\min} (E^T_{\Theta}W_SE_{\Theta}) }{ \sigma_{\star}(S) } \]

olarak hesaplanmaktadır.

NOI değeri büyüdükçe ilgili hareket segmentinin seçilen parametre grubuna daha fazla bağımsız bilgi sağladığı kabul edilmektedir.

Kalibrasyon yolunun görev dağılımı

Verianla Live: RINS/CNS tam parametre kalibrasyonunun segment mantığı

Aşağıdaki süreç kaynak Tablo 2 ve Tablo 3'ün bilimsel mantığını özetler. Tablo 2'de S4 süresine ilişkin kaynak içi zaman uyuşmazlığı bulunduğu için Live süreç tablosuna kesin segment süreleri eklenmemiştir.

SegmentHareket / temel amaçÖzellikle bilgilendirilen parametreler
S0Statik başlangıç ve opsiyonel ZUPTBaşlangıç kararlılığı, ivmeölçer ilişkili durumlar
S1β = 0° bazında iç eksen simetrik taramasıİç encoder, jiroskop bias/ölçek/kurulum, ivmeölçer parametreleri
S2β = +45° bazında iç eksen taramasıBasis switching, EOP, IMU eksen ayrıştırması
S3β = −45° bazında iç eksen taramasıBasis switching, EOP ve çapraz eksen parametreleri
S4Dış eksen üçgen hareketi ve frekans etiketlemeDış encoder δβ, dış eksen ilişkili ölçek ve non-orthogonality
S5Görüntü merkezinden kenar/köşe bölgelerine yıldız taramasıKamera IOP, skew ve Brown–Conrady distorsiyon
S6Statik son beklemeYakınsamanın ve son parametre platosunun kontrolü
 

Verianla Live: Görselleştirme görünür bilimsel süreç tablosundan oluşturulur; kaynakta bulunmayan bir hareket veya parametre eşlemesi eklenmemiştir.

Hangi segment hangi parametreyi gözlemlenebilir yapıyor?

Parametre grubuTemel uyarımÖne çıkan segmentler
Jiroskop biasYıldız tutum güncellemeleri + iç eksen taramasıS1–S3
Jiroskop X ölçeği+Xb ekseni etrafında dönmeS1
Jiroskop Y ölçeğiβ = ±45° ile eksen projeksiyonuS2–S3
Jiroskop Z ölçeğiDış eksen hareketiS4
İvmeölçer bias / scale / crossYerçekimi yönünün sensör eksenlerine göre değiştirilmesiS1–S3
İç encoder δαSimetrik iç eksen üçgen taramalarıS1–S3
Dış encoder δβDış eksen üçgen hareketi + demodülasyonS4
Gimbal non-orthogonalityÇift eksenli hareket + quasi-DC birikimS1–S4
Kamera EOPBasis switchingS2–S3
Kamera IOPGeniş görüntü alanı taramasıS5
Lens distorsiyonuBüyük görüntü yarıçapı ve azimut çeşitliliğiS5

Yöntem ve Bulgular

Saf simülasyon koşulları

ParametreKaynak değeri
Jiroskop sabit bias0,2°/h
Jiroskop ARW0,02°/√h
Jiroskop kurulum hatası40 arcsec
Jiroskop ölçek faktörü50 ppm
İvmeölçer sabit bias100 µg
İvmeölçer velocity random walk10 µg/√h
İvmeölçer kurulum hatası40 arcsec
İvmeölçer ölçek faktörü30 ppm
Gimbal non-orthogonality2 arcmin
Encoder hataları40 arcsec
Kamera EOP2 arcmin
Piksel gürültüsü σ0,3 px
IMU güncelleme oranı100 Hz
Piksel güncelleme oranı2 Hz

Simülasyon yakınsama davranışı

Kaynak Şekil 5–15'te parametre tahminleri zaman ekseni boyunca gösterilmektedir. Grafiklerin ortak özelliği, parametrelerin bütün 12 dakika boyunca aynı hızda yakınsamamasıdır.

Bir durum kendisini uyaran segmente ulaşmadan önce tahmin eğrisi büyük ölçüde düz kalmakta; ilgili hareket başladığında hızla gerçek değere yaklaşmaktadır.

Bu davranış, çalışmanın “segment-wise observability” iddiasının en önemli sayısal göstergesidir.

Jiroskop bias sonuçları

Üç jiroskop bias'ı ağırlıklı olarak S1–S3 içinde yakınsamıştır.

EksenGerçek (°/h)Son tahmin (°/h)Kaynakta verilen bağıl hata
X0,2000,209%4,50
Y0,2000,208%3,89
Z0,2000,192%3,87

Encoder sonuçları simülasyonda yolu doğruluyor

Kaynak Şekil 12'de iç encoder hatası \(\delta\alpha\) ilk iç eksen taramalarında hızla 40 arcsec civarındaki gerçek değere ulaşırken dış encoder hatası \(\delta\beta\) S4'e kadar düşük bilgiyle kalmaktadır.

S4 dış eksen üçgen hareketi başladığında \(\delta\beta\) tahmini keskin biçimde gerçek değere yakınsamaktadır.

Bu grafik, dış encoder hatasının yalnız iç eksen taramasıyla güvenilir biçimde tahmin edilemediğini doğrudan göstermektedir.

Kamera parametreleri neden S5'e kadar bekliyor?

Kaynak Şekil 13 ve 14'te kamera intrinsics ve Brown–Conrady distorsiyon tahminlerinin S5'e kadar sınırlı bilgi topladığı görülmektedir.

S5'te yıldız noktaları görüntü merkezinden görüntü kenarına/köşeye doğru taşındığında:

  • \(f_x,f_y\),
  • \(c_x,c_y\),
  • pixel skew,
  • \(k_1,k_2,p_1,p_2\)

parametreleri hızlı biçimde yakınsamaktadır.

Bu sonuç, optik kalibrasyonda yalnız çok sayıda yıldız görmenin değil, yıldızların görüntü düzlemindeki geometrik dağılımının da önemli olduğunu göstermektedir.

Normalize gözlemlenebilirlik grafiği ne söylüyor?

Kaynak sayfa 19'daki Şekil 16, S1–S6 boyunca farklı parametre bloklarının NOI değerlerini yan yana göstermektedir.

Grafikte:

  • iç encoder bilgisi S1–S3'te güçlü,
  • dış encoder bilgisi S4'te belirgin,
  • EOP bilgisi S2–S3'te güçlenmiş,
  • IOP/distorsiyon bilgisi S5'te belirgin biçimde yükselmiştir.

Bu nedenle kalibrasyon yolu yalnız sezgisel bir hareket dizisi olarak değil, Fisher bilgi matrisi ve SVD sonuçlarıyla ilişkilendirilmiştir.

Saf simülasyonun genel sonucu

Araştırmacılar 12 dakikalık çalışma sonunda IMU, gimbal ve kamera parametrelerinin genel olarak yaklaşık %5 sınırı içinde yakınsadığını bildirmektedir.

Ancak kamera intrinsiklerine ilişkin Tablo 5'in “Rel. %” sütunundaki bazı değerler truth–estimate oranıyla doğrudan uyuşmadığından Verianla makalesinde kamera için her satırdan yeniden türetilmiş bağıl yüzde sunulmamıştır.

HIL doğrulamasında gerçek olan ve sentetik olan nedir?

Kaynak sayfa 21'deki Şekil 17 fiziksel HIL test düzeneğini göstermektedir. Görselde üç eksenli turntable düzeni, kontrol bilgisayarı ve test ekipmanı görülmektedir.

HIL yapısında:

  • IMU: rate-table üzerinde gerçek sensör verileri, 100 Hz,
  • Gimbal: gerçek çift eksenli encoder akışı, 1 kHz,
  • Yıldız sensörü: SAO kataloğundan sentetik yıldız doğrultuları,
  • Piksel gürültüsü: yaklaşık 0,3 px,
  • Piksel güncellemesi: 2 Hz

kullanılmıştır.

Bu nedenle çalışma “hardware-in-the-loop” olarak adlandırılsa da optik sensör zinciri bütünüyle gerçek yıldız görüntüsü/centroid ölçümünden oluşmamaktadır.

HIL'de neden sonuçlar kötüleşiyor?

Gerçek sistemde yıldız görünürlüğü ve ölçüm kalitesi ideal değildir. Araştırmacılar HIL sürecinde yaklaşık %10–15 görüntü karesinin segment geçişlerindeki görünürlük değişimi veya centroid aykırı değerleri nedeniyle elendiğini bildirmektedir.

Bunun sonucunda yakınsama eğrileri saf simülasyona göre daha gürültülü ve daha yavaş hale gelmektedir.

Saf simülasyon ve HIL jiroskop bias karşılaştırması

Verianla Live: Jiroskop bias kalibrasyonunda saf simülasyon ve HIL farkı

Grafik yalnız kaynak tablolarında aritmetik olarak truth–estimate farkıyla uyumlu olan jiroskop bias satırlarını kullanır. Değerler kaynakta raporlanan bağıl hata büyüklükleridir.

Doğrulama / eksenBağıl hata (%)Kaynak
Saf simülasyon – Gyro X bias4,50Tablo 5
Saf simülasyon – Gyro Y bias3,89Tablo 5
Saf simülasyon – Gyro Z bias3,87Tablo 5
HIL – Gyro X bias12,87Tablo 6
HIL – Gyro Y bias11,15Tablo 6
HIL – Gyro Z bias13,36Tablo 6
 

Verianla Live: Görselleştirme görünür bilimsel tablodan oluşturulmuştur. Kaynakta aritmetik uyuşmazlık bulunan HIL encoder ve kamera satırları bu karşılaştırmaya dahil edilmemiştir.

HIL çalışmasının genel doğrulama sonucu

Araştırmacılar HIL deneylerinde bütün parametre gruplarının yine kendileri için bilgi taşıyan segmentlerde yakınsadığını ve S6 statik bölümünde geç dönem diverjans göstermeden kararlı kaldığını raporlamaktadır.

Genel hata düzeyi saf simülasyondaki yaklaşık ≤%5 seviyesinden pratik HIL koşullarında yaklaşık %6–16 bandına yükselmektedir.

Tablo 6'da bazı satırlarda %16'yı hafif aşan değerler de bulunmaktadır; örneğin accelerometer Y bias için %16,30 ve kamera EOP X için %16,34 yazılmıştır. Dolayısıyla “%16” ifadesi kesin matematiksel üst sınır yerine araştırmacıların yaklaşık performans özetidir.

Hesaplama yükü uçuş bilgisayarı açısından yüksek mi?

Kaynakta hata durum boyutu \(n=47\) ve tipik piksel ölçüm yığını yaklaşık \(m\approx40\) olarak verilmiştir.

Bir ölçüm güncellemesinin baskın karmaşıklığı:

\[ O(n^2m+m^3) \]

ve en kötü durumda kovaryans güncellemesi:

\[ O(n^3) \]

olarak değerlendirilmiştir.

Araştırmacılar 2 Hz piksel güncellemesi başına yaklaşık \(10^5\) kayan nokta işlemi mertebesinin güncel düşük güçlü uçuş işlemcileri için yönetilebilir olduğunu belirtmektedir.

Yöntemin güçlü yönleri

  • IMU, gimbal ve kamera parametreleri tek hata durum modelinde birleştirilmiştir.
  • Kamera yalnız boresight açısından değil iç parametreleri ve Brown–Conrady distorsiyonuyla birlikte kalibre edilmektedir.
  • Parametre ayrıştırması yalnız filtre ayarına bırakılmamış, fiziksel hareket rotası gözlemlenebilirlik için özel olarak tasarlanmıştır.
  • İç ve dış encoder hataları farklı simetri/frekans imzalarıyla ayrıştırılmaktadır.
  • PWCS, Fisher Information Matrix ve SVD ile hangi segmentin hangi parametreye bilgi sağladığı sayısal olarak incelenmiştir.
  • Saf simülasyona ek olarak gerçek IMU ve gimbal verileriyle HIL doğrulaması yapılmıştır.
  • Keşif kamerasının yıldız sensörü olarak yeniden kullanılması ayrı sensör yükünün azaltılmasına yönelik sistem seviyesi bir yaklaşım sunmaktadır.

Çalışmanın desteklediği sonuçlar

  • Çift eksenli hareket yolu altında IMU, gimbal ve kamera hata grupları aynı EKF içinde tahmin edilebilmektedir.
  • İç encoder hatası ağırlıklı olarak S1–S3 iç eksen taramalarında gözlemlenebilir hale gelmektedir.
  • Dış encoder hatası S4 dış eksen hareketi olmadan güçlü biçimde yakınsamamaktadır.
  • Kamera EOP'ları dış eksen basis switching ile S2–S3'te bilgi kazanmaktadır.
  • Kamera IOP ve Brown–Conrady distorsiyonu S5'teki merkez–kenar görüntü taramasından güçlü biçimde yararlanmaktadır.
  • Saf simülasyonda parametre grupları genel olarak yaklaşık %5 hata seviyesinin içinde yakınsamaktadır.
  • Gerçek IMU/gimbal verileri içeren HIL ortamında segment bazlı yakınsama davranışı korunmaktadır.
  • HIL'deki görünürlük kaybı ve ideal olmayan centroid dağılımları saf simülasyona kıyasla hata seviyelerini artırmaktadır.

Çalışmanın henüz kanıtlamadığı noktalar

  • Yöntem gerçek bir HALE İHA uçuşunda doğrulanmamıştır.
  • Gerçek yıldız görüntüleriyle kamera centroid extraction zinciri HIL çalışmasında kullanılmamıştır.
  • Atmosferik görüş, bulut, gün ışığı ve gerçek gökyüzü arka planı modellenmemiştir.
  • Kamera–IMU zaman senkronizasyon hatası açık biçimde kalibre edilmemiştir.
  • Rolling-shutter kamera etkisi modellenmemiştir.
  • Sıcaklık gradyanlarına bağlı IMU/gimbal parametre değişimleri mevcut kısa kalibrasyon penceresinde modellenmemiştir.
  • Uzun dönem sensör driftleri zamanla değişen parametre modeliyle doğrulanmamıştır.
  • Gimbal veya rate-table mekanik arızaları saf simülasyonda dikkate alınmamıştır.
  • Kalibrasyon sonrasında gerçek görev navigasyon hatasının metre veya derece cinsinden ne kadar iyileştiğine ilişkin gerçek uçuş sonucu verilmemiştir.
  • SAO katalog tabanlı sentetik optik doğrulama, gerçek yıldız sensörünün optik ve centroid hatalarının tamamını temsil etmez.

Gelecek çalışma

Araştırmacılar gelecekte sıcaklık ve zamana bağlı sensör hata modellerini kalibrasyona eklemeyi, görüş alanı ve mekanik hız limitleri altında kalibrasyon rotasını doğrudan gözlemlenebilirlik ölçütleriyle optimize etmeyi, hareketli UAV üzerinde kapalı çevrim uçuş testi gerçekleştirmeyi ve gerçek zamanlı Fisher bilgi miktarına göre kalibrasyon segmentlerini adaptif olarak uzatmayı planlamaktadır.

Kaynak ve Yöntem Notu

Tam özgün çalışma adı: A Full-Parameter Calibration Method for an RINS/CNS Integrated Navigation System in High-Altitude Drones

Yazarlar: Huanrui Zhang; Xiaoyue Zhang; Chunhua Cheng; Xinyi Lv; Chunxi Zhang.

Yazar sırası: Kaynak çalışmadaki sıra aynen korunmuştur.

Eş katkı/eş birinci yazar: Kaynakta eş birinci yazarlık veya eş katkı bildirimi bulunmamaktadır.

Sorumlu yazar: Chunhua Cheng.

Kurumlar:

  1. School of Instrument Science and Opto-Electronics Engineering, Beihang University, Beijing, China.
  2. Qingdao Branch, Naval Aeronautic University, Qingdao, China.

Dergi: Vehicles.

Yayınevi: MDPI, Basel, Switzerland.

Bibliyografik kayıt: Vehicles 2026, 8, 11.

DOI: 10.3390/vehicles8010011

Gönderim tarihi: 28 Ekim 2025.

Revizyon tarihi: 1 Aralık 2025.

Kabul tarihi: 27 Aralık 2025.

Yayın tarihi: 5 Ocak 2026.

Kaynak türü ve hakemlik: Hakemli araştırma makalesi. Ana doğrulama saf simülasyon ve yarı-fiziksel Hardware-in-the-Loop testlerinden oluşmaktadır.

Preprint durumu: İncelenen kaynakta ayrı bir preprint sürümüne ilişkin bilgi verilmemektedir; bu Verianla değerlendirmesinde yayımlanmış hakemli sürüm esas alınmıştır.

Resmî yayın bağlantısı:https://www.mdpi.com/2624-8921/8/1/11

DOI bağlantısı:https://doi.org/10.3390/vehicles8010011

Lisans: Creative Commons Attribution 4.0 International (CC BY 4.0). Çalışma açık erişimlidir.

Finansman: Araştırmacılar çalışmanın dış finansman almadığını bildirmiştir.

Veri erişilebilirliği: Kaynakta “No new data were created or analyzed in this study” ifadesi yer almaktadır.

Çıkar çatışması: Yazarlar çıkar çatışması bulunmadığını beyan etmektedir.

Yazar katkıları: Yöntem Huanrui Zhang; doğrulama Xiaoyue Zhang; biçimsel analiz Chunxi Zhang; araştırma Chunhua Cheng ve Xinyi Lv tarafından gerçekleştirilmiştir.

Deneysel doğrulama sınırı: HIL testinde fiziksel IMU/rate-table verileri ve yüksek hassasiyetli çift eksenli gimbal encoder akışları kullanılmıştır. Yıldız ölçümleri ise gerçek optik görüntülerden değil, SAO kataloğu ve gerçek yönelim/FOV kısıtları kullanılarak sentetik olarak üretilmiştir. Bu nedenle çalışma tam fiziksel uçuş veya gerçek yıldız sensörü görüntü testi değildir.

Kaynak içindeki önemli teknik ve editoryal tutarsızlıklar

  • Çalışma kalibrasyon rotasını toplam 12 dakika olarak tanımlamakta ve Şekil 4 yaklaşık 0–12 dakikalık segment zaman çizelgesi göstermektedir. Ancak Tablo 2'de S4 dış eksen üçgen hareketinin süresi 12 s olarak yazılmıştır. Tablo 2'deki bütün süreler bu değerle toplandığında yaklaşık 609 s, yani 10,15 dakika elde edilmektedir. S4 süresi 120 s olduğunda toplam yaklaşık 717 s veya 11,95 dakika olmaktadır. Bu nedenle “12 s” değeri güçlü biçimde yazım hatası izlenimi verse de Verianla bunu sessizce değiştirmemiştir.
  • Bölüm 4.3 kalibrasyon yolunu geçişlerle birlikte “11-segment decoupling path” olarak tanımlar. Sonuç bölümü ise S0–S6 dizisini “six rate-limited segments” olarak adlandırmaktadır. S0'dan S6'ya yedi S-segmenti bulunduğundan terminoloji kendi içinde tam tutarlı değildir.
  • Tablo 5'te kamera iç parametrelerinin “Rel. %” değerleri standart \(|estimate-truth|/|truth|\) hesabıyla uyuşmamaktadır. Ayrıca Bölüm 5.4, \(f_y\) için %4,70 en büyük bağıl hata değerini söylerken Tablo 5 \(f_y\) satırında %1,71 göstermektedir. Bu nedenle Verianla kamera intrinsikleri için kaynaktaki yüzde sütununu yeniden yorumlayarak yeni sonuç üretmemiştir.
  • HIL Tablo 6'da bazı parametrelerde verilen Truth, Estimate, Error ve Rel. % değerleri birbirleriyle aritmetik olarak uyuşmamaktadır. Özellikle bazı encoder, kamera intrinsics, distorsiyon ve EOP satırlarında doğrudan çıkarma işlemi tabloda yazılı “Error” değerini vermemektedir. Bu yüzden HIL'e ilişkin genel %6–16 performans ifadesi kaynak yazarlarının özeti olarak aktarılmış, sorunlu satırlardan bağımsız yeni yüzde türetilmemiştir.
  • Kaynak özeti HIL sonuçlarını yaklaşık %6–16 hata aralığında tanımlamaktadır. Tablo 6'da accelerometer Y bias için %16,30 ve EOP X için %16,34 gibi %16'yı hafif aşan değerler bulunmaktadır; dolayısıyla %16 kesin matematiksel üst sınır olarak yorumlanmamalıdır.

Bilimsel içerik notu: Bu Verianla makalesindeki sistem mimarisi, denklemler, hareket segmentleri, hata durumları, filtre ayarları, simülasyon sonuçları, HIL verileri ve sınırlılıklar yalnız incelenen kaynak çalışmaya dayanmaktadır. Dış kaynaklar bilimsel sonuç eklemek için değil, yalnız bibliyografik yayın kimliği ve hakemli yayın durumunu doğrulamak amacıyla kullanılmıştır.

Verianla Live notu: Süreç bileşeni kaynak Tablo 2–3 ve Bölüm 4'teki gözlemlenebilirlik mantığından oluşturulmuştur. Jiroskop bias karşılaştırmasında kullanılan saf simülasyon ve HIL bağıl hata değerleri kaynak Tablo 5 ve 6'daki aritmetik olarak tutarlı satırlardan aynen alınmıştır. Kaynakta bulunmayan veri veya interpolasyon üretilmemiştir.


Paylaş:

Yorumlar incelendikten sonra yayımlanır.Gönderdiğiniz yorum onay sürecine alınır ve uygun bulunduğunda görünür hâle gelir.

Bir yorum bırakın

E-posta adresiniz yayınlanmayacaktır. Gerekli alanlar * ile işaretlenmiştir

Your experience on this site will be improved by allowing cookies Cookie Policy