线性规划模型方案.docVIP

下载本文档

6
0
约6.18千字
约 9页
2017-07-08 发布于湖北
举报
版权申诉

线性规划模型方案.doc

1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
7、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

线性规划模型方案

线性规划模型 4.1 奶制品的生产与销售[1] 例2 奶制品的生产销售计划（P88~92） % plan.m c = [-24 -16 -44 -32 3 3] A = [4 3 0 0 4 3; 4 2 0 0 6 4;1 0 0 0 1 0] b = [600; 480; 100] aeq = [0 0 1 0 -0.8 0;0 0 0 1 0 -0.75] beq = zeros(2,1) xLB = zeros(6,1) xUB = inf * ones(6,1) [x,fval] = linprog(c,A,b,aeq,beq,xLB,xUB) 非线性规划模型 12.1 供应与选址[2] （1）编写M文件liaoch.m定义目标函数 % liaoch.m function f=liaoch(x) a=[1.25 8.75 0.5 5.75 3 7.25]; b=[1.25 0.75 4.75 5 6.5 7.75]; d=[3 5 4 7 6 11]; e=[20 20]; f1=0; for i=1:6 s(i)=sqrt((x(13)-a(i))^2+(x(14)-b(i))^2); f1=s(i)*x(i)+f1; end f2=0; for i=7:12 s(i)=sqrt((x(15)-a(i-6))^2+(x(16)-b(i-6))^2); f2=s(i)*x(i)+f2; end f = f1 + f2; （2）工地分布及需求量示意图 a=[1.25 8.75 0.5 5.75 3 7.25]; b=[1.25 0.75 4.75 5 6.5 7.75]; scatter(a,b) （3）编写主程序xuanzhi.m % xuanzhi.m A=[1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 1 1 1 1 0 0 0 0]; b=[20;20]; Aeq=[eye(6) eye(6) zeros(6,4)]; beq=[3 5 4 7 6 11] VLB=[zeros(12,1);-inf;-inf;-inf;-inf]; x0=[3 0 4 5 4 0 0 5 0 2 2 11 3 4 7 6.5]; [x,fval,exitflag]=fmincon(@liaoch, x0, A, b, Aeq, beq,VLB) （4）结果为 x = Columns 1 through 6 3.0000 0 4.0000 7.0000 6.0000 0 Columns 7 through 12 0 5.0000 0 0.0000 0 11.0000 Columns 13 through 16 3.2549 5.6523 7.2500 7.7500 fval = 85.2660 exitflag = 5 统计回归模型 10.1牙膏的销售量[1] x1 = [-0.05 0.25 0.60 0 0.25 0.20 0.15 0.05 -0.15 0.15 0.20 0.10 0.40 0.45 0.35 0.30 0.50 0.50 0.40 -0.05 -0.05 -0.10 0.20 0.10 0.50 0.60 -0.05 0 0.05 0.55]; y = [7.38 8.51 9.52 7.50 9.33 8.28 8.75 7.87 7.10 8.00 7.89 8.15 9.10 8.86 8.90 8.87 9.26 9.00 8.75 7.95 7.65 7.27 8.00 8.50 8.75 9.21 8.27 7.67 7.93 9.26]; scatter(x1,y), title(图1 y对x1的散点图) x2 = [5.50 6.75 7.25 5.50 7.00 6.50 6.75 5.25 5.25 6.00 6.50 6.25 7.00 6.90 6.80 6.80 7.10 7.00 6.80 6.50 6.25 6.00 6.50 7.00 6.80 6.80 6.50 5.75 5.80 6.80]; scatter(x2,y), title(图2 y对x2的散点图) x = [ones(size(x1));x1;x2;x2.^2]; X = x.; Y = y.; [b,bint,r,rint,stats] = regress(Y,X,0.05) b = 17.3244