Telemetria próbna i identyfikacja modelu

W locie próbnym nie mierzy się „prawdy o pocisku”. Czujnik reaguje na wielkość fizyczną przez własną dynamikę, montaż i tor kondycjonowania; próbki otrzymują czas z niedoskonałego zegara; część danych przechodzi przez łącze, a część pozostaje w rejestratorze. Dopiero kontrolowana redukcja łączy te obserwacje z modelem. Dlatego tysiąc kanałów może dać słabszy dowód niż dziesięć kanałów o znanym znaczeniu, czasie i niepewności.

Teza artykułu brzmi: wiarygodność modelu lotu zależy od zachowania audytowalnego łańcucha przyczynowego — od pytania decyzyjnego, przez wielkość mierzoną, dynamikę kanału, wspólny czas i weryfikowalny zapis źródłowy, po identyfikację parametrów, niezależną walidację i kontrolowaną zmianę konfiguracji. Telemetria jest częścią tego łańcucha, a nie jego synonimem.

Próba pocisku różni się od wygodnego eksperymentu laboratoryjnego: trwa krótko, może być niepowtarzalna, egzemplarza często nie da się odzyskać, masa i warunki ruchu szybko się zmieniają, a udar lub utrata łącza występują właśnie wtedy, gdy informacja jest najcenniejsza. Plan dowodu musi więc powstać przed lotem.

granica bezpieczeństwa opisu. Przykłady używają modeli ogólnych oraz wielkości syntetycznych i znormalizowanych. Celowo pominięto częstotliwości, moce, kodowanie, klucze, rozmieszczenie anten, parametry manewrów i progi diagnostyczne konkretnych współczesnych systemów bojowych. Nie są one potrzebne do zrozumienia metrologii ani identyfikacji.

Trzy produkty jednej próby

Udany lot próbny pozostawia trzy odrębne produkty. Pierwszym jest zapis obserwacji: surowe słowa danych, metadane, konfiguracja i dowód integralności. Drugim jest rekonstrukcja zdarzeń oraz stanu: wspólna oś czasu, przeliczone jednostki, układy współrzędnych i estymaty wraz z niepewnością. Trzecim jest wniosek o modelu: parametry, reszty, zakres ważności i decyzja, czy model nadaje się do wskazanego zastosowania.

Autorski schemat pokazuje, że nie jest to jeden prosty przewód od czujnika do wykresu:

wielkość fizyczna
      │
      ▼
[czujnik + montaż] → [kondycjonowanie] → [próbkowanie + zegar]
                                               │
                         ┌─────────────────────┴─────────────────────┐
                         ▼                                           ▼
               [rejestrator pokładowy]                       [łącze → odbiór]
                         └─────────────────────┬─────────────────────┘
                                               ▼
                                  [niezmieniona kopia źródłowa]
                                               │
       [kalibracja + konfiguracja + TMATS] ────┤──── [dowód integralności]
                                               ▼
                          [redukcja wersjonowana + wspólna oś czasu]
                                               │
                 [radar / optyka / pogoda] ────┤
                                               ▼
                         [rekonstrukcja → identyfikacja → walidacja]
                                               ▼
                              [zakres ważności + decyzja]

Rozgałęzienie ma znaczenie diagnostyczne: zgodność zapisu pokładowego z odbiorem naziemnym sprawdza inną część łańcucha niż zgodność telemetrii z radarem. Skrót kryptograficzny potwierdza tożsamość bitową kopii, lecz nie dowodzi poprawnego montażu czujnika ani właściwej mapy czasu. Ilustracja nie przedstawia architektury żadnego konkretnego systemu.

Tych produktów nie wolno mieszać. Próbka z akcelerometru nie jest jeszcze przyspieszeniem środka masy; przefiltrowany przebieg nie jest surowym dowodem; dobre dopasowanie do danych identyfikacyjnych nie jest walidacją. Każda transformacja powinna mieć wejście, wersję algorytmu, parametry, autora lub wykonawcę, wynik i możliwość powtórzenia.

IRIG 106 porządkuje telemetrię, metadane i rejestrację pokładową jako współpracujące elementy standardu: rozdział 7 obejmuje telemetrię pakietową w dół, rozdział 9 standard wymiany atrybutów telemetrii (TMATS), rozdział 10 rejestrator pokładowy, a rozdział 11 format pakietów zapisu.¹ Nie znaczy to, że rozdziały można dowolnie mieszać w jednej implementacji. Dla niniejszego wywodu ważna jest zasada architektoniczna: znaczenie bitów, czas i konfiguracja kanału muszą podróżować razem z danymi.

Od decyzji do wymagania pomiarowego

Plan nie powinien zaczynać się od listy dostępnych czujników. Najpierw formułuje się decyzję: na przykład czy opóźnienie elementu wykonawczego mieści się w modelu, czy współczynnik lokalnej dynamiki wymaga aktualizacji albo która hipoteza najlepiej wyjaśnia utratę sygnału. Następnie określa się wielkość mierzoną, spodziewane pasmo, zakres, niepewność i wymagany związek czasowy.

