翁云華，杜　娟，王雪嬌（成都理工大學管理科學學院，四川成都610059）

變分法應用條件的探索及其靈敏度分析

翁云華，杜娟，王雪嬌
（成都理工大學管理科學學院，四川成都610059）

考慮和探究變分法中的各種條件，通過對極值必要條件的探索和對不固定的自然條件的端點進行分類討論，以及對橫截性條件的端點問題進行探究，最終得到對變分法應用條件的總結(jié)以及各種情形的分析討論，使得變分法可以更好地拓展和運用到許多實際問題中.在此基礎(chǔ)上對變分法的各類討論情形進行了靈敏度分析，考察當參數(shù)變化時，原目標函數(shù)的最優(yōu)解變化情況，以及變化率的高低，推導并得出可直接求變分法的靈敏度的公式.由此簡化了獲取靈敏度的難度，及復雜性，并將所得的靈敏度公式應用于實踐.

歐拉拉格朗日方程；自然橫截性條件；變分法；橫截性條件；靈敏度分析

Weng YH,Du J,Wang XJ.The Exploration of the Condition of Variation Method and Sensityvity[J].Journal of Yibin University,2015,15（6）：117-120.

變分法是17世紀末葉開始發(fā)展起來的一個數(shù)學分支，它是為了解決實踐問題中的極值問題而產(chǎn)生的.不過它所討論的不是求普通函數(shù)的極值，而是泛函的極值問題.近年來，對變分問題各類條件的探索已經(jīng)成為人們研究的熱點.由老大中編寫的《變分法基礎(chǔ)》一書中對變分法的4類應用條件作了詳細的描述和總結(jié).本文在此基礎(chǔ)上，考慮變分的每一個變量都與參數(shù)p有關(guān)，觀察當參數(shù)改變時，目標函數(shù)相應的最優(yōu)解的變化，推導出直接計算目標函數(shù)靈敏度的公式.主要意義在于，對給定的目標函數(shù)，可以直接利用公式計算其靈敏度，從而大大減小了計算量，節(jié)約了計算成本.

1　準備知識

問題[1]：設(shè)函數(shù)F（x）在[a,b]上連續(xù)，任意函數(shù)u（x）在區(qū)間[a,b]上具有n階連續(xù)導數(shù)，考慮變分法在約束條件下的極值問題（泛函求極值問題）

約束條件為向量

其中向量

為了分析這個問題的平衡條件（目標函數(shù)取得極值時的條件），定義函數(shù)關(guān)于其參數(shù)變分，函數(shù)J與參數(shù)u,a,b有關(guān).變分δJ定義為：

2　條件探索

定理1[2]若J和H定義如上，設(shè)u*（t）為問題（1）-（2）的局部極值，設(shè)集合G={u（t）|：H（u（t））=0}非空，并且設(shè)J,H在u*的一個領(lǐng)域內(nèi)，J和Hi是弱連續(xù)的.如果對于G中的任一函數(shù)u（x）都有 δJ= ΔJ=J[u（x）]-J[u*（x）]=0，則泛函J[u（x）]在u*（x）出取得極值.

極值必要條件[3]：下述條件（I）和條件（II）必有一個成立.

條件（I）：對于所有增量 δu1,…,δum和增量δH1,…,δHm有：

條件（II）：對于u,a,b的任意增量δu,δa,δb，存在常數(shù)使得

若（5）不成立，則（6）可化為

其中 δu（t）、δa、δb為u（t）、a、b各自的增量，λ*= （…,），H=（H1,…,Hm），δa=δ（t）|a.

以上等式要求F和H對u和u′的偏導數(shù)存在.在下面的討論中，假設(shè)所有涉及到的函數(shù)充分光滑，在實際應用中通常會遇到以下四種情形：

情形1[4]：由于δu的任意性（7）式可化簡為滿足極值條件的歐拉－拉格朗日方程：

