钢管订购和运输问题一代码和结果.docxVIP

下载本文档

60
0
约3.04千字
约 7页
2017-09-28 发布于重庆
举报
版权申诉

钢管订购和运输问题一代码和结果.docx

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

钢管订购和运输问题一代码和结果

function f=result(t)%求解问题1tic;x0=zeros(8,15);vlb=zeros(8,15);m=zeros(1,7);s=[800 800 1000 2000 2000 2000 3000];s(t)=s(t)-50;N=[1 1 1 0 1 1 0]; %每公里钢管从Si到达Ai站点的最小费用 C=[330.7 320.3000 300.2000 258.6000 198.0000 180.5000 163.1000 181.2000 224.2000 252.0000 256.0000 266.0000 281.2000 288.0000 302.0000; 370.7 360.3000 345.2000 326.6000 266.0000 249.6000 241.0000 226.2000 269.2000 297.0000 301.0000 311.0000 326.2000 333.0000 347.0000; 385.7 375.3000 355.2000 336.6000 276.0000 260.5000 251.0000 241.2000 203.2000 237.0000 241.0000 251.0000 266.2000 273.0000 287.0000; 420.7 410.3000 395.2000 376.6000 316.0000 299.6000 291.0000 276.2000 244.2000 222.0000 211.0000 221.0000 236.2000 243.0000 257.0000; 410.7 400.3000 380.2000 361.6000 301.0000 285.5000 276.0000 266.2000 234.2000 212.0000 188.0000 206.0000 226.2000 228.0000 242.0000; 415.7 405.3000 385.2000 366.6000 306.0000 290.5000 281.0000 271.2000 234.2000 212.0000 201.0000 195.0000 176.2000 161.0000 178.0000; 435.7 425.3000 405.2000 386.6000 326.0000 310.5000 301.0000 291.2000 259.2000 236.0000 226.0000 216.0000 198.2000 186.0000 162.0000];options=optimset(LargeScale,off,Algorithm ,active-set,MaxFunEvals ,50000);%,Tolx,1.0000e-032);[x,f]=fmincon(myfun,x0,[],[],[],[],vlb,[],mycon,options,C,N,s);for i=1:7for j=1:15 m(i)=m(i)+N(i)*x(i,j);endendx,m,f;b=(f-1278600)/1278600*(s(t)+50)/50tocfunction f=myfun(XX,C,N,s)%问题1的目标函数x=XX(1:7,1:15);rl=XX(8,1:15);L=[104 301 750 606 194 205 201 680 480 300 220 210 420 500];f=0;for i=1:7for j=1:15 f=f+N(i)*x(i,j)*C(i,j);%运输费和成本费endendfor i=1:14 f=f+(rl(i)*(rl(i)+1)/2+(L(i)-rl(i))*(L(i)-rl(i)+1)/2)*0.1;%铺设时的运输费endffunction [c,ceq]=mycon(XX,C,N,s)%问题1的约束条件x=XX(1:7,1:15);rl=XX(8,1:15);L=[104 301 750 606 194 205 201 680 480 300 220 210 420 500];m=zeros(1,7);a=zeros(1,15);cc=0;for i=1:7for j=1:15 m(i)=m(i)+N(i)*x(i,j);end c(i)