NaveGo icon indicating copy to clipboard operation
NaveGo copied to clipboard

Question in generating the gyro about g_err_b

Open chengjunlong0524 opened this issue 3 years ago • 6 comments

g_err_b(i,:) = ( omega_in_b * gyro_b (i,:)' )' or gyro_err_b(i,:) = ( omega_in_b )';

chengjunlong0524 avatar Aug 16 '22 15:08 chengjunlong0524

omega_in_b is the projection of both Earth and transport rates in the body frame.

It is the inverse operation of:

wb = (wb - DCMbn' * (om_ie_n + om_en_n)); % Titterton, Eq. 3.29, p. 32

in att_update().

rodralez avatar Aug 16 '22 20:08 rodralez

邮件已收到!!!

chengjunlong0524 avatar Aug 16 '22 20:08 chengjunlong0524

Outputs from three ideal orthogonal gyroscopes are expressed as:

omega_in_b is the projection of both Earth and transport rates in the body frame, ie

What does the "g_err_b" means in code? The finall output is

wb_sim = gyro_b + g_err_b + g_wn + gb_sta + gb_dyn + g_rrw

gyro_b is same as omega_in_b? What does the code "g_err_b(i,:) = ( omega_in_b * gyro_b (i,:)' )'" means?

chengjunlong0524 avatar Aug 17 '22 03:08 chengjunlong0524

Any idea?

chengjunlong0524 avatar Aug 19 '22 04:08 chengjunlong0524

Hello,

I think we've got a bug here. Thank you very much for spot it.

I modified files gyro_gen.m and acc_gen.m in branch develop.

Please, check them and let me know your comments.

rodralez avatar Aug 22 '22 21:08 rodralez

Thanks for your reply. In noise_rrw.m, it probably should be

rrw_n (j, i) = rrw_n(j-1, i) + rrw(i) * sqrt(dt) .* b_noise(j-1);

See that the unit of rrw may be rad/s/root-Hz, rad/s^2/root-Hz, m/s^2/root-Hz or m/s^3/root-Hz. In another form, it is rad/sqrt(s), rad/s/sqrt(s), m/s/sqrt(s) or m/s^2/sqrt(s).

In addition, angle rate random walks(rad/s^2/root-Hz or rad/s/sqrt(s)) seems also reflects bias instability. Are gb_dyn and gyro_rrw repeated? This question also applys to acc_gen.

chengjunlong0524 avatar Aug 23 '22 03:08 chengjunlong0524

邮件已收到!!!

chengjunlong0524 avatar Oct 23 '22 15:10 chengjunlong0524