非线性规划教案.ppt

1.无约束优化问题的Matlab解法 求函数　　　　　　　　　　　　　 的极小值。 Matlab程序如下： function f=fun(x) 　f=100*(x(2)-x(1))^2+(1-x(1))^2; 　函数名Rosenbroken.m x0=[-1.9,2]; opt(6)=1; [x,opt]=fminu(fun,x0,opt); f=opt(8) 运行结果 x(1)=0.9999 x(2)= 0.9999 F(x)= 7.7199e-016 fmincon应用求解示例： 请问： 1、结合fmincon函数，需要提供哪些参数 第一步：编写一个M文件返回目标函数f在点x处的值函数程序 function f = myfun(x) f = -x(1) * x(2) * x(3); 第二步：为了调用MATLAB函数，必须将模型中的约束转化为如下形式(=)。 这里： A=[-1 -2 -2; 1 2 2 ]; b=[0 72]’; 第三步：提供一个搜索起点，然后调用相应函数，程序如下： % 给一个初始搜索点 x0 = [10; 10; 10]; [x,fval] = fmincon(myfun,x0,A,b) 主程序（整体）： A=[-1 -2 -2; 1 2 2 ]; b=[0 72]’; % 给一个初始搜索点 x0 = [10; 10; 10]; [x,fval] = fmincon(myfun,x0,A,b) 最后得到如下结果：? x = 24.0000 12.0000 12.0000 ? fval = -3.4560e+03 某装饰材料公司欲以每桶2元的价钱购进一批彩漆 一般来说随着彩漆售价的提高，预期销售量将减少，并对此进行了估算，见表1。为了尽快收回资金并获得较多的赢利，装饰材料公司打算做广告，投入一定的广告费后，销售量将有一个增长，可由销售增长因子来表示。根据经验，广告费与销售增长因子关系见 表2。现在的问题是装饰材料公司采取怎样的营销 战略预期的利润最大？ 表1 表2 符号说明及问题的分析 设x表示售价（单位：元），y表示预期销售量(单位：万桶），z表示广告费（单位：万元），k表示销售增长因子。投入广告费后，实际销售量记为s 获得的利润记为P（单位：元）。由表1易见预期 销售量 y 随着售价x 的增加而单调下降，而销售增长因子k在开始时随着广告费z的增加而增加，在广告费z等于50000元时达到最大值，然后在广告费增加时反而有所回落，为此可用Matlab画出散点图. 文件名：graph1.m x=(2.0:0.5:6.0); y=[4.1,3.8,3.4,3.2,2.9,2.8,2.5,2.2,2.0]; plot(x,y,r*) title(售价和预期销售量关系图) xlabel(售价（元）) ylabel(预期销售量（万桶）) 文件名：graph2.m x=(0:1:7); y=[1.0,1.4,1.7,1.85,1.95,2.0,1.95,1.80]; plot(x,y,rp) title(广告费和销售增长因子关系图) xlabel(广告费（万元）) ylabel(销售增长因子) ? 图-２ ? 从图1和图2易见，售价x与预期销售量y近似于 一条直线，广告费 z 与销售增长因子k近似于一条 二次曲线。为此可令： y=a+bx k=c+dz+ez2 系数a，b，c，d，e是待定参数。 模型的建立 投入广告费后，实际销售量s等于预期销售量y乘 以销售增长因子k，即s=ky。所获得的利润：  ? 由于目标函数不是线性函数，因此这一问题的数学 模型为有约束条件的非线性规划模型。在日常生活 中非线性规划问题要比线性规划问题普遍。 模型求解 首先利用Matlab计算（1）（2）中的参数a，b，c，d，e，并画出散点图和拟合曲线。 文件名：yihe1.m x=(2.0:0.5:6.0); y=[4.1,3.8,3.4,3.2,2.9,2.