车辆横向MPC控制中的前馈设计

该小节旨在通过设计控制律中的前馈项,解决上一节中MPC控制器无法达到零稳态横向误差的问题.

前馈信号基本形式

已知,从第一节得到的车辆横向动力学方程为
x˙=Ax+Bδ+Bcψ˙des\dot{x} = Ax + B\delta + B_c\dot{\psi}_{des} x˙=Ax+Bδ+Bc​ψ˙​des​
假设,系统已有一状态反馈控制器,其控制率为
δ=−Kxx\delta = -K_xx δ=−Kx​x
将以上控制律代入车辆横向动力学模型,则有
x˙=(A−BKx)x+Bcψ˙des\dot{x} = (A-BK_x)x + B_c\dot{\psi}_{des} x˙=(A−BKx​)x+Bc​ψ˙​des​
对闭环系统方程进行分析, (A−BKx)≺0(A-BK_x) \prec 0(A−BKx​)≺0 ,系统稳定,状态量 xxx 收敛.但是因 Bcψ˙desB_c\dot{\psi}_{des}Bc​ψ˙​des​ 项的存在,状态量 xxx 并不一定收敛于 0 ,即 xss=lim⁡t→∞x(t)≠0x_{ss} = \lim_{t \to \infty} x(t) \neq 0xss​=limt→∞​x(t)​=0

那么,我们是否可以通过在控制律中加入一个前馈项,使稳态状态量 xxx 收敛于 0 呢? 因此, 我们尝试把控制律改写为
δ=−Kxx+δff\delta = -K_xx + \delta_{ff} δ=−Kx​x+δff​
代入原动力学模型,可得新的闭环系统方程为
x˙=(A−BKx)X+Bδff+Bcψ˙des\dot{x} = (A-BK_x)X + B\delta_{ff} + B_c\dot{\psi}_{des} x˙=(A−BKx​)X+Bδff​+Bc​ψ˙​des​

接下来,对闭环系统稳态状态进行分析.为了简化问题, 我们针对稳态作出以下假设
1.车辆重心在 FLU坐标下的 x方向(车前进方向)速度保持恒定,即 Vx:=Constant2.车辆转弯半径 R保持不变,即 R:=Constant3.由(1),(2)可得, 车辆重心力量横向角加速度保持恒定3)ψ˙des=Vx2R:=Constant4.前馈项保持恒定,即 δff:=Constant5.初始状态量 x0=0\begin{aligned} &\text{1.车辆重心在 $FLU$ 坐标下的 $x$ 方向(车前进方向)速度保持恒定,即 $V_x := Constant$}\\ &\text{2.车辆转弯半径 $R$ 保持不变,即 $R := Constant$} \\ &\text{3.由(1),(2)可得, 车辆重心力量横向角加速度保持恒定} \\ &\qquad 3) \ \dot{\psi}_{des} = \frac{V_x^2}{R} := Constant \\ &\text{4.前馈项保持恒定,即 $\delta_{ff} := Constant$}\\ &\text{5.初始状态量 $x_0 = 0$} \end{aligned} ​1.车辆重心在 FLU 坐标下的 x 方向(车前进方向)速度保持恒定,即 Vx​:=Constant2.车辆转弯半径 R 保持不变,即 R:=Constant3.由(1),(2)可得, 车辆重心力量横向角加速度保持恒定3) ψ˙​des​=RVx2​​:=Constant4.前馈项保持恒定,即 δff​:=Constant5.初始状态量 x0​=0​

