Azzera filtri
Azzera filtri

noise calculation in the gps-simulator

7 visualizzazioni (ultimi 30 giorni)
Siyu
Siyu il 22 Giu 2023
Commentato: Ryan Salvo il 26 Giu 2023
Hi,
I am trying to understand the principle of gps sensor simulator in Matlab. According to the documentation, the horizontal position accuracy specifies the standard deviation of the noise in the horizontal position measurement and the vertical position accuracy specifies the standard deviation of the noise in the vertical position measurement. From my understanding, the measurement noise matrix be
. Where I assumed that the the standard deviation in x and y directions are the same.
But after reading the source code of matlab implementation, I am wondering where comes from the factor 2. The relevant source code is shown in the below.
dt = 1 ./ obj.SampleRate;
% in case of white noise process the decayfactor is 0.
decayFactor = obj.DecayFactor;
% tau = dt for white noise process
tau = dt / (1-decayFactor);
horzSigma = obj.HorizontalPositionAccuracy;
vertSigma = obj.VerticalPositionAccuracy;
sigmas = [horzSigma horzSigma vertSigma];
% updated standard deviation
% it equals to obj.pSigmaScaled = sigmas .* sqrt(2);
obj.pSigmaScaled = sigmas .* sqrt(2.*dt./tau);
obj.pPositionErrorFilterNum = 1;
obj.pPositionErrorFilterDen = [1 -decayFactor];
According to the implementation, the measurement noise is:
.
I don't know where I got it wrong. Is there someone could help me about it?
Thanks for any help and hints in advance.
Best regards

Risposte (2)

Ryan Salvo
Ryan Salvo il 23 Giu 2023
Hi Siyu,
For the noise calculation, the standard deviation is scaled by a bandwidth equal to 1/2 of the sampling rate. That results in the multiplication of the standard deviation by 2*dt.
Thanks,
Ryan

Siyu
Siyu il 26 Giu 2023
Hi Ryan,
thanks for your reply.
I can understand that acoording to the Nyquist Sampling Theorem the sampling rate is 2 times of the noise bandwidth. But I am still not clear what is the relationship between the noise bandwidth and the rescaling. Further, why should the term be divided to calculate the scaled sigma? Maybe I am missing some background information. Could you provide me more detailed information about it?
Thanks for your help in adavance.
Best
Siyu
  1 Commento
Ryan Salvo
Ryan Salvo il 26 Giu 2023
Hi Siyu,
We use the generic sensor noise model detailed in the TSim manual to model this random noise in the GPS sensor. See the Derivation of Modified Rate Random Walk Noise Model section for the complete details.
Thanks,
Ryan

Accedi per commentare.

Prodotti


Release

R2023a

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by