2019年五一杯数学建模B题木板最优切割方案解题全过程文档及程序

2019年五一杯数学建模

B题 木板最优切割方案

原题再现

  徐州某家具厂新进一批木板如表 1 所示,在家具加工的过程中,需要使用切割工具生产表 2所示的产品。假设:木板厚度和割缝宽度忽略不计。
2019年五一杯数学建模B题木板最优切割方案解题全过程文档及程序_第1张图片
  请为该家具厂给出如下问题的木板最优切割方案。
  1. 在一块木板上切割 P1 产品,建立数学模型,给出木板利用率最高(即剩余木板面积最小)的切割方案,并将最优方案的结果填入表 3。
2019年五一杯数学建模B题木板最优切割方案解题全过程文档及程序_第2张图片
  2. 在一块木板上切割 P1 和 P3 产品,建立数学模型,给出按照木板利用率由高到低排序的前 3 种切割方案,并将结果填入表 4。
2019年五一杯数学建模B题木板最优切割方案解题全过程文档及程序_第3张图片
  3. 需要完成表 2 中 P1 和 P3 产品的生产任务,建立数学模型,给出木板总利用率最高的切割方案,并将结果填入表 5。
2019年五一杯数学建模B题木板最优切割方案解题全过程文档及程序_第4张图片
  4. 需要完成表 2 中 P1、P2、P3、P4 产品的生产任务,建立数学模型,给出木板总利用率最高的切割方案,并将结果填入表 6。
2019年五一杯数学建模B题木板最优切割方案解题全过程文档及程序_第5张图片
  5. 不考虑产品P1,P2,P3,P4的需求数量,给定100张S1木板,按照表2中给出的利润,建立数学模型,给出总利润最大的切割方案,并将结果填入表7。
2019年五一杯数学建模B题木板最优切割方案解题全过程文档及程序_第6张图片

整体求解过程概述(摘要)

  本文主要研究了在给定条件下如何对木板最优切割问题,重点研究分析了不同情况下的切割方案,通过建立线性规划模型和混合整数规划模型,并利用 MATLAB和 LINGO 对提出的模型进行求解,具体结果如下:
  对于问题一:考虑到只切割 P1 产品,若想木板利用率最高,只需切割产品的总面积最大,根据两个产品不能重叠,引入了坐标的思想,建立线性规划模型,结合 MATLAB 里的 fmincon 函数求解,得出:在只切割 P1 产品时,一块木板上能够切割 59 块 P1 产品,最高木板利用率为 98.30%。
  对于问题二:在问题一模型的基础上,添加变量表示产品种类,建立新的线性规划模型,结合遗传算法求解,在两种产品数量可以为 0 的情况下,木板利用率由高到底排序的前三种切割方案为:
  方案1:在木板全部切割P3产品,能够切割48 块P3,此时木板利用率为99.17%;
  方案2:在木板上切割47块P3产品和 1块P1 产品,此时木板利用率为98.77%;
  方案3:在木板上切割46块P3产品和 2块P1 产品,此时木板利用率为98.37%。
  对于问题三:根据问题一、二中提出的线性规划模型,对问题 3 分析并建立新的整数规划模型,得到最优切割方案如下:用 33 个木板全部切割 P3,得到 1584 个P3,剩余 39 个,利用率为 99.17%;用 13 个木板全部切割 P1,得到 767 个 P1,剩余 7 个 P1,利用率为 98.30%;再将剩余的 7 个 P1 和 39 个 P3 用 1 个木板切割,利用率为 92.23%;所需木板总数量为 47,木板总利用率为 98.78%。
  对于问题四:首先考虑切割面积较大的产品,每次切割按面积大小依次切割,由于选取方案较多,只枚举部分切割方案,以切割的木板数量最少为目标函数,建立混合整数规划模型。运用 LINGO 求解,得到结果:需要木板 S1 的总数量为 138,木板总利用率为 98.47%,具体方案见表 6。
  对于问题五:同问题 4 的思想相同,在分割时,首先考虑尽可能的切割单位面积利润较大的产品,即首先切割 P1 产品,建立线性规划模型,通过 LINGO 求解,得到结果:在 100 块木板上只切割 P1,总数量为 5900,其余产品数量均为 0,总利润为 117410 元,木板总利用率为 98.30%
  最后,深入研究了论文中提出的模型并加以分析,对模型的优化提供了一些改进思路。

模型假设:

  由于产品要求,需要对木板进行一定的打磨去除粗边的加工,切割后的木板也需要进行一定的加工处理才可投入使用。切割时需要考虑切割损失等影响因素,因此,对于本题,我们做出如下假设:
  (1) 不考虑切割时木板产生的损耗。
  (2) 木板可以沿着任意方向切割。
  (3) 木板可以任意摆放。
  (4) 用材不能重叠,不超过原材料的大小。
  (5) 原板材密度厚度均匀,形状为规则矩形。
  (6) 木板厚度和割缝宽度忽略不计。

问题分析:

  此问题为板材下料问题,只切割一种产品时,切割方案比较容易,对于切割多种产品时,需要尽可能的求解板材切割方案,确定不同的方案使得木板得到更多的利用。
  对于问题 1,求解出只切割 P1 产品的木板利用率最高的切割方案,若想木板利用率最高,只需切割产品的总面积最大,即在只切割同一种产品的情况下,只需切割产品的数量最多,问题可以转化为:经过对木板切割,使得获得的 P1 产品数量最多,建立线性规划模型,结合 MATLAB 里的 fmincon 函数求解。
  对于问题 2,求解出在一块木板上切割 P1 和 P3 产品,木板利用率由高到低排序前 3 的切割方案,对问题进行分析,可以将此问题转化为:在 P1 和 P3 产品个数可以为 0 的情况下,在木板上切割 P1 和 P3 产品,可将问题转化为两种产品的总面积之和最大的问题。同问题 1 一样,用线性规划模型,结合 MATLAB 里的 fmincon函数求解。
  对于问题 3,要求出木板总利用率最高的切割方案,经过对问题 3 的分析,可将问题转化为所需木板数量最少的问题,即两种木板的总面积除以所需木板的总面积最大。同问题 1、2 一样,用线性规划建立模型,结合 MATLAB 里的 fmincon 函数求解。
  对于问题 4,在完成 P1,P2,P3,P4 产品的生产任务的条件下,找出木板总利用率最高的切割方案。通过建立相关模型,我们首先找出利用率较高的 12 种切割方案,再建立线性规划模型对这 12 种切割方案进行优化组合,通过 LINGO 函数进行求解。
  对于问题 5,在给定 100 张 S1 木板,不考虑产品 P1,P2,P3,P4 的需求数量的条件下,找出总利润最大切割方案,经过对问题 5 的分析,可以将其转化为找出单张木板的切割利润最大的问题。类似问题 4,建立线性规划模型,通过 LINGO 函数进行求解。

模型的建立与求解整体论文缩略图

2019年五一杯数学建模B题木板最优切割方案解题全过程文档及程序_第7张图片

全部论文请见下方“ 只会建模 QQ名片” 点击QQ名片即可

程序代码:(代码和文档not free)