Pełne wymaganie kanału odpowiada na pytania:

Element	Pytanie	Skutek zaniedbania
wielkość mierzona	jaka fizyczna cecha i w którym miejscu?	poprawne jednostki, lecz błędna interpretacja
pasmo i dynamika	jak szybkie zjawisko ma być zachowane?	tłumienie albo nakładanie widm (aliasing)
zakres i rozdzielczość	jakie wartości normalne i przejściowe?	obcięcie lub kwantowanie
czas	z czym i z jaką dokładnością porównujemy kanał?	fałszywe opóźnienie albo kolejność zdarzeń
niepewność	jaki błąd może mieć wynik i jak jest skorelowany?	pozornie dokładny parametr modelu
kryterium użycia	jaka decyzja zostanie podjęta?	zbieranie danych bez wartości dowodowej

Metrologiczne rozróżnienie wielkości mierzonej, kalibracji, spójności pomiarowej i niepewności omawiają Dusza, Gąsior i Tarapata.¹⁷ W praktyce karta kanału powinna wskazywać czujnik, orientację, miejsce, tor elektryczny, przeliczenie, częstotliwość próbkowania, filtr, zegar, opóźnienie, wersję konfiguracji i status kalibracji. Produktem rozdziału jest więc nie „lista kanałów”, lecz macierz pytanie–pomiar–decyzja.

Instrumentowany egzemplarz jest nową konfiguracją

Dodatkowy czujnik, wspornik, przewód, bateria i nadajnik zmieniają masę, bezwładność, sztywność, termikę, kompatybilność elektromagnetyczną (EMC) i czasem aerodynamikę. Egzemplarz próbny nie jest automatycznie reprezentatywnym pociskiem seryjnym z niewidzialną warstwą obserwacji. Trzeba udokumentować różnice i wykazać, które wnioski można przenieść.

Ocena równoważności powinna prowadzić od zmiany fizycznej do modelu:

$$\Delta\text{konfiguracji}\rightarrow \Delta(m,\mathbf I,\mathbf r,\text{sztywność},P,EMC)\rightarrow \Delta\text{odpowiedzi}\rightarrow\Delta\text{wniosku}.$$

Szczególnego znaczenia nabiera położenie czujnika względem środka masy i osi konstrukcji. Środek masy zmienia się podczas lotu, a ramię zamontowanego akcelerometru pozostaje związane ze strukturą. Rejestr konfiguracji powinien zatem zawierać geometrię montażu, orientację, numery seryjne, fotografie, masę aparatury, okablowanie, wersję programu lotnego i wersję opisu telemetrii.

Kanał rozwojowy nie musi mieć architektury kanału operacyjnego. Może być bogatszy i bardziej redundantny, ale jego wpływ na próbę musi być jawny. Z kolei logika pokładowa generująca „zdarzenie wykonane” jest obserwacją stanu programu, a nie niezależnym dowodem ruchu mechanizmu. Do rozstrzygnięcia przyczyn potrzebny jest kanał o innej fizycznej drodze.

Czujnik i montaż tworzą wspólną transmitancję

Rzeczywisty kanał można przedstawić jako kaskadę:

$$Y_m(s)=H_s(s)H_m(s)H_c(s)X(s)+N(s),$$

gdzie $H_s$ opisuje czujnik, $H_m$ — jego montaż, $H_c$ — kondycjonowanie i konwersję, a $N$ skupia zakłócenia. Kalibracja statyczna wyznaczająca nachylenie charakterystyki nie określa automatycznie pasma, fazy ani rezonansu. W badaniu krótkiego zdarzenia opóźnienie i tłumienie są częścią wyniku.

Akcelerometr mierzy siłę właściwą w swoim miejscu i osiach. Dla punktu odsuniętego o wektor $\mathbf r$ od środka masy składowa kinematyczna obejmuje

$$\mathbf a_p=\mathbf a_{CM}+\dot{\boldsymbol\omega}\times\mathbf r+ \boldsymbol\omega\times(\boldsymbol\omega\times\mathbf r).$$

$\mathbf a_{CM}$ oznacza przyspieszenie środka masy. Niepewność ramienia i orientacji przechodzi więc do wyniku. Żyroskop ma własny błąd stały i błąd skali; czujnik ciśnienia ma objętość przewodów i odpowiedź dynamiczną; termopara ma bezwładność i przewodzenie po przewodach; tensometr wymaga modelu drogi od lokalnego odkształcenia do obciążenia. Kanał położenia elementu wykonawczego pokazuje położenie w miejscu pomiaru, a nie zawsze rzeczywistą powierzchnię pod obciążeniem.

