主頁 > 軟體設計 > 數值最優化—優化問題的解(二)

數值最優化—優化問題的解(二)

2020-09-25 19:05:09 軟體設計

一、定理

區域最小值點一階必要條件: ? 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)+21?t2pT?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={zz?x?<r} ,我們假設 p p p 為任意滿足 ∥ p ∥ < r \|p\|<r p<r 的向量,那么我們可以得到 x ? + p ∈ D x^*+p\in \mathcal{D} x?+pD
(這里做一點簡要的解釋,在 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?)+21?pT?2f(z)p=f(x?)+21?pTf(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

轉載請註明出處,本文鏈接:https://www.uj5u.com/ruanti/127585.html

標籤:其他

上一篇:香港大學計算機院面試 for MSc in Computer Science (19/20 Round)

下一篇:基于立體視覺的三維模型重建系統設計

標籤雲
其他(157675) Python(38076) JavaScript(25376) Java(17977) C(15215) 區塊鏈(8255) C#(7972) AI(7469) 爪哇(7425) MySQL(7132) html(6777) 基礎類(6313) sql(6102) 熊猫(6058) PHP(5869) 数组(5741) R(5409) Linux(5327) 反应(5209) 腳本語言(PerlPython)(5129) 非技術區(4971) Android(4554) 数据框(4311) css(4259) 节点.js(4032) C語言(3288) json(3245) 列表(3129) 扑(3119) C++語言(3117) 安卓(2998) 打字稿(2995) VBA(2789) Java相關(2746) 疑難問題(2699) 细绳(2522) 單片機工控(2479) iOS(2429) ASP.NET(2402) MongoDB(2323) 麻木的(2285) 正则表达式(2254) 字典(2211) 循环(2198) 迅速(2185) 擅长(2169) 镖(2155) 功能(1967) .NET技术(1958) Web開發(1951) python-3.x(1918) HtmlCss(1915) 弹簧靴(1913) C++(1909) xml(1889) PostgreSQL(1872) .NETCore(1853) 谷歌表格(1846) Unity3D(1843) for循环(1842)

熱門瀏覽
  • 面試突擊第一季,第二季,第三季

    第一季必考 https://www.bilibili.com/video/BV1FE411y79Y?from=search&seid=15921726601957489746 第二季分布式 https://www.bilibili.com/video/BV13f4y127ee/?spm_id_fro ......

    uj5u.com 2020-09-10 05:35:24 more
  • 第三單元作業總結

    1.前言 這應該是本學期最后一次寫作業總結了吧。總體來說,對作業的節奏也差不多掌握了,作業做起來的效率也更高了。雖然和之前的作業一樣,作業中都要用到新的知識,但是相比之前,更加懂得了如何利用工具以及資料。雖然之間卡過殼,但總體而言,這幾次作業還算完成的比較好。 2.作業程序總結 相比前兩個單元,此單 ......

    uj5u.com 2020-09-10 05:35:41 more
  • 北航OO(2020)第四單元博客作業暨課程總結博客

    北航OO(2020)第四單元博客作業暨課程總結博客 本單元作業的架構設計 在本單元中,由于UML圖具有比較清晰的樹形結構,因此我對其中需要進行查詢操作的元素進行了包裝,在樹的父節點中存盤所有孩子的參考。考慮到性能問題,我采用了快取機制,一次查詢后盡可能快取已經遍歷過的資訊,以減少遍歷次數。 本單元我 ......

    uj5u.com 2020-09-10 05:35:48 more
  • BUAA_OO_第四單元

    一、UML決議器設計 ? 先看下題目:第四單元實作一個基于JDK 8帶有效性檢查的UML(Unified Modeling Language)類圖,順序圖,狀態圖分析器 MyUmlInteraction,實際上我們要建立一個有向圖模型,UML中的物件(元素)可能與同級元素連接,也可與低級元素相連形成 ......

    uj5u.com 2020-09-10 05:35:54 more
  • 6.1邏輯運算子

    邏輯運算子 1. && 短路與 運算式1 && 運算式2 01.運算式1為true并且運算式2也為true 整體回傳為true 02.運算式1為false,將不會執行運算式2 整體回傳為false 03.只要有一個運算式為false 整體回傳為false 2. || 短路或 運算式1 || 運算式2 ......

    uj5u.com 2020-09-10 05:35:56 more
  • BUAAOO 第四單元 & 課程總結

    1. 第四單元:StarUml檔案決議 本單元采用了圖模型決議UML。 UML檔案可以抽象為圖、子圖、邊的邏輯結構。 在實作中,圖的節點包括類、介面、屬性,子圖包括狀態圖、順序圖等。 采用了三次遍歷UML元素的方法建圖,第一遍遍歷建點,第二、三次遍歷設定屬性、連邊,實作圖物件的初始化。這里借鑒了一些 ......

    uj5u.com 2020-09-10 05:36:06 more
  • 談談我對C# 多型的理解

    面向物件三要素:封裝、繼承、多型。 封裝和繼承,這兩個比較好理解,但要理解多型的話,可就稍微有點難度了。今天,我們就來講講多型的理解。 我們應該經常會看到面試題目:請談談對多型的理解。 其實呢,多型非常簡單,就一句話:呼叫同一種方法產生了不同的結果。 具體實作方式有三種。 一、多載 多載很簡單。 p ......

    uj5u.com 2020-09-10 05:36:09 more
  • Python 資料驅動工具:DDT

    背景 python 的unittest 沒有自帶資料驅動功能。 所以如果使用unittest,同時又想使用資料驅動,那么就可以使用DDT來完成。 DDT是 “Data-Driven Tests”的縮寫。 資料:http://ddt.readthedocs.io/en/latest/ 使用方法 dd. ......

    uj5u.com 2020-09-10 05:36:13 more
  • Python里面的xlrd模塊詳解

    那我就一下面積個問題對xlrd模塊進行學習一下: 1.什么是xlrd模塊? 2.為什么使用xlrd模塊? 3.怎樣使用xlrd模塊? 1.什么是xlrd模塊? ?python操作excel主要用到xlrd和xlwt這兩個庫,即xlrd是讀excel,xlwt是寫excel的庫。 今天就先來說一下xl ......

    uj5u.com 2020-09-10 05:36:28 more
  • 當我們創建HashMap時,底層到底做了什么?

    jdk1.7中的底層實作程序(底層基于陣列+鏈表) 在我們new HashMap()時,底層創建了默認長度為16的一維陣列Entry[ ] table。當我們呼叫map.put(key1,value1)方法向HashMap里添加資料的時候: 首先,呼叫key1所在類的hashCode()計算key1 ......

    uj5u.com 2020-09-10 05:36:38 more
最新发布
  • 【中介者設計模式詳解】C/Java/JS/Go/Python/TS不同語言實作

    * 中介者模式是一種行為型設計模式,它可以用來減少類之間的直接依賴關系,
    * 將物件之間的通信封裝到一個中介者物件中,從而使得各個物件之間的關系更加松散。
    * 在中介者模式中,物件之間不再直接相互互動,而是通過中介者來中轉訊息。 ......

    uj5u.com 2023-04-20 08:20:47 more
  • 露天煤礦現場調研和交流案例分享

    他們集團的資訊化公司及研究院在一個礦區正在做智能礦山的統一平臺的 試點,專案投資大概1億,包括了礦山的各方面的內容,顯示得我們這次交流有點多余。他們2年前開始做智能礦山的規劃,有很多煤礦行業專家的加持,他們的描述是非常完美,但是去年底應該上線的平臺,現在還沒有看到影子。他們確實有很多場景需求,但是被... ......

    uj5u.com 2023-04-20 08:20:25 more
  • 《社區人員管理》實戰案例設計&個人案例分享

    設計是一個讓人夢想成真程序,開始編碼、測驗、除錯之前進行需求分析和架構設計,才能保證關鍵方面都做正確 ......

    uj5u.com 2023-04-20 08:20:17 more
  • 軟體架構生態化-多角色交付的探索實踐

    作為一個技術架構師,不僅僅要緊跟行業技術趨勢,還要結合研發團隊現狀及痛點,探索新的交付方案。在日常中,你是否遇到如下問題 “ 業務需求排期長研發是瓶頸;非研發角色感受不到研發技改提效的變化;引入ISV 團隊又擔心質量和安全,培訓周期長“等等,基于此我們探索了一種新的技術體系及交付方案來解決如上問題。 ......

    uj5u.com 2023-04-20 08:20:10 more
  • 【中介者設計模式詳解】C/Java/JS/Go/Python/TS不同語言實作

    * 中介者模式是一種行為型設計模式,它可以用來減少類之間的直接依賴關系,
    * 將物件之間的通信封裝到一個中介者物件中,從而使得各個物件之間的關系更加松散。
    * 在中介者模式中,物件之間不再直接相互互動,而是通過中介者來中轉訊息。 ......

    uj5u.com 2023-04-20 08:19:44 more
  • 露天煤礦現場調研和交流案例分享

    他們集團的資訊化公司及研究院在一個礦區正在做智能礦山的統一平臺的 試點,專案投資大概1億,包括了礦山的各方面的內容,顯示得我們這次交流有點多余。他們2年前開始做智能礦山的規劃,有很多煤礦行業專家的加持,他們的描述是非常完美,但是去年底應該上線的平臺,現在還沒有看到影子。他們確實有很多場景需求,但是被... ......

    uj5u.com 2023-04-20 08:19:07 more
  • 《社區人員管理》實戰案例設計&個人案例分享

    設計是一個讓人夢想成真程序,開始編碼、測驗、除錯之前進行需求分析和架構設計,才能保證關鍵方面都做正確 ......

    uj5u.com 2023-04-20 08:18:57 more
  • 軟體架構生態化-多角色交付的探索實踐

    作為一個技術架構師,不僅僅要緊跟行業技術趨勢,還要結合研發團隊現狀及痛點,探索新的交付方案。在日常中,你是否遇到如下問題 “ 業務需求排期長研發是瓶頸;非研發角色感受不到研發技改提效的變化;引入ISV 團隊又擔心質量和安全,培訓周期長“等等,基于此我們探索了一種新的技術體系及交付方案來解決如上問題。 ......

    uj5u.com 2023-04-20 08:18:49 more
  • 05單件模式

    #經典的單件模式 public class Singleton { private static Singleton uniqueInstance; //一個靜態變數持有Singleton類的唯一實體。 // 其他有用的實體變數寫在這里 //構造器宣告為私有,只有Singleton可以實體化這個類! ......

    uj5u.com 2023-04-19 08:42:51 more
  • 【架構與設計】常見微服務分層架構的區別和落地實踐

    軟體工程的方方面面都遵循一個最基本的道理:沒有銀彈,架構分層模型更是如此,每一種都有各自優缺點,所以請根據不同的業務場景,并遵循簡單、可演進這兩個重要的架構原則選擇合適的架構分層模型即可。 ......

    uj5u.com 2023-04-19 08:42:41 more