clc,clear
x1=0;
x2=0; 
n2=0; 
n=0;
for n1=0:8              % 长为3000的边全部放 横放 最多放 8  块
%     n2=0;             % 若将 n2 定义成局部变量 后面的 n2 清零操作可以不需要,因为每次循环结束 局部变量的值自动释放
    while n2*201+n1*373<=3000 
        n2=n2+1;
    end
    if n2*201+n1*373>3000 
        n2=n2-1;        % 最后一次循环n2的值使 n2*201+n1*373>3000  超出 边界30 ,所以 n2 需要减 1
    end
    n=n2*4+n1*7;    % 竖放 最多放 4 行, 横放 最多放 7q=((n2*4+n1*7)*(373*201))/(3000*1500);
    fprintf('竖放%2d 列,横放%2d 列,总数%d,利用率%f\n',n2,n1,n,q); 
%     n=0;        % 清零 进行下次循环   重新计算下次循环总数
    n2=0;       % 清零 进行下次循环   为 "竖放" 计数
end
clc;clear; x=3000; y=1500;
x1=0; 
x2=0; 
i=0;
%同行不同列
fprintf('横放\n');
fprintf('(1)同行不同列\n');
for x1=1:7
    while x1*373+x2*406<=3000       % P1、P3按最长边横放在长为 3000 的边上  且 P1 、P3至少横放一个
        x2=x2+1;
    end
    while x1*373+x2*406>3000 
        x2=x2-1;
    end
    m1=x1*7;    % P1 最多放 715/2.01=7.4627
    m2=x2*6;    % P3 最多放 615/2.29=6.5502
    q=((m1*373*201)+(m2*406*229))/(3000*1500);
    fprintf('P1:%2d 列 ,	P1:%2d 个 ,	P3:%2d 列 ,	P3:%2d 个 , 利 用率%f\n',x1,m1,x2,m2,q);
    x2=0;
end
 
%同列不同行
y2=0;
fprintf('(2)同列不同行\n'); 
for y1=1:7
    while 201*y1+229*y2<=1500       % P1、P3 最短边 放在 长 1500的边
        y2=y2+1;
    end
    while 201*y1+229*y2>1500 
        y2=y2-1;
    end
    m1=y1*8; 
    m2=y2*7;
    q=((m1*373*201)+(m2*406*229))/(3000*1500);
    fprintf('P1:%2d 行 ,	P1:%2d 个 ,	P3:%2d 行 ,	P3:%2d 个 , 利 用率%f\n',y1,m1,y2,m2,q);
    y2=0;
end

fprintf('竖放\n');
fprintf('(3)同行不同列\n'); 
x2=0;
for x1=1:14
    while x1*201+x2*229<=3000       % P1、P3 最短边 放在 长 3000的边
        x2=x2+1;
    end
    while x1*201+x2*229>3000 
        x2=x2-1;
    end
    m1=x1*4; 
    m2=x2*3;
    q=((m1*373*201)+(m2*406*229))/(3000*1500);
    fprintf('P1:%2d 列 ,	P1:%2d 个 ,	P3:%2d 列 ,	P3:%2d 个 , 利 用率%f\n',x1,m1,x2,m2,q);
    x2=0;
end

fprintf('(4)同列不同行\n'); 
y2=0;
for y1=1:3
    while y1*373+y2*406<=1500   % P1、P3最长边放在  长1500的边
        y2=y2+1;
    end
    while y1*373+y2*406>1500
        y2=y2-1;
    end
    m1=y1*14; 
    m2=y2*13;
    q=((m1*373*201)+(m2*406*229))/(3000*1500);
    fprintf('P1:%2d 行 ,	P1:%2d 个 ,	P3:%2d 行 ,	P3:%2d 个 , 利 用率%f\n',y1,m1,y2,m2,q);
    y2=0;
end

%交叉
fprintf('(5)p1 竖放,p3 横放\n');
fprintf('同行不同列\n'); 
x2=0;
y2=0;
for x1=1:13
    while x1*201+x2*406<=3000       % P1最短边、P3 最长边 放在 长 3000的边
        x2=x2+1;
    end
    while x1*201+x2*406>3000 
        x2=x2-1;
    end
    m1=x1*4; 
    m2=x2*6;
    q=((m1*373*201)+(m2*406*229))/(3000*1500);
    fprintf('P1:%2d 列 ,	P1:%2d 个 ,	P3:%2d 列 ,	P3:%2d 个 , 利 用率%f\n',x1,m1,x2,m2,q);
    x2=0;
end
y2=0;
fprintf('(6)p1 竖放,p3 横放,不同行\n'); 
for y1=1:4
    while y1*373+y2*229<=1500       % P1最长边、P3 最短边 放在 长 1500的边
        y2=y2+1;
    end
    while y1*373+y2*229>1500 
        y2=y2-1;
    end
    m1=y1*14; 
    m2=y2*7;
    q=((m1*373*201)+(m2*406*229))/(3000*1500);
    fprintf('P1:%2d 列 ,	P1:%2d 个 ,	P3:%2d 列 ,	P3:%2d 个 , 利 用率%f\n',y1,m1,y2,m2,q);
    x2=0;
end
 
fprintf('(7)p1 横放,p3 竖放,同行\n'); 
for x1=1:7
    while x1*373+x2*229<=3000       % P1最长边、P3 最短边 放在 长 3000的边
        x2=x2+1;
    end
    while x1*373+x2*229>3000 
        x2=x2-1;
    end
    m1=x1*7; 
    m2=x2*3;
    q=((m1*373*201)+(m2*406*229))/(3000*1500);
    fprintf('P1:%2d 列 ,	P1:%2d 个 ,	P3:%2d 列 ,	P3:%2d 个 , 利 用率%f\n',x1,m1,x2,m2,q);
    x2=0;
end

y2=0;
fprintf('(8)p1 横放,p3 竖放,不同行\n'); 
for y1=1:5
    while y1*201+y2*406<=1500       % P1最短边、P3 最长边 放在 长 1500的边
        y2=y2+1;
    end
    while y1*201+y2*406>1500 
        y2=y2-1;
    end
    m1=y1*8; 
    m2=y2*13;
    q=((m1*373*201)+(m2*406*229))/(3000*1500);
    fprintf('P1:%2d 列 ,	P1:%2d 个 ,	P3:%2d 列 ,	P3:%2d 个 , 利 用率%f\n',y1,m1,y2,m2,q);
    x2=0;
end
clc;clear; 
x=3000; 
y=1500;
x1=0; 
x2=0; 
i=0; 
s=0;
qq=0;

%同行不同列
 
fprintf('横放\n');
fprintf('(1)同行不同列\n'); 
for x1=1:6
    while x1*373+x2*406<=3000 
        x2=x2+1;
    end
    while x1*373+x2*406>3000 
        x2=x2-1;
    end
    m1=x1*7; 
    m2=x2*6;
    q=((m1*373*201)+(m2*406*229))/(3000*1500);
    fprintf('P1:%2d 列 ,	P1:%2d 个 ,	P3:%2d 列 ,	P3:%2d 个 , 利 用率%f',x1,m1,x2,m2,q);
    while s*m1<=774||s*m2<=1623 
        s=s+1;
    end
    qq=(373*201*774+406*229*1623)/(s*3000*1500);
    fprintf('\tS1 的块数%d,总的利用率%f\n',s,qq); 
    qq=0;
    s=0; 
    x2=0;
end
 

%同列不同行
y2=0;
fprintf('(2)同列不同行\n'); 
for y1=1:6
    while 201*y1+229*y2<=1500 
        y2=y2+1;
    end
    while 201*y1+229*y2>1500 
        y2=y2-1;
    end
    m1=y1*8; 
    m2=y2*7;
    q=((m1*373*201)+(m2*406*229))/(3000*1500);
    fprintf('P1:%2d 行 ,	P1:%2d 个 ,	P3:%2d 行 ,	P3:%2d 个 , 利 用率%f',y1,m1,y2,m2,q);
    y2=0;
    while s*m1<=774||s*m2<=1623 
        s=s+1;
    end
    qq=(373*201*774+406*229*1623)/(s*3000*1500);
    fprintf('\tS1 的块数%d,总的利用率%f\n',s,qq); 
    qq=0;
    s=0;
end
fprintf('竖放\n');
fprintf('(3)同行不同列\n'); 
x2=0;
for x1=1:13
    while x1*201+x2*229<=3000 
        x2=x2+1;
    end
    while x1*201+x2*229>3000 
        x2=x2-1;
    end
    m1=x1*4; 
    m2=x2*3;
    q=((m1*373*201)+(m2*406*229))/(3000*1500);
    fprintf('P1:%2d 列 ,	P1:%2d 个 ,	P3:%2d 列 ,	P3:%2d 个 , 利 用率%f',x1,m1,x2,m2,q);
    x2=0;
    while s*m1<=774||s*m2<=1623 
        s=s+1;
    end
    qq=(373*201*774+406*229*1623)/(s*3000*1500);
    fprintf('\tS1 的块数%d,总的利用率%f\n',s,qq); 
    qq=0;
    s=0;
end

fprintf('(4)同列不同行\n'); 
y2=0;
for y1=1:2
    while y1*373+y2*406<=1500 
        y2=y2+1;
    end
    while y1*373+y2*406>1500 
        y2=y2-1;
    end
    m1=y1*14; 
    m2=y2*13;
    q=((m1*373*201)+(m2*406*229))/(3000*1500);
    fprintf('P1:%2d 行 ,	P1:%2d 个 ,	P3:%2d 行 ,	P3:%2d 个 , 利 用率%f',y1,m1,y2,m2,q);
    y2=0;
    while s*m1<=774||s*m2<=1623 
        s=s+1;
    end
    qq=(373*201*774+406*229*1623)/(s*3000*1500);
    fprintf('\tS1 的块数%d,总的利用率%f\n',s,qq); 
    s=0;
    qq=0;
end
 

%交叉
fprintf('(5)p1 竖放,p3 横放,同行\n');
fprintf('同行不同列\n'); 
x2=0;
y2=0;
for x1=1:12
    while x1*201+x2*406<=3000 
        x2=x2+1;
    end
    while x1*201+x2*406>3000 
        x2=x2-1;
    end
    m1=x1*4; m2=x2*6;
    q=((m1*373*201)+(m2*406*229))/(3000*1500);
    fprintf('P1:%2d 列 ,	P1:%2d 个 ,	P3:%2d 列 ,	P3:%2d 个 , 利 用率%f',x1,m1,x2,m2,q);
    x2=0;
    while s*m1<=774||s*m2<=1623 
        s=s+1;
    end
    qq=(373*201*774+406*229*1623)/(s*3000*1500);
    fprintf('\tS1 的块数%d,总的利用率%f\n',s,qq); 
    s=0;
    qq=0;
end
 
y2=0;
fprintf('(6)p1 竖放,p3 横放,不同行\n'); 
for y1=1:3
    while y1*373+y2*229<=1500 
        y2=y2+1;
    end
    while y1*373+y2*229>1500
        y2=y2-1;
    end
    m1=y1*14; 
    m2=y2*7;
    q=((m1*373*201)+(m2*406*229))/(3000*1500);
    fprintf('P1:%2d 列 ,	P1:%2d 个 ,	P3:%2d 列 ,	P3:%2d 个 , 利 用率%f',y1,m1,y2,m2,q);
    x2=0;
    while s*m1<=774||s*m2<=1623 
        s=s+1;
    end
    qq=(373*201*774+406*229*1623)/(s*3000*1500);
    fprintf('\tS1 的块数%d,总的利用率%f\n',s,qq); 
    s=0;
    qq=0;
end

fprintf('(7)p1 横放,p3 竖放,同行\n'); 
for x1=1:7
    while x1*373+x2*229<=3000 
        x2=x2+1;
    end
    while x1*373+x2*229>3000 
        x2=x2-1;
    end
    m1=x1*7; 
    m2=x2*3;
    q=((m1*373*201)+(m2*406*229))/(3000*1500);
    fprintf('P1:%2d 列 ,	P1:%2d 个 ,	P3:%2d 列 ,	P3:%2d 个 , 利 用率%f',x1,m1,x2,m2,q);
    x2=0;
    while s*m1<=774||s*m2<=1623 
        s=s+1;
    end
    qq=(373*201*774+406*229*1623)/(s*3000*1500);
    fprintf('\tS1 的块数%d,总的利用率%f\n',s,qq); 
    s=0;
    qq=0;
end

y2=0;
fprintf('(8)p1 横放,p3 竖放,不同行\n'); 
for y1=1:5
    while y1*201+y2*406<=1500 
        y2=y2+1;
    end
    while y1*201+y2*406>1500 
        y2=y2-1;
    end
    m1=y1*8; 
    m2=y2*13;
    q=((m1*373*201)+(m2*406*229))/(3000*1500);
    fprintf('P1:%2d 列 ,	P1:%2d 个 ,	P3:%2d 列 ,	P3:%2d 个 , 利 用率%f',y1,m1,y2,m2,q);
    x2=0;
    while s*m1<=774||s*m2<=1623 
        s=s+1;
    end
    qq=(373*201*774+406*229*1623)/(s*3000*1500);
    fprintf('\tS1 的块数%d,总的利用率%f\n',s,qq); 
    s=0;
    qq=0;
end
clc;clear; 
x1=0;y1=1;m1=0; 
x2=0;y2=1;m2=0; 
x3=0;y3=1;m3=0; 
x4=0;y4=1;m4=0; 
q0=3000*1500; 
q1=373*201; 
q2=477*282; 
q3=406*229; 
q4=311*225; 
s=0;
i=1;
fprintf('(1)同列不同行 横放\n');%1 
for y1=1:10
    for y2=1:10
        for y3=1:10
            for y4=1:10
                if y1*201+y2*282+y3*229+y4*225<=1500&&y1*201+y2*282+y3*229+(y4+1)*225>1500
                    m1=y1*8;m2=y2*6;m3=y3*7;m4=y4*9; 
                    Q0=(m1*q1+m2*q2+m3*q3+m4*q4)/q0;
                    while s*m1<774||s*m2<2153||s*m3<1623||s*m4<1614 
                        s=s+1;
                    end
                    Q1=(q1*774+2153*q2+1623*q3+1614*q4)/(q0*s); 
                    if Q1>=0.7
                        fprintf('P1:%d 行、%2d 个,P2:%d 行、%2d 个,P3:%d 行 、 %2d 个 , P4:%d 行 、 %2d 个 , 单 块 利 用 率 %f,S1 的 块 数 %d, 总 利 用率%f\n',y1,m1,y2,m2,y3,m3,y4,m4,Q0,s,Q1);
                    end
                    s=0;
                end
            end
        end
    end
end