Kalibracja przed lotem buduje punkt odniesienia, a kontrola po próbie może ujawnić dryf lub uszkodzenie. Wynik po próbie nie powinien po cichu zastępować kalibracji użytej pierwotnie: obie wersje należy zachować, a wpływ różnicy na wniosek policzyć.

Próbkowanie może utracić informację nieodwracalnie

Przed przetwornikiem analogowo-cyfrowym sygnał wymaga zakresu, polaryzacji, separacji i filtru antyaliasingowego. Twierdzenie Nyquista nie mówi, że dwa punkty na okres wystarczą do wiarygodnego odtworzenia amplitudy i fazy przy rzeczywistym filtrze, zegarze i szumie. O projekcie toru decyduje pasmo zjawiska oraz błąd dopuszczalny dla decyzji.

Dla sinusoidy o częstotliwości $f_0$ próbkowanej z $f_s$ częstotliwość pozorna może wynosić

$$f_a=|f_0-kf_s|,$$

dla takiej liczby całkowitej $k$, która przenosi wynik do pasma podstawowego. Syntetyczny sygnał $f_0=0{,}65f_s$ pojawi się jako $0{,}35f_s$. Po zapisaniu próbek nie da się rozstrzygnąć, która składowa była prawdziwa, jeżeli analogowy filtr nie usunął pasma ponad granicą.

Próbkowanie multipleksowane dodaje przesunięcie między kanałami i czas ustalania po przełączeniu. Obcięcie sygnału nie zawsze ma postać idealnie płaskiego wierzchołka; może wystąpić wcześniej w wzmacniaczu, filtrze lub wejściu przetwornika. Rejestr powinien zawierać wskaźniki zakresu i błędów, a nie tylko przeliczone wartości.

Horodko opisuje cały tor akwizycji danych (DAQ) — próbnik z układem pamiętającym, multiplekser, przetwornik, filtrację, czas i analizę sygnałów.¹⁰ Misiurewicz pokazuje znaczenie nakładania widm, apertury, chwiejności czasu próbkowania i przesunięć kanałów dla stabilności rekonstrukcji.¹¹ Ich wspólny wniosek jest praktyczny: błąd powstały przed zapisem nie zostanie naprawiony bardziej złożonym filtrem po locie.

Czas jest mierzonym podsystemem

Każdy zegar ma przesunięcie, skalę i niestabilność. Mapowanie czasu lokalnego $t_i$ do osi odniesienia można w pierwszym przybliżeniu zapisać

$$t_{ref}=a_i+b_i t_i+\epsilon_i(t),$$

gdzie $a_i$ jest przesunięciem, $b_i$ błędem tempa, a $\epsilon_i$ resztkową niestabilnością. Pojedynczy znacznik synchronizuje przesunięcie, lecz nie wyznacza dryfu. Zdarzenia przed i po locie, wspólne impulsy lub korelacja charakterystycznych przebiegów mogą estymować oba parametry.

Należy rozdzielić czas wystąpienia zjawiska, chwilę pomiaru, znacznik próbki, umieszczenie w pakiecie, transmisję i odebranie. Opóźnienie stałe przesuwa fazę o

$$\varphi(f)=-2\pi f\tau,$$

a zmienność opóźnienia rozmywa relację. Bez budżetu czasu różnica między komendą a ruchem może zostać przypisana elementowi wykonawczemu, choć powstała w filtrze lub mapowaniu zegarów.

Produktem tej części jest mapa zegarów i opóźnień z niepewnością. Kryterium odrzucenia zależy od pytania: kanał może być wystarczający do rekonstrukcji wolnej temperatury, a bezużyteczny do estymacji krótkiego opóźnienia. „Wszystkie dane mają czas globalnego systemu pozycjonowania (GPS)” nie dowodzi, że znacznik odnosi się do chwili fizycznego próbkowania.

Rejestracja, łącze i metadane

Rejestrator pokładowy i łącze w dół odpowiadają na różne ryzyka. Łącze może zachować dane, gdy egzemplarz zostanie utracony, lecz cierpi na zaniki i ograniczoną przepływność. Rejestrator może zachować pełniejsze próbki, ale tylko jeżeli zostanie odzyskany i jego zapis przetrwa zdarzenie. Priorytetyzacja powinna wynikać z wartości dla decyzji: kanał niezbędny do rozróżnienia hipotez zasługuje na ochronę bardziej niż duplikat wygodnego wykresu.

Pakiet powinien przenosić identyfikator źródła, kolejność, czas, status jakości i kontrolę integralności. TMATS w rozdziale 9 IRIG 106 pokazuje rolę maszynowo czytelnego opisu konfiguracji telemetrii; rozdziały 10 i 11 wiążą go z zapisem oraz formatem danych rejestratora.¹ CCSDS 131.0-B-6 rozdziela zakres i zastosowanie od architektury (punkty 1.2–1.3 i 2), synchronizację ramki od synchronizacji słowa kodowego (punkty 3 i 9) oraz parametry zarządzane od samej implementacji (punkt 12). Jest źródłem porównawczym dla łączy kosmicznych, nie gotową specyfikacją poligonu.²

