最優(yōu)化-馬昌鳳-第五章作業(yè)-全文預(yù)覽
【正文】 un39。 選取初始點為，調(diào)用函數(shù)程序，得出最小極值點為，極小值為，在界面框中輸入的命令如下：x0=[3 2]39。x0=x。*Bk*sk)+(yk*yk39。yk=feval(gfun,x,varargin{:})gk。break。 while(m20)%用Armijo搜索求步長 newf=feval(fun,x0+rho^m*dk,varargin{:})。 if(norm(gk)epsilon),break。k=0。 用算法求的極小點，選取初始點為解： 仿照書上編寫程序具體如下：function[x,val,k]=bfgs(fun,gfun,x0,varargin)%功能：用BFGS算法求解無約束問題：min f(x)；%輸入：x0是初始點，fun,gfun分別是目標(biāo)函數(shù)及梯度；%varargin是輸入的可變參數(shù)變量，簡單調(diào)用bfgs時可以忽略它；%但若其他程序循環(huán)調(diào)用該程序時將會發(fā)生重要的作用%輸出：x,val分別為近似最優(yōu)點和最優(yōu)值，k是迭代次數(shù)maxk=500。fun39。 然后仿照書上建立兩個目標(biāo)函數(shù)和梯度的M文件：function f=fun(x)f=x(1)^2+3*x(2)^2。*yk)。*yk0) Hk=Hk(Hk*yk*yk39。 end x=x0+rho^mk*dk。 while(m20)%用Armijo搜索求步長 if(feval(fun,x0+rho^m*dk)feval(fun,x0)+sigma*rho^m*gk39。%計算梯度 if(norm(gk)epsilon),break。 n=length(x0)。%給出最大迭代次數(shù)rho=。 解： 仿照書上編寫程序，將程序矩陣變?yōu)榫唧w如下：function[x,val,k]=dfp(fun,gfun,x0)%功能：用DFP算法求解吳宇舒問題：min f(x)%輸入：x0是初始點，fun,gfun分別是目標(biāo)函數(shù)及其梯度%輸出：x,va
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1