fprintf('(2)同行不同列 全部竖放\n');%1 
for x1=1:15
    for x2=1:15
        for x3=1:15
            for x4=1:15 
                if x1*201+x2*282+x3*229+x4*225<=3000&&x1*201+x2*282+x3*229+(x4+1)*225>3000
                    m1=x1*4;m2=x2*3;m3=x3*3;m4=x4; 
                    Q0=(m1*q1+m2*q2+m3*q3+m4*q4)/q0;
                    while s*m1<774||s*m2<2153||s*m3<1623||s*m4<1614 
                        s=s+1;
                    end
                    Q1=(q1*774+2153*q2+1623*q3+1614*q4)/(q0*s); 
                    if Q1>=0.7
                        fprintf('P1:%d 列、%2d 个,P2:%d 列、%2d 个,P3:%d 列 、 %2d 个 , P4:%d 列 、 %2d 个 , 单 块 利 用 率 %f,S1 的 块 数 %d, 总 利 用率%f\n',x1,m1,x2,m2,x3,m3,x4,m4,Q0,s,Q1);
                    end
                    s=0;
                end
            end
        end
    end
end
 

fprintf('(3)同行不同列 P1 竖放,P3 竖放,P3 竖放,P2 横放\n');%1 
for y1=1:15
    for x2=1:15
        for y3=1:15
            for y4=1:15 
                if y1*201+x2*477+y3*229+y4*225<=3000&&y1*201+x2*477+y3*229+(y4+1)*225>3000
                    m1=y1*4;m2=x2*5;m3=y3*3;m4=y4*4; 
                    Q0=(m1*q1+m2*q2+m3*q3+m4*q4)/q0;
                    while s*m1<774||s*m2<2153||s*m3<1623||s*m4<1614 
                        s=s+1;
                    end
                    Q1=(q1*774+2153*q2+1623*q3+1614*q4)/(q0*s); 
                    if Q1>=0.7
                        fprintf('P1:%d 行、%2d 个,P2:%d 列、%2d 个,P3:%d 行 、 %2d 个 , P4:%d 行 、 %2d 个 , 单 块 利 用 率 %f,S1 的 块 数 %d, 总 利 用率%f\n',y1,m1,y2,m2,y3,m3,y4,m4,Q0,s,Q1);
                    end
                    s=0;
                end
            end
        end
    end
end
 

fprintf('(4)同行不同列 P2 竖放,P3 竖放,P4 竖放,P1 横放\n');%1 
for x1=1:15
    for x2=1:15
        for x3=1:15
            for x4=1:15 
                if x1*373+x2*282+x3*229+x4*225<=3000&&x1*373+x2*282+x3*229+(x4+1)*225>3000
                    m1=x1*7;m2=x2*3;m3=x3*3;m4=x4*4; 
                    Q0=(m1*q1+m2*q2+m3*q3+m4*q4)/q0;
                    while s*m1<774||s*m2<2153||s*m3<1623||s*m4<1614 
                        s=s+1;
                    end
                    Q1=(q1*774+2153*q2+1623*q3+1614*q4)/(q0*s); 
                    if Q1>=0.7
                        fprintf('P1:%d 列、%2d 个,P2:%d 列、%2d 个,P3:%d 列 、 %2d 个 , P4:%d 列 、 %2d 个 , 单 块 利 用 率 %f,S1 的 块 数 %d, 总 利 用率%f\n',x1,m1,x2,m2,x3,m3,x4,m4,Q0,s,Q1);
                    end
                    s=0;
                end
            end
        end
    end
end

fprintf('(5)同行不同列 P2 竖放,P1 竖放,P4 竖放,P3 横放\n');%1 
for x1=1:15
    for x2=1:15
        for x3=1:15
            for x4=1:15 
                if x1*201+x2*282+x3*406+x4*225<=3000&&x1*201+x2*282+x3*406+(x4+1)*225>3000
                    m1=x1*4;m2=x2*3;m3=x3*6;m4=x4*4; 
                    Q0=(m1*q1+m2*q2+m3*q3+m4*q4)/q0;
                    while s*m1<774||s*m2<2153||s*m3<1623||s*m4<1614 
                        s=s+1;
                    end
                    Q1=(q1*774+2153*q2+1623*q3+1614*q4)/(q0*s); 
                    if Q1>=0.7
                        fprintf('P1:%d 列、%2d 个,P2:%d 列、%2d 个,P3:%d 列 、 %2d 个 , P4:%d 列 、 %2d 个 , 单 块 利 用 率 %f,S1 的 块 数 %d, 总 利 用率%f\n',x1,m1,x2,m2,x3,m3,x4,m4,Q0,s,Q1);
                    end
                    s=0;
                end
            end
        end
    end
end
 

fprintf('(6)同行不同列 P2 竖放,P1 竖放,P3 竖放,P4 横放\n');%1 
for x1=1:15
    for x2=1:15
        for x3=1:15
            for x4=1:15 
                if x1*201+x2*282+x3*229+x4*311<=3000&&x1*201+x2*282+x3*229+(x4+1)*311>3000
                    m1=x1*4;m2=x2*3;m3=x3*3;m4=x4*6; 
                    Q0=(m1*q1+m2*q2+m3*q3+m4*q4)/q0;
                    while s*m1<774||s*m2<2153||s*m3<1623||s*m4<1614 
                        s=s+1;
                    end
                    Q1=(q1*774+2153*q2+1623*q3+1614*q4)/(q0*s); 
                    if Q1>=0.7
                        fprintf('P1:%d 列、%2d 个,P2:%d 列、%2d 个,P3:%d 列 、 %2d 个 , P4:%d 列 、 %2d 个 , 单 块 利 用 率 %f,S1 的 块 数 %d, 总 利 用率%f\n',x1,m1,x2,m2,x3,m3,x4,m4,Q0,s,Q1);
                    end
                    s=0;
                end
            end
        end
    end
end

fprintf('(7)同行不同列 P3 竖放,P4 竖放,P2 横放,P1 横放\n');%1 
for x1=1:15
    for x2=1:15
        for x3=1:15
            for x4=1:15 
                if x1*373+x2*477+x3*229+x4*225<=3000&&x1*373+x2*477+x3*229+(x4+1)*225>3000
                    m1=x1*7;m2=x2*5;m3=x3*3;m4=x4*4; 
                    Q0=(m1*q1+m2*q2+m3*q3+m4*q4)/q0;
                    while s*m1<774||s*m2<2153||s*m3<1623||s*m4<1614 
                        s=s+1;
                    end
                    Q1=(q1*774+2153*q2+1623*q3+1614*q4)/(q0*s); 
                    if Q1>=0.7
                        fprintf('P1:%d 列、%2d 个,P2:%d 列、%2d 个,P3:%d 列 、 %2d 个 , P4:%d 列 、 %2d 个 , 单 块 利 用 率 %f,S1 的 块 数 %d, 总 利 用率%f\n',x1,m1,x2,m2,x3,m3,x4,m4,Q0,s,Q1);
                    end
                    s=0;
                end
            end
        end
    end
end
 