基于上述假设, 对闭环系统方程进行拉普拉斯变换,可得
sX(s)=(A−BKx)X(s)+Bδff1s+Bcψ˙des1s\begin{aligned} sX(s) &= (A-BK_x)X(s) + B\delta_{ff}\frac{1}{s} + B_c\dot{\psi}_{des}\frac{1}{s}\\ \end{aligned} sX(s)​=(A−BKx​)X(s)+Bδff​s1​+Bc​ψ˙​des​s1​​
进一步的
X(s)=(sI−A+BKx)−1[Bδff+Bcψ˙des]1sX(s) = (sI - A + BK_x)^{-1}[B\delta_{ff} + B_c\dot{\psi}_{des}]\frac{1}{s}\\ X(s)=(sI−A+BKx​)−1[Bδff​+Bc​ψ˙​des​]s1​
由终值定理可得
xss=lim⁡t→∞x(t)=lim⁡s→0sX(s)=(−A+BKx)−1[Bδff+Bcψ˙des]\begin{aligned} x_{ss} &= \lim_{t \to \infty} x(t) = \lim_{s \to 0}sX(s) \\ &= (- A + BK_x)^{-1}[B\delta_{ff} + B_c\dot{\psi}_{des}] \end{aligned} xss​​=t→∞lim​x(t)=s→0lim​sX(s)=(−A+BKx​)−1[Bδff​+Bc​ψ˙​des​]​
令稳态误差为 0 ,反推前馈项 δff\delta_{ff}δff​ 的形式, 则有
xss=(−A+BKx)−1[Bδff+Bcψ˙des]=0\begin{aligned} x_{ss} = (- A + BK_x)^{-1}[B\delta_{ff} + B_c\dot{\psi}_{des}] = 0 \end{aligned} xss​=(−A+BKx​)−1[Bδff​+Bc​ψ˙​des​]=0​
将在第一节中得到的横向动力学模型中的系数矩阵代入上式
[ess1e˙ss1ess2e˙ss2]={[0−10002(Cαf+Cαr)mVx−2(Cαf+Cαr)m2(Cαflf−Cαrlr)mVx000−102(Cαflf−Cαrlr)IzVx−2(Cαflf−Cαrlr)Iz2(Cαflf2+Cαrlr2)IzVx]+[02Cαfm02CαflfIz][kx1kx2kx3kx4]}−1{[02Cαfm02CαflfIz]δff+[02(−Cαflf+Cαrlr)mVx−Vx0−2(Cαflf2+Cαrlr2)IzVx]ψ˙des}\begin{bmatrix} e_{ss1} \\ \dot{e}_{ss1} \\ e_{ss2} \\ \dot{e}_{ss2} \\ \end{bmatrix}= \left\{ \begin{bmatrix} 0 & -1 & 0 & 0 \\ 0 & \frac{2(C_{\alpha_f}+C_{\alpha_r})}{mV_x} & \frac{-2(C_{\alpha_f}+C_{\alpha_r})}{m} & \frac{2(C_{\alpha_f}l_f-C_{\alpha_r}l_r)}{mV_x}\\ 0 & 0 & 0 & -1 \\ 0 & \frac{2(C_{\alpha_f}l_f-C_{\alpha_r}l_r)}{I_zV_x} & \frac{-2(C_{\alpha_f}l_f-C_{\alpha_r}l_r)}{I_z} & \frac{2(C_{\alpha_f}l^2_f+C_{\alpha_r}l^2_r)}{I_zV_x} \\ \end{bmatrix} + \begin{bmatrix} 0 \\ \frac{2C_{\alpha_f}}{m} \\ 0 \\ \frac{2C_{\alpha_f}l_f}{I_z} \\ \end{bmatrix} \begin{bmatrix} k_{x1} & k_{x2} & k_{x3} & k_{x4} \end{bmatrix} \right\} ^ {-1} \left \{ \begin{bmatrix} 0 \\ \frac{2C_{\alpha_f}}{m} \\ 0 \\ \frac{2C_{\alpha_f}l_f}{I_z} \\ \end{bmatrix}\delta_{ff}+ \begin{bmatrix} 0 \\ \frac{2(-C_{\alpha_f}l_f+C_{\alpha_r}l_r)}{mV_x}-V_x \\ 0 \\ \frac{-2(C_{\alpha_f}l^2_f+C_{\alpha_r}l^2_r)}{I_zV_x} \\ \end{bmatrix} \dot{\psi}_{des} \right \} ⎣⎢⎢⎡​ess1​e˙ss1​ess2​e˙ss2​​⎦⎥⎥⎤​=⎩⎪⎪⎪⎨⎪⎪⎪⎧​⎣⎢⎢⎢⎡​0000​−1mVx​2(Cαf​​+Cαr​​)​0Iz​Vx​2(Cαf​​lf​−Cαr​​lr​)​​0m−2(Cαf​​+Cαr​​)​0Iz​−2(Cαf​​lf​−Cαr​​lr​)​​0mVx​2(Cαf​​lf​−Cαr​​lr​)​−1Iz​Vx​2(Cαf​​lf2​+Cαr​​lr2​)​​⎦⎥⎥⎥⎤​+⎣⎢⎢⎢⎡​0m2Cαf​​​0Iz​2Cαf​​lf​​​⎦⎥⎥⎥⎤​[kx1​​kx2​​kx3​​kx4​​]⎭⎪⎪⎪⎬⎪⎪⎪⎫​−1⎩⎪⎪⎪⎨⎪⎪⎪⎧​⎣⎢⎢⎢⎡​0m2Cαf​​​0Iz​2Cαf​​lf​​​⎦⎥⎥⎥⎤​δff​+⎣⎢⎢⎢⎡​0mVx​2(−Cαf​​lf​+Cαr​​lr​)​−Vx​0Iz​Vx​−2(Cαf​​lf2​+Cαr​​lr2​)​​⎦⎥⎥⎥⎤​ψ˙​des​⎭⎪⎪⎪⎬⎪⎪⎪⎫​
经过一系列的矩阵运算可得
[ess1e˙ss1ess2e˙ss2]=[δffkx1−mv2kx1R(lf+lr)(lr2Cαf+(kx3−1)lf2Cαr)−1kx1R(lf+(1−kx3)lr)0−lrR+mv2RLlf2Cαr0]\begin{bmatrix} e_{ss1} \\ \dot{e}_{ss1} \\ e_{ss2} \\ \dot{e}_{ss2} \\ \end{bmatrix}= \begin{bmatrix} \frac{\delta_{ff}}{k_{x1}} - \frac{mv^2}{k_{x1}R(l_f+l_r)}(\frac{l_r}{2C_{\alpha_f}}+(k_{x3}-1)\frac{l_f}{2C_{\alpha_r}})-\frac{1}{k_{x1}R}(l_f + (1-k_{x3})l_r) \\ 0 \\ -\frac{l_r}{R} + \frac{mv^2}{RL}\frac{l_f}{2C_{\alpha_r}} \\ 0 \end{bmatrix} ⎣⎢⎢⎡​ess1​e˙ss1​ess2​e˙ss2​​⎦⎥⎥⎤​=⎣⎢⎢⎢⎡​kx1​δff​​−kx1​R(lf​+lr​)mv2​(2Cαf​​lr​​+(kx3​−1)2Cαr​​lf​​)−kx1​R1​(lf​+(1−kx3​)lr​)0−Rlr​​+RLmv2​2Cαr​​lf​​0​⎦⎥⎥⎥⎤​
观察可知,合理的前馈项 δff\delta_{ff}δff​ ,可使稳态横向误差 ess1e_{ss1}ess1​ 为 0.但前馈项无法影响到稳态航向误差 ess2e_{ss2}ess2​.

