拉格朗日傳 六 拉格朗日乘數法

{

\"code\": 200,

\"title\": \"\",

\"content\": \"語言定義\\n\\n設給定二元函數z=?(x,y)和附加條件φ(x,y)=0，為尋找z=?(x,y)在附加條件下的極值點，先做拉格朗日函數L(x,y)=?(x,y) λφ(x,y)，其中λ為參數。求L(x,y)對x和y的一階偏導數，令它們等於零，並與附加條件聯立，即\\n\\nL＇x(x,y)=?＇x(x,y) λφ＇x(x,y)=0，\\n\\nL＇y(x,y)=?＇y(x,y) λφ＇y(x,y)=0，\\n\\nφ(x,y)=0\\n\\n由上述方程組解出x，y及λ，如此求得的(x,y)，就是函數z=?(x,y)在附加條件φ(x,y)=0下的可能極值點。\\n\\n用\\\"拉格朗日乘數法\\\"求極值\\n\\n求函數f(x,y,z)在條件φ(x,y,z)=0下的極值\\n\\n方法（步驟）是：\\n\\n1.做拉格朗日函數L=f(x,y,z) λφ(x,y,z),λ稱拉格朗日乘數\\n\\n2.求L分彆對x,y,z,λ求偏導,得方程組,求出駐點P(x,y,z)\\n\\n如果這個實際問題的最大或最小值存在,一般說來駐點唯一,於是最值 可求.\\n\\n條件極值問題也可以化為無條件極值求解,但有些條件關係比較複雜,代換和運算很繁,而相對來說,\\\"拉格朗日乘數法\\\"不需代換,運算簡單一點.這就是優勢.\\n\\n條件φ(x,y,z)一定是個等式，不妨設為φ(x,y,z)=m\\n\\n則再建一個函數g(x,y,z)=φ(x,y,z)-m\\n\\ng(x,y,z)=0,以g(x,y,z)代替φ(x,y,z)\\n\\n在許多極值問題中，函數的自變量往往要受到一些條件的限製，比如，要設計一個容積為 V的長方體形開口水箱，確定長、寬和高, 使水箱的表麵積最小. 設水箱的長、寬、高分彆為 x,y,z， 則水箱容積V=xyz\\n\\n焊製水箱用去的鋼板麵積為 S=2xz 2yz xy\\n\\n這實際上是求函數 S 在 V 限製下的最小值問題。\\n\\n這類附有條件限製的極值問題稱為條件極值問題，其一般形式是在條件\\n\\n限製下，求函數F的極值\\n\\n條件極值與無條件極值的區彆\\n\\n條件極值是限製在一個子流形上的極值，條件極值存在時無條件極值不一定存在，即使存在二者也不一定相等。\\n\\n例如，求馬鞍麵 z=x.^2-y.^2 1 被平麵X·Z 平麵所截的曲線上的最低點。\\n\\n從其幾何圖形可以看出整個馬鞍麵冇有極值點，但限製在馬鞍麵被平麵 平麵所截的曲線上，有極小值 1，這個極小值就稱為條件極值。\\n\\n條件極值點的必要條件\\n\\n設在約束條件 之下求函數 的極值。滿足約束條件的點 是函數 的條件極值點 , 且在該點函數 滿足隱函數存在條件時, 由方程 定隱函數 , 於是點就是一元函數 的極限點, 有\\n\\n代入 , 就有\\n\\n( 以下 、 、 、 均表示相應偏導數在點 的值 . )\\n\\n即 - , 亦即 ( , ) , ) .\\n\\n可見向量( , )與向量 , )正交. 注意到向量 , )也與向量 , )正交, 即得向量( , )與向量 , )線性相關, 即存在實數 , 使\\n\\n( , ) , ) .\\n\\n亦即\\n\\nLagrange乘數法 :\\n\\n由上述討論可見 , 函數 在約束條件 之下的條件極值點應是方程組\\n\\n的解.\\n\\n引進所謂Lagrange函數\\n\\n, ( 稱其中的實數 為Lagrange乘數 )\\n\\n則上述方程組即為方程組\\n\\n因此，解決條件極值通常有兩種方法\\n\\n1）直接的方法是從方程組（１）中解出 並將其表示為 代入 消去 成為變量為 的函數將問題化為函數 的無條件極值問題；\\n\\n2）在一般情形下，要從方程組（1）中解出 來是困難的，甚至是不可能的，因此上麵求解方法往往是行不通的。通常采用的拉格朗日乘數法，是免去解方程組（1）的困難，將求 的條件極值問題化為求下麵拉格朗日函數的穩定點問題，然後根據所討論的實際問題的特性判斷出哪些穩定點是所求的極值的。\\n\\n解應用問題舉例\\n\\n例1\\n\\n拋物麵 被平麵 截成一個橢圓. 求該橢圓到座標\\n\\n原點的最長和最短距離.\\n\\n例3求函數 在條件\\n\\n下的極小值. 並證明不等式 , 其中 為任意正常數 .\\n\\n現在就以上麵水箱設計為例，看一看拉格朗日乘數法求解條件極值的過程\\n\\n解： 這個問題的實質是求函數\\n\\n在條件下的最小值問題， 應用拉格朗日乘法，令\\n\\nL='2*(x*z y*z) x*y v*(x*y*z-V)';\\n\\ndLdx=diff(L,'x')\\n\\ndLdy=diff(L,'y')\\n\\ndLdz=diff(L,'z')\\n\\ndLdv=diff(L,'v')\\n\\ndLdx =2*z y v*y*z\\n\\ndLdy =2*z x v*x*z\\n\\ndLdz =2*x 2*y v*x*y\\n\\ndLdv =x*y*z-V\\n\\n令 Ｌ 的各偏導等零，解方程組求穩定點\\n\\ns1='2*z y v*y*z';\\n\\ns2='2*z x v*x*z';\\n\\ns3='2*x 2*y v*x*y';\\n\\ns4='x*y*z-V';\\n\\n[v,x0,y0,z0]=s·lve(s1,s2,s3,s4)\\n\\nv =\\n\\n[ -2*2^(2\\/3)\\/V^(1\\/3)]\\n\\n[ -8*(-1\\/4*2^(1\\/3)*V^(1\\/3) 1\\/4*i*3^(1\\/2)*2^(1\\/3)*V^(1\\/3))^2\\/V]\\n\\n[ -8*(-1\\/4*2^(1\\/3)*V^(1\\/3)-1\\/4*i*3^(1\\/2)*2^(1\\/3)*V^(1\\/3))^2\\/V]\\n\\nx0 =[ 2^(1\\/3)*V^(1\\/3)]\\n\\ny0 =[ 2^(1\\/3)*V^(1\\/3)]\\n\\nz0 =[ 1\\/2*2^(1\\/3)*V^(1\\/3)]\\n\\n這裡顯然隻有實數解纔有意義，所以 Ｌ 的穩定點隻有下麵一個\\n\\n又已知所求的問題確實存在最小值，從而解出的穩定點就是最小值點，即水箱長寬與為高的２倍時用鋼板最省。\\n\\n例2\\n\\n再看一個條件極值求解問題\\n\\n拋物麵 被平麵 截成一個橢圓，求這個橢圓到座標原點的最長最短距離。(x73)\\n\\n解 這個問題的實質是求函數\\n\\n在條件 與 下的最大、最小值問題，應用拉格朗日乘法，令\\n\\nL='x^2 y^2 z^2 v*(x^2 y^2-z) h*(x y z-1)';\\n\\ndLdx=diff(L,'x')\\n\\ndLdy=diff(L,'y')\\n\\ndLdz=diff(L,'z')\\n\\ndLdv=diff(L,'v')\\n\\ndLdh=diff(L,'h')\\n\\ndLdx =2*x 2*v*x h\\n\\ndLdy =2*y 2*v*y h\\n\\ndLdz =2*z-v h\\n\\ndLdv =x^2 y^2-z\\n\\ndLdh =x y z-1\\n\\ns1='2*x 2*v*x h';\\n\\ns2='2*y 2*v*y h';\\n\\ns3='2*z-v h';\\n\\ns4='x^2 y^2-z';\\n\\ns5='x y z-1';\\n\\n[h,v,x0,y0,z0]=s·lve(s1,s2,s3,s4,s5);\\n\\nx0,y0,z0\\n\\nx0 =\\n\\n[ 3\\/4-1\\/4*i*13^(1\\/2)]\\n\\n[ 3\\/4 1\\/4*i*13^(1\\/2)]\\n\\n[ -1\\/2 1\\/2*3^(1\\/2)]\\n\\n[ -1\\/2-1\\/2*3^(1\\/2)]\\n\\ny0 =\\n\\n[ 3\\/4 1\\/4*i*13^(1\\/2)]\\n\\n[ 3\\/4-1\\/4*i*13^(1\\/2)]\\n\\n[ -1\\/2 1\\/2*3^(1\\/2)]\\n\\n[ -1\\/2-1\\/2*3^(1\\/2)]\\n\\nz0 = -1\\/2， -1\\/2， 2-3^(1\\/2)， 2 3^(1\\/2)\\n\\n即 的穩定點有兩個\\n\\n因為函數 在有界閉集 上連續，必有最大值和最小值，而求得的穩定點又恰是兩個，所以它們一個是最大點，另一個是最小，其最大\\n\\n最小值為。(x73)\\n\\nx1=-1\\/2 1\\/2*3^(1\\/2);\\n\\nx2=-1\\/2-1\\/2*3^(1\\/2);\\n\\ny1=-1\\/2 1\\/2*3^(1\\/2);\\n\\ny2=-1\\/2-1\\/2*3^(1\\/2);\\n\\nz1=2-3^(1\\/2);\\n\\nz2=2 3^(1\\/2);\\n\\nf1=(x1^2 y1^2 z1^2)^(1\\/2)\\n\\nf2=(x2^2 y2^2 z2^2)^(1\\/2)\\n\\nf1 = 0.5829 ; f2 = 4.2024\\n\\n拉格朗日乘數法舉例\\n\\n我們知道, 對於\\\"限製條件為等式,x值均為正值\\\"的最大化問題, 滿足最大化的x組合一定滿足: F(i)(x*)-Σλj Gj(i)(x*)=0, i=1,2,3,.....n, j=1,2,...m. 從這裡我們看到,如果限製條件 Gj(x*)=cj 中的 cj 變化 dcj , 如果全部作用於x(i),那麼引起的dx(i)=dcj\\/Gj(i)(x*),從而導致目標方程取值變化dF=F(i)(x*)dcj\\/Gj(i)(x*)=λj*dcj [注意:對於同一個限製條件j,我們由上一節已經知道必然有: F(i)(x*)\\/Gj(i)(x*)=F(i')(x*)\\/Gj(i')(x*)=λj (i不等於i')]. 那麼我們得到:λj=dF\\/dcj.也就是說,拉格朗日乘數其實代表的是cj對最大化目標函數F的邊際影響. 雖然這裡考慮的是僅僅cj發生變化,我們可以對此加以推廣,比如整體的c向量發生變化到 c dc, dc是一個m-維向量, 那麼F的總變化量dF就是Σλj dCj, j=1,2,...m.\\n\\n舉一個具體的實例: 假如一個計劃經濟體係下,政府實施如前所述的最大化問題(在有限資源如勞動力,自然礦產,人力資本等的限製下使社會整體效用\\/福利最大化),並已經找到了滿足最大化條件的x組合. 現在假設萬能的上帝允許該國的勞動力資源可以額外增加dc1, 那麼根據拉格朗日乘數的經濟學含義我們知道給整個社會帶來的福利將是λ1*dc1. 但是上帝說:要獲得這個額外的勞動力資源,你們必須以一定數量的其他資源比如土地來跟我交換,以示公平.那麼我們人類政府該拿多少土地來跟上帝換呢?指定該土地數量為dx2,那麼由此減少的社會福利是λ2*dc2. 如果λ1*dc1>λ2*dc2,上帝不會答應,如果反之我們不會答應.所以必然有λ1*dc1=λ2*dc2,也就是dc2=(λ1\\/λ2)dc1. 學過初級微觀的朋友馬上可以看出,這跟微觀經濟學中相對價格的概念十分相似.相對價格反映物與物之間的交換價值,即人們願意怎麼樣進行物與物的交換.不同的是,這裡的價格不是以錢來計算,而是以社會福利來衡量;這裡的相對價格λ1\\/λ2中的λ1和λ2是基於解決社會福利最大化問題而計算出來的,不同於市場中的價格P1,P2. 由於這個原因,我們把λ叫做\\\"影子價格\\\"(shad·w price). 如果我們偶爾發現某個市場經濟下市場價格之比恰恰等於影子價格之比,我們稱這個市場被一雙看不見的手所指引,因為該市場居然可以自發調整解決社會福利的最大化問題.\\n\\n再來考慮\\\"限製條件為非等式\\\"的情況. 我們知道市場價格通常都不可能為零或負數.但是影子價格確不同,它描述的是限製方程右方cj對整體目標函數值的邊際影響.在限製條件Gj為非等式的情況下,增加額外的cj不一定就意味著目標函數值的增加.比如, 限製條件為\\\"社會某消費產品不得高於cj\\\",目標函數為投資量.如果cj提高,那麼消費該產品增加,導致投資量減少,目標函數值減少.這時影子價格就是一個負值.再比如,目標函數為產量,限製條件為\\\"同時參加勞動的工人數量不得高於cj\\\".如果cj增加,那麼同時勞動的工人數量增加,可能導致勞動力邊際產量遞減效應的發生,這時總產量可能不增反降.這時我們情願不增加工人;換句話說,我們情願把一些資源放在一旁不予利用(free dsip·sal).這時候再增加這些勞動力資源,對總產量已經冇有作用了,所以影子價格為零. 事實上,根據前一節所述的庫恩-塔克定理,這一點是很明顯的.庫恩-塔克定理說,滿足最大化問題解的x一定使得下麵的條件滿足:\\n\\nLλ(x, λ)>=0, λ>=0, 互補鬆散\\n\\n就是說,如果Lλ(x, λ)=c-G(x*)>0, 那麼說明有資源餘缺閒置,這時λ=0.如果Lλ(x, λ)=c-G(x*)=0,那麼說明資源全部被使用,其邊際效用λ>0.\\n\\n注意:這裡我們通過對拉格朗日乘數的解釋考查了cj的微小變動dcj對目標函數最大值的變化的影響,這就是開篇所說的比較靜態研究--研究參數θ的變化對最大值的影響.所以我們在進行比較靜態研究的時候必須把目標函數看成是同時關於x和參數θ的函數. 基於這一點,我們從另一個角度來看λ的確定. 考察參數cj,如果cj變化一點點到cj dcj,那麼相應地最佳組合x*變動到x* dx*,最大目標值也由F(x*)變化為F(x* dx*).由泰勒一階展開我們得到: dF=F(x* dx*)-F(x*)=Fx(x*)dx* Fcj(x*)dcj.根據拉格朗日乘數法一階必要條件,我們有:Fx(x*)=λj Gx(x*),所以dF=λj Gx(x*)dx* Fcj(x*)dcj=λj Gx(x*)dx*,我們又知道根據限製條件方程G(x*)=cj,在cj變化到cj dcj的過程中,Gx(x*)dx*=dcj,所以dF=λj dcj.同樣推導出了λ的定義式. 更一般地, 如果F和G都是關於x和參數θ的函數,如果參數θ變動到θ dθ,x隨之變動到x dx,那麼:\\n\\ndF=F(x dx,θ dθ)-F(x,θ)=Fx(x,θ)dx Fθ(x,θ)dθ=λGx(x,θ)dx Fθ(x,θ)dθ...(1)\\n\\n由於G是關於x和θ的函數G(x,θ)=c,所以在θ變化的過程中始終有\\n\\nGx(x,θ)dx Gθ(x,θ)dθ=dc...................................................(2)\\n\\n代入(1)式,我們得到:\\n\\ndF=λdc -λGθ(x,θ)dθ Fθ(x,θ)dθ=Lθ(x,λ,θ)dθ λdc....................(3)\\n\\n這就是最一般化的比較靜態公式.我們在研究影子價格λ的時候,冇有考慮任何參數θ的變化,所以公式(3)的第一項為零,這樣dF=λdc.反之,我們在某些情況下不考慮c的變化,而側重於參數θ的變化,這時公式(3)變化為: dF=Lθ(x,λ,θ)dθ. 如果隻有函數F跟θ有關,而G跟θ無關,那麼公式(3)簡化為dF=Fθ(x,θ)dθ. 注意:1. 在參數θ變化的過程中,θ-->θ dθ,x-->x dx,但是對目標函數值的影響卻隻要考慮拉格朗日函數對θ的偏微分,而且該偏微分在原來最優點x處取值.這是我們用泰勒一階展開應該得到的結論. 2.這裡的x雖然冇有標上星號*,但不言自明的是它們都應該是最優組合,而且它們也都是關於參數θ的函數x(θ). 如果我們把最大化了的F定義成一個新函數最優目標方程V(θ),那麼由剛剛推導出來的公式(3): dF=Fθ(x,θ)dθ 我們有 Vθ(θ)=Fθ(x(θ),θ). 再次提醒注意,這裡的x(θ)是滿足最大化條件的最優點.如果我們再定義一個普通目標函數F(x',θ),但是這裡的x'是任意值,不一定是最優點x(θ).假設對應這個x'的能使 F 函數值最大的θ是θ'.那麼V(θ)在θ'點處的斜率為:Vθ(θ')=Fθ(x(θ'),θ'). 但我們知道,x(θ')=x'.所以Fθ(x(θ'),θ')=Fθ(x',θ').而後者就是函數F(x',θ)在點θ'的斜率.這就是說,函數V(θ)和函數F(x',θ)在點(x',θ')處的斜率相等. 這個結論對於x'取任意一個固定值都是成立的,所以從幾何圖形上來看,見圖示,最優目標函數V(θ)把普通目標函數曲線族緊緊包圍住.因此,dF=Fθ(x,θ)dθ 往往又稱為\\\"包絡定理\\\"(envel·pe the·rem). 微觀經濟學裡麵的短期成本和長期成本之間的關係就是符合信封定理的,因為這裡的成本都是滿足了成本最小化之後的成本。\\n\\n拉格朗日乘數法在消費者均衡原則中的應用\\n\\n微觀經濟學研究消費者行為時，所要闡述的核心問題是消費者均衡的原則。所謂消費者均衡指的是一個有理性的消費者所采取的均衡購買行為。進一步說，它是指保證消費者實現效用最大化的均衡購買行為。\\n\\n但人的需要或**是無限的，而滿足需要的手段是有限的。所以微觀經濟學所說的效用最大化隻能是一種有限製的效用最大化。而這種限製的因素就是各種商品的價格和消費者的貨幣收入水平。\\n\\n首先，我們先引入一些名詞解釋:\\n\\n總效用(TU):消費者在一定時間內消費一定數量某種商品或商品組合所得到的總的滿足。\\n\\n邊際效用(MU):消費者在所有其它商品的消費水平保持不變時，增加消費一單位某種商品所帶來的滿足程度的增加，也就是說指增加一單位某種商品所引起的總效用的增加。\\n\\n商品數量(Q)，商品價格(P), 收入(I)\\n\\n邊際效用的公式表達為:MU=?TU\\/?Q\\n\\n那麼如何才能實現在製約條件下效用最大化的商品組合呢？\\n\\n就是當消費者把全部收入用於購買各種商品時，他從所購買的每一種商品所得到的邊際效用與其價格的比例都相同，這樣的商品組合就是最佳的或均衡的商品組合。\\n\\n假設當消費者選擇兩種商品x,y時，消費者均衡原則的公式表達為:\\n\\nMUx ? MUy\\n\\nPx Py\\n\\n製約條件的公式表達式為:I=Px?Qx Py?Qy。那麼這一結論是如何推導出來的呢？解決這一問題最直接的方法就是拉格朗日乘數法。\\n\\n設效用函數U(Qx,Qy),為使它在製約條件下取得極值，首先建立拉格朗日函數:L=U(Qx,Qy) λ(I-Px?Qx-Py?Qy)，λ為參數。求L(x,y)對x和y的一階偏導數，令它們等於零，並與附加條件連立。\\n\\n即\\n\\n?L\\/?Qx=?U\\/?Qx-λPx=0 ⑴\\n\\n?L\\/?Qy=?U\\/?Qy-λPy=0 ⑵\\n\\nI-Px?Qx-Py?Qy=0 ⑶\\n\\n將方程⑴除以方程⑵，得:\\n\\n?U\\/?Qx ? Px 即 MUx ? MUy\\n\\n?U\\/?Qy Py PX Py\\n\\n所以，消費者要實現兩種商品的效用最大化，邊際效用的比率應該等於價格比率。\\n\\n以上是關於x和y兩種商品所說的，是否同樣適用於多種商品呢？答案是肯定的。如果消費者在n種商品中做出選擇，則消費者均衡的原則可表達為:\\n\\nMU1 ? MU2 ? MU3 ? … ? MUn\\n\\nP1 P2 P3 Pn\\n\\n這一結論同樣可用拉格朗日乘數法證明。\\n\\n拉格朗日乘數法可推廣到求n元函數?(x1,x2,…,xn)在m個附加條件φ(x1,x2,…,xn)下的條件極值。\\n\\n方法如下:\\n\\nm\\n\\n⑴做拉格朗日函數L(x1,x2,…,xn)=?(x1,x2,…,xn) ∑λiφi(x1,…x2);\\n\\ni=1\\n\\n⑵求L(x1,…xn)關於x1,…xn的偏導數，令它們等於零，並與附加條件聯立，即\\n\\nm\\n\\nL＇xi==?＇xi ∑λiφ'i=0 ,i=1,2,…,n\\n\\ni=1\\n\\nφk(x1,x2,…,xn)=0 ,k=1,2,…,n\\n\\n求解此方程組，可得到極值點。\\n\\n現在回到我們的問題中，設效用函數U(Qx1,Qx2,…Qxn),為使它在製約條件下取得極值,首先建立拉格朗日函數:\\n\\nL=U(Qx1,Qx2,…Qxn ) λ(I-Px1?Qx1-P2?Qy2-…-Pxn?Qxn)，λ為參數。求L(x1,x2,…xn)對x1,…,xn的一階偏導數，令它們等於零，並與附加條件聯立。\\n\\n即\\n\\n?L\\/?Qx1=?U\\/?Qx1-λPx1=0 （1）\\n\\n?L\\/?Qx2=?U\\/?Qx2-λPx2=0 （2）\\n\\n…… …\\n\\n?L\\/?Qxn=?U\\/?Qxn-λPxn=0 （n）\\n\\nI-Px1?Qx1-P2?Qy2-…-Pxn?Qxn\\n\\n將方程⑴到(n)相除，即得，\\n\\nMUx1 ? MUx2 ? … ? MUxn\\n\\nPx1 Px2 Pn\\n\\n所以，消費者要實現n種商品的效用最大化，邊際效用的比率應該等於價格比率。\\n\\n\"

}