CN113110045B - 一种基于计算图的模型预测控制实时优化并行计算方法 - Google Patents
一种基于计算图的模型预测控制实时优化并行计算方法 Download PDFInfo
- Publication number
- CN113110045B CN113110045B CN202110344736.0A CN202110344736A CN113110045B CN 113110045 B CN113110045 B CN 113110045B CN 202110344736 A CN202110344736 A CN 202110344736A CN 113110045 B CN113110045 B CN 113110045B
- Authority
- CN
- China
- Prior art keywords
- control
- model
- prediction
- objective function
- parallel computing
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000005457 optimization Methods 0.000 title claims abstract description 37
- 238000004364 calculation method Methods 0.000 title claims abstract description 36
- 238000000034 method Methods 0.000 claims abstract description 26
- 238000011478 gradient descent method Methods 0.000 claims abstract description 11
- 239000011159 matrix material Substances 0.000 claims description 27
- 230000008569 process Effects 0.000 claims description 17
- 238000003860 storage Methods 0.000 claims description 8
- 238000000354 decomposition reaction Methods 0.000 claims description 3
- 238000006073 displacement reaction Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 230000001133 acceleration Effects 0.000 description 2
- 230000008878 coupling Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- 238000009795 derivation Methods 0.000 description 2
- 238000004904 shortening Methods 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 238000011217 control strategy Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000005096 rolling process Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000001960 triggered effect Effects 0.000 description 1
- 239000013598 vector Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
- G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
- G05B13/04—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
- G05B13/042—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
- G06N3/084—Backpropagation, e.g. using gradient descent
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N5/00—Computing arrangements using knowledge-based models
- G06N5/02—Knowledge representation; Symbolic representation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/06—Physical realisation, i.e. hardware implementation of neural networks, neurons or parts of neurons
- G06N3/063—Physical realisation, i.e. hardware implementation of neural networks, neurons or parts of neurons using electronic means
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T10/00—Road transport of goods or passengers
- Y02T10/10—Internal combustion engine [ICE] based vehicles
- Y02T10/40—Engine management systems
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Software Systems (AREA)
- Artificial Intelligence (AREA)
- Evolutionary Computation (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Computational Linguistics (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Physics (AREA)
- Computing Systems (AREA)
- Health & Medical Sciences (AREA)
- Biomedical Technology (AREA)
- Molecular Biology (AREA)
- General Health & Medical Sciences (AREA)
- Biophysics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Medical Informatics (AREA)
- Automation & Control Theory (AREA)
- Feedback Control In General (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及一种基于计算图的模型预测控制实时优化并行计算方法,包括以下步骤:搭建***状态量的预测模型,并构建***的目标函数;建立预测模型和目标函数的模型预测控制并行计算架构,所述的并行计算架构采用触发式并行计算方法同步计算预测模型和目标函数;采用反向传播的方式求解计算梯度,利用梯度下降法优化***的控制量,实现***的实时优化控制。与现有技术相比,本发明大幅提升了运算效率,保证模型预测控制器的实时性,拓展了模型预测控制的应用领域。
Description
技术领域
本发明涉及模型预测控制的实时优化技术领域,尤其是涉及一种基于计算图的模型预测控制实时优化并行计算方法。
背景技术
由于模型预测控制具有滚动优化、反馈校正、显式考虑***约束等特点,因此越来越多的应用领域特别是快速动态***(例如电力电子、机电工程、汽车电子等)迫切需要模型预测控制来处理其复杂约束优化问题和提高控制性能。然而模型预测控制的当前控制动作是在每一个采样时刻通过模型求解一个有限时域开环最优控制问题而获得,涉及很大的计算量和计算时间,因此,模型预测控制的在线优化计算量大是制约其应用的主要瓶颈。
近年来,预测控制的快速计算研究取得了许多有价值的成果。在控制策略方面,通过控制器结构优化设计简化预测控制求解过程,有效降低计算复杂度,但现有方法多数是时域内串行计算策略,提速空间有限;在优化问题求解方面,现有方法大多采用标准或者改进规划算法进行迭代计算求解,但是直接求解非线性优化问题涉及大量复杂的梯度计算和矩阵求逆运算,计算量很大,同时一般优化迭代逻辑较复杂,并且多数属于串行迭代算法,这一特点决定了它们并不具备多大的并行加速空间,只能是尽可能的加快每次迭代的速度,这一点很大程度上要靠处理器的速度提升,不利于硬件的并行实现。
发明内容
本发明的目的就是为了克服上述现有技术存在的缺陷而提供一种基于计算图的模型预测控制实时优化并行计算方法。
本发明的目的可以通过以下技术方案来实现:
一种基于计算图的模型预测控制实时优化并行计算方法,包括以下步骤:
S1:搭建***状态量的预测模型,并构建***的目标函数;
S2:建立预测模型和目标函数的模型预测控制并行计算架构,所述的并行计算架构采用触发式并行计算方法同步计算预测模型和目标函数;
S3:采用反向传播的方式求解计算梯度,利用梯度下降法优化***的控制量,实现***的实时优化控制。
优选地,所述的步骤S2的模型预测控制并行计算架构中,用当前步序中的预测模型和目标函数的求解均完成的标志作为下一步序预测计算开始的标志,实现预测模型、目标函数的并行计算。
优选地,所述的预测模型与目标函数的递推关系为:
其中,J为目标函数,f为***的预测模型,xk+i|k为第k时刻第i步的***状态量,uk+i|k为k时刻第i步的***控制量,
优选地,所述的步骤S4的具体步骤包括:
S41:搭建多个计算节点,每个计算节点设置一个存储单元,所述的存储单元内存储相关计算参数;
S42:根据多个计算节点中的计算参数,基于反向传播获取目标函数对输入量的梯度;
S43:利用梯度下降法优化***的控制量,获取最优控制序列,实现***的并行预测控制。
优选地,所述步骤S43的具体包括:
利用梯度下降法优化控制量:
uk|k,uk+1|k…uk+N-1|k,
其中,uk|k,uk+1|k…uk+N-1|k分别为第k时刻内第0、1…N-1步的控制量,
当满足优化条件其中之一时,完成优化过程,获取最优控制序列Uk *:
重复以上步骤,完成第k+1时刻的预测、优化过程。
优选地,所述的最优控制序列的计算公式为:
优选地,所述的优化条件为当前迭代的步长的目标函数与上一步的目标函数的差值小于设定值或达到限定优化次数或目标函数变化量为0。
优选地,所述的***为车辆模型,所述的预测模型为车辆路径与速度跟踪模型。
优选地,所述的车辆路径与速度跟踪模型的控制目标是:
快速准确跟踪车辆纵向速度vx和横向位移Y,设控制时域与预测时域均为N,目标函数为:
其中,xi为当前时刻第i步的状态量,ui为当前时刻第i步的***控制量,Q为正定加权矩阵,P为正定终端惩罚矩阵,R为正定控制量惩罚矩阵。
优选地,所述的目标函数进行计算时分解为:
Ji=h(yi,Ji-1)=Ji-1+(yi-yr)TQ(yi-yr)
=(xji-xjr)TQ(xji-xjr),j=1,6,i=1,2,…,N
其中,g、h分别为分解计算函数。
与现有技术相比,本发明基于多指令、多数据的并行计算思想提出了采用前向求解与目标函数求解并行计算的架构,并考虑并行任务数据的耦合性,结合触发式并行计算方式,保证数据计算数的顺序性,最终达到缩短计算时间的目的。同时,该架构基于计算图的思想,将一次正向传播求解预测状态的过程作为一个节点,并且该节点包含输入量、输出量以及输入量对输出量的偏导数。通过反向传播的方式,将节点中保存的输入与输出的偏导数依次取出,并将它们相乘可计算梯度。一般的,模型预测控制中的控制序列矩阵维数要高于目标函数矩阵的维数,因此,相对于前向计算,反向计算可以减少了运算次数,进一步大幅提升了运算效率,保证模型预测控制器的实时性,拓展了模型预测控制的应用领域,具有实用性和应用价值。
附图说明
图1为本发明的流程图;
图2为本发明的预测***未来状态递推过程的计算模型图;
图3为本发明的并行计算架构图;
图4为前向传播的计算方式示意图;
图5为反向传播的计算方式示意图;
图6为梯度下降法部分计算图;
图7为参考速度曲线图;
图8为路径跟踪曲线图;
图9为速度跟踪曲线图;
图10为位移跟踪误差曲线图;
图11为速度跟踪误差曲线图。
具体实施方式
下面结合附图和具体实施例对本发明进行详细说明。注意,以下的实施方式的说明只是实质上的例示,本发明并不意在对其适用物或其用途进行限定,且本发明并不限定于以下的实施方式。
实施例
一种基于计算图的模型预测控制实时优化并行计算方法,如图1所示,包括以下步骤:
S1:搭建***状态量的预测模型,并构建***的目标函数。
本实施例中,考虑如下控制***:
***离散动力学方程为xk+1=f(xk,uk);
***输出量等于***状态量,即y=x;
***期望输出为零矩阵,即yr=0;
控制目标是预测误差的平方和与控制序列变化率的平方和同时达到最小;
控制时域与预测时域相等,均为N。
依据以上***,结合图模型思想,可以得到如下预测***未来状态递推过程的计算模型如图2所示,其中xk|k为k时刻的状态量,uk|k为k时刻的输入量,f为递推函数(预测模型),J为待优化的目标函数,P为终端的加权矩阵,Q为状态加权矩阵R为控制量的加权矩阵。
预测模型与目标函数的递推关系为:
其中,J为目标函数,f为***的预测模型,xk+i|k为第k时刻第i步的***状态量,uk+i|k为k时刻第i步的***控制量。
S2:建立预测模型和目标函数的模型预测控制并行计算架构,所述的并行计算架构采用触发式并行计算方法同步计算预测模型和目标函数。步骤S2的模型预测控制并行计算架构中,用当前步序中的预测模型和目标函数的求解均完成的标志作为下一步序预测计算开始的标志,实现预测模型、目标函数的并行计算。
依据以上前向计算过程,编写计算程序。搭建的每一节点对应单片机的一个存储单元,包含以下5个部分,以xk+N-1|k处的节点为例,当完成前向计算时,每一个节点的输出以及该节点对应的偏导数都存储在对应的存储空间内。
xk+N-1|k
uk+N-1|k
xk+N|k=f(xk+N-1|k,uk+N-1|k)
在进行前向递推的过程中,目标函数的计算将同步进行,如图3所示。由于并行任务中数据不是完全独立,而是存在耦合。但是针对每一步的预测过程,数据是是独立的,为此本发明结合触发式并行计算方式,即用第N步中的前向推导(即并行计算内容一)和目标函数的求解(即并行计算内容二)均完成的标志作为第N+1步预测计算开始的标志,以保证数据计算数的顺序性,最终达到缩短计算时间的目的。
S3:采用反向传播的方式求解计算梯度,利用梯度下降法优化***的控制量,实现***的实时优化控制。
通过反向传播计算梯度,即把所需节点对应存储的偏导数依次取出并相乘,获得目标函数对输入量的梯度。由于整个控制时域内存在N个控制量,1个输出变量的控制模型,正向求解需要遍历全部节点N次才能计算出目标函数对各输入的偏导数,而反向求解只需要遍历全部节点一次便可计算出目标函数对各输入的偏导数,其求解过程分别如图4和图5所示。
步骤S3的具体步骤包括:
S31:搭建多个计算节点,每个计算节点设置一个存储单元,所述的存储单元内存储相关计算参数;
S32:根据多个计算节点中的计算参数,基于反向传播获取目标函数对输入量的梯度;
S33:利用梯度下降法优化***的控制量,获取最优控制序列,实现***的并行预测控制。
步骤S33的具体包括:
利用梯度下降法优化控制量:
uk|k,uk+1|k…uk+N-1|k,
其中,uk|k,uk+1|k…uk+N-1|k分别为第k时刻内第0、1…N-1步的控制量,
当满足优化条件其中之一时,完成优化过程,获取最优控制序列Uk *:
重复以上步骤,完成第k+1时刻的预测、优化过程。
最优控制序列的计算公式为:
本实施例中,***为车辆模型,所述的预测模型为车辆路径与速度跟踪模型,进行并行计算,考虑如下三自由度车辆非线性模型:
因此车辆连续非线性模型可以改写为
***状态预测(前向传播)
首先将连续***模型进行离散化,为了提高离散化精度使用三阶三段龙格-库塔公式进行离散化得到***离散化模型:
然后对控制输入赋初值,控制时域与预测时域均为N,由于是双输入***,设初值为2N行1列的列向量。
U=[u(0)T u(1)T … u(N-1)T]T
初始状态为x(0),状态预测根据公式(7)进行计算。
车辆路径与速度跟踪模型的控制目标是:
快速准确跟踪车辆纵向速度vx和横向位移Y,设控制时域与预测时域均为N,目标函数为:
其中,xi为当前时刻第i步的状态量,ui为当前时刻第i步的***控制量,Q为正定加权矩阵,P为正定终端惩罚矩阵,R为正定控制量惩罚矩阵。
目标函数进行计算时分解为:
Ji=h(yi,Ji-1)=Ji-1+(yi-yr)TQ(yi-yr)
=(xji-xjr)TQ(xji-xjr),j=1,6,i=1,2,…,N
其中,g、h分别为分解计算函数。
根据图2计算图模型,将其分为N层,前向预测和反向求导时每次计算一层。
其中Ac,k为f对x在(xk,uk)处的雅可比矩阵,Bc,k为f对u在(xk,uk)处的雅可比矩阵。
根据图2计算图模型,可以依次计算出每一层的局部偏导数如下:
对于第N-1层:
得到局部的偏导数之后,根据图6以及链式法则可以求:
以及
综上,输入加权矩阵R,Q和P,以及***的雅可比矩阵的公式,则可根据龙格-库塔公式预测状态,根Ak、Bk求解公式,求解未来状态xi对应的Ai和Bi,根据递推公式即可计算得到目标函数对优化变量的梯度:
每次迭代沿梯度相反方向更新控制输入
仿真实验结果
取预测时域和控制时域N=10,s=0.01,加权矩阵R=S=0,P=Q=diag(0.5,0,0.2,0.2,0,0.5)。所要跟踪的参考路径yref为:
其中,Xpre为横向预测位移。
车辆初速度为15m/s,参考速度vx,ref为:
参考速度曲线如图7所示,减速阶段加速度设置为-1m/s2。车辆路径和速度跟踪曲线结果如图8和9所示,位移和速度跟踪误差曲线如图10和11所示。
上述实施方式仅为例举,不表示对本发明范围的限定。这些实施方式还能以其它各种方式来实施,且能在不脱离本发明技术思想的范围内作各种省略、置换、变更。
Claims (4)
1.一种基于计算图的模型预测控制实时优化并行计算方法,其特征在于,包括以下步骤:
S1:搭建***状态量的预测模型,并构建***的目标函数,所述的***为车辆模型,所述的预测模型为车辆路径与速度跟踪模型;
S2:建立预测模型和目标函数的模型预测控制并行计算架构,所述的并行计算架构采用触发式并行计算方法同步计算预测模型和目标函数;
S3:采用反向传播的方式求解计算梯度,利用梯度下降法优化***的控制量,实现***的实时优化控制;
所述的车辆路径与速度跟踪模型的控制目标是:
快速准确跟踪车辆纵向速度vx和横向位移Y,设控制时域与预测时域均为N,目标函数为:
其中,xi为当前时刻第i步的状态量,ui为当前时刻第i步的***控制量,Q为正定加权矩阵,P为正定终端惩罚矩阵,R为正定控制量惩罚矩阵;
所述的步骤S2的模型预测控制并行计算架构中,用当前步序中的预测模型和目标函数的求解均完成的标志作为下一步序预测计算开始的标志,实现预测模型、目标函数的并行计算;
所述的步骤S3的具体步骤包括:
S31:搭建多个计算节点,每个计算节点设置一个存储单元,所述的存储单元内存储相关计算参数;
S32:根据多个计算节点中的计算参数,基于反向传播获取目标函数对输入量的梯度;
S33:利用梯度下降法优化***的控制量,获取最优控制序列,实现***的并行预测控制;
所述步骤S33的具体包括:
利用梯度下降法优化控制量:
uk|k,uk+1|k…uk+N-1|k,
其中,uk|k,uk+1|k…uk+N-1|k分别为第k时刻内第0、1…N-1步的控制量,
当满足优化条件其中之一时,完成优化过程,获取最优控制序列Uk *:
重复以上步骤,完成第k+1时刻的预测、优化过程;
所述的优化条件为当前迭代的步长的目标函数与上一步的目标函数的差值小于设定值或达到限定优化次数或目标函数变化量为0。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110344736.0A CN113110045B (zh) | 2021-03-31 | 2021-03-31 | 一种基于计算图的模型预测控制实时优化并行计算方法 |
US17/562,001 US11763166B2 (en) | 2021-03-31 | 2021-12-27 | Method for real time optimization and parallel computing of model prediction control based on computing chart |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110344736.0A CN113110045B (zh) | 2021-03-31 | 2021-03-31 | 一种基于计算图的模型预测控制实时优化并行计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113110045A CN113110045A (zh) | 2021-07-13 |
CN113110045B true CN113110045B (zh) | 2022-10-25 |
Family
ID=76712938
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110344736.0A Active CN113110045B (zh) | 2021-03-31 | 2021-03-31 | 一种基于计算图的模型预测控制实时优化并行计算方法 |
Country Status (2)
Country | Link |
---|---|
US (1) | US11763166B2 (zh) |
CN (1) | CN113110045B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115712243B (zh) * | 2022-11-04 | 2024-07-09 | 北京和利时工业软件有限公司 | 一种校正方法、装置、设备及介质 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109885883A (zh) * | 2019-01-21 | 2019-06-14 | 江苏大学 | 一种基于gk聚类算法模型预测的无人车横向运动的控制方法 |
CN111267834A (zh) * | 2020-02-17 | 2020-06-12 | 北京理工大学 | 一种车辆横摆稳定预测控制方法及*** |
Family Cites Families (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102981408B (zh) * | 2012-12-10 | 2015-05-27 | 华东交通大学 | 一种动车组运行过程建模与自适应控制方法 |
US20170139423A1 (en) * | 2015-11-12 | 2017-05-18 | King Fahd University Of Petroleum And Minerals | Control system and method for multi-vehicle systems |
CN106021364B (zh) * | 2016-05-10 | 2017-12-12 | 百度在线网络技术(北京)有限公司 | 图片搜索相关性预测模型的建立、图片搜索方法和装置 |
CN108038580A (zh) * | 2017-12-30 | 2018-05-15 | 国网江苏省电力公司无锡供电公司 | 基于同步挤压小波变换的光伏功率多模型综合预测方法 |
CN108900128B (zh) * | 2018-09-06 | 2021-09-28 | 吉林大学 | 基于模型预测控制的永磁同步电机直接转矩控制方法 |
CN110362084B (zh) * | 2019-07-22 | 2022-08-26 | 东南大学 | 一种基于蚁群算法的车辆路径追踪控制方法 |
CN110989625B (zh) * | 2019-12-25 | 2020-11-27 | 湖南大学 | 一种车辆路径跟踪控制方法 |
US11794775B2 (en) * | 2020-03-03 | 2023-10-24 | Motional Ad Llc | Control architectures for autonomous vehicles |
CN111665861A (zh) * | 2020-05-19 | 2020-09-15 | 中国农业大学 | 一种轨迹跟踪控制方法、装置、设备和存储介质 |
CN111626219B (zh) * | 2020-05-28 | 2023-06-09 | 深圳地平线机器人科技有限公司 | 轨迹预测模型生成方法、装置、可读存储介质及电子设备 |
DE102021203353A1 (de) * | 2021-04-01 | 2022-10-06 | Volkswagen Aktiengesellschaft | Verfahren zum Bestimmen einer Trajektorie eines zumindest teilweise assistiert betriebenen Kraftfahrzeugs, Computerprogramm sowie Assistenzsystem |
-
2021
- 2021-03-31 CN CN202110344736.0A patent/CN113110045B/zh active Active
- 2021-12-27 US US17/562,001 patent/US11763166B2/en active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109885883A (zh) * | 2019-01-21 | 2019-06-14 | 江苏大学 | 一种基于gk聚类算法模型预测的无人车横向运动的控制方法 |
CN111267834A (zh) * | 2020-02-17 | 2020-06-12 | 北京理工大学 | 一种车辆横摆稳定预测控制方法及*** |
Also Published As
Publication number | Publication date |
---|---|
US20220327388A1 (en) | 2022-10-13 |
US11763166B2 (en) | 2023-09-19 |
CN113110045A (zh) | 2021-07-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107688849B (zh) | 一种动态策略定点化训练方法及装置 | |
CN109472418B (zh) | 基于卡尔曼滤波的机动目标状态预测优化方法 | |
CN106767780B (zh) | 基于Chebyshev多项式插值逼近的扩展椭球集员滤波方法 | |
CN108287471B (zh) | Mimo偏格式无模型控制器基于***误差的参数自整定方法 | |
Deng et al. | A parallel Newton-type method for nonlinear model predictive control | |
CN105955031B (zh) | 非线性预测控制的fpga硬件加速控制器及其加速实现方法 | |
CN111158237B (zh) | 基于神经网络的工业炉炉温多步预测控制方法 | |
JP2004280792A (ja) | リアルタイムモデル予測制御における二次のプログラミングに対する加速されたアクティブセット探索のシステムおよび方法 | |
CN108170029B (zh) | Mimo全格式无模型控制器基于偏导信息的参数自整定方法 | |
CN113110045B (zh) | 一种基于计算图的模型预测控制实时优化并行计算方法 | |
CN108345213B (zh) | Mimo紧格式无模型控制器基于***误差的参数自整定方法 | |
CN109782586B (zh) | 参数自整定的miso异因子紧格式无模型控制方法 | |
CN108181809B (zh) | Miso紧格式无模型控制器基于***误差的参数自整定方法 | |
Wilamowski et al. | NBN algorithm | |
CN108287470B (zh) | Mimo偏格式无模型控制器基于偏导信息的参数自整定方法 | |
CN107942654B (zh) | Siso偏格式无模型控制器基于偏导信息的参数自整定方法 | |
CN108154231B (zh) | Miso全格式无模型控制器基于***误差的参数自整定方法 | |
JPH07200512A (ja) | 最適化問題解決装置 | |
CN108132600B (zh) | Mimo紧格式无模型控制器基于偏导信息的参数自整定方法 | |
CN114777762B (zh) | 一种基于贝叶斯nas的惯性导航方法 | |
CN113432608A (zh) | 适于ins/cns组合导航***的基于最大相关熵的广义高阶ckf算法 | |
CN113076596B (zh) | 融合车辆模型和轮胎力曲线的轨迹跟随并行预测方法 | |
CN108008634B (zh) | Miso偏格式无模型控制器基于偏导信息的参数自整定方法 | |
CN114967472A (zh) | 一种无人机轨迹跟踪状态补偿深度确定性策略梯度控制方法 | |
CN108107715B (zh) | Miso全格式无模型控制器基于偏导信息的参数自整定方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |