矩形排料问题,组合优化问题 - 范文中心

矩形排料问题,组合优化问题

03/30

《二维矩形条带装箱问题的底部左齐择优匹配算法_蒋兴波》matlab的实现，不包括遗传算法部分。 function area =

PackingAlgorithm(length,width,length1,width1,length2,width2,length3,width3,restrict1,restrict2)

area = 0;

frameCount = 1;

count1 = 0;

count2 = 0;

runLLABF;

function runLLABF

rectBig.length = length;

rectBig.width = width;

rectSmall(1).length = length1;

rectSmall(1).width = width1;

rectSmall(1).color = 'r';

rectSmall(2).length = length2;

rectSmall(2).width = width2;

rectSmall(2).color = 'b';

rectSmall(3).length = length3;

rectSmall(3).width = width3;

rectSmall(3).color = 'g';

edges(1).x = 0;

edges(1).y = 0;

edges(1).length = rectBig.length;

edges(2).x = -100;

edges(2).y = 10000;

edges(2).length = 0;

edges(3).x = rectBig.length+100;

edges(3).y = 10000;

edges(3).length = 0;

while(1)

flag = -1;

if(flag

[sortedEdges,lowestEdge,id] = edgesSort(edges);

[edges,flag] = FullFitFirst(sortedEdges,lowestEdge,id,rectSmall); if(flag

[sortedEdges,lowestEdge,id] = edgesSort(edges);

[edges,flag] =

WidthFitFirst(sortedEdges,lowestEdge,id,rectSmall);

end

if(flag

[sortedEdges,lowestEdge,id] = edgesSort(edges);

[edges,flag] =

HeightFitFirst(sortedEdges,lowestEdge,id,rectSmall);

end

if(flag

[sortedEdges,lowestEdge,id] = edgesSort(edges);

[edges,flag] =

PlaceabelFirst(sortedEdges,lowestEdge,id,rectSmall);

end

if(flag

[sortedEdges,lowestEdge,id] = edgesSort(edges);

[edges,flag] =

cannotPalce(sortedEdges,lowestEdge,id,rectSmall,flag)

end

end

if count1 >= restrict1

rectSmall(1).length = 100000;

rectSmall(1).width = 100000;

end

if count2 >= restrict2

rectSmall(2).length = 100000;

rectSmall(2).width = 100000;

end

sortRect = sort([rectSmall(1).length,rectSmall(1).width,... rectSmall(2).length,rectSmall(2).width,...

rectSmall(3).length,rectSmall(3).width]);

minRect = sortRect(1);

minRect2 = sortRect(2);

[sortedEdges,lowestEdge,id] = edgesSort(edges);

[~,h] = size(sortedEdges);

for i = 1:h

if(sortedEdges(i).y+minRect

break;

end

end

if i == h

break;

end

if i == h-1 && lowestEdge.x + minRect2 > length

break

end

if frameCount > 300

break;

end

end

end

function initial

rectBig.length = 30;

rectBig.width = 20;

rectSmall(1).length = 4;

rectSmall(1).width = 3;

rectSmall(2).length = 3;

rectSmall(2).width = 3;

rectSmall(3).length = 4;

rectSmall(3).width = 1;

edges(1).x = 0;

edges(1).y = 0;

edges(1).length = rectBig.length;

% edges(1).x = 12;

% edges(1).y = 10;

% edges(1).length = 2;

%

% edges(2).x = 3;

% edges(2).y = 8;

% edges(2).length = 2;

%

% edges(3).x = 6;

% edges(3).y = 4;

% edges(3).length = 1;

%

% edges(4).x = 1;

% edges(4).y = 8;

% edges(4).length = 2;

end

function [sortedEdges,lowestEdge,id] = edgesSort(edges)

sortedEdges = edges;

[~,m] = size(sortedEdges);

for j = 1:m

for i = j:m

if(sortedEdges(i).x

tmpedge = sortedEdges(j);

sortedEdges(j) = sortedEdges(i);

sortedEdges(i) = tmpedge;

end

end

end

[~,m] = size(sortedEdges);

disp(m)

if(m>=2)

i = 2;

while(1)

if( sortedEdges(i-1).y == sortedEdges(i).y )

sortedEdges(i-1).length = sortedEdges(i-1).length + sortedEdges(i).length;

for j = i:(m-1)

sortedEdges(j) = sortedEdges(j+1);

end

sortedEdges(m) = [];

[~,n] = size(sortedEdges);

m = n;

continue;

end

[~,n] = size(sortedEdges);

m = n;

if i == n

break;

end

i = i+1;

end

else

lowestEdge = sortedEdges(1);

end

lowestEdges = sortedEdges;

[~,n] = size(lowestEdges);

y = lowestEdges(1).y;

for i = 2:n

if(y>lowestEdges(i).y)

y = lowestEdges(i).y;

end

end

for i = 1:n

if(lowestEdges(i).y == y )

lowestEdge = lowestEdges(i);

id = i;

break;

end

end

end

function [Edges,flag] = FullFitFirst(Edges,lEdge,lEdgeId,rectSmall) for i = 1:3

if (( rectSmall(i).length == lEdge.length )...

&& ( (lEdge.y+rectSmall(i).width == Edges(lEdgeId-1).y) ||

( lEdge.y+rectSmall(i).width == Edges(lEdgeId+1).y )))

if( lEdge.y+rectSmall(i).width

Edges(lEdgeId).y = lEdge.y+rectSmall(i).width;

flag = 1;

figurePlot(lEdge,rectSmall(i),rectSmall(i).color);

area = area+rectSmall(i).width*rectSmall(i).length;

if i == 1

count1 = count1+1;

end

if i == 2

count2 = count2+1;

end

break;

else

flag = -1;

end

else

flag = -1;

end

if (( rectSmall(i).width == lEdge.length )...

&& ( (lEdge.y+rectSmall(i).length == Edges(lEdgeId-1).y) ||

( lEdge.y+rectSmall(i).length == Edges(lEdgeId+1).y )))

if( lEdge.y+rectSmall(i).length

Edges(lEdgeId).y = lEdge.y+rectSmall(i).length;

flag = 1;

figurePlotRotation(lEdge,rectSmall(i),rectSmall(i).color); area = area+rectSmall(i).width*rectSmall(i).length;

if i == 1

count1 = count1+1;

end

if i == 2

count2 = count2+1;

end

break;

else

flag = -1;

end

else

flag = -1;

end

end

end

function [Edges,flag] = WidthFitFirst(Edges,lEdge,lEdgeId,rectSmall) count = 1;

% selected = zeros(1,3);

for i = 1:3

if ( rectSmall(i).length == lEdge.length && rectSmall(i).width + lEdge.y

selected(count).index = i;

selected(count).area = rectSmall(i).length * rectSmall(i).width; selected(count).rotation = 0;

count = count + 1;

flag = 1;

continue;

else

flag = -1;

end

if ( rectSmall(i).width == lEdge.length && rectSmall(i).length + lEdge.y

selected(count).index = i;

selected(count).area = rectSmall(i).length * rectSmall(i).width; selected(count).rotation = 1;

count = count + 1;

flag = 1;

else

flag = -1;

end

end

if(flag == 1)

[~,n] = size(selected);

for i =1:n

for j = i:n

if(selected(i).area>selected(j).area)

tmpSelected = selected(i);

selected(i) = selected(j);

selected(j) = tmpSelected;

end

end

end

index = selected(n).index;

if(selected(n).rotation == 0)

if lEdge.y+rectSmall(index).width

Edges(lEdgeId).y = lEdge.y+rectSmall(index).width;

flag = 1;

figurePlot(lEdge,rectSmall(index),rectSmall(index).color); area = area+rectSmall(index).width*rectSmall(index).length; if index == 1

count1 = count1+1;

end

if index == 2

count2 = count2+1;

end

else

flag = -1;

end

end

if(selected(n).rotation == 1)

if lEdge.y+rectSmall(index).length

Edges(lEdgeId).y = lEdge.y+rectSmall(index).length; flag = 1;

figurePlotRotation(lEdge,rectSmall(index),rectSmall(index).color); area =

area+rectSmall(index).width*rectSmall(index).length;

if index == 1

count1 = count1+1;

end

if index == 2

count2 = count2+1;

end

else

flag = -1;

end

end

end

end

function [Edges,flag] = HeightFitFirst(Edges,lEdge,lEdgeId,rectSmall)

[~,n] = size(Edges);

for i = 1:3

if ( rectSmall(i).length

Edges(n+1).x = Edges(lEdgeId).x+rectSmall(i).length; Edges(n+1).y = Edges(lEdgeId).y;

Edges(n+1).length =

Edges(lEdgeId).length-rectSmall(i).length;

Edges(lEdgeId).y = lEdge.y+rectSmall(i).width;

Edges(lEdgeId).length = rectSmall(i).length;

flag = 1;

figurePlot(lEdge,rectSmall(i),rectSmall(i).color);

area = area+rectSmall(i).width*rectSmall(i).length; if i == 1

count1 = count1+1;

end

if i == 2

count2 = count2+1;

end

else

flag = -1;

end

if(flag == 1)

break;

end

if( rectSmall(i).width

Edges(n+1).x = Edges(lEdgeId).x+rectSmall(i).width; Edges(n+1).y = Edges(lEdgeId).y;

Edges(n+1).length =

Edges(lEdgeId).length-rectSmall(i).width;

Edges(lEdgeId).y = lEdge.y+rectSmall(i).length;

Edges(lEdgeId).length = rectSmall(i).width;

flag = 1;

figurePlotRotation(lEdge,rectSmall(i),rectSmall(i).color); area = area+rectSmall(i).width*rectSmall(i).length; if i == 1

count1 = count1+1;

end

if i == 2

count2 = count2+1;

end

else

flag = -1;

end

if(flag == 1)

break;

end

end

end

function [Edges,flag] = PlaceabelFirst(Edges,lEdge,lEdgeId,rectSmall) count = 1;

[~,m] = size(Edges);

selected(1).index = 1;

selected(1).area = 0;

selected(1).rotation = 0;

selected(2).index = 1;

selected(2).area = 0;

selected(2).rotation = 0;

selected(3).index = 1;

selected(3).area = 0;

selected(3).rotation = 0;

for i = 1:3

if ( rectSmall(i).length

selected(count).index = i;

selected(count).area = rectSmall(i).length * rectSmall(i).width; selected(count).rotation = 0;

count = count + 1;

flag = 1;

continue;

else

flag = -1;

end

if ( rectSmall(i).width

selected(count).index = i;

selected(count).area = rectSmall(i).length * rectSmall(i).width; selected(count).rotation = 1;

count = count + 1;

flag = 1;

else

flag = -1;

end

end

if flag == 1

n = count -1;

for i =1:n

for j = i:n

if(selected(i).area>selected(j).area)

tmpSelected = selected(i);

selected(i) = selected(j);

selected(j) = tmpSelected;

end

end

end

index = selected(n).index;

if(selected(n).rotation == 0)

Edges(m+1).x = lEdge.x+rectSmall(index).length;

Edges(m+1).y = lEdge.y;

Edges(m+1).length = lEdge.length-rectSmall(index).length; Edges(lEdgeId).y = lEdge.y+rectSmall(index).width;

Edges(lEdgeId).length = rectSmall(index).length;

figurePlot(lEdge,rectSmall(index),rectSmall(index).color); area = area+rectSmall(index).width*rectSmall(index).length; if index == 1

count1 = count1+1;

end

if index == 2

count2 = count2+1;

end

end

if(selected(n).rotation == 1)

Edges(m+1).x = lEdge.x+rectSmall(index).width;

Edges(m+1).y = lEdge.y;

Edges(m+1).length = lEdge.length-rectSmall(index).width; Edges(lEdgeId).y = lEdge.y+rectSmall(index).length;

Edges(lEdgeId).length = rectSmall(index).width;

figurePlotRotation(lEdge,rectSmall(index),rectSmall(index).color); area = area+rectSmall(index).width*rectSmall(index).length; if index == 1

count1 = count1+1;

end

if index == 2

count2 = count2+1;

end

end

end

end

function [Edges,flag] = cannotPalce(Edges,lEdge,lEdgeId,rectSmall,flag2) count = 0;

for i = 1:3

if (rectSmall(i).width > lEdge.length) && (rectSmall(i).length > lEdge.length) || (flag2 == -1)

count = count + 1;

end

end

if count == 3

flag = 1;

if Edges(lEdgeId-1).y

Edges(lEdgeId).y = Edges(lEdgeId-1).y;

else

Edges(lEdgeId).y = Edges(lEdgeId+1).y;

end

end

flag = 1;

end

function figurePlot(lEdge,rect,color)

x1 = lEdge.x;

y1 = lEdge.y;

x2 = x1+rect.length;

y2 = y1;

x3 = x2;

y3 = y2 + rect.width;

x4 = x1;

y4 = y3;

x = [x1,x2,x3,x4];

y = [y1,y2,y3,y4];

patch(x,y,color,'facealpha',0.5);

axis([-1 length+1 -1 width+1])

% set(gca,'XTick',0:2:length);

% set(gca,'YTick',0:2:width);

% m(frameCount) = getframe;

m = getframe(gcf);

im = frame2im(m);

[I,map] = rgb2ind(im,256);

if frameCount == 1

imwrite(I,map,'out.gif','gif','loopcount',inf,'Delaytime',0.2) else

imwrite(I,map,'out.gif','gif','writemode','append','Delaytime',0.2) end

frameCount = frameCount+1;

% hold on

% pause(0.1);

axis equal;

end

function figurePlotRotation(lEdge,rect,color)

x1 = lEdge.x;

y1 = lEdge.y;

x2 = x1+rect.width;

y2 = y1;

x3 = x2;

y3 = y2 + rect.length;

x4 = x1;

y4 = y3;

axis([-1 length+1 -1 width+1])

x = [x1,x2,x3,x4];

y = [y1,y2,y3,y4];

patch(x,y,color,'facealpha',0.5);

% set(gca,'XTick',0:2:length);

% set(gca,'YTick',0:2:width);

% m(frameCount) = getframe;

m = getframe(gcf);

im = frame2im(m);

[I,map] = rgb2ind(im,256);

if frameCount == 1

imwrite(I,map,'out.gif','gif','loopcount',inf,'Delaytime',0.2) else

imwrite(I,map,'out.gif','gif','writemode','append','Delaytime',0.2) end

frameCount = frameCount+1;

% hold on

% pause(0.1);

axis equal;

end

end


相关内容

  • 假期走访一
    我很庆幸能在放寒假的当天跟学长.学姐们一起去苏州的一个手机厂工作.在手机厂刚开始做清屏工作,后来转成看外观.刚进无尘室的第一感觉觉得一个小小的手机屏一定不会难倒我.可是并没想象的简单.进到无尘室有一个小师傅带我,她说清屏需要很大的耐心,我说 ...
  • 上海市黄浦江卢浦大桥设计.pdf
    第 38 卷第 1 期 2 0 0 5 年 1 月 土 木 工 程 学 报 CHINA CIVIL ENGINEERING JOURNAL Vol138 No11 Jan1 2 0 0 5 上海市黄浦江卢浦大桥设计 林元培 章曾焕 马 周 ...
  • 整体包装设计系统
    整体包装设计系统 ✓ 软件版本:V3.2 ✓ 软件功能介绍: 1. 木包装制品设计 采用目前最新的三位设计技术,自动生成包装产品模型 采用有限元仿真分析技术,对构件进行强度分析.包装结构设计优化,从设计角度降低包装成本 快速生成包装设计的总 ...
  • 均匀分布的随机数
    随机数的产生 摘要 本文研究了连续型随机数列的产生,先给出了均匀分布的随机数的产生算法,在通过均匀分布的随机数变换得到其他连续型随机数的产生算法.在vc环境下,我们给出了产生均匀分布随机数的算法,然后探讨了同余法的理论原理.通过均匀随机数产 ...
  • 服装QC岗位职责
    qc岗位责任制 岗位名称:qc检验员 报告对象:qc检验室负责人 岗位目标: 1 完成分管范围的产品检验工作,保证检验结果的准确.真实. 2 遵守公司质量管理方面及检验室的各项规定,保证检验工作的安全及检验室的清洁. 岗位职责: 1 qc检 ...
  • 大跨度建筑结构体系分析探讨
    1 前言 为了满足实际生产建设的需要和人们的审美观念需求,各色各样的建筑类型不断出现,其结构类型也多种多样.随着科学技术的不断发展,一些新技术应用的到建筑结构上,改造者传统的大跨度建筑结构,使大跨度的建筑结构越来越具有良好的抗压.抗震.稳定 ...
  • 矩形断面明渠均匀流水力计算的直接计算公式
    第36卷 第9期2008年9月西北农林科技大学学报(自然科学版) Jour nal of N or thwest A &F U niversity(Nat. Sci. Ed. ) V ol. 36No. 9Sep. 2008 矩形断 ...
  • 一种电子式电流互感器的研制
    # xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 一种电子式电流互感器的研制 申 烛!王士敏!罗承沐 清华大学电机系!北京市#"$$$%&a ...
  • 桥梁设计理论第五讲
    第五讲 薄壁箱梁的自由扭转 第一节 基本假定 在材料力学中,我们曾经讨论过圆截面杆的扭转问题,那时,我们假定杆件变形后截面保持为平面,只是相对地转动了一个角度,而截面的大小和形状都保持不变.这个假定对于圆截面杆来说是比较符合实际情况的,那么 ...
  • 通风管道加工制作流程
    通风管道的作用就是为我们工作生活的大楼提供清新空气的流通.宁夏凯润沣洋通风管道公司的工程师指出,由于是直接和空气接触,空气中含有大量的水分会对金属材质的通风管道造成腐蚀,所以在风管加工中需要注意选好制作材料.此外风管加工过程中对尺寸也有严格 ...

热门内容