运筹学第一次实验.docVIP

熟练利用0.618法，共轭梯度法以及非二次函数的共轭梯度法求解相关问题。 问题描述 0.618法： min f(x)=，[]=[-1,1],精度L 共轭梯度法(CG算法)： min f(x)=+2 非二次的共轭梯度法（CG算法推广）： min f(x)= (三) 算法介绍 (1)0.618法 Step1：对区间[a,b]=[]中取两点：=， Step2：若,停止运算，输出结果，否则计算并比较。 若，则，，，) 若，则，，，) Step3:置k：=k+1，返回Step2 (2)CG算法原理及步骤： n元正定二次函数f(x)=1/2,给定任一初始点,计算，置k：=0 step1：, step2：判断，若成立则终止，输出。否则转step3 step3：计算下一次的搜索方向，其中 step4：置k:=k+1,转step1进行下一次迭代。 (3)CG算法的推广： 对改用直接的e.l.s，求得步长因子。 对 （四）程序代码及运行结果： （1）0.618法源程序代码 a=-1;%区间端点 b=1; %区间端点 fprintf(初始区间为); interval=[a,b] e=0.16;%误差 q=b-a; u1=a+0.382*(b-a); u2=a+0.618*(b-a); n=0; while qe n=n+1; f1=2*u1^2-u1-1; f2=