为了达到零稳态横向误差,可得
δff=mv2R(lf+lr)(lr2Cαf+(kx3−1)lf2Cαr)+1R(lf+(1−kx3)lr)=mv2RL(lr2Cαf+(kx3−1)lf2Cαr)+LR−lrRkx3=mv2RL(lr2Cαf−lf2Cαr)+kx3(−lrR+mv2RLlf2Cαr)+LR=mv2RL(lr2Cαf−lf2Cαr)+kx3ess2+LR=v2R(mlr2CαfL−mlf2CαrL)+kx3ess2+LR\begin{aligned} \delta_{ff} &= \frac{mv^2}{R(l_f+l_r)}(\frac{l_r}{2C_{\alpha_f}}+(k_{x3}-1)\frac{l_f}{2C_{\alpha_r}})+\frac{1}{R}(l_f + (1-k_{x3})l_r)\\ &= \frac{mv^2}{RL}(\frac{l_r}{2C_{\alpha_f}}+(k_{x3}-1)\frac{l_f}{2C_{\alpha_r}})+\frac{L}{R} - \frac{l_r}{R}k_{x3}\\ &= \frac{mv^2}{RL}(\frac{l_r}{2C_{\alpha_f}}-\frac{l_f}{2C_{\alpha_r}})+k_{x3}(-\frac{l_r}{R} + \frac{mv^2}{RL}\frac{l_f}{2C_{\alpha_r}})+\frac{L}{R}\\ &= \frac{mv^2}{RL}(\frac{l_r}{2C_{\alpha_f}}-\frac{l_f}{2C_{\alpha_r}})+k_{x3}e_{ss2}+\frac{L}{R}\\ &= \frac{v^2}{R}(\frac{ml_r}{2C_{\alpha_f}L}-\frac{ml_f}{2C_{\alpha_r}L})+k_{x3}e_{ss2}+\frac{L}{R}\\ \end{aligned} δff​​=R(lf​+lr​)mv2​(2Cαf​​lr​​+(kx3​−1)2Cαr​​lf​​)+R1​(lf​+(1−kx3​)lr​)=RLmv2​(2Cαf​​lr​​+(kx3​−1)2Cαr​​lf​​)+RL​−Rlr​​kx3​=RLmv2​(2Cαf​​lr​​−2Cαr​​lf​​)+kx3​(−Rlr​​+RLmv2​2Cαr​​lf​​)+RL​=RLmv2​(2Cαf​​lr​​−2Cαr​​lf​​)+kx3​ess2​+RL​=Rv2​(2Cαf​​Lmlr​​−2Cαr​​Lmlf​​)+kx3​ess2​+RL​​
定义
kv=mlr2CαfL−mlf2CαrLk_v = \frac{ml_r}{2C_{\alpha_f}L}-\frac{ml_f}{2C_{\alpha_r}L} kv​=2Cαf​​Lmlr​​−2Cαr​​Lmlf​​
则前馈项可表示为
δff=LR+kvv2R+kx3ess2\delta_{ff} = \frac{L}{R} + k_v \frac{v^2}{R} + k_{x3}e_{ss2} δff​=RL​+kv​Rv2​+kx3​ess2​
因此,理想稳态控制量 δss\delta_{ss}δss​ 为
δss=−Kxxss+δff=−kx3ess2+δff=LR+kvv2R\begin{aligned} \delta_{ss} &= -K_xx_{ss} + \delta_{ff} \\ &= -k_{x3}e_{ss2} + \delta_{ff} \\ &= \frac{L}{R} + k_v \frac{v^2}{R} \end{aligned} δss​​=−Kx​xss​+δff​=−kx3​ess2​+δff​=RL​+kv​Rv2​​
从这个角度上去理解,前馈项中最后一项 kx3ess2k_{x3}e_{ss2}kx3​ess2​ 的作用就是抵消稳态状态反馈项.

