資源描述:
《數(shù)學(xué)實(shí)驗(yàn)教程_實(shí)驗(yàn)17(拉格朗日乘數(shù)法)》由會員上傳分享,免費(fèi)在線閱讀,更多相關(guān)內(nèi)容在行業(yè)資料-天天文庫。
1、實(shí)驗(yàn)17拉格朗日乘數(shù)法-99-實(shí)驗(yàn)17拉格朗日乘數(shù)法實(shí)驗(yàn)?zāi)康耐ㄟ^二維情形Lagrange乘數(shù)法的幾何觀察,幫助學(xué)生理解Lagrange乘數(shù)法,并掌握利用數(shù)學(xué)軟件求解帶等式約束條件的極值問題。預(yù)備知識條件極值、Lagrange乘子法實(shí)驗(yàn)內(nèi)容Lagrange乘數(shù)法是解有約束最優(yōu)化問題的一種方法。問題的形式為:最大化(或最小化);約束條件。【步驟】:為了具體,我們考察如下問題:最大化(或最小化);約束條件?!維tep1】:畫出約束曲線和的一些等值線,觀察可能的最優(yōu)值:圖17-1等值線與可行域圖從圖中觀察,要得到最優(yōu)值,需要尋找約束曲線和等值
2、曲線相切時(shí)的等值曲線值和它們的交點(diǎn)?!境绦颉浚篗athematica程序constrainpic=Plot[{Sqrt[1-x^2],-Sqrt[1-x^2]},{x,-1,1},AspectRatio?實(shí)驗(yàn)17拉格朗日乘數(shù)法-99-Automatic,Axes?True,PlotStyle?{{Thickness[0.01],RGBColor[0,0,1]}}];f[x_,y_]:=y-x^2;fpic=ContourPlot[f[x,y],{x,-1.4,1.4},{y,-1.4,1.4},ContourShading?False
3、];Show[{constrainpic,fpic}]【Step2】:尋找曲線之間相切的點(diǎn)由于的梯度是垂直于約束的(因?yàn)槭堑牡戎登€),而且的梯度是垂直于的等值曲線的,因此當(dāng)兩個(gè)梯度互為倍數(shù)關(guān)系時(shí),這兩條曲線就相切了。我們要做的就是解關(guān)于變量的方程組:我們可以利用Mathematica來解這個(gè)方程組?!境绦颉浚篗athematica程序g[x_,y_]:=x^2+y^2-1;gradf[x_,y_]={D[f[x,y],x],D[f[x,y],y]};gradg[x_,y_]={D[g[x,y],x],D[g[x,y],y]};Sol
4、ve[{gradf[x,y]lamdagradg[x,y],g[x,y]==0},{x,y,lamda}]注:并不是所有的這種方程都可由Solve求解,可以使用FindRoot命令求出某點(diǎn)附近的根。FindRoot[{g[x,y]0,gradf[x,y][[1]]lambdagradg[x,y][[1]],gradf[x,y][[2]]lambdagradg[x,y][[2]]},{x,1},{y,-0.6},{lambda,0.8}]{x?0.866025,y?-0.5,lambda?-1.}【Step3】:確定最大值或者最小值。實(shí)
5、驗(yàn)17拉格朗日乘數(shù)法-99-f[-Sqrt[3]/2,-1/2]f[Sqrt[3]/2,-1/2]f[0,-1]-1f[0,1]1最小值為:,最大值為:1。實(shí)驗(yàn)練習(xí)求位于橢球面上的最大值和最小值,并對得到的結(jié)果作出幾何解釋。