12345qiupeng
/
gtsam
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Make two acceleration scenarios

release/4.3a0
Frank Dellaert 2015-12-23 15:04:36 -08:00
parent 31335608a8
commit f79a9b8d3a
1 changed files with 30 additions and 7 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

27
gtsam/navigation/tests/testScenarioRunner.cpp
View File

@ -76,13 +76,17 @@ TEST(ScenarioRunner, Loop) {
} }
/* ************************************************************************* */ /* ************************************************************************* */
TEST(ScenarioRunner, Accelerating) { namespace initial {
// Set up body pointing towards y axis, and start at 10,20,0 with velocity // Set up body pointing towards y axis, and start at 10,20,0 with velocity
// going in X. The body itself has Z axis pointing down // going in X. The body itself has Z axis pointing down
const Rot3 nRb(Point3(0, 1, 0), Point3(1, 0, 0), Point3(0, 0, -1)); const Rot3 nRb(Point3(0, 1, 0), Point3(1, 0, 0), Point3(0, 0, -1));
const Point3 P0(10, 20, 0); const Point3 P0(10, 20, 0);
const Vector3 V0(50, 0, 0); const Vector3 V0(50, 0, 0);
}
/* ************************************************************************* */
TEST(ScenarioRunner, Accelerating) {
using namespace initial;
const double a = 0.2; // m/s^2 const double a = 0.2; // m/s^2
const Vector3 A(0, a, 0); const Vector3 A(0, a, 0);
const AcceleratingScenario scenario(nRb, P0, V0, A); const AcceleratingScenario scenario(nRb, P0, V0, A);
@ -93,7 +97,26 @@ TEST(ScenarioRunner, Accelerating) {
ImuFactor::PreintegratedMeasurements pim = runner.integrate(T); ImuFactor::PreintegratedMeasurements pim = runner.integrate(T);
EXPECT(assert_equal(scenario.pose(T), runner.predict(pim).pose, 1e-9)); EXPECT(assert_equal(scenario.pose(T), runner.predict(pim).pose, 1e-9));
Matrix6 estimatedCov = runner.estimatePoseCovariance(T); Matrix6 estimatedCov = runner.estimatePoseCovariance(T,10000);
EXPECT(assert_equal(estimatedCov, runner.poseCovariance(pim), 0.1));
cout << estimatedCov << endl << endl;
cout << runner.poseCovariance(pim) << endl;
}
/* ************************************************************************* */
TEST(ScenarioRunner, AcceleratingAndRotating) {
using namespace initial;
const double a = 0.2; // m/s^2
const Vector3 A(0, a, 0), W(0, 0.1, 0);
const AcceleratingScenario scenario(nRb, P0, V0, A, W);
const double T = 3; // seconds
ScenarioRunner runner(&scenario, T / 10, kGyroSigma, kAccelerometerSigma);
ImuFactor::PreintegratedMeasurements pim = runner.integrate(T);
EXPECT(assert_equal(scenario.pose(T), runner.predict(pim).pose, 1e-9));
Matrix6 estimatedCov = runner.estimatePoseCovariance(T,10000);
EXPECT(assert_equal(estimatedCov, runner.poseCovariance(pim), 0.1)); EXPECT(assert_equal(estimatedCov, runner.poseCovariance(pim), 0.1));
} }