1. 程式人生 > >變分法與最速降線問題

變分法與最速降線問題

在諸多優化問題中,優化的變數是一個函式,也就是說需要找到一個函式,使得某個特定的目標最小。

固體力學中的最小勢能原理

處於平衡狀態的彈性體,真實位移場使得系統總勢能取最小值。

基於勢能原理,也發展出了固體力學中的變分解法,也就是求出使得總勢能最小的位移場,該位移就是真實的位移解答。

同樣變分法也可以用於最優控制中,求取控制律,使得系統系統的效能指標最優,而控制律通常是時間的函式。

1. 泛函

泛函其實是挺難理解的概念,首先看微積分中函式的概念。

定義在實數域RR上的一元函式y=f(x)y=f(x),值域也是實數集RR,實際上是定義了一種對映關係fff:RRf:~R \to R

而泛函的定義域並不是數集,而是函式集。

假設函式集合CCyC\forall y \in C,yy的定義域為[x0,x1][x_0,x_1]yy都為可微函式,且

y(x0)=y0,y(x1)=y1y(x_0)=y_0,~y(x_1)=y_1(1)

也就是說函式幾何CC中的元素均為平面上從點A(x0,y0)A(x_0,~y_0)到點B(x1,y1)~B(x_1,~y_1)的光滑曲線。

這裡討論泛函J(y)J(y)的定義域就是函式集合CC,其值域是實數集RR。泛函J(y)J(y)

(y)同樣定義了一種對映關係J:CRJ:~C \to R。一般泛函可以用積分式表示:

J[y(x)]=abF(x,y,y)dxJ[y(x)]=\int_{a}^{b}F(x,~y,~y')\text{d}x(2)

2. 泛函變分

函式y=f(x)y=f(x)的微分dydy指的是由於自變數xx+dxx \to x+dx所引起的因變數yy的變化量dy=f(x+dx)f(x)dy=f(x+dx)-f(x)

而泛函J[y(x)]J[y(x)]的變分δ

J\delta J指的是由於自變數yy+δyy \to y+\delta y所引起的因變數JJ的變化δJ=J[y+δy]J[y]\delta J=J[y+\delta y]-J[y]

yy+δyy \to y+\delta y時,利用多元函式的一階Taylor展開

F(x,y+δy,y+δy)=F(x,y,y)+Fyδy+FyδyF(x,~y+\delta y,~y'+\delta y')=F(x,~y,~y')+\frac{\partial F}{\partial y}\delta y+\frac{\partial F}{\partial y'}\delta y'(3)

因此,泛函δJ\delta J可以表為

δJ=ab[Fyδy+Fyδy]dx\delta J=\int_{a}^{b} \left[ \frac{\partial F}{\partial y}\delta y+\frac{\partial F}{\partial y'}\delta y' \right]\text{d}x(4)

3. 泛函極值

函式y=f(x)y=f(x)取在x=x0x=x_0處取極值的必要條件是

dydxx=x0=0\frac{dy}{dx}\Big|_{x=x_0}=0

泛函J[y(x)]J[y(x)]取極值的必要條件是

δJ=0\delta J = 0(5)

Fyddx(Fy)=0\frac{\partial F}{\partial y}-\frac{\text{d}}{\text{d}x}\left(\frac{\partial F}{\partial y'}\right)=0(6)

下面證明上述必要條件,證明過程很長,並且有點不嚴謹。

理論上δy\delta y是任意小的函式,可以取δy=εφ(x)\delta y=\varepsilon \varphi(x)ε\varepsilon是任意小的實數,φ(x)\varphi(x)為可微函式,並滿足

φ(a)=0,φ(b)=0\varphi(a)=0,~\varphi(b)=0(7)

這樣函式J[y(x)+εφ(x)]J[y(x)+\varepsilon\varphi(x)]ε\varepsilon的函式,在ε=0\varepsilon=0處取得極值。因此有

ddε(J[y(x)+εφ(x)])ε=0=0\frac{\text{d}}{\text{d}\varepsilon}\left(J[y(x)+\varepsilon\varphi(x)]\right)\Big|_{\varepsilon=0}=0(8)

參考式(3)有

J[y(x)+εφ(x)]=J[y(x)]+ab[Fyεφ(x)+Fyεφ(x)]dxJ[y(x)+\varepsilon\varphi(x)]=J[y(x)]+\int_{a}^{b}\left[\frac{\partial F}{\partial y}\varepsilon\varphi(x)+\frac{\partial F}{\partial y'}\varepsilon\varphi'(x)\right]\text{d}x(9)

將式(9)帶入式(8)中,

ddε(J[y(x)+εφ(x)])ε=0=ab[Fyφ(x)+Fyφ(x)]dx=0\frac{\text{d}}{\text{d}\varepsilon}\left(J[y(x)+\varepsilon\varphi(x)]\right)\Big|_{\varepsilon=0}=\int_{a}^{b}\left[\frac{\partial F}{\partial y}\varphi(x)+\frac{\partial F}{\partial y'}\varphi'(x)\right]\text{d}x=0(10)

再以ε\varepsilon乘以上式有

ab[F