大连理工大学优化方法上机大作业.docxVIP

大连理工大学优化方法上机大作业精选文档

TTMSsystemofficeroom【TTMS16H-TTMS2A-TTMS8Q8-TTMSHHJ8】

2016年大连理工大学优化方法上机大作业

学院：

专业：

班级：

学号：

姓名：

上机大作业1：

1.最速下降法：

g=zeros(2,1);

functionf=fun(x)

f=(1-x(1))^2+100*(x(2)-x(1)^2)^2;

end

functiong=grad(x)

g(1)=2*(x(1)-1)+400*x(1)*(x(1)^2-x(2));

g(2)=200*(x(2)-x(1)^2);

end

functionx_star=steepest(x0,eps)

gk=grad(x0);

res=norm(gk);

k=0;

whileresepsk=1000

dk=-gk;

ak=1;f0=fun(x0);

f1=fun(x0+ak*dk);

slope=dot(gk,dk);

whilef1f0+*ak*slope

ak=ak/4;

xk=x0+ak*dk;

f1=fun(xk);

end

k=k+1;

x0=xk;

gk=grad(xk);

res=norm(gk);

fprintf(--The%d-thiter,theresidualis%f\n,k,res);

end

x_star=xk;

end

clear

x0=[0,0];

eps=1e-4;

x=steepest(x0,eps)

2.牛顿法：

g=zeros(2,2);

functionf=fun(x)

f=(1-x(1))^2+100*(x(2)-x(1)^2)^2;

end

functiong=grad2(x)

res=norm(gk);

g(1,1)=2+400*(3*x(1)^2-x(2));

g(1,2)=-400*x(1);

g(2,1)=-400*x(1);

g(2,2)=200;

end

functiong=grad(x)

g=zeros(2,1);

g(1)=2*(x(1)-1)+400*x(1)*(x(1)^2-x(2));

g(2)=200*(x(2)-x(1)^2);

end

functionx_star=newton(x0,eps)

gk=grad(x0);

bk=[grad2(x0)]^(-1);

x0=[0,0];

k=0;

whileresepsk=1000

dk=-bk*gk;

xk=x0+dk;

k=k+1;

x0=xk;

gk=grad(xk);

bk=[grad2(xk)]^(-1);

res=norm(gk);

fprintf(--The%d-thiter,theresidualis%f\n,k,res);

end

x_star=xk;

end

clear

eps=1e-4;

x1=newton(x0,eps)

--The1-thiter,theresidualis

--The2-thiter,theresidualis

x1=

法：

g=zeros(2,1);

functionf=fun(x)

f=(1-x(1))^2+100*(x(2)-x(1)^2)^2;

end

functiong=grad(x)

g(1)=2*(x(1)-1)+400*x(1)*(x(1)^2-x(2));

g(2)=200*(x(2)-x(1)^2);

end

functionx_star=bfgs(x0,eps)

g0=grad(x0);

gk=g0;

res=norm(gk);

Hk=eye(2);

k=0;

whileresepsk=1000

dk=-Hk*gk;

ak=1;f0=fun(x0);

f1=fun(x0+ak*dk);

slope=dot(gk,dk);

whilef1f0+*ak*slope

x_star=xk;

ak=ak/4;

xk=x0+ak*dk;

f1=fun(xk);

end

k=k+1;

fa0=xk-x0;

x0=xk;

go=gk;

gk=grad(xk);

y0=gk-g0;

Hk=((eye(2)-fa0*(y0))/((fa0)*(y0)))*((eye(2)-

(y0)*(fa0))/((fa