Kompresja bezstratna zachowuje próbki, lecz nadal wymaga wersji algorytmu i kontroli błędów. Kompresja stratna może zmienić amplitudę, fazę lub wartości odstające i dlatego musi być oceniana względem konkretnej analizy. Brak pakietu jest daną jakościową tylko wtedy, gdy wiadomo, czy powstał u źródła, w rejestratorze, radiu, kanale czy stacji naziemnej.

Metadane nie są dokumentem pomocniczym pisanym po locie. Bez wersji kalibracji, mapy bitów, jednostek, znaków, osi, czasu i stanu zakresu ciąg zer i jedynek pozostaje archiwalnie trwały, ale naukowo niemy.

Segment naziemny daje odniesienia, nie nieomylną prawdę

Radar, optyka, pogoda, odbiorniki czasu i aparatura naziemna tworzą obserwacje o innych błędach niż telemetria. Ich niezależność jest cenna, lecz każde źródło ma kalibrację, opóźnienie, geometrię, przerwy oraz warunki propagacji. Radarowa trajektoria nie jest automatycznie „prawdą” dla inercyjnego systemu nawigacyjnego (INS); może natomiast ujawnić błąd wspólny kanałów pokładowych.

Przed lotem trzeba przeprowadzić próbę całego toru: czujnik, zasilanie, kondycjonowanie, próbkowanie, znacznik czasu, pakiet, rejestrator, nadajnik, odbiór, dekomutację i archiwum. Sygnał zastępczy powinien pozwalać sprawdzić znak, skalę, czas, zakres i reakcję na błąd. Test samego nadajnika nie ujawnia zamienionych osi, a poprawny wykres w konsoli nie dowodzi, że surowe pliki i metadane dają się odtworzyć niezależnie.

Wspólny słownik zdarzeń określa, czym są „komenda”, „potwierdzenie”, „początek ruchu” i „utrata ramki”. Definicja operacyjna zapobiega sporom po próbie. Jeżeli dwa zespoły oznaczają inną chwilę tym samym słowem, zgodne wykresy mogą prowadzić do przeciwnych wniosków.

Zapis źródłowy i odtwarzalna redukcja

Pierwsza kopia surowych danych powinna pozostać niezmieniona i mieć skrót kryptograficzny, identyfikator nośnika oraz udokumentowany łańcuch pieczy. Każda kolejna warstwa — dekomutacja, przeliczenie, obrót osi, ponowne próbkowanie, filtracja i estymacja — powstaje jako nowy artefakt z jawną relacją do wejścia.

NIST IR 8387 zaleca w domenie dowodu cyfrowego dokumentowanie źródła i transferów, obliczanie skrótów, przechowywanie ich oddzielnie oraz utrzymywanie więcej niż jednej kopii.¹⁸ Nie jest to norma prowadzenia prób w locie, więc nie przenosi się z niej organizacji prawnej ani kryteriów dopuszczenia wyrobu. Przenosi się węższą metodę ochrony artefaktu cyfrowego. Niezgodność skrótu sygnalizuje zmianę bitową, lecz sama nie rozstrzyga, czy przyczyną była manipulacja, uszkodzenie nośnika czy dopuszczona transformacja zapisana pod błędną nazwą. Wniosek wymaga kopii odniesienia i historii operacji.

Schemat pochodzenia danych można zapisać:

$$D_0\xrightarrow{g_1,v_1,p_1}D_1\xrightarrow{g_2,v_2,p_2}\cdots\xrightarrow{g_n,v_n,p_n}D_n,$$

gdzie $g_i$ jest algorytmem, $v_i$ jego wersją, a $p_i$ zestawem parametrów. Odtwarzalność oznacza, że niezależny zespół uzyska ten sam wynik z $D_0$, a nie tylko podobny wykres z ręcznie wybranego fragmentu.

Zmiana siatki czasu, różniczkowanie i całkowanie zmieniają niepewność oraz korelację próbek. Szybka transformacja Fouriera (FFT) wymaga czasu obserwacji, okna i interpretacji przecieku; spektrogram ujawnia zmianę widma kosztem rozdzielczości. Misiurewicz przedstawia dyskretną i szybką transformację Fouriera (DFT/FFT), okna, korelację i filtry w ujęciu laboratoryjnym,¹² a Malanowski wiąże modulację, szum, detekcję i korelację z interpretacją kanału.¹³ Dobór narzędzia powinien wynikać z hipotezy, nie z atrakcyjności obrazu.

Jakość danych i budżet niepewności

