apis/examples_matlab_control/CutstockModel_8m_source.html

classdef CutstockModel < 핸들


속성(SetAccess =비공개)

너비

raw_width

수요

폭

pat_rep

ws

cutstock_data

정보

선택

직업

끝


메소드


함수 무료 슬롯 게임 = CutstockModel(ws)


무료 슬롯 게임ws = ws;

무료 슬롯 게임opt = ws.addOptions();


무료 슬롯 게임cutstock_data = 무료 슬롯 게임ws.addDatabase('gdxincname');


무료 슬롯 게임opt.defines('gdxincname', 'gdxincname');

무료 슬롯 게임opt.defines('dbOut1', 'dbOut1');

무료 슬롯 게임opt.solveLink = gams.control.options.SolveLink.LoadLibrary;


무료 슬롯 게임widths = 무료 슬롯 게임cutstock_data.addSet('i', 'widths');

무료 슬롯 게임raw_width = 무료 슬롯 게임cutstock_data.addParameter('r', '원시 너비');

무료 슬롯 게임demand = 무료 슬롯 게임cutstock_data.addParameter('d', '수요', 무료 슬롯 게임widths);

무료 슬롯 게임width = 무료 슬롯 게임cutstock_data.addParameter('w', 'width', 무료 슬롯 게임widths);


무료 슬롯 게임job = 무료 슬롯 게임ws.addJobFromString(무료 슬롯 게임modelSource());

끝


함수 실행(무료 슬롯 게임, 출력)

if ~무료 슬롯 게임cutstock_data.checkDomains()

error('컷스톡 데이터베이스의 도메인 오류');

끝


비어 있는 경우(출력)

무료 슬롯 게임job.run(무료 슬롯 게임opt, false, 무료 슬롯 게임cutstock_data);

그밖에

무료 슬롯 게임job.run(무료 슬롯 게임opt, 출력, false, 무료 슬롯 게임cutstock_data);

끝

무료 슬롯 게임dbout = 무료 슬롯 게임ws.addDatabaseFromGDX([무료 슬롯 게임opt.getDefinitionOf('dbOut1'), '.gdx']);

무료 슬롯 게임pat_rep = 무료 슬롯 게임dbout.getParameter('patrep');

끝


함수 모델 = modelSource(무료 슬롯 게임)

