92 lines
2.4 KiB
C++
92 lines
2.4 KiB
C++
/**
|
|
* @file Mechanization_bRn2.h
|
|
* @date Jan 25, 2012
|
|
* @author Chris Beall
|
|
* @author Frank Dellaert
|
|
*/
|
|
|
|
#pragma once
|
|
|
|
#include <gtsam/geometry/Rot3.h>
|
|
#include <gtsam/base/Vector.h>
|
|
#include <gtsam_unstable/dllexport.h>
|
|
#include <list>
|
|
|
|
namespace gtsam {
|
|
|
|
class GTSAM_UNSTABLE_EXPORT Mechanization_bRn2 {
|
|
|
|
private:
|
|
Rot3 bRn_; ///< rotation from nav to body
|
|
Vector3 x_g_; ///< gyroscope bias
|
|
Vector3 x_a_; ///< accelerometer bias
|
|
|
|
public:
|
|
|
|
/**
|
|
* Constructor
|
|
* @param initial_bRn initial rotation from nav to body frame
|
|
* @param initial_x_g initial gyro bias
|
|
* @param r3 Z-axis of rotated frame
|
|
*/
|
|
Mechanization_bRn2(const Rot3& initial_bRn = Rot3(),
|
|
const Vector3& initial_x_g = Z_3x1, const Vector3& initial_x_a = Z_3x1) :
|
|
bRn_(initial_bRn), x_g_(initial_x_g), x_a_(initial_x_a) {
|
|
}
|
|
|
|
/// Copy constructor
|
|
Mechanization_bRn2(const Mechanization_bRn2& other) :
|
|
bRn_(other.bRn_), x_g_(other.x_g_), x_a_(other.x_a_) {
|
|
}
|
|
|
|
/// gravity in the body frame
|
|
Vector3 b_g(double g_e) const {
|
|
Vector3 n_g(0, 0, g_e);
|
|
return bRn_ * n_g;
|
|
}
|
|
|
|
const Rot3& bRn() const {return bRn_; }
|
|
const Vector3& x_g() const { return x_g_; }
|
|
const Vector3& x_a() const { return x_a_; }
|
|
|
|
/**
|
|
* Initialize the first Mechanization state
|
|
* @param U a list of gyro measurement vectors
|
|
* @param F a list of accelerometer measurement vectors
|
|
* @param g_e gravity magnitude
|
|
* @param flat flag saying whether this is a flat trim init
|
|
*/
|
|
static Mechanization_bRn2 initializeVector(const std::list<Vector>& U,
|
|
const std::list<Vector>& F, const double g_e = 0, bool flat=false);
|
|
|
|
/// Matrix version of initialize
|
|
static Mechanization_bRn2 initialize(const Matrix& U,
|
|
const Matrix& F, const double g_e = 0, bool flat=false);
|
|
|
|
/**
|
|
* Correct AHRS full state given error state dx
|
|
* @param obj The current state
|
|
* @param dx The error used to correct and return a new state
|
|
*/
|
|
Mechanization_bRn2 correct(const Vector9& dx) const;
|
|
|
|
/**
|
|
* Integrate to get new state
|
|
* @param obj The current state
|
|
* @param u gyro measurement to integrate
|
|
* @param dt time elapsed since previous state in seconds
|
|
*/
|
|
Mechanization_bRn2 integrate(const Vector3& u, const double dt) const;
|
|
|
|
/// print
|
|
void print(const std::string& s = "") const {
|
|
bRn_.print(s + ".R");
|
|
|
|
std::cout << s + ".x_g" << x_g_ << std::endl;
|
|
std::cout << s + ".x_a" << x_a_ << std::endl;
|
|
}
|
|
|
|
};
|
|
|
|
} // namespace gtsam
|