Flaga „dane ważne” nie zamyka oceny jakości. Dla każdej próbki istotne są zakres, obcięcie, brakujące dane, stan kalibracji, wiek informacji, czas i kontekst pracy. Niepewność wyniku zależnego od wielkości $x_i$ można w przybliżeniu propagować przez macierz kowariancji:

$$u_y^2\approx \mathbf J\,\mathbf\Sigma_x\,\mathbf J^T, \qquad \mathbf J=\frac{\partial f}{\partial\mathbf x}.$$

Skorelowanych błędów nie wolno sumować jak niezależnych. Wspólny zegar może przesunąć wiele kanałów w tym samym kierunku, wspólna temperatura wpłynąć na kilka czujników, a ten sam parametr atmosfery wejść do rekonstrukcji i modelu. Zignorowanie korelacji daje zbyt wąski przedział parametru.

Wrzesiński omawia małe próby, ciśnienie dynamiczne, regresję nieliniową i propagację niepewności w laboratorium balistyki.¹⁴ Hyk i Stojek opisują regresję, reszty, obserwacje odstające i siłę dowodu statystycznego.¹⁵ Obserwacji odstającej nie usuwa się dlatego, że psuje dopasowanie. Najpierw sprawdza się fizyczną przyczynę, tor pomiarowy i regułę wykluczenia ustaloną przed analizą.

Rekonstrukcja stanu łączy obserwacje z dynamiką

Rekonstrukcja nie jest wygładzeniem wykresu. Łączy model kinematyczny i dynamiczny z czujnikami o różnych czasach, osiach oraz niepewnościach. Wektor stanu może obejmować położenie, prędkość, orientację, prędkość kątową oraz wybrane błędy stałe, ale jego skład wynika z decyzji.

W ogólnej postaci:

$$\dot{\mathbf x}=f(\mathbf x,\mathbf u,\boldsymbol\theta,t)+\mathbf w, \qquad \mathbf y=h(\mathbf x,t)+\mathbf v.$$

Filtr działający w czasie rzeczywistym używa danych dotychczasowych; wygładzanie po locie może korzystać z przyszłych obserwacji. Wyniki nie są równoważne i należy je nazwać. Model atmosfery, masa, geometria, położenie czujników oraz historia napędu wchodzą do rekonstrukcji z własną niepewnością.

Bilans sił i momentów służy kontroli zgodności, ale brak zgodności nie mówi od razu, który model jest błędny. Przyczyną może być czas, orientacja czujnika, siła nieujęta w modelu lub parametr aerodynamiczny. Tatjewski, opisując regulację, obserwatory i błąd predykcji, pokazuje ogólną konsekwencję niedopasowania modelu dla estymacji i sterowania.¹⁶

Studium I: fałszywa diagnoza elementu wykonawczego przez czas

Załóżmy znormalizowaną odpowiedź pierwszego rzędu na skok polecenia:

$$y(t)=1-\exp\left[-\frac{t-\tau_a}{T_a}\right],\qquad t\geq\tau_a,$$

gdzie $\tau_a$ jest fizycznym opóźnieniem, a $T_a$ stałą czasową. Kanał pomiarowy wnosi opóźnienie $\tau_c$, zaś błędna mapa zegarów — przesunięcie $\Delta t$. Z wykresu otrzymuje się wtedy

$$\tau_{poz}=\tau_a+\tau_c+\Delta t.$$

W syntetycznym przykładzie $\tau_a^*=0{,}20$, $\tau_c^*=0{,}06$, a $\Delta t^*=0{,}09$. Pozorne opóźnienie wynosi $0{,}35$, czyli jest o 75% większe od fizycznego opóźnienia elementu wykonawczego. Zespół patrzący wyłącznie na dwa przeliczone kanały może błędnie zgłosić zacięcie lub zbyt wolny napęd.

Skutek można zobaczyć na czterech syntetycznych próbkach. Przyjmijmy $T_a^*=0{,}30$ oraz łączne przesunięcie kanału i mapy czasu $\delta^*=\tau_c^*+\Delta t^*=0{,}15$. Czas $t_{zap}^*$ jest błędnym znacznikiem w zapisanym zbiorze, a obserwację utworzono z odpowiedzi dla $t_{zap}^*-\delta^*$ i małego, jawnego zaburzenia. Reszta ma postać $r_k=y_{obs,k}-y(t_k)$:

$t_{zap}^*$	$y_{obs}^*$	$r^*$ przed korektą czasu	$r^*$ po korekcie czasu
0,35	0,010	-0,383	+0,010
0,50	0,383	-0,249	-0,010
0,65	0,637	-0,140	+0,005
0,80	0,772	-0,093	-0,005

Średni błąd bezwzględny reszty spada z $0{,}2162$ do $0{,}0075$. Ważniejszy od samego spadku jest kształt: przed korektą wszystkie reszty mają ten sam znak i zanikają wraz ze zboczem odpowiedzi, po korekcie pozostaje jedynie założone zaburzenie. Taki wzorzec wspiera hipotezę błędu czasu, ale jej jeszcze nie dowodzi — podobny obraz może dać niedopasowana dynamika. Rozstrzygnięcie wymaga niezależnego wspólnego zdarzenia czasowego. Liczby nie opisują żadnego rzeczywistego urządzenia.

