1 projection introduction

本项目希望大家根据Stanley算法实现车辆的横向控制，结合上一章的PID算法实现轨迹跟踪。基本的系统框架已经给出，仅需要完成stanley_control.cpp中todo部分。

给定的系统框架如下

#include "stanley_control.h"#include <algorithm>
#include <iomanip>
#include <utility>
#include <vector>
#include "Eigen/LU"
#include <math.h>using namespace std;namespace shenlan {
namespace control {double atan2_to_PI(const double atan2) {return atan2 * M_PI / 180;
}double PointDistanceSquare(const TrajectoryPoint &point, const double x,const double y) {const double dx = point.x - x;const double dy = point.y - y;return dx * dx + dy * dy;
}void StanleyController::LoadControlConf() {k_y_ = 0.5;
}// /** to-do **/ 计算需要的控制命令, 实现对应的stanley模型,并将获得的控制命令传递给汽车
// 提示，在该函数中你需要调用计算误差
void StanleyController::ComputeControlCmd(const VehicleState &vehicle_state,const TrajectoryData &planning_published_trajectory, ControlCmd &cmd) {}// /** to-do **/ 计算需要的误差，包括横向误差，纵向误差void StanleyController::ComputeLateralErrors(const double x, const double y,const double theta, double &e_y,double &e_theta) {}TrajectoryPoint StanleyController::QueryNearestPointByPosition(const double x,const double y) {double d_min = PointDistanceSquare(trajectory_points_.front(), x, y);size_t index_min = 0;for (size_t i = 1; i < trajectory_points_.size(); ++i) {double d_temp = PointDistanceSquare(trajectory_points_[i], x, y);if (d_temp < d_min) {d_min = d_temp;index_min = i;}}// cout << " index_min: " << index_min << endl;//cout << "tarjectory.heading: " << trajectory_points_[index_min].heading << endl;theta_ref_ = trajectory_points_[index_min].heading;return trajectory_points_[index_min];
}}  // namespace control
}  // namespace shenlan

2 思路提示

从掌握stanley算法到用C++实现还需要思考不少问题

2.1 ComputeControlCmd

ComputeControlCmd是计算控制命令的函数，因此调用计算误差的函数，得到航向误差与横向误差，通过如下公式即可计算出车辆前轮转角

但是这里还需要考虑车辆转向角有范围，当超过范围则需要取边值

//计算需要的控制命令, 实现对应的stanley模型,并将获得的控制命令传递给汽车
void StanleyController::ComputeControlCmd(const VehicleState &vehicle_state,const TrajectoryData &planning_published_trajectory, ControlCmd &cmd) {
trajectory_points_=planning_published_trajectory.trajectory_points;
double e_y=0.0;
double e_theta=0.0;
cout<< "vehicle_state.heading:" <<vehicle_state.heading <<endl;
ComputeLateralErrors(vehicle_state.x,vehicle_state.y,vehicle_state.heading,e_y,e_theta);
cout<< "e_y:" <<e_y <<endl;
cout<< "e_theta:" <<e_theta <<endl;
cmd.steer_target=e_theta+std::atan2(k_y_*e_y,vehicle_state.velocity);
if (cmd.steer_target>M_PI/4.0){
cmd.steer_target=M_PI/4.0;
}
else if(cmd.steer_target<-M_PI/4.0){
cmd.steer_target=-M_PI/4.0;
}
}

我还把 vehicle_state.heading、e_y、e_theta 输出，便于观察误差是否能快速收敛到0并且稳态误差较小，同时设定最大转角为45度（根据实际车辆可以调整该参数）

2.2 ComputeLateralErrors

这一部分需要计算航向误差与横向误差

航向误差：离车辆前轮中心最近点的切线与车辆纵轴之间的夹角，但是我们需要小心此处的正负

横向误差：离车辆前轮中心最近点与车辆前轮中心之间的距离，同时这里也会涉及到正负

在理解算法时可能没有关注到上述两个量的正负，但是用代码实现时，这是非常重要的问题

横向误差：因为当一辆车在参考轨迹的右侧和左侧，转向角必须做出相反的反应，这里通过向量叉乘判断车在参考轨迹的哪一侧

航向误差：会涉及到的问题就是需要把两角之差限制到-pi~pi，因为原始角度会出现差2*k*pi是代表一样的角度，导致把多出来的2*k*pi附加到转向角中

//计算需要的误差，包括横向误差，纵向误差
void StanleyController::ComputeLateralErrors(const double x, const double y,const double theta, double &e_y,double &e_theta) {
TrajectoryPoint target_point=QueryNearestPointByPosition(x,y);e_y=sqrt((target_point.x-x)*(target_point.x-x)+(target_point.y-y)*(target_point.y-y));
double judge=(target_point.y-y)*cos(target_point.heading)-(target_point.x-x)*sin(target_point.heading);
if (judge>0){
e_y=-e_y;
}
e_theta=theta-target_point.heading;
if(e_theta>M_PI){
e_theta-=2.0*M_PI;
}
else if(e_theta<-M_PI){
e_theta+=2.0*M_PI;
}
return;
}