fprintf('(8)同行不同列 P2 竖放,P4 竖放,P3 横放,P1 横放\n');%1 
for x1=1:15
    for x2=1:15
        for x3=1:15
            for x4=1:15 
                if x1*373+x2*282+x3*406+x4*225<=3000&&x1*373+x2*282+x3*406+(x4+1)*225>3000
                    m1=x1*7;m2=x2*3;m3=x3*6;m4=x4*4; 
                    Q0=(m1*q1+m2*q2+m3*q3+m4*q4)/q0;
                    while s*m1<774||s*m2<2153||s*m3<1623||s*m4<1614
                        s=s+1;
                    end
                    Q1=(q1*774+2153*q2+1623*q3+1614*q4)/(q0*s); 
                    if Q1>=0.7
                        fprintf('P1:%d 列、%2d 个,P2:%d 列、%2d 个,P3:%d 列 、 %2d 个 , P4:%d 列 、 %2d 个 , 单 块 利 用 率 %f,S1 的 块 数 %d, 总 利 用率%f\n',x1,m1,x2,m2,x3,m3,x4,m4,Q0,s,Q1);
                    end
                    s=0;
                end
            end
        end
    end
end

fprintf('(9)同行不同列 P2 竖放,P3 竖放,P4 横放,P1 横放\n');%1 
for x1=1:15
    for x2=1:15
        for x3=1:15
            for x4=1:15 
                if x1*373+x2*282+x3*229+x4*311<=3000&&x1*373+x2*282+x3*229+(x4+1)*311>3000
                    m1=x1*7;m2=x2*3;m3=x3*3;m4=x4*6; 
                    Q0=(m1*q1+m2*q2+m3*q3+m4*q4)/q0;
                    while s*m1<774||s*m2<2153||s*m3<1623||s*m4<1614 
                        s=s+1;
                    end
                    Q1=(q1*774+2153*q2+1623*q3+1614*q4)/(q0*s); 
                    if Q1>=0.7
                        fprintf('P1:%d 列、%2d 个,P2:%d 列、%2d 个,P3:%d 列 、 %2d 个 , P4:%d 列 、 %2d 个 , 单 块 利 用 率 %f,S1 的 块 数 %d, 总 利 用率%f\n',x1,m1,x2,m2,x3,m3,x4,m4,Q0,s,Q1);
                    end
                    s=0;
                end
            end
        end
    end
end 

fprintf('(10)同行不同列 P1 竖放,P4 竖放,P2 横放,P3 横放\n');%1 
for x1=1:15
    for x2=1:15
        for x3=1:15        
            for x4=1:15 
                if x1*201+x2*477+x3*406+x4*225<=3000&&x1*201+x2*477+x3*406+(x4+1)*225>3000
                    m1=x1*4;m2=x2*5;m3=x3*6;m4=x4*4; 
                    Q0=(m1*q1+m2*q2+m3*q3+m4*q4)/q0;
                    while s*m1<774||s*m2<2153||s*m3<1623||s*m4<1614 
                        s=s+1;
                    end
                    Q1=(q1*774+2153*q2+1623*q3+1614*q4)/(q0*s); 
                    if Q1>=0.7
                        fprintf('P1:%d 列、%2d 个,P2:%d 列、%2d 个,P3:%d 列 、 %2d 个 , P4:%d 列 、 %2d 个 , 单 块 利 用 率 %f,S1 的 块 数 %d, 总 利 用率%f\n',x1,m1,x2,m2,x3,m3,x4,m4,Q0,s,Q1);
                    end
                    s=0;
                end
            end
        end
    end
end
 

fprintf('(11)同行不同列 P1 竖放,P3 竖放,P2 横放,P4 横放\n');%1 
for x1=1:15
    for x2=1:15
        for x3=1:15
            for x4=1:15 
                if x1*201+x2*477+x3*229+x4*311<=3000&&x1*201+x2*477+x3*229+(x4+1)*311>3000
                    m1=x1*4;m2=x2*5;m3=x3*3;m4=x4*6; 
                    Q0=(m1*q1+m2*q2+m3*q3+m4*q4)/q0;
                    while s*m1<774||s*m2<2153||s*m3<1623||s*m4<1614 
                        s=s+1;
                    end
                    Q1=(q1*774+2153*q2+1623*q3+1614*q4)/(q0*s); 
                    if Q1>=0.7
                        fprintf('P1:%d 列、%2d 个,P2:%d 列、%2d 个,P3:%d 列 、 %2d 个 , P4:%d 列 、 %2d 个 , 单 块 利 用 率 %f,S1 的 块 数 %d, 总 利 用率%f\n',x1,m1,x2,m2,x3,m3,x4,m4,Q0,s,Q1);
                    end
                    s=0;
                end
            end
        end
    end
end

fprintf('(12)同行不同列 P1 竖放,P2 竖放,P3 横放,P4 横放\n');%1 
for x1=1:15
    for x2=1:15
        for x3=1:15
            for x4=1:15 
                if x1*201+x2*282+x3*406+x4*311<=3000&&x1*201+x2*282+x3*406+(x4+1)*311>3000
                    m1=x1*4;m2=x2*3;m3=x3*6;m4=x4*6; 
                    Q0=(m1*q1+m2*q2+m3*q3+m4*q4)/q0;
                    while s*m1<774||s*m2<2153||s*m3<1623||s*m4<1614 
                        s=s+1;
                    end
                    Q1=(q1*774+2153*q2+1623*q3+1614*q4)/(q0*s); 
                    if Q1>=0.7
                        fprintf('P1:%d 列、%2d 个,P2:%d 列、%2d 个,P3:%d 列 、 %2d 个 , P4:%d 列 、 %2d 个 , 单 块 利 用 率 %f,S1 的 块 数 %d, 总 利 用率%f\n',x1,m1,x2,m2,x3,m3,x4,m4,Q0,s,Q1);
                    end
                    s=0;
                end
            end
        end
    end
end
 

fprintf('(13)同行不同列 P1 竖放,P2 横放,P3 横放,P4 横放\n');%1 
for x1=1:15
    for x2=1:15
        for x3=1:15
            for x4=1:15 
                if x1*201+x2*477+x3*406+x4*311<=3000&&x1*201+x2*477+x3*406+(x4+1)*311>3000
                    m1=x1*4;m2=x2*5;m3=x3*6;m4=x4*6; 
                    Q0=(m1*q1+m2*q2+m3*q3+m4*q4)/q0;
                    while s*m1<774||s*m2<2153||s*m3<1623||s*m4<1614 
                        s=s+1;
                    end
                    Q1=(q1*774+2153*q2+1623*q3+1614*q4)/(q0*s); 
                    if Q1>=0.7
                        fprintf('P1:%d 列、%2d 个,P2:%d 列、%2d 个,P3:%d 列 、 %2d 个 , P4:%d 列 、 %2d 个 , 单 块 利 用 率 %f,S1 的 块 数 %d, 总 利 用率%f\n',x1,m1,x2,m2,x3,m3,x4,m4,Q0,s,Q1);
                    end
                    s=0;
                end
            end
        end
    end
end 