求解线性MPC稳态输出项 −Kxxss-K_xx_{ss}−Kx​xss​

前馈推导式中仍有一未知项,即稳态航向角误差增益 kx3k_{x3}kx3​.接下来,我们针对MPC问题,对其进行求解.开始之前,我们作出以下几点前提假设
1.稳态状态反馈项 −Kxxss:=−kx3ess2, 来源于上一节设计的 MPC 控制器的稳态输出2.建立 MPC 问题时, 考虑动力学方程中 Bcψ˙des项的影响3.接近稳态时,可忽略 MPC 的约束条件. 即将线性MPC简化为一个有限时域LQR问题4.在有限时域 LQR 的最终时刻 t=N时,系统达到稳态\begin{aligned} &\text{1.稳态状态反馈项 $-K_xx_{ss} := -k_{x3}e_{ss2}$, 来源于上一节设计的 MPC 控制器的稳态输出}\\ &\text{2.建立 MPC 问题时, 考虑动力学方程中 $B_c\dot{\psi}_{des}$ 项的影响} \\ &\text{3.接近稳态时,可忽略 MPC 的约束条件. 即将线性MPC简化为一个有限时域LQR问题} \\ &\text{4.在有限时域 LQR 的最终时刻 $t = N$ 时,系统达到稳态} \\ \end{aligned} ​1.稳态状态反馈项 −Kx​xss​:=−kx3​ess2​, 来源于上一节设计的 MPC 控制器的稳态输出2.建立 MPC 问题时, 考虑动力学方程中 Bc​ψ˙​des​ 项的影响3.接近稳态时,可忽略 MPC 的约束条件. 即将线性MPC简化为一个有限时域LQR问题4.在有限时域 LQR 的最终时刻 t=N 时,系统达到稳态​

