优化设功计机械工程陈宇.docxVIP

机械优化设计姓 名： 陈 宇学 号：152080200010 班 级：机械工程15班黄金分割:用0.618法求一元函数f（t）=t^2+2*t在区间[-3,5]中的极小点，要求计算到区间的长度小于0.05。function q = line_s(a,b)N=10000;r=0.005;for k=0:N; v=a+0.382*(b-a); u=a+0.618*(b-a); fv=v^2+2*v; fu=u^2+2*u; if fvfu if b-v=r q=u; break; else a=v;v=u; u=a+0.618*(b-a); end else if u-a=r q=v; break; else b=u;u=v; v=a+0.382*(b-a); end end k=k+1endreturn q=line_s(-3,5)maltlab 运行结果牛顿法用牛顿法求函数f = x1^2 - 2*x2*x1 + 3/2*x2^2 + x1 - 2x2的极小点。e = 1;x1 = 0;x2 = 0; while e0.00001 f_d = function_c1(x1,x2); f = subs(f_d, x1,x2, [x1,x2]); H = function_c2(x1,x2); H_inv = inv(H); x = [x1;x2] - H_inv * f; e = abs(x(1)-x1) / x(1); x1 = x(1); x2 = x(2); end disp(x)function [diff_x] = function_c1(x1,x2) syms x1 x2 f = x1^2 - 2*x2*x1 + 3/2 * x2^2 + x1 - 2*x2; diff_x1 = diff(f,x1);diff_x2 = diff(f,x2); diff_x = [diff_x1; diff_x2];returnendfunction [H] = function_c2(x1,x2) syms x1 x2 f = x1^2 - 2*x2*x1 + 3/2 * x2^2 + x1 - 2*x2; diff_x1 = diff(f,x1);diff_x2 = diff(f,x2); diff_x1_x1 = diff(diff_x1,x1);diff_x2_x2 = diff(diff_x2,x2); diff_x1_x2 = diff(diff_x1,x2);diff_x2_x1 = diff_x1_x2; H = [diff_x1_x1 diff_x1_x2;... diff_x2_x1 diff_x2_x2]; return endmaltlab 运行结果鲍威尔法 用Powell法求函数f=t^2+s^2-t*s-10*t-4*s+60的最优点X*=[t,s]T。计算收敛精度eps=0.0001。function [y]=Funval(f,~,x0) y = subs(f,findsym(f),x0);function [x,minf]=minHJ(f,a,b,eps)format long;if nargin==3 eps=1.0e-6;end l=a+0.382*(b-a);u=a+0.618*(b-a);k=1;tol=b-a; while toleps k100000 fl=subs(f,findsym(f),l); fu=subs(f,findsym(f),u); if flfu a=l; 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);endif k==100000 disp; x=NaN; minf=NaN; return;endx=(a+b)/2;minf=subs(f,findsym(f),x);minf = double(minf);format short;function[minx,maxx]=minJT(f,x0,h0,eps)format long;if nargin==3 eps=1.0e-6;end x1=x0; k=0; h=h0; while 1 x4=x1+h; k=k+1; f4=subs(f,findsym(f),x4); f1=subs(f,findsym(f),x1); if f4f1 x2=x1; x1=x4; f2=f1; f1=f4; h=2*h; else if k==1 h=-h; x2=x4; f2=f4; else x3=x2; x2=x1; x1=x4; break; end end end minx=min(x1