fprintf('(14)同行不同列 P1 横放,P2 竖放,P3 横放,P4 横放\n');%1 
for x1=1:15
    for x2=1:15
        for x3=1:15
            for x4=1:15 
                if x1*373+x2*282+x3*406+x4*311<=3000&&x1*373+x2*282+x3*406+(x4+1)*311>3000
                    m1=x1*7;m2=x2*3;m3=x3*6;m4=x4*6; 
                    Q0=(m1*q1+m2*q2+m3*q3+m4*q4)/q0;
                    while s*m1<774||s*m2<2153||s*m3<1623||s*m4<1614 
                        s=s+1;
                    end
                    Q1=(q1*774+2153*q2+1623*q3+1614*q4)/(q0*s); 
                    if Q1>=0.7
                        fprintf('P1:%d 列、%2d 个,P2:%d 列、%2d 个,P3:%d 列 、 %2d 个 , P4:%d 列 、 %2d 个 , 单 块 利 用 率 %f,S1 的 块 数 %d, 总 利 用率%f\n',x1,m1,x2,m2,x3,m3,x4,m4,Q0,s,Q1);
                    end
                    s=0;
                end
            end
        end
    end
end 

fprintf('(15)同行不同列 P1 横放,P2 横放,P3 竖放,P4 横放\n');%1 
for x1=1:15
    for x2=1:15
        for x3=1:15
            for x4=1:15 
                if x1*373+x2*477+x3*229+x4*311<=3000&&x1*373+x2*477+x3*229+(x4+1)*311>3000
                    m1=x1*7;m2=x2*5;m3=x3*3;m4=x4*6; 
                    Q0=(m1*q1+m2*q2+m3*q3+m4*q4)/q0;
                    while s*m1<774||s*m2<2153||s*m3<1623||s*m4<1614 
                        s=s+1;
                    end
                    Q1=(q1*774+2153*q2+1623*q3+1614*q4)/(q0*s); 
                    if Q1>=0.7
                        fprintf('P1:%d 列、%2d 个,P2:%d 列、%2d 个,P3:%d 列 、 %2d 个 , P4:%d 列 、 %2d 个 , 单 块 利 用 率 %f,S1 的 块 数 %d, 总 利 用率%f\n',x1,m1,x2,m2,x3,m3,x4,m4,Q0,s,Q1);
                    end
                    s=0;
                end
            end
        end
    end
end
 

fprintf('(16)同行不同列 P1 横放,P2 横放,P3 横放,P4 竖放\n');%1 
for x1=1:15
    for x2=1:15
        for x3=1:15
            for x4=1:15 
                if x1*373+x2*477+x3*406+x4*225<=3000&&x1*373+x2*477+x3*406+(x4+1)*225>3000
                    m1=x1*7;m2=x2*5;m3=x3*6;m4=x4*4; 
                    Q0=(m1*q1+m2*q2+m3*q3+m4*q4)/q0;
                    while s*m1<774||s*m2<2153||s*m3<1623||s*m4<1614 
                        s=s+1;
                    end
                    Q1=(q1*774+2153*q2+1623*q3+1614*q4)/(q0*s); 
                    if Q1>=0.7
                        fprintf('P1:%d 列、%2d 个,P2:%d 列、%2d 个,P3:%d 列 、 %2d 个 , P4:%d 列 、 %2d 个 , 单 块 利 用 率 %f,S1 的 块 数 %d, 总 利 用率%f\n',x1,m1,x2,m2,x3,m3,x4,m4,Q0,s,Q1);
                    end
                    s=0;
                end
            end
        end
    end
end
 

fprintf('(17)同行不同列 全部横放\n');%1 
for x1=1:15
    for x2=1:15
        for x3=1:15
            for x4=1:15
                if x1*373+x2*477+x3*406+x4*311<=3000&&x1*373+x2*477+x3*406+(x4+1)*311>3000
                    m1=x1*7;m2=x2*5;m3=x3*6;m4=x4*6; 
                    Q0=(m1*q1+m2*q2+m3*q3+m4*q4)/q0;
                    while s*m1<774||s*m2<2153||s*m3<1623||s*m4<1614 
                        s=s+1;
                    end
                    Q1=(q1*774+2153*q2+1623*q3+1614*q4)/(q0*s); 
                    if Q1>=0.7
                        fprintf('P1:%d 列、%2d 个,P2:%d 列、%2d 个,P3:%d 列 、 %2d 个 , P4:%d 列 、 %2d 个 , 单 块 利 用 率 %f,S1 的 块 数 %d, 总 利 用率%f\n',x1,m1,x2,m2,x3,m3,x4,m4,Q0,s,Q1);
                    end
                    s=0;
                end
            end
        end
    end
end
 

fprintf('(18)同列不同行 全部竖放\n');%1 
for x1=1:15
    for x2=1:15
        for x3=1:15
            for x4=1:15
                if x1*373+x2*477+x3*406+x4*311<=1500&&x1*373+x2*477+x3*406+(x4+1)*311>1500
                    m1=x1*14;m2=x2*10;m3=x3*13;m4=x4*13;
                    Q0=(m1*q1+m2*q2+m3*q3+m4*q4)/q0;
                    while s*m1<774||s*m2<2153||s*m3<1623||s*m4<1614
                        s=s+1;
                    end
                    Q1=(q1*774+2153*q2+1623*q3+1614*q4)/(q0*s);
                    if Q1>=0.7
                        fprintf('P1:%d 列、%2d 个,P2:%d 列、%2d 个,P3:%d 列 、 %2d 个 , P4:%d 列 、 %2d 个 , 单 块 利 用 率 %f,S1 的 块 数 %d, 总 利 用率%f\n',x1,m1,x2,m2,x3,m3,x4,m4,Q0,s,Q1);
                    end
                    s=0;
                end
            end
        end
    end
end
 
%1500 为基准
fprintf('(19)同列不同行 P2 竖放,P3 竖放,P4 竖放,P1 横放\n');%1 
for x1=1:15
    for x2=1:15
        for x3=1:15
            for x4=1:15 
                if x1*201+x2*477+x3*406+x4*311<=1500&&x1*201+x2*477+x3*406+(x4+1)*311>1500
                    m1=x1*8;m2=x2*10;m3=x3*13;m4=x4*13; 
                    Q0=(m1*q1+m2*q2+m3*q3+m4*q4)/q0;
                    while s*m1<774||s*m2<2153||s*m3<1623||s*m4<1614 
                        s=s+1;
                    end
                    Q1=(q1*774+2153*q2+1623*q3+1614*q4)/(q0*s); 
                    if Q1>=0.7
                        fprintf('P1:%d 列、%2d 个,P2:%d 列、%2d 个,P3:%d 列 、 %2d 个 , P4:%d 列 、 %2d 个 , 单 块 利 用 率 %f,S1 的 块 数 %d, 总 利 用率%f\n',x1,m1,x2,m2,x3,m3,x4,m4,Q0,s,Q1);
                    end
                    s=0;
                end
            end
        end
    end
end
 

fprintf('(20)同列不同行 P1 竖放,P3 竖放,P4 竖放,P2 横放\n');%1 
for x1=1:15
    for x2=1:15
        for x3=1:15
            for x4=1:15 
                if x1*373+x2*282+x3*406+x4*311<=1500&&x1*373+x2*282+x3*406+(x4+1)*311>1500
                    m1=x1*14;m2=x2*6;m3=x3*13;m4=x4*13; 
                    Q0=(m1*q1+m2*q2+m3*q3+m4*q4)/q0;
                    while s*m1<774||s*m2<2153||s*m3<1623||s*m4<1614 
                        s=s+1;
                    end
                    Q1=(q1*774+2153*q2+1623*q3+1614*q4)/(q0*s); 
                    if Q1>=0.7
                        fprintf('P1:%d 列、%2d 个,P2:%d 列、%2d 个,P3:%d 列 、 %2d 个 , P4:%d 列 、 %2d 个 , 单 块 利 用 率 %f,S1 的 块 数 %d, 总 利 用率%f\n',x1,m1,x2,m2,x3,m3,x4,m4,Q0,s,Q1);
                    end
                    s=0;
                end
            end
        end
    end
