運筹学实验报告共轭梯度内外点法.docVIP

运筹与优化 实 验 报 告） 试验程序： function [x,minf]=minHJ(f,a,b,eps) format long; eps=0.16; syms x; f=2*x^2-x-1; a=-1; b=1; l=a+0.382*(b-a); u=a+0.618*(b-a); k=1; tol=b-a; while toleps f1=2*l^2-l-1; fu=2*u^2-u-1; if f1fu a=1; l=u; u=a+0.618*(b-a); else b=u; u=l; l=a+0.382*(b-a); end k=k+1; tol=abs(b-a); end k if k==10000 disp(找不到最小值!); x=NuN; minf =NuN; return end x=(a+b)/2; min=subs(f,findsym(f),x); format short; 运行结果： 黄金分割法的结果： K=7,ans=-1.0341 实验二 实验名称：MATLAB实现共轭梯度法 实验要求：用MATLAB实现共轭梯度法求解实例（） 试验程序： function f=conjugate(x0,t) x=x0; syms xi,yi,a; f=xi^2-xi*yi+yi^2+2*xi-4*yi; fx=diff(f,xi); fy=diff(f,yi); fx=subs(fx,{xi,yi},x0); fy=subs(fy,{xi,yi},x0); fi=[fx,fy]; count=0; while double(sqrt(fx^2+fy^2))t s=-fi; if count=0 s=-fi else x=x+a*s; f=subs(f,{xi,yi},x); f1=diff(f); f1=solve(f); if f1~=0 ai=double(f1); else break x.f=subs(f,{xi,yi},x),count end x=subs(x,a,ai); f=xi^2-xi*yi+yi^2+2*xi-4*yi; fxi=diff(f,xi); fyi=diff(f,yi); fxi=subs(fxi,{xi,yi},x); fyi=subs(fyi,{xi,yi},x); fii=[fxi,fyi]; 运行结果： 共轭梯度法：conjugate([0,0],0.000001) ,f=-4,count=2,ans=-4 实验三 实验名称：MATLAB实现外点法 实验要求：用MATLAB实现外点法和内点法求解实例: min s.t., 实验程序： m=zeros(1,50);a=zeros(1,50);b=zeros(1,50);f0=zeros(1,50); syms d x1 x2 e; m(1)=1;c=10;a(1)=0;b(1)=0; f=x1^2+x2^2+e*(1-x1)^2; f0(1)=1; fx1=diff(f,x1); fx2=diff(f,x2); for k=1:100 x1=a(k);x2=b(k);e=m(k); for n=1:100 f1=subs(fx1); f2=subs(fx2); if(double(sqrt(f1^2+f2^2))=0.002)