Christiaan C. J. M. Tiberius
Current Organizations / Institutes
The organizations/institutes from the scholar's latest article
Research Subjects
Christiaan C. J. M. Tiberius has a total of 102 co-authors. Here are the top co-authors by number of shared publications:
- Peter J. G. Teunissen
Delft University of Technology
16 shared publications Gerard J. M. Janssen
Signal Processing Systems group, Delft University of Technology, Delft, The Netherlands
14 shared publicationsJunlin Yan
Intel Mobile Communications GmbH, Dusseldorfer Landstr. 401, D-47259 Duisburg, Germany
13 shared publicationsGiovanni Bellusci
Xsens, Pantheon 6a, 7521 PR, Enschede, Netherlands
12 shared publications- Hans van der Marel
Department of Geoscience and Remote Sensing, Delft University of Technology, Delft, The Netherlands
6 shared publications - Gerard J. M. Janssen
Delft University of Technology, Delft, The Netherlands
5 shared publications
Latest Publications
- Open Access
- 6 May 2026
- Published by Institute of Electrical and Electronics Engineers (IEEE) in IEEE Access
- Vol. 14, 70581-70599
On OFDM Ranging Performance Degradation in Multipath Scenarios: Bias and Misspecified Cramér–Rao Bounds
Multipath propagation represents a dominant error source limiting the accuracy of radio-based ranging and positioning in urban environments. Conventional ranging estimators typically estimate only the range of the first arriving path while neglecting the multipath components. This introduces a bias and causes the estimator’s variance to diverge from the Cramér-Rao Bound (CRB). While extensive research has established performance bounds for estimators that jointly estimate all paths, the specific impact on the variance of ignoring these components remains largely unexplored. This paper derives bounds on the bias, variance, and mean square error (MSE) for ranging in multipath channels for the purpose of positioning. We consider Orthogonal Frequency Division Multiplexing (OFDM) signals in multipath channels with deterministic path delays and gains. The paper focuses on time-delay estimation (TDE) in the mid-tohigh signal-to-noise ratio (SNR) regime under two scenarios: (i) when the receiver jointly estimates all propagation paths, and (ii) when the receiver underestimates the number of paths. Specifically, for the latter case, we derive a semi closed-form expression for the variance of the time-delay estimator that considers a single-path when in reality there are L paths. This derivation is based on the misspecified Cramér-Rao bound (MCRB). Although unmodelled multipath has been traditionally viewed as detrimental to time-delay estimation, we reveal that in some cases the estimation variance improves. For a two-path channel we show that the variance depends on the relative gain, carrier phase, and separation of the paths. Additionally, the estimation bias is upper-bounded by constructive and destructive path interference. Finally, we empirically process a six-path simulated channel frequency response with fixed path delays and gains, to demonstrate that the derived MSE bound is tight for mid-to-high SNRs. This work characterizes the impact of multipath propagation on the variance of time-delay estimation, which is essential for designing accurate ranging signals and estimators in urban scenarios.
- 6 May 2026
- Published by Institute of Electrical and Electronics Engineers (IEEE) in IEEE Access
- Vol. 14, 70581-70599
Open AccessOn OFDM Ranging Performance Degradation in Multipath Scenarios: Bias and Misspecified Cramér–Rao Bounds
Multipath propagation represents a dominant error source limiting the accuracy of radio-based ranging and positioning in urban environments. Conventional ranging estimators typically estimate only the range of the first arriving path while neglecting the multipath components. This introduces a bias and causes the estimator’s variance to diverge from the Cramér-Rao Bound (CRB). While extensive research has established performance bounds for estimators that jointly estimate all paths, the specific impact on the variance of ignoring these components remains largely unexplored. This paper derives bounds on the bias, variance, and mean square error (MSE) for ranging in multipath channels for the purpose of positioning. We consider Orthogonal Frequency Division Multiplexing (OFDM) signals in multipath channels with deterministic path delays and gains. The paper focuses on time-delay estimation (TDE) in the mid-tohigh signal-to-noise ratio (SNR) regime under two scenarios: (i) when the receiver jointly estimates all propagation paths, and (ii) when the receiver underestimates the number of paths. Specifically, for the latter case, we derive a semi closed-form expression for the variance of the time-delay estimator that considers a single-path when in reality there are L paths. This derivation is based on the misspecified Cramér-Rao bound (MCRB). Although unmodelled multipath has been traditionally viewed as detrimental to time-delay estimation, we reveal that in some cases the estimation variance improves. For a two-path channel we show that the variance depends on the relative gain, carrier phase, and separation of the paths. Additionally, the estimation bias is upper-bounded by constructive and destructive path interference. Finally, we empirically process a six-path simulated channel frequency response with fixed path delays and gains, to demonstrate that the derived MSE bound is tight for mid-to-high SNRs. This work characterizes the impact of multipath propagation on the variance of time-delay estimation, which is essential for designing accurate ranging signals and estimators in urban scenarios.
- Open Access
- 20 April 2026
- Published by TU Delft OPEN Publishing
Engineering Signal Analysis: from Fourier to filtering: exercises
This book, Exercises, contains hundreds of exercises, including answers and worked examples, for studying and practicing the theory. Python scripts illustrate how to perform spectral signal analysis on a computer.
- 20 April 2026
- Published by TU Delft OPEN Publishing
Open AccessEngineering Signal Analysis: from Fourier to filtering: exercises
This book, Exercises, contains hundreds of exercises, including answers and worked examples, for studying and practicing the theory. Python scripts illustrate how to perform spectral signal analysis on a computer.
- Open Access
- 9 March 2026
- Published by Springer Nature in GPS Solutions
- Vol. 30 (2) , 87
DIA-Estimator and Multidimensional Model Misspecifications: GNSS-based Positioning Safety Analysis for UAVs
The Detection, Identification, and Adaptation (DIA)-estimator integrates parameter estimation and hypothesis testing for model misspecifications. This contribution presents a positioning safety analysis approach grounded in the DIA-estimator framework, with a particular emphasis on multidimensional model misspecifications, such as simultaneous outliers in the observations. While recent work has focused on the performance of the detection and identification of multidimensional model misspecifications, we turn our attention to how they affect the probability density function (PDF) of the DIA-estimator and, consequently, the probability of positioning failure–an indicator relevant for safety-of-life applications (e.g., automotive, aviation, rail, maritime). This work formulates and quantifies the probability of positioning failure and its conditional components. A representative simulation-based study is presented for a UAV equipped with a GPS receiver configured to achieve performance comparable to Technical Standard Order (TSO)-certified receivers. The analysis is carried out for two scenarios: a fixed GPS satellite geometry at a single time snapshot, and for a varying GPS satellite geometry over a 24-hour period over an authorized UAV airspace region in the Netherlands using real satellite ephemeris data. Together, these scenarios provide insights into the structure of the DIA-estimator’s PDF, such as multimodality and orientation with respect to the chosen positioning safety region, and support comprehensive evaluation of positioning safety. Although the current focus is on GPS-based positioning, the presented approach is general and can be extended to include multisensor configurations, additional GNSS constellations, and applied to other safety-critical applications, which are subjects of future work.
- 9 March 2026
- Published by Springer Nature in GPS Solutions
- Vol. 30 (2) , 87
Open AccessDIA-Estimator and Multidimensional Model Misspecifications: GNSS-based Positioning Safety Analysis for UAVs
The Detection, Identification, and Adaptation (DIA)-estimator integrates parameter estimation and hypothesis testing for model misspecifications. This contribution presents a positioning safety analysis approach grounded in the DIA-estimator framework, with a particular emphasis on multidimensional model misspecifications, such as simultaneous outliers in the observations. While recent work has focused on the performance of the detection and identification of multidimensional model misspecifications, we turn our attention to how they affect the probability density function (PDF) of the DIA-estimator and, consequently, the probability of positioning failure–an indicator relevant for safety-of-life applications (e.g., automotive, aviation, rail, maritime). This work formulates and quantifies the probability of positioning failure and its conditional components. A representative simulation-based study is presented for a UAV equipped with a GPS receiver configured to achieve performance comparable to Technical Standard Order (TSO)-certified receivers. The analysis is carried out for two scenarios: a fixed GPS satellite geometry at a single time snapshot, and for a varying GPS satellite geometry over a 24-hour period over an authorized UAV airspace region in the Netherlands using real satellite ephemeris data. Together, these scenarios provide insights into the structure of the DIA-estimator’s PDF, such as multimodality and orientation with respect to the chosen positioning safety region, and support comprehensive evaluation of positioning safety. Although the current focus is on GPS-based positioning, the presented approach is general and can be extended to include multisensor configurations, additional GNSS constellations, and applied to other safety-critical applications, which are subjects of future work.
- Open Access
- 1 February 2026
- Published by TU Delft OPEN Publishing
Engineering Signal Analysis: from Fourier to filtering: theory
Engineering Signal Analysis is an introductory textbook on the analysis of signals in time and frequency domain. It discusses how to characterize, model, analyze, interpret and operate on signals in time and frequency. In the first part of the book the basic theoretical concepts are introduced in continuous time, covering the Fourier Series and Fourier Transform. The second part introduces discrete-time signals, addressing sampling and finite signal duration, and their implications on spectral analysis. The third part elaborates on spectral estimation, covering the basic periodogram and more advanced methods. The fourth part concludes the book with an introduction to linear system theory, and addresses in particular signal filtering.
- 1 February 2026
- Published by TU Delft OPEN Publishing
Open AccessEngineering Signal Analysis: from Fourier to filtering: theory
Engineering Signal Analysis is an introductory textbook on the analysis of signals in time and frequency domain. It discusses how to characterize, model, analyze, interpret and operate on signals in time and frequency. In the first part of the book the basic theoretical concepts are introduced in continuous time, covering the Fourier Series and Fourier Transform. The second part introduces discrete-time signals, addressing sampling and finite signal duration, and their implications on spectral analysis. The third part elaborates on spectral estimation, covering the basic periodogram and more advanced methods. The fourth part concludes the book with an introduction to linear system theory, and addresses in particular signal filtering.
- 1 January 2026
- Published by Institute of Electrical and Electronics Engineers (IEEE) in IEEE Transactions on Vehicular Technology
- Vol. 75 (9) , 1-15
Probability of Positioning Failure for UAVs in Multiple Authorized European Airspace Regions
- 1 January 2026
- Published by Institute of Electrical and Electronics Engineers (IEEE) in IEEE Transactions on Vehicular Technology
- Vol. 75 (9) , 1-15
Probability of Positioning Failure for UAVs in Multiple Authorized European Airspace Regions
- 18 November 2025
- Published by Springer Nature
- p. 45-59
Parameter Estimation in Terrestrial Networked Positioning and Timing
- 18 November 2025
- Published by Springer Nature
- p. 45-59
Parameter Estimation in Terrestrial Networked Positioning and Timing
- Open Access
- 12 November 2025
- Published by Institute of Electrical and Electronics Engineers (IEEE) in IEEE Transactions on Aerospace and Electronic Systems
- Vol. 62 (00189251) , 15-25
GNSS Ambiguity-Resolved Detector: Implementation With a Lookup Table
The ambiguity-resolved detector (ARD) is developed for the validation of high-precision Global Navigation Satellite Systems (GNSS) mixed-integer observation models. The ARD releases the high success rate restriction of resolving the float ambiguities to integers, enabling ambiguity resolution to contribute to model validation even if the resolution success rate is not close to one. However, there are no closed-form expressions for the distribution of the ARD test statistic. The ARD critical value can only be obtained by Monte Carlo simulation, which requires heavy computations. To reduce the effort of applying the ARD, we provide a lookup table for the critical value of the ARD that utilizes the integer least-squares (ILS) estimator to resolve the ambiguities ($\mathbf {ARD_{ILS}}$). We obtain 263 408 $\mathbf {ARD_{ILS}}$ critical values corresponding to GNSS observation models that vary in the precision of observables, satellite geometry, constellation, and the number of signal frequencies. We treat critical values as functions of the integer bootstrapping success rate and conduct curve fitting with second-order polynomials. The polynomial parameters are provided in the lookup table with which the $\mathbf {ARD_{ILS}}$ critical values can be simply computed rather than simulated. Numerical experiments demonstrate that the critical values computed with the lookup table provide significance levels close to the specified ones. The $\mathbf {ARD_{ILS}}$ can be implemented easily and efficiently with the lookup table. The lookup table with the polynomial coefficients to compute the ARD critical values can be downloaded from (Yin, 2024).- 12 November 2025
- Published by Institute of Electrical and Electronics Engineers (IEEE) in IEEE Transactions on Aerospace and Electronic Systems
- Vol. 62 (00189251) , 15-25
Open AccessGNSS Ambiguity-Resolved Detector: Implementation With a Lookup Table
The ambiguity-resolved detector (ARD) is developed for the validation of high-precision Global Navigation Satellite Systems (GNSS) mixed-integer observation models. The ARD releases the high success rate restriction of resolving the float ambiguities to integers, enabling ambiguity resolution to contribute to model validation even if the resolution success rate is not close to one. However, there are no closed-form expressions for the distribution of the ARD test statistic. The ARD critical value can only be obtained by Monte Carlo simulation, which requires heavy computations. To reduce the effort of applying the ARD, we provide a lookup table for the critical value of the ARD that utilizes the integer least-squares (ILS) estimator to resolve the ambiguities ($\mathbf {ARD_{ILS}}$). We obtain 263 408 $\mathbf {ARD_{ILS}}$ critical values corresponding to GNSS observation models that vary in the precision of observables, satellite geometry, constellation, and the number of signal frequencies. We treat critical values as functions of the integer bootstrapping success rate and conduct curve fitting with second-order polynomials. The polynomial parameters are provided in the lookup table with which the $\mathbf {ARD_{ILS}}$ critical values can be simply computed rather than simulated. Numerical experiments demonstrate that the critical values computed with the lookup table provide significance levels close to the specified ones. The $\mathbf {ARD_{ILS}}$ can be implemented easily and efficiently with the lookup table. The lookup table with the polynomial coefficients to compute the ARD critical values can be downloaded from (Yin, 2024). - Open Access
- 30 May 2025
- Published by Institute of Electrical and Electronics Engineers (IEEE) in IEEE Transactions on Vehicular Technology
- Vol. 74 (10) , 15238-15253
A Method to Compute the Probability of Positioning Failure for Vehicles in the Context of Dependence Between Parameter Estimation and Statistical Hypothesis Testing
Positioning technologies are widely used in automotive, aviation, rail, and maritime safety-critical applications. Therefore, the computation of the probability of positioning failure for vehicles, which is the probability that the position estimator is outside a safety region, is of interest for positioning safety analyses. Since parameter estimation and statistical hypothesis testing for model misspecifications are commonly employed in positioning algorithms, the resulting position estimator is conditioned on the statistical hypothesis testing outcome. Hence, the probability density function (PDF) of the vehicle position estimator that accounts for the dependence between the two inference concepts should be used in the computations. In this contribution, we propose a method to compute the probability of positioning failure using the PDF of the vehicle position estimator, which accounts for the aforementioned dependence and is based on rare event simulation techniques, specifically Importance Sampling and the Cross-Entropy method. We apply the proposed method to a satellite-based positioning scenario, in decimeter precision, of an automated vehicle. The results show that the proposed method enables extensive positioning safety-analyses giving insights that can be used in the development of positioning algorithms and deciding whether safety targets and/or requirements are met. Finally, we discuss some limitations of the method and propose several further improvements.
- 30 May 2025
- Published by Institute of Electrical and Electronics Engineers (IEEE) in IEEE Transactions on Vehicular Technology
- Vol. 74 (10) , 15238-15253
Open AccessA Method to Compute the Probability of Positioning Failure for Vehicles in the Context of Dependence Between Parameter Estimation and Statistical Hypothesis Testing
Positioning technologies are widely used in automotive, aviation, rail, and maritime safety-critical applications. Therefore, the computation of the probability of positioning failure for vehicles, which is the probability that the position estimator is outside a safety region, is of interest for positioning safety analyses. Since parameter estimation and statistical hypothesis testing for model misspecifications are commonly employed in positioning algorithms, the resulting position estimator is conditioned on the statistical hypothesis testing outcome. Hence, the probability density function (PDF) of the vehicle position estimator that accounts for the dependence between the two inference concepts should be used in the computations. In this contribution, we propose a method to compute the probability of positioning failure using the PDF of the vehicle position estimator, which accounts for the aforementioned dependence and is based on rare event simulation techniques, specifically Importance Sampling and the Cross-Entropy method. We apply the proposed method to a satellite-based positioning scenario, in decimeter precision, of an automated vehicle. The results show that the proposed method enables extensive positioning safety-analyses giving insights that can be used in the development of positioning algorithms and deciding whether safety targets and/or requirements are met. Finally, we discuss some limitations of the method and propose several further improvements.
- Open Access
- 8 March 2025
- Published by Springer Nature in GPS Solutions
- Vol. 29 (2) , 1-17
Performance of ambiguity-resolved detector for GNSS mixed-integer model
Teunissen (J Geod 98(83):1–16, 2024) proposed the ambiguity-resolved (AR) detection theory for GNSS mixed-integer model validation. In this contribution, we study the performance of the AR detector through analysis and simulation experiments and compare it with the ambiguity-float (AF) and ambiguity-known (AK) detectors. We describe how the detectors can be implemented and how to evaluate their performance by computing the power as functions of the model misspecifications’ size. We present two simulation experiments with single- and dual-frequency GPS models and demonstrate that the AR detector can provide a larger detection power than the AF detector, even if the success rate is not close to one. Then, we obtain power functions over 25 user locations with five observation models and 72 satellite geometries per location per model. We find that the AR detector increases the detection probability of ionosphere and troposphere delays by 47% and 60% on average when the success rate is larger than 97.5% and the level of significance is 0.01. We also find the AR detection power to be larger than that of the AF detector in case of multi-dimensional misspecifications.
- 8 March 2025
- Published by Springer Nature in GPS Solutions
- Vol. 29 (2) , 1-17
Open AccessPerformance of ambiguity-resolved detector for GNSS mixed-integer model
Teunissen (J Geod 98(83):1–16, 2024) proposed the ambiguity-resolved (AR) detection theory for GNSS mixed-integer model validation. In this contribution, we study the performance of the AR detector through analysis and simulation experiments and compare it with the ambiguity-float (AF) and ambiguity-known (AK) detectors. We describe how the detectors can be implemented and how to evaluate their performance by computing the power as functions of the model misspecifications’ size. We present two simulation experiments with single- and dual-frequency GPS models and demonstrate that the AR detector can provide a larger detection power than the AF detector, even if the success rate is not close to one. Then, we obtain power functions over 25 user locations with five observation models and 72 satellite geometries per location per model. We find that the AR detector increases the detection probability of ionosphere and troposphere delays by 47% and 60% on average when the success rate is larger than 97.5% and the level of significance is 0.01. We also find the AR detection power to be larger than that of the AF detector in case of multi-dimensional misspecifications.
- Open Access
- 24 January 2025
- Published by Institute of Electrical and Electronics Engineers (IEEE) in IEEE Transactions on Intelligent Transportation Systems
- Vol. 26 (4) , 5509-5521
Dependence Between Parameter Estimation and Statistical Hypothesis Testing: Positioning Safety Analysis for Automated/Autonomous Vehicles
The analysis of positioning safety often employs a probability-based formulation. This approach quantifies the probability of positioning failure, which is the probability of the position estimator being outside a safety-region, and compares it against an application specific requirement. The design of positioning algorithms for safety-critical applications, such as automated/autonomous vehicles, should consider the dependence between parameter or state estimation and statistical hypothesis testing for model misspecifications in the evaluation of positioning safety. If this dependence is not considered, as this article shows, the conclusions drawn from the positioning safety analysis might be overly-optimistic. Therefore, this article focuses on the aforementioned dependence through a vehicle positioning scenario based on an Extended Kalman Filter (EKF) and the Detection, Identification, and Adaptation (DIA) method for misspecifications in the motion and measurement models. Grounded in the distributional theory for the DIA method, our positioning safety analysis utilizes the conditional probability density functions (PDFs) of the combined EKF and DIA position error, which are generally nonnormal. We compute the probability of vehicle positioning failure in two cases 1) when the dependence is considered and 2) when it is not, to quantify the over-optimism introduced by ignoring this dependence. Finally, we present our conclusions and recommendations.
- 24 January 2025
- Published by Institute of Electrical and Electronics Engineers (IEEE) in IEEE Transactions on Intelligent Transportation Systems
- Vol. 26 (4) , 5509-5521
Open AccessDependence Between Parameter Estimation and Statistical Hypothesis Testing: Positioning Safety Analysis for Automated/Autonomous Vehicles
The analysis of positioning safety often employs a probability-based formulation. This approach quantifies the probability of positioning failure, which is the probability of the position estimator being outside a safety-region, and compares it against an application specific requirement. The design of positioning algorithms for safety-critical applications, such as automated/autonomous vehicles, should consider the dependence between parameter or state estimation and statistical hypothesis testing for model misspecifications in the evaluation of positioning safety. If this dependence is not considered, as this article shows, the conclusions drawn from the positioning safety analysis might be overly-optimistic. Therefore, this article focuses on the aforementioned dependence through a vehicle positioning scenario based on an Extended Kalman Filter (EKF) and the Detection, Identification, and Adaptation (DIA) method for misspecifications in the motion and measurement models. Grounded in the distributional theory for the DIA method, our positioning safety analysis utilizes the conditional probability density functions (PDFs) of the combined EKF and DIA position error, which are generally nonnormal. We compute the probability of vehicle positioning failure in two cases 1) when the dependence is considered and 2) when it is not, to quantify the over-optimism introduced by ignoring this dependence. Finally, we present our conclusions and recommendations.