end
 

fprintf('(21)同列不同行 P1 竖放,P2 竖放,P4 竖放,P3 横放\n');%1 
for x1=1:15
    for x2=1:15
        for x3=1:15
            for x4=1:15 
                if x1*373+x2*477+x3*229+x4*311<=1500&&x1*373+x2*477+x3*229+(x4+1)*311>1500
                    m1=x1*14;m2=x2*10;m3=x3*7;m4=x4*13; 
                    Q0=(m1*q1+m2*q2+m3*q3+m4*q4)/q0;
                    while s*m1<774||s*m2<2153||s*m3<1623||s*m4<1614 
                        s=s+1;
                    end
                    Q1=(q1*774+2153*q2+1623*q3+1614*q4)/(q0*s); 
                    if Q1>=0.7
                        fprintf('P1:%d 列、%2d 个,P2:%d 列、%2d 个,P3:%d 列 、 %2d 个 , P4:%d 列 、 %2d 个 , 单 块 利 用 率 %f,S1 的 块 数 %d, 总 利 用率%f\n',x1,m1,x2,m2,x3,m3,x4,m4,Q0,s,Q1);
                    end
                    s=0;
                end
            end
        end
    end
end
 

fprintf('(22)同列不同行 P1 竖放,P2 竖放,P3 竖放,P4 横放\n');%1 
for x1=1:15
    for x2=1:15
        for x3=1:15
            for x4=1:15 
                if x1*373+x2*477+x3*406+x4*225<=1500&&x1*373+x2*477+x3*406+(x4+1)*225>1500
                    m1=x1*14;m2=x2*10;m3=x3*13;m4=x4*9; 
                    Q0=(m1*q1+m2*q2+m3*q3+m4*q4)/q0;
                    while s*m1<774||s*m2<2153||s*m3<1623||s*m4<1614 
                        s=s+1;
                    end
                    Q1=(q1*774+2153*q2+1623*q3+1614*q4)/(q0*s); 
                    if Q1>=0.7
                        fprintf('P1:%d 列、%2d 个,P2:%d 列、%2d 个,P3:%d 列 、 %2d 个 , P4:%d 列 、 %2d 个 , 单 块 利 用 率 %f,S1 的 块 数 %d, 总 利 用率%f\n',x1,m1,x2,m2,x3,m3,x4,m4,Q0,s,Q1);
                    end
                    s=0;
                end
            end
        end
    end
end
 

fprintf('(23)同列不同行 P1 竖放,P2 竖放,P3 横放,P4 横放\n');%1 
for x1=1:15
    for x2=1:15
        for x3=1:15
            for x4=1:15 
                if x1*373+x2*477+x3*229+x4*225<=1500&&x1*373+x2*477+x3*229+(x4+1)*225>1500
                    m1=x1*14;m2=x2*10;m3=x3*7;m4=x4*9; 
                    Q0=(m1*q1+m2*q2+m3*q3+m4*q4)/q0;
                    while s*m1<774||s*m2<2153||s*m3<1623||s*m4<1614 
                        s=s+1;
                    end
                    Q1=(q1*774+2153*q2+1623*q3+1614*q4)/(q0*s); 
                    if Q1>=0.7
                        fprintf('P1:%d 列、%2d 个,P2:%d 列、%2d 个,P3:%d 列 、 %2d 个 , P4:%d 列 、 %2d 个 , 单 块 利 用 率 %f,S1 的 块 数 %d, 总 利 用率%f\n',x1,m1,x2,m2,x3,m3,x4,m4,Q0,s,Q1);
                    end
                    s=0;
                end
            end
        end
    end
end
 

fprintf('(24)同列不同行 P1 竖放,P2 横放,P3 竖放,P4 横放\n');%1 
for x1=1:15
    for x2=1:15
        for x3=1:15
            for x4=1:15
                if x1*373+x2*282+x3*406+x4*225<=1500&&x1*373+x2*282+x3*406+(x4+1)*225>1500
                    m1=x1*14;m2=x2*6;m3=x3*13;m4=x4*9; 
                    Q0=(m1*q1+m2*q2+m3*q3+m4*q4)/q0;
                    while s*m1<774||s*m2<2153||s*m3<1623||s*m4<1614 
                        s=s+1;
                    end
                    Q1=(q1*774+2153*q2+1623*q3+1614*q4)/(q0*s); 
                    if Q1>=0.7
                        fprintf('P1:%d 列、%2d 个,P2:%d 列、%2d 个,P3:%d 列 、 %2d 个 , P4:%d 列 、 %2d 个 , 单 块 利 用 率 %f,S1 的 块 数 %d, 总 利 用率%f\n',x1,m1,x2,m2,x3,m3,x4,m4,Q0,s,Q1);
                    end
                    s=0;
                end
            end
        end
    end
end
 

fprintf('(25)同列不同行 P1 竖放,P3 横放,P4 竖放,P2 横放\n');%1 
for x1=1:15
    for x2=1:15
        for x3=1:15
            for x4=1:15 
                if x1*373+x2*282+x3*229+x4*311<=1500&&x1*373+x2*282+x3*229+(x4+1)*311>1500
                    m1=x1*14;m2=x2*6;m3=x3*7;m4=x4*13; 
                    Q0=(m1*q1+m2*q2+m3*q3+m4*q4)/q0;
                    while s*m1<774||s*m2<2153||s*m3<1623||s*m4<1614 
                        s=s+1;
                    end
                    Q1=(q1*774+2153*q2+1623*q3+1614*q4)/(q0*s); 
                    if Q1>=0.7
                        fprintf('P1:%d 列、%2d 个,P2:%d 列、%2d 个,P3:%d 列 、 %2d 个 , P4:%d 列 、 %2d 个 , 单 块 利 用 率 %f,S1 的 块 数 %d, 总 利 用率%f\n',x1,m1,x2,m2,x3,m3,x4,m4,Q0,s,Q1);
                    end
                    s=0;
                end
            end
        end
    end
end
 
fprintf('(26)同列不同行 P1 横放,P2 竖放,P3 竖放,P4 横放\n');%1 
for x1=1:15
    for x2=1:15
        for x3=1:15
            for x4=1:15 
                if x1*201+x2*477+x3*406+x4*225<=1500&&x1*201+x2*477+x3*406+(x4+1)*225>1500
                    m1=x1*8;m2=x2*10;m3=x3*13;m4=x4*9; 
                    Q0=(m1*q1+m2*q2+m3*q3+m4*q4)/q0;
                    while s*m1<774||s*m2<2153||s*m3<1623||s*m4<1614 
                        s=s+1;
                    end
                    Q1=(q1*774+2153*q2+1623*q3+1614*q4)/(q0*s); 
                    if Q1>=0.7
                        fprintf('P1:%d 列、%2d 个,P2:%d 列、%2d 个,P3:%d 列 、 %2d 个 , P4:%d 列 、 %2d 个 , 单 块 利 用 率 %f,S1 的 块 数 %d, 总 利 用率%f\n',x1,m1,x2,m2,x3,m3,x4,m4,Q0,s,Q1);
                    end
                    s=0;
                end
            end
        end
    end