把這樣定義的泛函變分是由拉格朗日泛函變分法得到的，稱之為拉格朗日泛函變分，它與前面定義的變分是等價的，而且也便于計算泛函的變分.

為方便書寫，引進符號：

所以歐拉－拉格朗日方程可簡記為：

情形2：如果u（a）或u（b）非定值，由δu（a）的任意性，可將（7）式化簡得到關(guān)于u的邊值不固定的自然邊界條件：

例1[5]：求泛函J[y]=[p（x）y′2+q（x）y2+2f（x）y]d x極值問題的自然邊界條件，其中x0和x1均為自由邊界，p（x）,q（x）和 f（x）均為已知函數(shù)，且p（x）≠0.

解：因為x0和x1均為自由邊界，根據(jù)定理1及情形1、2的自然邊界條件為

情形3：如果端點a或b為自由端點.則自由端點a或b滿足泛函求駐值的必要條件，這種條件稱之為端點不固定自然邊界條件，同時把固定邊邊界條件稱之為強制邊界條件或本質(zhì)邊界條件，由于δ（a）或δ（b）的任意性，有：

情形4[6]：如果右端點固定，而端點a在曲線Ψ（t）上滑動，由y=u（x）與已知曲線y=Ψ（x）交點處的u′與Ψ′的斜率關(guān)系，有橫截性條件（斜截條件）：

解：由以上橫截性條件可得：

化簡得：

3　靈敏度分析

考慮帶參數(shù)組的變分問題，分析當參數(shù)變化時相應的最優(yōu)解的變化，為了方便起見，首先考慮有限個參數(shù)的情形，考慮問題[8]

約束條件為

其中

且0∈Rm,p=（p1,…,pm）∈Rk是參數(shù)向量.

這個問題能夠取到最優(yōu)解u*的必要條件已經(jīng)通過歐拉拉格朗日方程（10）給出，現(xiàn)在假設(shè)約束條件（15）以及依賴于這個問題的等式（11）-（13）都與參數(shù)p有關(guān).

為獲得方程組的靈敏度，計算所有參數(shù)的變分.由（14）和（16）得（17）：

其中算子ν定義如下：

如果邊界是固定的，由變分的邊界條件得

在端點a處如果u（a）是不固定的，則有變分自然橫截性條件[9]：

同樣地，如果a是不固定的，則有變分的自然邊界條件為：

如果在端點a處，若端點在曲線?（t；p）（對比（13）式）上滑動，有變分的橫截性條件：

其中：

對于每一個λ,p,δλ,δp都在一個非齊次線性二階積分方程中，和（24）一起構(gòu)成邊值問題，它的解結(jié)合（17）、（18）可獲得關(guān)于J,u,λ及其相應的靈敏度.通常用單位矩陣來代替δp.需要注意的是（17）、（18）、（24）給出了關(guān)于δJ,δp,δλ的線性方程組.

所以原問題（14）-（16）的目標函數(shù)的關(guān)于參數(shù)p的靈敏度向量可給出如下：

這是關(guān)于參數(shù) p的拉格朗日函數(shù)在最優(yōu)解u*,λ*處的梯度.

雖然原問題和雙變量的靈敏度很難獲得，以上實際結(jié)論卻為提供了計算靈敏度的直接公式.

對于有無限個參數(shù)的目標函數(shù)的靈敏度也可以類似得到.下面簡單看一下求靈敏度的實際應用

例3：考慮下面關(guān)于參數(shù)p的問題

約束條件為

由于已經(jīng)給出靈敏度的公式，所以可以直接進行計算.

為了計算J,λ,u（t；p）的偏導數(shù)，首先解決關(guān)于p的參數(shù)問題，問題的歐拉方程為：

其解為：

在（29）中應用邊界條件，得到最優(yōu)解

則u（t,p）,λ,J的偏導數(shù)為

由此便得到了目標函數(shù)的靈敏度.

4　結(jié)語