首先,定义有限时域 LQR 的价值函数为 Vt(xt)V_t(x_t)Vt​(xt​) ,其中 t=0,...,Nt = 0,...,Nt=0,...,N, 因此
Vt(xt)=xNTQfxN+min⁡δt,...,δN−1∑τ=tN−1(xτTQxτ+δτTRδτ)V_t(x_t) = x_N^TQ_fx_N+\min_{\delta_t,...,\delta_{N-1}}\sum_{\tau=t}^{N-1}(x_{\tau}^TQx_{\tau}+\delta_{\tau}^TR\delta_{\tau}) Vt​(xt​)=xNT​Qf​xN​+δt​,...,δN−1​min​τ=t∑N−1​(xτT​Qxτ​+δτT​Rδτ​)
其中
1.Vt给出的是 t时刻, 有限时域LQR的最小剩余代价.2.Vt可被构造为二次型,即 Vt(xt)=xtTPtxt,Pt=PtT⪰03.VN(xN)=xNTQfxN,PN=Qf\begin{aligned} &\text{1.$V_t$给出的是 $t$ 时刻, 有限时域LQR的最小剩余代价.}\\ &\text{2.$V_t$可被构造为二次型,即 $V_t(x_t) = x_t^TP_tx_t, \ P_t = P_t^T \succeq 0$ } \\ &\text{3.$V_N(x_N) = x_N^TQ_fx_N, \ P_N = Q_f$} \end{aligned} ​1.Vt​给出的是 t 时刻, 有限时域LQR的最小剩余代价.2.Vt​可被构造为二次型,即 Vt​(xt​)=xtT​Pt​xt​, Pt​=PtT​⪰0 3.VN​(xN​)=xNT​Qf​xN​, PN​=Qf​​
下一步,我们需要问题的假设,明确MPC输出 δt\delta_tδt​ 的结构形式.假设,我们已知 Vt+1(xt+1)V_{t+1}(x_{t+1})Vt+1​(xt+1​), 需求 ttt 时刻的最优控制量 δt\delta_tδt​.因此
Vt(xt)=min⁡δt(xtTQxt+δtTRδt+Vt+1(xt+1))=xtTQxt+min⁡δt(δtTRδt+Vt+1(Adxt+Bdδt+Bcdψ˙des))=xtTQxt+min⁡δt(δtTRδt+(Adxt+Bdδt+Bcdψ˙des)TPt+1(Adxt+Bdδt+Bcdψ˙des))\begin{aligned} V_t(x_t) &= \min_{\delta_t}(x_t^TQx_t+\delta_t^TR\delta_t+V_{t+1}(x_{t+1}))\\ &= x_t^TQx_t+\min_{\delta_t}(\delta_t^TR\delta_t+V_{t+1}(A_dx_t + B_d\delta_t + B_{cd}\dot{\psi}_{des}))\\ &= x_t^TQx_t+\min_{\delta_t}(\delta_t^TR\delta_t+(A_dx_t + B_d\delta_t + B_{cd}\dot{\psi}_{des})^TP_{t+1}(A_dx_t + B_d\delta_t + B_{cd}\dot{\psi}_{des})) \end{aligned} Vt​(xt​)​=δt​min​(xtT​Qxt​+δtT​Rδt​+Vt+1​(xt+1​))=xtT​Qxt​+δt​min​(δtT​Rδt​+Vt+1​(Ad​xt​+Bd​δt​+Bcd​ψ˙​des​))=xtT​Qxt​+δt​min​(δtT​Rδt​+(Ad​xt​+Bd​δt​+Bcd​ψ˙​des​)TPt+1​(Ad​xt​+Bd​δt​+Bcd​ψ˙​des​))​
对上式求导可得
∂Vt∂t=2Rδt∗+2BdTPt+1(Adxt+Bdδt∗+Bcdψ˙des)=0\frac{\partial V_t}{\partial t} = 2R\delta_t^*+2B_d^TP_{t+1}(A_dx_t+B_d\delta_t^*+B_{cd}\dot{\psi}_{des}) = 0 ∂t∂Vt​​=2Rδt∗​+2BdT​Pt+1​(Ad​xt​+Bd​δt∗​+Bcd​ψ˙​des​)=0
由此可得线性MPC输出项形式为
δt∗=−(R+BdTPt+1Bd)−1BdTPt+1(Adxt+Bcdψ˙des)=−(R+BdTPt+1Bd)−1BdTPt+1Adxt−(R+BdTPt+1Bd)−1BdTPt+1Bcdψ˙des=−Klqrxt−(R+BdTPt+1Bd)−1BdTPt+1Bcdψ˙des=−Klqrxt+ϕcomp\begin{aligned} \delta_t^* &= -(R+B_d^TP_{t+1}B_d)^{-1}B_d^TP_{t+1}(A_dx_t+B_{cd}\dot{\psi}_{des}) \\ &= -(R+B_d^TP_{t+1}B_d)^{-1}B_d^TP_{t+1}A_dx_t - (R+B_d^TP_{t+1}B_d)^{-1}B_d^TP_{t+1}B_{cd}\dot{\psi}_{des} \\ &= - K_{lqr} x_t - (R+B_d^TP_{t+1}B_d)^{-1}B_d^TP_{t+1}B_{cd}\dot{\psi}_{des} \\ &= - K_{lqr} x_t + \phi_{comp} \end{aligned} δt∗​​=−(R+BdT​Pt+1​Bd​)−1BdT​Pt+1​(Ad​xt​+Bcd​ψ˙​des​)=−(R+BdT​Pt+1​Bd​)−1BdT​Pt+1​Ad​xt​−(R+BdT​Pt+1​Bd​)−1BdT​Pt+1​Bcd​ψ˙​des​=−Klqr​xt​−(R+BdT​Pt+1​Bd​)−1BdT​Pt+1​Bcd​ψ˙​des​=−Klqr​xt​+ϕcomp​​
进一步的,根据前提假设4与价值函数性质3,可得MPC稳态输出为
δssmpc=δN−1=−(R+BdTPNBd)−1BdTPN(Adxss+Bcdψ˙des)=−(R+BdTQfBd)−1BdTQfAdxss−(R+BdTQf(Bd)−1BdTQfBcdψ˙des=−Klqrssxss+ϕcompss=−klqrss3ess2+ϕcompss\begin{aligned} \delta_{ss_{mpc}} = \delta_{N-1} &= -(R+B_d^TP_{N}B_d)^{-1}B_d^TP_{N}(A_dx_{ss}+B_{cd}\dot{\psi}_{des}) \\ &= -(R+B_d^TQ_fB_d)^{-1}B_d^TQ_fA_dx_{ss} - (R+B_d^TQ_f(B_d)^{-1}B_d^TQ_fB_{cd}\dot{\psi}_{des} \\ &= - K_{lqr_{ss}} x_{ss} + \phi_{comp_{ss}} \\ &= - k_{lqr_{ss3}} e_{ss2} + \phi_{comp_{ss}} \end{aligned} δssmpc​​=δN−1​​=−(R+BdT​PN​Bd​)−1BdT​PN​(Ad​xss​+Bcd​ψ˙​des​)=−(R+BdT​Qf​Bd​)−1BdT​Qf​Ad​xss​−(R+BdT​Qf​(Bd​)−1BdT​Qf​Bcd​ψ˙​des​=−Klqrss​​xss​+ϕcompss​​=−klqrss3​​ess2​+ϕcompss​​​
为了在稳态时,得到零横向误差,需使
δss=δssmpc+δff=LR+kvv2R\delta_{ss} = \delta_{ss_{mpc}} + \delta_{ff} = \frac{L}{R} + k_v \frac{v^2}{R} δss​=δssmpc​​+δff​=RL​+kv​Rv2​
因此,最终的前馈项为
δff=LR+kvv2R+klqrss3ess2−ϕcompss\delta_{ff} = \frac{L}{R} + k_v \frac{v^2}{R} + k_{lqr_{ss3}} e_{ss2} - \phi_{comp_{ss}} δff​=RL​+kv​Rv2​+klqrss3​​ess2​−ϕcompss​​

模型预测控制器(MPC)系列: 3.车辆横向MPC控制中的前馈设计相关推荐

  1. 黄大叔学MPC系列之Lyapunov-based Nonlinear MPC

    黄大叔学MPC系列之Lyapunov-based Nonlinear MPC Lyapunov-based Nonlinear MPC 仿真对象 主函数 模型函数 目标函数 Lyapunov 约束函数 ...

  2. 模型预测控制器(MPC)系列: 2.求解车辆横向控制中的MPC

    求解MPC: 在滚动时间窗内建立并求解QP问题 该小节旨在根据上一节得到的离散车辆横向动力学模型,在滚动时间窗内建立并求解QP问题,实现MPC横向控制. 在开始求解之前,我们需要回答以下两个问题: 1 ...

  3. 模型预测控制器(MPC)系列: 1.建立车辆横向动力学模型

    勘误 Update 02/23/2021 之前的文章中有不严谨的地方,这里做一个勘误.错误就在下面描述坐标系的图中.<更正后的图已覆盖到坐标系小节下> 在这个图中,我指出ENU坐标下,车自 ...

  4. 模型预测控制算法基础与车辆纵向控制仿真分析

    模型预测控制算法基础与车辆纵向控制仿真分析 第三章 模型预测控制算法基础与控制仿真分析 模型预测控制算法基础 模型预测控制的基本思想就是利用已有的模型.系统当前的状态和未来的控制量去预测系统未来的输出 ...

  5. 无人驾驶模型预测控制(五)-基于动力学模型的无人驾驶车辆主动转向控制

    理论基础 各种概念 非线性系统:状态变量之间有相乘,无法表达为x点=Ax+Bu形式 根据是否有时间t作为参数,还分为非线性时变和非线性定常系统. 线性系统,可以表达为线性形式. 线性时变:参数随时间变 ...

  6. 预测电流FCS-MPC模型预测电流控制

    预测电流FCS-MPC模型预测电流控制 级联(链式):5个单元 下载地址[https://download.csdn.net/download/a_zxswer/20009112](https://d ...

  7. 滚动时域控制 matlab,在 Simulink 中设计神经网络预测控制器

    在 Simulink 中设计神经网络预测控制器 在 Deep Learning Toolbox™ 软件中实现的神经网络预测控制器使用非线性被控对象的神经网络模型来预测被控对象将来的性能.然后,控制器计 ...

  8. 用Python搭建机器学习模型预测房租价格

    毫无疑问,机器学习是当前大数据分析中最热门的话题.这也是一些最令人兴奋的技术领域的基本概念,例如自动驾驶汽车和预测分析.百度上的机器学习搜索在2019年4月创历史新高,自此以来兴趣一直没有下降. 但是 ...

  9. R语言ARIMA集成模型预测时间序列分析

    全文链接:http://tecdat.cn/?p=18493 本文我们使用4个时间序列模型对每周的温度序列建模.第一个是通过auto.arima获得的,然后两个是SARIMA模型,最后一个是Buys- ...

  10. python用ARIMA模型预测CO2浓度时间序列实现

    全文下载链接:http://tecdat.cn/?p=20424 时间序列为预测未来数据提供了方法.根据先前的值,时间序列可用于预测经济,天气的趋势.时间序列数据的特定属性意味着通常需要专门的统计方法 ...

最新文章

  1. 【运维】使用FC命令辅助查杀DLL木马
  2. java三大范_Java深度学习系列——数据库的三大范式
  3. FFmpeg使用遇到问题记录
  4. 9012教你如何使用gulp4开发项目脚手架
  5. [转]解决IE下CSS背景图片闪烁的Bug
  6. 雷军:小米推动山寨机灭亡;苹果回应禁止第三方换电池;Javascript 引擎 V8 7.7 发布 | 极客头条...
  7. 我的Android进阶之旅------百度地图学习:BDLocation.getLocType ( )值分析
  8. oracle 11g安装自己理解
  9. Servlet 生命周期、工作原理(转)
  10. 【机器学习系列】MCMC第一讲:蒙特卡罗方法初认识
  11. 腾讯大讲堂 微信红包系统设计 优化
  12. python编程实战:暴力破解WIFI密码!亲测运行有效!
  13. 加密WebSphere Application Server系统密码
  14. Android 10.0 关机界面全屏显示(UI全屏显示)
  15. 你又知不知道,日有所思,夜有所梦。
  16. QQ音乐关键字搜索并生成下载url
  17. python类中的魔方方法
  18. 九秀直播助手下载|九秀直播助手下载
  19. MongoDB多层嵌套查询
  20. GAU : Transformer Quality in Linear Time(new attention+)

热门文章

  1. 微信朋友圈html5广告,微信朋友圈原生广告页是什么
  2. YUV422_UYVY图像格式转RGB565
  3. 数据线CE测试标准 准备资料
  4. 小米路由器青春版(R1CL)刷高恪软路由系统
  5. ADB连接手机和简单命令
  6. 6.2 数据通路的建立
  7. NRF52840 NRF52833 nRF5 SDK 开发
  8. metabase开源BI
  9. FreeBSD和JunOS安装手册
  10. win10下网易云网络异常问题解决办法