end
 

fprintf('(27)同列不同行 P1 横放,P2 竖放,P4 竖放,P3 横放\n');%1 
for x1=1:15
    for x2=1:15
        for x3=1:15
            for x4=1:15 
                if x1*201+x2*477+x3*229+x4*311<=1500&&x1*201+x2*477+x3*229+(x4+1)*311>1500
                    m1=x1*8;m2=x2*10;m3=x3*7;m4=x4*13; 
                    Q0=(m1*q1+m2*q2+m3*q3+m4*q4)/q0;
                    while s*m1<774||s*m2<2153||s*m3<1623||s*m4<1614 
                        s=s+1;
                    end
                    Q1=(q1*774+2153*q2+1623*q3+1614*q4)/(q0*s); 
                    if Q1>=0.7
                        fprintf('P1:%d 列、%2d 个,P2:%d 列、%2d 个,P3:%d 列 、 %2d 个 , P4:%d 列 、 %2d 个 , 单 块 利 用 率 %f,S1 的 块 数 %d, 总 利 用率%f\n',x1,m1,x2,m2,x3,m3,x4,m4,Q0,s,Q1);
                    end
                    s=0;
                end
            end
        end
    end
end 

fprintf('(28)同列不同行 P1 横放,P3 竖放,P4 竖放,P2 横放\n');%1 
for x1=1:15
    for x2=1:15
        for x3=1:15
            for x4=1:15 
                if x1*201+x2*282+x3*406+x4*311<=1500&&x1*201+x2*282+x3*406+(x4+1)*311>1500
                    m1=x1*8;m2=x2*6;m3=x3*13;m4=x4*13; 
                    Q0=(m1*q1+m2*q2+m3*q3+m4*q4)/q0;
                    while s*m1<774||s*m2<2153||s*m3<1623||s*m4<1614 
                        s=s+1;
                    end
                    Q1=(q1*774+2153*q2+1623*q3+1614*q4)/(q0*s); 
                    if Q1>=0.7
                        fprintf('P1:%d 列、%2d 个,P2:%d 列、%2d 个,P3:%d 列 、 %2d 个 , P4:%d 列 、 %2d 个 , 单 块 利 用 率 %f,S1 的 块 数 %d, 总 利 用率%f\n',x1,m1,x2,m2,x3,m3,x4,m4,Q0,s,Q1);
                    end
                    s=0;
                end
            end
        end
    end
end
 

fprintf('(29)同列不同行 P1 竖放,P2 横放,P3 横放,P4 横放\n');%1 
for x1=1:15
    for x2=1:15
        for x3=1:15
            for x4=1:15 
                if x1*373+x2*282+x3*229+x4*225<=1500&&x1*373+x2*282+x3*229+(x4+1)*225>1500
                    m1=x1*14;m2=x2*6;m3=x3*7;m4=x4*9; 
                    Q0=(m1*q1+m2*q2+m3*q3+m4*q4)/q0;
                    while s*m1<774||s*m2<2153||s*m3<1623||s*m4<1614 
                        s=s+1;
                    end
                    Q1=(q1*774+2153*q2+1623*q3+1614*q4)/(q0*s); 
                    if Q1>=0.7
                        fprintf('P1:%d 列、%2d 个,P2:%d 列、%2d 个,P3:%d 列 、 %2d 个 , P4:%d 列 、 %2d 个 , 单 块 利 用 率 %f,S1 的 块 数 %d, 总 利 用率%f\n',x1,m1,x2,m2,x3,m3,x4,m4,Q0,s,Q1);
                    end
                    s=0;
                end
            end
        end
    end
end 

fprintf('(30)同列不同行 P1 横放,P2 竖放,P4 横放,P3 横放\n');%1 
for x1=1:15
    for x2=1:15
        for x3=1:15
            for x4=1:15 
                if x1*201+x2*477+x3*229+x4*225<=1500&&x1*201+x2*477+x3*229+(x4+1)*225>1500
                    m1=x1*8;m2=x2*10;m3=x3*7;m4=x4*9; 
                    Q0=(m1*q1+m2*q2+m3*q3+m4*q4)/q0;
                    while s*m1<774||s*m2<2153||s*m3<1623||s*m4<1614 
                        s=s+1;
                    end
                    Q1=(q1*774+2153*q2+1623*q3+1614*q4)/(q0*s); 
                    if Q1>=0.7
                        fprintf('P1:%d 列、%2d 个,P2:%d 列、%2d 个,P3:%d 列 、 %2d 个 , P4:%d 列 、 %2d 个 , 单 块 利 用 率 %f,S1 的 块 数 %d, 总 利 用率%f\n',x1,m1,x2,m2,x3,m3,x4,m4,Q0,s,Q1);
                    end
                    s=0;
                end
            end
        end
    end
end
 

fprintf('(31)同列不同行 P1 横放,P3 竖放,P4 横放,P2 横放\n');%1 
for x1=1:15
    for x2=1:15
        for x3=1:15
            for x4=1:15 
                if x1*201+x2*282+x3*406+x4*225<=1500&&x1*201+x2*282+x3*406+(x4+1)*225>1500
                    m1=x1*8;m2=x2*6;m3=x3*13;m4=x4*9; 
                    Q0=(m1*q1+m2*q2+m3*q3+m4*q4)/q0;
                    while s*m1<774||s*m2<2153||s*m3<1623||s*m4<1614
                        s=s+1;
                    end
                    Q1=(q1*774+2153*q2+1623*q3+1614*q4)/(q0*s); 
                    if Q1>=0.7
                        fprintf('P1:%d 列、%2d 个,P2:%d 列、%2d 个,P3:%d 列 、 %2d 个 , P4:%d 列 、 %2d 个 , 单 块 利 用 率 %f,S1 的 块 数 %d, 总 利 用率%f\n',x1,m1,x2,m2,x3,m3,x4,m4,Q0,s,Q1);
                    end
                    s=0;
                end
            end
        end
    end
end
 

fprintf('(32)同列不同行 P1 横放,P3 横放,P4 竖放,P2 横放\n');%1 
for x1=1:15
    for x2=1:15
        for x3=1:15
            for x4=1:15 
                if x1*201+x2*282+x3*229+x4*311<=1500&&x1*201+x2*282+x3*229+(x4+1)*311>1500
                    m1=x1*8;m2=x2*6;m3=x3*7;m4=x4*13; 
                    Q0=(m1*q1+m2*q2+m3*q3+m4*q4)/q0;
                    while s*m1<774||s*m2<2153||s*m3<1623||s*m4<1614 
                        s=s+1;
                    end
                    Q1=(q1*774+2153*q2+1623*q3+1614*q4)/(q0*s); 
                    if Q1>=0.7
                        fprintf('P1:%d 列、%2d 个,P2:%d 列、%2d 个,P3:%d 列 、 %2d 个 , P4:%d 列 、 %2d 个 , 单 块 利 用 率 %f,S1 的 块 数 %d, 总 利 用率%f\n',x1,m1,x2,m2,x3,m3,x4,m4,Q0,s,Q1);
                    end
                    s=0;
                end
            end
        end
    end
end
全部论文请见下方“ 只会建模 QQ名片” 点击QQ名片即可

你可能感兴趣的:(数学建模,数据分析,五一杯,数学建模,五一杯数学建模,数据分析,数学建模数据分析)