모델 =

            '$Title 절단 스톡 - 열 생성 접근 방식(CUTSTOCK,SEQ=294) '

            '                                                                           '

            '$ontext'

            '과제는 '에서 다양한 크기의 종이 제품을 잘라내는 것입니다.

            '고객을 만나기 위한 대형 원시 종이 롤'님의 명령입니다. 목표 '

            '필요한 용지 롤 수를 최소화하는 것입니다.                         '

            '                                                                           '

            'P. C. Gilmore 및 R. E. Gomory, '에 대한 선형 프로그래밍 접근 방식

            '재고 문제, 파트 I, Operations Research 9(1961), 849-859.      '

            '                                                                           '

            'P. C. Gilmore 및 R. E. Gomory, '에 대한 선형 프로그래밍 접근 방식

            '재고 문제, 파트 II, Operations Research 11(1963), 863-888.    '

            '$offtext '

            '                                                                           '

            'i 너비 설정'

            '매개변수'

            '    r 원시 너비 '

            '    w(i) 너비 '

            '    d(i) 수요;                                                          '

            '                                                                           '

            '$gdxincname $abort가 설정되지 않은 경우 ''제공된 데이터 파일에 대한 포함 파일 이름이 없습니다'''

            '$gdxin %gdxincname% '

            '$load r i w d '

            '$gdxin '

            '                                                                           '

            '* 길모어-고모리 열 생성 알고리즘 '

            '                                                                           '

            'p개의 가능한 패턴 설정 /p1*p1000/ '

            '     pp(p) p '의 동적 하위 집합

            '매개변수'

            '    aip(i,p) p에서 성장하는 패턴의 너비 i 수;                    '

            '                                                                           '

            '* 마스터 모델 '

            '사용된 가변 xp(p) 패턴'

            '         z 목적 변수 '

            '정수 변수 xp; xp.up(p) = sum(i, d(i));                              '

            '                                                                           '

            '사용된 패턴의 방정식 숫자 '

            '         수요(i) 수요를 충족합니다.                                            '

            '                                                                           '

            'numpat.. z =e= sum(pp, xp(pp));                                        '

            '수요(i).. sum(pp, aip(i,pp)*xp(pp)) =g= d(i);                           '

            '                                                                           '

            '모델 마스터 /numpat, 수요/;                                             '

            '                                                                           '

            '* 가격 문제 - 배낭 모델 '

            '변수 y(i) 새로운 패턴;                                                '

            '정수 변수 y; y.up(i) = ceil(r/w(i));                                '

            '                                                                           '

            '방정식 defobj'

            '         배낭 배낭 제약;                                     '

            '                                                                           '

            'defobj.. z =e= 1 - sum(i, 수요.m(i)*y(i));                           '

            '배낭.. sum(i, w(i)*y(i)) =l= r;                                      '

            '                                                                           '

            '모델 가격 /defobj, 배낭/;                                          '

            '                                                                           '

            '* 초기화 - 초기 패턴은 단일 너비를 갖습니다. '

            'pp(p) = ord(p)<=카드(i);                                                   '

            'aip(i,pp(p))$(ord(i)=ord(p)) = Floor(r/w(i));                              '

            '*디스플레이 aip;                                                              '

            '                                                                           '

            '스칼라 완료 루프 표시기 /0/ '

            '마지막 패턴의 pi(p) 세트를 설정합니다; 파이(p) = ord(p)=카드(pp)+1;           '

            '                                                                           '

            '옵션 optcr=0,limrow=0,limcol=0,solprint=off;                             '

            '                                                                           '

            '그동안(완료되지 않았고 카드(pp)<카드(p), '

            '    z를 최소화하는 rmip를 사용하여 마스터를 해결합니다.                                  '

            '    z를 최소화하는 mip를 사용하여 가격 책정을 해결합니다.                                  '

            '                                                                           '

            '* 마스터 모델을 개선할 수 있는 패턴을 찾았습니까?                       '

            '  if(z.l < -0.001, '

            '     aip(i,pi) = round(y.l(i));                                            '

            '     pp(pi) = 예; 파이(p) = 파이(p-1);                                        '

            '  그렇지 않으면 '

            '     완료 = 1;                                                             '

            '  );                                                                       '

            ');                                                                         '

            '디스플레이''롤 수의 하한값'', master.objval;                '

            '                                                                           '

            '옵션 solprint=on;                                                        '

            'z를 최소화하는 mip를 사용하여 마스터를 해결합니다.                                       '

            '                                                                           '

            '매개변수 patrep 솔루션 패턴 보고서'

            '        demrep 솔루션 수요 공급 보고서;                              '

            '                                                                           '

            '패트렙(''# 생산됨'',p) = round(xp.l(p));                                 '

            'patrep(i,p)$patrep(''# 생산됨'',p) = aip(i,p);                           '

            'patrep(i,''합계'') = 합계(p, patrep(i,p));                                 '

            '패트렙(''# 생산됨'',''합계'') = 합계(p, 패트렙(''# 생산됨'',p));       '

            '                                                                           '

            'demrep(i,''생산됨'') = 합계(p,patrep(i,p)*patrep(''# 생산됨'',p));      '

            'demrep(i,''수요'') = d(i);                                               '

            'demrep(i,''끝'') = demrep(i,''생산됨'') - demrep(i,''수요'');        '

            '                                                                           '

            '표시 patrep, demrep;                                                    '

            '$dbOut1 $abort를 설정하지 않은 경우 ''제공된 데이터베이스 외부 1 파일에 대한 파일 이름이 없습니다'''

            'execute_unload ''%dbOut1%'', 패트렙;                                       '

            '                                                                           '};

모델 = 스프린트f('%s\n', 모델:);

끝


끝


끝