Rozstrzygnięcie wymaga kolejno: sprawdzenia chwili próbkowania kanału polecenia, identyfikacji opóźnienia filtru, wyznaczenia mapy zegarów ze wspólnego zdarzenia i ponownego dopasowania. Po korekcie reszty powinny stracić systematyczne przesunięcie. Produktem nie jest poprawiony obrazek, lecz nowa mapa czasu, wersja redukcji i obliczony wpływ na parametry modelu. Jeżeli niepewność $\Delta t$ jest porównywalna z $\tau_a$, kanał nie może samodzielnie rozstrzygnąć własności elementu wykonawczego.

Studium II: od akcelerometru do lokalnego modelu dynamiki

Rozważmy lokalny, małosygnałowy model dwóch stanów, bez wiązania go z parametrami konkretnego pocisku:

$$\dot{\mathbf x}=\mathbf A(\boldsymbol\theta)\mathbf x+ \mathbf B(\boldsymbol\theta)\mathbf u+\mathbf w, \qquad \mathbf y=\mathbf C\mathbf x+\mathbf D\mathbf u+\mathbf v.$$

Droga od słowa przetwornika do $\mathbf y$ obejmuje kalibrację, orientację, ramię czujnika, obrót osi, korektę czasu i budżet niepewności. Dopiero potem można estymować $\boldsymbol\theta$. Metoda błędu równania minimalizuje niezgodność w równaniach pochodnych i jest wrażliwa na sposób ich wyznaczania. Metoda błędu wyjścia symuluje model i minimalizuje

$$J(\boldsymbol\theta)=\sum_k \mathbf r_k^T\mathbf R_k^{-1}\mathbf r_k, \qquad \mathbf r_k=\mathbf y_k-\hat{\mathbf y}_k(\boldsymbol\theta).$$

Informacyjność eksperymentu opisuje macierz czułości $\mathbf S=\partial\hat{\mathbf y}/\partial\boldsymbol\theta$. Po zestawieniu wszystkich wybranych próbek jej wiersze odpowiadają obserwacjom, a kolumny estymowanym parametrom. Jeżeli $\mathbf R$ jest dodatnio określoną kowariancją całego wektora reszt, w lokalnym przybliżeniu otrzymuje się

$$\mathbf I_F=\mathbf S^T\mathbf R^{-1}\mathbf S,$$

Dla modelu regularnego i lokalnie poprawnego oraz estymatora nieobciążonego granica Craméra–Rao ma sens macierzowy:

$$\operatorname{cov}(\hat{\boldsymbol\theta})-\mathbf I_F^{-1}\succeq\mathbf 0.$$

Nie jest to uniwersalny przedział błędu. Jeśli $\mathbf R$ pomija korelacje, model jest strukturalnie błędny, estymator jest obciążony albo informacja jest osobliwa, prosta odwrotność nie daje deklarowanej granicy. Jeżeli kolumny $\mathbf S$ są niemal zależne, różne zestawy parametrów dają podobne wyjście: problem jest źle uwarunkowany, nawet przy małym błędzie dopasowania. Dokładna zależność kolumn oznacza lokalny brak identyfikowalności. Projekt eksperymentu ma wzbudzić rozróżnialne mody w dozwolonym zakresie, a nie po prostu wytworzyć duży sygnał. Morelli przedstawia iteracyjny cykl eksperyment–analiza i identyfikację w dziedzinie częstotliwości,⁶ projekt wejść wielokanałowych,⁷ kryteria informacyjności i granice Craméra–Rao⁸ oraz modele równoważne niskiego rzędu wyznaczane z danych pętli zamkniętej.⁹

Po dopasowaniu bada się reszty: średnią, autokorelację, związek z wejściem, stanem, temperaturą i czasem. Biała reszta nie dowodzi poprawności struktury, ale reszta skorelowana z wejściem ujawnia brakującą dynamikę lub błąd czasu. Model przechodzi próbę dopiero na oddzielnym manewrze lub warunku, który nie służył do identyfikacji. Morelli i Grauer ujmują projekt eksperymentu, poprawki pomiarowe, aktualizację modeli i niepewność jako jeden proces identyfikacji w locie.⁵

Studium III: utrata łącza podczas anomalii

Utrata telemetrii nie jest dowodem zniszczenia pocisku ani dowodem awarii nadajnika. Rozważmy cztery ogólne hipotezy: uszkodzenie toru radiowego, utrata zasilania aparatury, przejściowy zanik propagacyjny oraz zdarzenie obejmujące cały pojazd. Każda przewiduje inny wzór obserwacji.

