实验2解无约束非线性规划.docxVIP

PAGE PAGE 1 实 验 报 告 解无约束非线性规划实验（运筹学与最优化方法，4 学时） — 实验目的 1．掌握迭代算法的思想。 2．掌握 0.618 法。 3．掌握最速下降法、共轭梯度法和拟牛顿法。 二 实验内容 用 0.618 法解下列问题： min f (x) ? 2x2 ? x ?1 初始区间为：[a , b ] ? [?1,1], ? ? 0.16 。 1 1 用最速下降法、共轭梯度法和拟牛顿法分别解下列问题： min f (x , x ) ? x2 ? 4x2 取(x(0) , x(0) ) ? (1,1) 1 2 1 2 1 2 三 实验步骤（算法）与结果 解:首先建立一个黄金分割计算最优值的函数文件并保存为HJFG， 可供调用： x(1)=input(a=); y(1)=input(b=); k=input(k=); n=1;m=1; p(1)=x(1)+0.382*(y(1)-x(1)); q(1)=x(1)+0.618*(y(1)-x(1)); while y(n)-x(n)=k if g(p(n))g(q(n)) n=n+1; x(n)=p(n-1);y(n)=y(n-1);p(n)=q(n-1); q(n)=x(n)+0.618*(y(n)-x(n)); elseif g(p(n))=g(q(n)) n=n+1; x(n)=x(n-1);y(n)=q(n-1);q(n)=p(n-1); p(n)=x(n)+0.382*(y(n)-x(n)); end end x(n),y(n),min=1/2*(x(n)+y(n)) 参数说明：x(1)：计算区间的左极限； y(1)：计算区间的右极限； k： 计 算 要 求 达 到 的 精 度 ； p(n)，q(n)：黄金分割的计算公式； g(x)： 输 入 的 计 算 函 数 ； min=1/2*(x(n)+y(n))：满足条件的最优解. 实际使用： 首先建立函数文件并保存： function G=g(x) G=2*x^2-x-1; 然后调用上面的函数 输出结果为： ans =0.1672 ans =0.2787 min =0.2229 也就是满足条件的解的存在区间之一是[0.1672,0.2787]，取平均值 0.2229 作为近似最优解. 解： 先建立所求目标的函数文件； function y=fun2(a,b) y=a.^2+4*b.^2; 以及他们的偏导数： function dy=dx1_fun(a) dy=2*a; function dy=dx2_fun(b) dy=8*b; 然后是主程序： 最速下降法： function f=zs(x1,x2,eps) x(1,1)=x1;x(1,2)=x2; k=1; while 1 c(1)=dx1_fun(x(k,1));c(2)=dx2_fun(x(k,2)); if norm(c,2)=eps d=-c; a=sym(a); g=diff(fun2(x(k,1)+a*d(1),x(k,2)+a*d(2)),a); a=eval(solve(g,a)); x(k+1,:)=x(k,:)+a*d; k=k+1; else break end end x 运行： zs(1,1,0.01) 结果是： x = 1.0000 1.0000 0.7385 -0.0462 0.1108 0.1108 0.0818 -0.0051 0.0123 0.0123 0.0091 -0.0006 0.0014 0.0014 0.0010 -0.0001 共轭梯度法： function f2=GE(x1,x2,eps) x(1,1)=x1;x(1,2)=x2; k=1; d=0; q=[1,0;0,4]; c(1)=dx1_fun(x(k,1));c(2)=dx2_fun(x(k,2)); if norm(c,2)=eps d=-c; a=sym(a); f=diff(fun2(x(k,1)+a*d(1),x(k,2)+a*d(2)),a); a=eval(solve(f,a)); x(k+1,:)=x(k,:)+a*d; end while 1 k=k+1; c(1)=dx1_fun(x(k,1));c(2)=dx2_fun(x(k,2)); if norm(c,2)eps r=(d*q*c)/(d*q*d); d=-c+r*d; a=sym(a); f=diff(fun2(x(k,1)+a*d(1),x(k,2)+a*d(2)),a); a=eval(solve(f,a)); x(k+1,:)=x(k,:)+a*d; else break end end x 运行： GE(1,1,0.