最优化思路第三步
Steps: k=0, 1, 2,…
F(α)=f(xk+α.dk)
xk,dk is fixed
∴其实就是一维问题,单变量的最优化问题。
我们要做的就是,告诉计算机两件事情:
- xk点的确定,方向在哪里
- 什么地方停下来
不管有多少变量,多少问题。最优化问题都能很快解决
Golden-section search
1980年,发明了一种算法:来来回回移动,套住一个区间,停在这里.
In matlab, inexact line search 中的 inex_lsearch.m
Ex:
alpha = inex.lsearch(四个输入)
四个输入分别为
- xk
- dk
- 函数程序 ftest.m
- 梯度程序 gtest.m
→小的interval
控制问题
一维的控制
不仅拉直,而且用到的能量最小
x(t)⎣⎢⎢⎢⎡Θ1(t)Θ2(t)Θ1′(t).Θ2′(t).⎦⎥⎥⎥⎤
矩阵前两个为偏角,后两个为角速度
u(0),u(Δ),u(2Δ),...,u(nΔ)
平方和最小 <—> 控制最小,能量最小
约束条件 四维向量 = 0
x12+x2=11
x1+x22=7
[x1x2]=[32]
x2=11−x12
x1+(11−x12)2=7
x14−22x12+x1+114=0 一元四次方程
根据高斯代数积分定理
MatLab: x1=roots([1 0 -22 1 114])
x1=⎣⎢⎢⎡3.58443−3.7793−2.8051⎦⎥⎥⎤
x2=11-x^2 enter =⎣⎢⎢⎡−1.84812−3.28323.1313⎦⎥⎥⎤
f1(x)=x12+x2−11
f2(x)=x1+x22−7
f1(x)=0;f2(x)=0
目标函数: f(x)=(x12+x2−11)2+(x1+x22−7)2 求最小值
▽f=[4(x12+x2−11)x1+2(x1+x22−7)2(x12+x2−11)+4(x1+x22−7)x2]