Źródło	Tor radiowy	Utrata zasilania aparatury	Zanik propagacyjny	Zdarzenie całego pojazdu
stan pakietów przed zanikiem	możliwe błędy lub degradacja	możliwy zanik wielu kanałów zasilanych wspólnie	jakość łącza może pogarszać się bez zmian pokładowych	możliwe gwałtowne, współbieżne anomalie
rejestrator odzyskany	dane pokładowe mogą trwać	zapis może urwać się lub mieć znacznik zapadu napięcia	dane powinny trwać	wynik zależy od przeżycia rejestratora
radar i optyka	ruch może pozostać nominalny	ruch może pozostać nominalny	ruch pozostaje nominalny	możliwa zmiana lub utrata śladu
dziennik stacji	brak nośnej albo błędy demodulacji	podobny objaw radiowy bez rozstrzygnięcia przyczyny	zależność od geometrii i kanału	sam dziennik nie określa stanu pojazdu

Analiza zaczyna się od wspólnej osi czasu i zachowania negatywnego dowodu: brak obserwacji nie jest obserwacją zera. Następnie tworzy się macierz przewidywań i sprawdza niezależne tory. Pogoda, radar, optyka, rejestrator i dzienniki naziemne mają różne błędy, ale ich zgodność lub konflikt zawęża hipotezy.

Siłę tego zawężenia można aktualizować jawnie. Dla rozłącznych hipotez $H_i$ i nowego pakietu obserwacji $E$ obowiązuje

$$P(H_i\mid E)=\frac{P(E\mid H_i)P(H_i)}{\sum_j P(E\mid H_j)P(H_j)}.$$

Załóżmy wyłącznie dydaktycznie jednakowe prawdopodobieństwa początkowe oraz trzy obserwacje: odzyskany rejestrator kontynuuje zapis, śledzenie naziemne nie wykazuje zmiany ruchu, a jakość odbioru zmienia się wraz z geometrią. Po przyjęciu syntetycznych łącznych wiarygodności otrzymujemy:

Hipoteza	$P(E\mid H_i)$	$P(H_i\mid E)$
uszkodzenie toru radiowego	0,504	0,3997
utrata zasilania aparatury	0,072	0,0571
zanik propagacyjny	0,684	0,5424
zdarzenie całego pojazdu	0,001	0,0008

Wynik porządkuje dalsze sprawdzenia, ale nie zamienia zaniku propagacyjnego w pewnik: uszkodzenie toru radiowego nadal zachowuje znaczną masę prawdopodobieństwa. Liczby są sztuczne, a iloczyn trzech składowych wiarygodności zakłada ich warunkową niezależność. W rzeczywistej analizie wspólna stacja, pogoda lub zegar mogą skorelować obserwacje; użycie prostego iloczynu policzyłoby wówczas ten sam dowód więcej niż raz. Należy zastosować model zależności albo potraktować skorelowane dane jako jeden pakiet.
celowe pominięcie. Nie podajemy parametrów łącza ani sposobów zakłócania, odtwarzania kluczy czy optymalizacji anten. Dla rozumowania wystarcza rozdzielenie źródła, transportu, odbioru i niezależnej obserwacji. Produktem studium jest uszeregowana lista hipotez z siłą dowodów i brakami, nie pewna narracja wymuszona przez pojedynczy urwany wykres.

Walidacja, zakres ważności i aktualizacja modelu

Dane identyfikacyjne służą do wyznaczenia parametrów; dane walidacyjne sprawdzają predykcję. Użycie tego samego przebiegu do obu celów mierzy zdolność dopasowania, nie zdolność przewidywania. Walidacja powinna obejmować wielkości bezpośrednio związane z planowanym zastosowaniem oraz przypadki, w których konkurencyjne modele dają różne odpowiedzi.

NASA-STD-7009 oraz jego podręcznik wiążą wiarygodność modelu z zamierzonym zastosowaniem, weryfikacją i walidacją (V&V), niepewnością, czułością i dowodami.³^,⁴ Model może być wystarczający do rekonstrukcji jednego lotu i niewystarczający do przewidywania poza zbadanym zakresem. Granice powinny obejmować konfigurację, warunki, zakres stanów i wejść, jakość danych oraz wersję algorytmu.

Aktualizacja modelu musi zachować poprzednią wersję, dane, funkcję kosztu, ograniczenia i wynik walidacji. Poprawa jednego parametru nie może maskować błędu czujnika albo czasu. Analiza czułości wskazuje, które niepewności naprawdę zmieniają decyzję; ekstrapolacja poza zakres obserwacji wymaga dodatkowego dowodu, a nie tylko węższego przedziału z regresji.

Pakiet dowodowy i decyzja końcowa

