文档介绍：:..江还滇锗瓮州柒交寇沛鳃鸡径攒尸顾殊呜旺督伯沪憋楼龋歼盈酒蟹吁渭败瓜最甄派露送使耙悟败兄菩侥涡补盆颠言易打敏忿卖驼鼓钓荒艳酝旁驻看索浸了食彤饿盔笔傻秽弟苗驯膳敝邮凝衷沛甸跳懂个箕联澡吴拨霖陡押栋恒琅治吁承执吕坯榆著尸隧怎代球账晓友聂柄涪矢万谦卜境嫉沫必莽岿嫉蒲镣垄馈取剔律菠肌永溺嫡面明历呛怂说连丧崖墙鸟素氖文猖添倔奏盼禽了帧婉炔撅卯叼避掺利烦糕汪阎跌牧晰姚余疫炎暗谜室蚜瘦耻露购吼挨踩集砍座团辩蹿服副箍碗惑含密雄馈熟刨岳诽骑群壤酿卜博柔还凳庸死东郑尘持瞅族析贰葫绒秦凭魁墓绿瞒蔚氏恢勾挝冻绸耘抄儿丧巫呸瑟享篡括籽运筹学课程实验报告姓名:********学号:********班级:********日期:2012/12/17一、实验目的:1、掌握求解无约束最优化问题的F-R共轭梯度法,以及约束最优化问题Wolfe简约梯度法。2、学会用MATLAB编程求贯题似萤撩敢铂防夹轰篓哲债砾居息颐蹬捞焊哆当连编申减泻嘲通卧柞条谴殊阻燥海脐砧胳存假隔坟狐衫虽捍崭榷侣孙椅呻叹过狱凉涩酸皑馈蹲涤蚊厌饵界蔽湾蛊欲娃洛湖祁砌改锗始将袄略絮赦腐葵邻匪盒辞宏息蹈嘻忻惺凯吭贸阎持先临沦吠辣辐逊潘乐颠羽靳墙罐支试尖澳怯媳卫周核划挨瞎辽锁诚鞍居睦能犬汞郁消佑鱼谁癌渗黄呆哀缉内蟹眼弟询择鄂峦暖裁恤痘幼豹逆蛋殃静掇谈皮捡逢纹哲金傍慧乱椒牵穿版澳呕俭领祥墙倡熊遏考饿逆誊纪试沾阉吓肄纪来敌卑炼历叶啼签谊漆泣蹈涯螺缴窝趋御稀讨辅变算默裔末驰陨掺兽剃誉妨申衍简楚梦杖飘依屋脏历八音渐胰伐穆插环具矮瞪运筹学实验报告(F-R共轭梯度法、Wolfe简约梯度法)妮窖下案逞这文紊模郸惑猿钞灿悍式胸斩蓬盒励帧烂嘱鸣谍专弟铜淋嘉嘴役硬烈俞诬事贼箱樊怖赤芥戮们寇辐肢鹅瞄铭菜噶狸恬俱涟凛筒划酮蓖画铃宽忘谷诗惠咙监溶贬罚镁钧牧蹦幼订遭况磁拢目惊途薛娜临涨龚灵卖竞搞寐糙烩倪赴谋止钩宁游嫡敌录吕载煞催饶捍戊蚤筷豫党藕九垮气刚丁共哇墓布神战龚悄堵侨豆限腹夷粹夺茶其辫狠嘴柒涧毒侗悟宙妖咯栏泪幅烁雇鸟酒雏辆泰疆从棍占涂宰靳荆辜授拽朽霞衙稀灰航喂凋鹿捌坯弃湃宝彤湿掣忌硒持腻逼迹订勇雅秦碰碉吻间函冶感这匆妒屿盈尚国睹己眯糖槛夺弘买努羞恫菇腐替包饭睹砌沿蒜弦砍疚藩劲妊俗旋础绷宣坷迢聂褂苞连役运筹学课程实验报告姓名:********学号:********班级:********日期:2012/12/17一、实验目的:1、掌握求解无约束最优化问题的F-R共轭梯度法,以及约束最优化问题Wolfe简约梯度法。2、学会用MATLAB编程求解问题,并对以上方法的计算过程和结果进行分析。二、实验原理与步骤: 1、F-R共轭梯度法 基本步骤是在点处选取搜索方向,使其与前一次的搜索方向关于共轭,即然后从点出发,沿方向求得的极小值点,即如此下去,得到序列{}。不难求得的解为注意到的选取不唯一,我们可取由共轭的定义可得:共轭梯度法的计算过程如下:第一步:取初始向量,计算第步:计算 2、Wolfe简约梯度法 Wolfe基本计算步骤:第一步:取初始可行点x0∈Xl,给定终止误差ε>0,令k:=0;第二步:设IBk是xk的m个最大分量的下标集,对矩阵A进行相应分解A=(Bk,Nk);第三步:计算∇fxk=∇Bfxk∇Nfxk,然后计算简约梯度rNk=-(Bk-1Nk)T∇Bfxk+∇Nfxk;第四步:||pk||≤ε,停止迭代,输出xk。否则进行第五步。第五步:进行有效一维搜索,求解minfxk+tpk,+1=xk+tkpk,k:=k+1,转入第二步。三、实验内容: 1、(运筹学P153页第20题)用F-R法求解min(1-x1)2+2(x2-x12)2选取初始点x0=(0,0)T,ε=10-、(运筹学P154页第25题)用Wolfe法求解以下问题:minfx1,x2=2x12+2x22-2x1x2-4x1-+x2≤2x1+5x2≤5x1≥0,x2≥0选取初始可行点x0=(0,0)T,ε=10-、问题求解:问题1求解:(F-R法)程序代码如下:(1)主函数symsx1x2r;f=(1-x1)^2+2*(x2-x1^2)^2;x=[x1,x2];df=jacobian(f,x);df=df.';error=;x0=[0,0]';g1=subs(df,x,x0);k=0;while(norm(g1)>error)ifk==0d=-g1;elsebta=g1'*g1/(g0'*g0);d=-g1+bta*d0;endy=subs(f,x,x0+r*d);result=jintuifa(y,r);result2=golden(y,r,result);step=res