Steps: k=0, 1, 2,…
F(α)=f(xk+α.dk)F(\alpha) =f(x_k+\alpha.d_k)F(α)=f(xk+α.dk)
xk,dkx_k, d_kxk,dk is fixed
∴\therefore∴其实就是一维问题,单变量的最优化问题。
我们要做的就是,告诉计算机两件事情:
不管有多少变量,多少问题。最优化问题都能很快解决
Golden-section search1980年,发明了一种算法:来来回回移动,套住一个区间,停在这里.
In matlab, inexact line search 中的 inex_lsearch.m
Ex:alpha = inex.lsearch(四个输入)
四个输入分别为
→\rightarrow→小的interval
控制问题一维的控制
不仅拉直,而且用到的能量最小
矩阵前两个为偏角,后两个为角速度
u(0),u(Δ),u(2Δ),...,u(nΔ)u(0), u(\Delta), u(2\Delta),..., u(n\Delta)u(0),u(Δ),u(2Δ),...,u(nΔ)
平方和最小 控制最小,能量最小
约束条件 四维向量 = 0
x12+x2=11x_1^2+x_2=11x12+x2=11
x1+x22=7x_1+x_2^2=7x1+x22=7
[x1x2]=[32]\begin{bmatrix} x_1\\ x_2\\ \end{bmatrix}=\begin{bmatrix} 3\\ 2\\ \end{bmatrix}[x1x2]=[32]
x2=11−x12x_2=11-x_1^2x2=11−x12
x1+(11−x12)2=7x_1+(11-x_1^2)^2=7x1+(11−x12)2=7
x14−22x12+x1+114=0x_1^4-22x_1^2+x_1+114=0x14−22x12+x1+114=0 一元四次方程
根据高斯代数积分定理
MatLab: x1x_1x1=roots([1 0 -22 1 114])
x1=[3.58443−3.7793−2.8051]x_1=\begin{bmatrix} 3.5844\\ 3\\ -3.7793\\ -2.8051\\ \end{bmatrix}x1=⎣⎢⎢⎡3.58443−3.7793−2.8051⎦⎥⎥⎤
x2x_2x2=11-x^2 enter =[−1.84812−3.28323.1313]\begin{bmatrix} -1.8481\\ 2\\ -3.2832\\ 3.1313\\ \end{bmatrix}⎣⎢⎢⎡−1.84812−3.28323.1313⎦⎥⎥⎤
f1(x)=x12+x2−11f_1(x)=x_1^2+x_2-11f1(x)=x12+x2−11
f2(x)=x1+x22−7f_2(x)=x_1+x_2^2-7f2(x)=x1+x22−7
f1(x)=0;f2(x)=0f_1(x)=0;f_2(x)=0f1(x)=0;f2(x)=0