matalb最优化的实现程序.docVIP

第1方程求根 在命令窗口输入程序： x=fzero(inline(2*x+1),0) 结果： x = -0.5000 第2方程组求解 建立fun002. m文件 function y=fun002(x) y(1)=x(1)^2+x(2)^2-10; y(2)=x(1)^2-x(2)^2-8; 在命令窗口输入程序： x0=[1 1];x=fsolve(fun002,x0) 结果： x = 3.0000 1.0000 第3有界单变量优化 在命令窗口输入程序： v1=1;v2=5;f=inline(2*x-1);x=fminbnd(f,v1,v2) 结果： x = 1.0001 第4无约束极小（非线性规划） 建立fun004. m文件 function f=fun004(x) f=(x(1)^2-x(2)-1)^2+((x(1)-2)^2+(x(2)-1/2)^2-1)^2; 在命令窗口输入程序： x0=[0 0];x=fminunc(fun004,x0) 结果： x = 1.0673 0.1392 第5 约束极小（非线性规划） 建立目标函数fun0501.m文件： function f=fun0501(x) f=exp(x(1))*(4*x(1)^2+2*x(2)^2+4*x(1)*x(2)+1); 建立非线性约束fun0502.m文件： function [c,ceq]=fun0502(x) c=[x(1)*x(2)-x(1)-x(2)+1.5;-x(1)*x(2)-10]; ceq=x(1)^2+x(2)-1; 在命令窗口输入程序： x0=[1 -1];x=fmincon(fun0501,x0,[],[],[],[],[],[],fun0502) 结果： x = 1.2972 -0.6826 第7 二次规划 建立fun07. m文件： H=[4 -5;-5 2]; c=[-2 3]; A1=[-2 1;1 -3]; b1=[1 5]; A2=[1 2]; b2=4; v1=[1 -Inf]; v2=[Inf 0]; x=quadprog(H,c,A1,b1,A2,b2,v1,v2) 结果： x = 4.0000 0 第8 极小极大 建立fun0801.m文件 function f=fun0801(x) f(1)=x(1)^2-x(1)*x(2)+3; f(2)=x(1)+x(2)^2; f(3)=x(1)^2+x(2); f(4)=x(1)+x(2); 建立非线性约束fun0802.m文件 function [c,ceq]=fun0802(x) c=[x(1)*x(2)-x(1)-x(2)+1.5;-x(1)*x(2)-10]; ceq=x(1)^2+x(2)-1; 在命令窗口输入程序： x0=[1 1];x=fminimax(fun0801,x0,[],[],[],[],[],[],fun0802) 结果： x = -0.1250 1.3283 第9 非负线性最小二乘 在命令窗口输入程序： C=[1 -1];d=-1;x=lsqnonneg(C,d) 结果： x = 0 1 第10 约束线性最小二乘 在命令窗口输入程序： C=[1-1];d=-1;A1=[1 -2];b1=3;A2=[1 1];b2=5;x=lsqlin(C,d,A1,b1,A2,b2,[],[]) 结果： x = 2.0000 3.0000 第11非线性最小二乘 建立fun11. m文件 function f=fun11(x,t,c) f=x(1)*exp(x(2)*t)-c; 在命令窗口输入程序： x0=[10 0.5];t=[0.25 0.5 1 1.5 2 3 4 6 8];c=[19.21 18.15 15.36 14.10 12.89 9.32 7.45 5.24 3.01];x=lsqnonlin(fun11,x0,[],[],[],t,c) 结果： x = 20.2413 -0.2420 第12 非线性拟合 建立fun12. m文件 function f=fun12(x,t,c) f=x(1)*exp(x(2)*t); 在命??窗口输入程序： x0=[10 0.5];t=[0.25 0.5 1 1.5 2 3 4 6 8];c=[19.21 18.15 15.36 14.10 12.89 9.32 7.45 5.24 3.01];x=lsqcurvefit(fun12,x0,t,c) 结果： x = 20.2413 -0.2420 第13多目标规划 第14 0-1规划 在命令窗口输入程序： c=[1