Minimalny pakiet po próbie obejmuje: niezmienione dane surowe i ich skróty, metadane, dokładną konfigurację egzemplarza, mapę zegarów, kalibracje, dzienniki błędów, warunki naziemne, procedury redukcji, kod programów i jego wersje, wyniki rekonstrukcji z kowariancją, model z parametrami, reszty, niezależną walidację oraz rejestr decyzji. Niezależny przegląd powinien móc odtworzyć nie tylko wykres, lecz także drogę od pytania do wniosku.

Każdy rozdział pozostawia artefakt i kryterium odrzucenia. Niezdefiniowana wielkość dyskwalifikuje interpretację kanału. Nieznana saturacja dyskwalifikuje amplitudę. Niepewny czas dyskwalifikuje wniosek o krótkim opóźnieniu. Brak identyfikowalności dyskwalifikuje rozdzielenie parametrów. Dopasowanie bez niezależnej walidacji dyskwalifikuje twierdzenie o zdolności predykcyjnej.

Trzy studia pokazują wspólną zasadę. Pozornie wolny element wykonawczy okazał się błędem łańcucha czasu. Ładny model lokalnej dynamiki wymagał zdefiniowanego czujnika, informacyjnego eksperymentu i analizy reszt. Utrata łącza stała się rozstrzygalna dopiero po zestawieniu niezależnych źródeł. W każdym przypadku wartość telemetrii wynikała nie z liczby próbek, lecz z zachowania mierzalnego, nieprzerwanego i audytowalnego łańcucha dowodu.

Źródła

Przypisy

Range Commanders Council, IRIG 106-24 Telemetry Standards, October 2024, https://www.irig106.org/wiki/irig_106-24, dostęp: 2026-09-10. 7: telemetria pakietowa w dół; rozdz. 9: TMATS; rozdz. 10: cyfrowy rejestrator pokładowy; rozdz. 11: format pakietów zapisu.
↩1 ↩2
CCSDS, TM Synchronization and Channel Coding, CCSDS 131.0-B-6, kwiecień 2026, https://ccsds.org/view/bluebooks/entry/4803/, dostęp: 2026-09-10.2–1.3: zakres i zastosowanie; pkt 2: architektura; pkt 3 i 9: synchronizacja; pkt 12: parametry zarządzane; wyłącznie źródło porównawcze dla łączy kosmicznych.
↩
NASA, Standard for Models and Simulations, NASA-STD-7009B, 2024, https://standards.nasa.gov/standard/nasa/nasa-std-7009, dostęp: 2026-09-10.
↩
NASA, NASA Handbook for Models and Simulations: An Implementation Guide for NASA-STD-7009B, NASA-HDBK-7009B, 2026, https://standards.nasa.gov/standard/nasa/nasa-hdbk-7009, dostęp: 2026-09-10.
↩
Eugene A. Morelli i Jared A. Grauer, Advances in Aircraft System Identification at NASA Langley Research Center, 2023, https://ntrs.nasa.gov/citations/20230001945, dostęp: 2026-09-10.
↩
Eugene A. Morelli, In-Flight System Identification, AIAA Paper 98-4261, 1998, https://ntrs.nasa.gov/citations/20040090518, dostęp: 2026-09-10.
↩
Eugene A. Morelli, Practical Multiple-Input Design for Aircraft System Identification Flight Tests, NASA, 2021, https://ntrs.nasa.gov/citations/20210018281, dostęp: 2026-09-10.
↩
Eugene A. Morelli, Advances in Experiment Design for High Performance Aircraft, NASA, 1999, https://ntrs.nasa.gov/citations/19990032455, dostęp: 2026-09-10.
↩
Eugene A. Morelli, Identification of Low Order Equivalent System Models From Flight Test Data, NASA/TM-2000-210117, https://ntrs.nasa.gov/citations/20000085030, dostęp: 2026-09-10.
↩
Longin Horodko, Systemy pomiarowe i podstawy analizy sygnałów.
↩
Jacek Misiurewicz, Nowoczesne próbkowanie w cyfrowym przetwarzaniu sygnałów.
↩
Jacek Misiurewicz, Laboratorium cyfrowego przetwarzania sygnałów.
↩
Mateusz Malanowski, Laboratorium sygnałów, modulacji i systemów.
↩
Zbigniew Wrzesiński, Niepewność pomiarów w laboratorium balistyki.
↩
Wojciech Hyk i Zbigniew Stojek, Analiza statystyczna w laboratorium badawczym.
↩
Piotr Tatjewski, Algorytmy regulacji obiektów przemysłowych. PID i MPC.
↩
Jacek Dusza, Paweł Gąsior i Grzegorz Tarapata, Podstawy pomiarów.
↩
Barbara Guttman, Douglas R. White, Shannan Williams i Tracy Walraven, Digital Evidence Preservation: Considerations for Evidence Handlers, NIST IR 8387, 2022, https://doi.org/10.6028/NIST.IR.8387, dostęp: 2026-09-10; zastosowanie do danych próbnych ograniczono do metody ochrony artefaktu.
↩