方程（10）-（13）是變分法求極值的著名的必要條件，對泛函求極值提供了事半功倍的作用.在自然科學和工程技術(shù)中所遇到的變分問題，有時要求極值函數(shù)除滿足給定的邊界條件外，還要滿足一定的附約束加條件，這些附加條件也可直接在變分法的應用條件基礎(chǔ)上再上附加條件，就可以直接應用各類條件即可，這就是泛函的求條件極值問題.泛函的條件極值的計算方法與函數(shù)的條件極值的計算方法類似，很多問題也可用拉格朗日乘數(shù)法來實現(xiàn).

變分方法已經(jīng)被越來越多的拓展和運用到許多實際問題中.（25）式給出了目標函數(shù)直接求靈敏度的公式，為在應用實踐中求靈敏度提供了很大的方便，同時對于變分法的靈敏度分析方法也可適用于（線性和非線性）最優(yōu)化問題，以及優(yōu)化控制問題，其靈敏度的獲取方法與求變分法的靈敏度相似.

[1]Bonnans JF,Shapiro A.Perturbation analysis ofoptimization prob?lems[M].NW：Springer,2000.

[2]Caillau JB,Noailles J.Sensitivity analysis for time optimal orbit transfer[J].Optimization：A Journal of Mathematical Programming and OperationsResearch,2001,49（4）：327-350.

[3]Castillo E,Conejo A,Pedregal P,etal.Building and solvingmathe?matical programming models in engineering and science[M].Pure ApplMath,2001.

[4]Conejo A,Castillo E,Mínguez R,et al.Decomposition techniques inmathematical programming engineering and science applications [M].NW：Springer,2011.

[5]Gelfand IM,Fomin SV.Calculus of Variations[M].NJ：Englewood Cliffs,1963.

[6]老大中.變分法基礎(chǔ)[M].北京：國防工業(yè)出版社,2004.

[7]Castillo E,Conejo A J,Castillo C,et al.Perturbation approach to sensitivity analysis inmathematical programming[J].JOptim Theo?ry and Applications,（2006）[2015-02-14].http：//citeseerx.ist.psu.edu/viewdoc/summary？doi=10.1.1.378.3590.

[8]Castillo E,Conejo A,Castillo C.A closed formula for localsensitivi?ty analysis inmathematical programming[J].Engineering Optimiza?tion,2006（38）：93-112.

[9]Gal T,Greenberg H J.Advances in sensitivity analysis and para?metric programming[M].NW：Springer,1997.

（編校：許潔）

The Exploration of the Condition of Variation Method and Sensitivity Analysis

WENGYunhua,DU Juan,WANGXuejiao
（College ofManagementScience,Chengdu University ofTechnology,Chengdu,Sichuan 610059,China）

The various conditionsofvariationalmethod wasexplored through theexploration ofnecessary conditionsofex?treme and classified discussion on the endpoints of unfixed natural conditions.Endpoint issues of transversality condi?tionswere explored so as to eventually get the summary ofapplication conditions of variationalmethod,ensuring the bet?ter developmentand further application of variationalmethod in solving practical problems.Sensitivity analysiswere then carried outon the basis of the above research of variationalmethod to observe the shiftof optimal solution of original ob?jective function with parameters change,aswellas the rate of change.And a formula of directly getting the sensitivity rate ofvariationalmethod were derived.Thissimplifies complexity in accessing sensitivity issues.

Euler-Lagrange equations；natural and transversality conditions；calculus of variations；transversality condi?tions；sensitivity analysis

O29

A

1671-5365（2015）06-00117-04

2015-03-24修回：2015-04-08

翁云華（1991-），男，碩士研究生，研究方向為應用泛函分析

網(wǎng)絡出版時間：2015-04-09 09：59網(wǎng)絡出版地址：http：//www.cnki.net/kcms/detail/51.1630.Z.20150409.0959.001.html

引用格式：翁云華，杜娟，王雪嬌.變分法應用條件的探索及其靈敏度分析[J].宜賓學院學報,2015,15（6）：117-120.