數值最優化—優化問題的解(二)
一、定理
區域性最小值點一階必要條件:
∇
f
(
x
∗
)
=
0
\nabla f(x^*)=0
∇f(x∗)=0
區域性最小值點二階必要條件:
∇
f
(
x
∗
)
=
0
且
∇
2
f
(
x
∗
)
\nabla f(x^*)=0 且 \nabla^2 f(x^*)
∇f(x∗)=0且∇2f(x∗) 正定。
區域性最小值點二階充分條件:
∇
2
f
(
x
)
\nabla^2 f(x)
∇2f(x) 在
x
∗
x^*
x∗的開鄰域內連續,
∇
f
(
x
∗
)
=
0
\nabla f(x^*)=0
∇f(x∗)=0並且
∇
2
f
(
x
∗
)
\nabla^2 f(x^*)
∇2f(x∗) 正定,那麼
x
∗
x^*
x∗是
f
(
x
)
f(x)
f(x) 的嚴格區域性最小值點。
二、證明
Proof 1
反證法
假設當
x
∗
x^*
x∗ 為區域性最小值時
∇
f
(
x
∗
)
≠
0
\nabla f(x^*)\neq0
∇f(x∗)=0 ,那麼定義
p
=
−
∇
f
(
x
∗
)
p = -\nabla f(x^*)
p=−∇f(x∗) 。這樣我們可以得到:
p
T
∇
f
(
x
∗
)
=
−
∇
f
(
x
∗
)
T
∇
f
(
x
∗
)
=
−
∥
∇
f
(
x
∗
)
∥
<
0
p^T\nabla f(x^*)=-\nabla f(x^*)^T\nabla f(x^*)=-\|\nabla f(x^*)\|<0
pT∇f(x∗)=−∇f(x∗)T∇f(x∗)=−∥∇f(x∗)∥<0 (1)
根據函式的連續性,必然可以得到,存在一個足夠小的正數
T
T
T ,使得:
p
T
∇
f
(
x
∗
+
t
p
)
<
0
p^T\nabla f(x^*+tp)<0
pT∇f(x∗+tp)<0 for all
t
∈
[
0
,
T
]
t\in[0,T]
t∈[0,T] (2)
考慮 \overline{t} \in(0,T] ,並對 f(x)$ 在
x
∗
x^*
x∗ 點處做泰勒展開,我們可以得到:
f
(
x
∗
+
t
‾
p
)
=
f
(
x
∗
)
+
t
‾
p
T
∇
f
(
x
∗
+
t
p
)
f
o
r
s
o
m
e
t
∈
(
0
,
t
‾
)
f(x^*+\overline{t}p) = f(x^*)+\overline{t}p^T\nabla f(x^*+tp) for some t\in (0,\overline{t})
f(x∗+tp)=f(x∗)+tpT∇f(x∗+tp)forsomet∈(0,t) (3)
觀察公式(3)的第二項,
t
‾
\overline{t}
t 是一個正實數,根據公式(2)我們又可以知道
p
T
∇
f
(
x
∗
+
t
p
)
p^T\nabla f(x^*+tp)
pT∇f(x∗+tp) 是一個負數。那麼我們可以我們一定可以找到一個正數
β
\beta
β ,使得:
f
(
x
∗
+
t
‾
p
)
=
f
(
x
∗
)
−
β
f(x^*+\overline{t}p) = f(x^*)-\beta
f(x∗+tp)=f(x∗)−β
也就是:
f
(
x
∗
+
t
‾
p
)
<
f
(
x
∗
)
f(x^*+\overline{t}p) < f(x^*)
f(x∗+tp)<f(x∗)
那麼 x^* 不是區域性最小值點,與假設相反,原結論成立。
Q.E.D
Proof 2:
由Proof 1我們已經證明了前半部分,現在證明後半部分,也就是
∇
2
f
(
x
∗
)
\nabla^2 f(x^*)
∇2f(x∗) 正定。
反證法
假設
∇
2
f
(
x
∗
)
\nabla^2 f(x^*)
∇2f(x∗) 不是正定的,那麼肯定可以找到一個向量
p
p
p 使得
p
T
∇
2
f
(
x
∗
)
p
<
0
p^T\nabla^2 f(x^*)p<0
pT∇2f(x∗)p<0 。同樣由連續性,可以得到,存在一個正數
T
T
T 使得:
p
T
∇
2
f
(
x
∗
+
t
p
)
p
<
0
p^T\nabla^2f(x^*+tp)p<0
pT∇2f(x∗+tp)p<0 for all
t
∈
[
0
,
T
]
t\in[0,T]
t∈[0,T]
考慮
t
‾
∈
(
0
,
T
]
\overline{t} \in(0,T]
t∈(0,T] ,並對
f
(
x
)
f(x)
f(x) 在
x
∗
x^*
x∗ 點處做二階泰勒展開,我們可以得到:
f
(
x
∗
+
t
‾
p
)
=
f
(
x
∗
)
+
t
‾
p
T
∇
f
(
x
∗
+
t
p
)
+
1
2
t
‾
2
p
T
∇
2
f
(
x
∗
+
t
p
)
p
<
f
(
x
∗
)
f(x^*+\overline{t}p) = f(x^*)+\overline{t}p^T\nabla f(x^*+tp)+\frac{1}{2}\overline{t}^2p^T\nabla^2f(x^*+tp)p<f(x^*)
f(x∗+tp)=f(x∗)+tpT∇f(x∗+tp)+21t2pT∇2f(x∗+tp)p<f(x∗)
同理我們可以得到假設錯誤,原結論成立。
Q.E.D
Proof 3
由於
∇
2
f
(
x
)
\nabla^2f(x)
∇2f(x) 的連續性,我們可以知道,在
x
∗
x^*
x∗ 的一個足夠小的鄰域內,
∇
2
f
(
x
)
\nabla^2f(x)
∇2f(x) 可以保持正定性。假設這個鄰域的半徑為
r
r
r ,那麼這個鄰域內所有點的集合 為
D
=
{
z
∣
∥
z
−
x
∗
∥
<
r
}
\mathcal{D} =\{z| \|z-x^*\|<r\}
D={z∣∥z−x∗∥<r} 。我們假設
p
p
p 為任意滿足
∥
p
∥
<
r
\|p\|<r
∥p∥<r 的向量。那麼我們可以得到
x
∗
+
p
∈
D
x^*+p\in \mathcal{D}
x∗+p∈D 。
(這裡做一點簡要的解釋,在
x
∗
x^*
x∗ 附近滿足正定性條件的點不一定都在集合
D
\mathcal{D}
D 裡面,之所以引入集合
D
\mathcal{D}
D 的概念主要是為了簡化描述。對於單一變數而言,
D
\mathcal{D}
D是一個以
x
∗
x^*
x∗為中點的線段,這個線段中任何一點都能使二階導正定。這個線段可能不是最長的線段,但是他必須是以
x
∗
x^*
x∗為中點的,二維的話就是以
x
∗
x^*
x∗ 為圓心的圓,三維就是以
x
∗
x^*
x∗ 為球心的球。)
同樣我們做泰勒展開:
f
(
x
∗
+
p
)
=
f
(
x
∗
)
+
p
T
∇
f
(
x
∗
)
+
1
2
p
T
∇
2
f
(
z
)
p
=
f
(
x
∗
)
+
1
2
p
T
f
(
z
)
p
f(x^*+p) = f(x^*) +p^T \nabla f(x^*) + \frac{1}{2}p^T\nabla^2 f(z) p\\=f(x^*) + \frac{1}{2} p^T f(z) p
f(x∗+p)=f(x∗)+pT∇f(x∗)+21pT∇2f(z)p=f(x∗)+21pTf(z)p
這裡需要解釋的有兩點。1、第二個等式是因為在
x
∗
x^*
x∗ 處
∇
f
(
x
∗
)
=
0
\nabla f(x^*) = 0
∇f(x∗)=0 。2、根據泰勒定理,這裡的
z
=
x
∗
+
t
p
z=x^*+tp
z=x∗+tp for some
t
∈
(
0
,
1
)
t \in (0,1)
t∈(0,1) 。對於最後一個等式的最後一項,根據正定的定義我們可以知道
p
T
f
(
z
)
p
p^T f(z) p
pTf(z)p 是一個正數。
我們可以得到:
f
(
x
∗
+
p
)
>
f
(
x
∗
)
f(x^*+p) > f(x^*)
f(x∗+p)>f(x∗)
由於
x
∗
+
p
x^*+p
x∗+p 可以認為是一個在超平面上以
x
∗
x^*
x∗ 為球心的開球域。這也就是說,在
x
∗
x^*
x∗ 附近所有的函式值都比
x
∗
x^*
x∗ 點處的大。
Q.E.D
相關文章
- 【機器學習】--迴歸問題的數值優化機器學習優化
- 凸優化問題優化
- 斜率優化(凸包優化)DP問題acm優化ACM
- matlab最優化問題的函式(fminbnd),fmincon,globalsearch,multistart(全域性區域性最優)Matlab優化函式
- go的編譯優化問題Go編譯優化
- Android效能優化——圖片優化(二)Android優化
- Oracle效能優化-SQL優化(案例二)Oracle優化SQL
- 效能優化(二) UI 繪製優化優化UI
- 03-凸優化問題優化
- 04 最優化方法優化
- [20211210]優化遇到的奇怪問題.txt優化
- MySQL之SQL優化詳解(二)MySql優化
- 【離散優化】覆蓋問題優化
- 記一個效能優化問題優化
- N皇后問題(各種優化)優化
- Android效能優化——效能優化的難題總結Android優化
- 優化-mysql子查詢索引失效問題解決優化MySql索引
- MySQL優化之系統變數優化MySql優化變數
- 關於vue的webpack打包優化問題VueWeb優化
- 優化GAN的分佈的梯度問題,WGAN優化梯度
- [20181119]使用sql profile優化問題.txtSQL優化
- webpack dll打包重複問題優化Web優化
- iOS問題整理08----效能優化iOS優化
- 資料庫sql的優化問題的面試題資料庫SQL優化面試題
- Shiro效能優化:解決Session頻繁讀寫問題優化Session
- hoj 1070 圖的m可著色優化問題優化
- GitHub 訪問優化Github優化
- 效能優化之關於畫素管道及優化(二)優化
- MySQL優化篇系列文章(二)——MyISAM表鎖與InnoDB鎖問題MySql優化
- 運籌優化(十七)--儲存論基礎及其最優化求解優化
- 運籌優化(十八)--對策論基礎及其最優化求解優化
- 運籌優化(十九)--決策論基礎及其最優化求解優化
- 運籌優化(十六)--排隊論基礎及其最優化求解優化
- Python小白的數學建模課-19.網路流優化問題Python優化
- 首屏優化系列(二)優化
- Mysql高階優化(二)MySql優化
- MySQL(二) MySql常用優化MySql優化
- sql優化專題SQL優化