apis/examples_matlab_control/cutstock_8m_source.html

함수컷스톡(varargin)


% 인수에서 작업공간 정보 확인

    if나긴 > 0

wsInfo = 피망 슬롯control.WorkspaceInfo();

wsInfo.systemDirectory = varargin1;

ws = 피망 슬롯control.Workspace(wsInfo);

    그밖에

ws = 피망 슬롯control.Workspace();

끝


% 옵션 인스턴스화 및 매개변수 정의

opt = ws.addOptions();

cutstockData = ws.addDatabase('csdata');

opt.setAllModelTypes('복합체');

opt.optCR = 0; % 최적성 해결

최대 패턴 = 35;

opt.defines('pmax', 스프린트('%d', 최대 패턴));

opt.defines('solveMasterAs', 'RMIP');


%가 입력 데이터를 정의

d = 컨테이너.맵('i1', 'i2', 'i3', 'i4'}, {97, 610, 395, 211});

w = 컨테이너.맵('i1', 'i2', 'i3', 'i4'}, {45, 36, 31, 14});

r = 100;


너비 = cutstockData.addSet('나', 1,'너비');

rawWidth = cutstockData.addParameter('r', 0,'원시 너비');

수요 = cutstockData.addParameter('d', 1,'수요');

너비 = cutstockData.addParameter('w', 1,'너비');


    for키 = 키(d)

widths.addRecord(키1);

rec =demand.addRecord(키1);

rec.value = d(키1);

끝

    for키 = 키(w)

rec = width.addRecord(키1);

rec.value = w(키1);

끝

rec = rawWidth.addRecord();

rec.값 = 100;


마스터모델 =

    '$제목 절단 재고 - 마스터 문제 '

    '                                                    '

    'i 너비 설정'

    '매개변수'

    ' w(i) 너비 '

    ' d(i) 수요 '

    '스칼라'

    ' r 원시 너비;                                    '

    '$gdxin csdata '

    '$load i w d r '

    '                                                    '

    '$pmax를 설정하지 않은 경우 $set pmax 1000 '

    'p개의 가능한 패턴 설정 /1*%pmax%/ '

    '     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, 수요/;                      '

    '                                                    '};

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


하위 모델 =

    '$제목 절단 재고 - 가격 문제는 배낭 모델입니다. '

    '                                                           '

    'i 너비 설정'

    '매개변수'

    ' w(i) 너비;                                               '

    '스칼라'

    ' r 원시 너비;                                           '

    '                                                           '

    '$gdxin csdata '

    '$load i w r '

    '                                                           '

    '매개변수'

    ' demdual(i) 마스터 수요 제약의 이중 /#i eps/;    '

    '                                                           '

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

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

    '                                                           '

    '방정식 defobj '

    '         배낭;                                         '

    '                                                           '

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

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

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

    '                                                           '};

하위 모델 = sprintf('%s\n', 하위 모델:);


% 초기 체크포인트 생성

masterCP = ws.addCheckpoint();

masterInitJob = ws.addJobFromString(masterModel);

masterInitJob.run(opt, masterCP, cutstockData);


masterJob = ws.addJobFromString('execute_load ''csdata'', aip, pp; %solveMasterAs%를 사용하여 마스터 최소 z를 해결합니다.', 마스터CP);


패턴 = cutstockData.addSet('pp', 1,'패턴 색인');

patternData = cutstockData.addParameter('아이피', 2,'패턴 데이터');


% 초기 패턴: 패턴 i는 너비 i를 유지합니다.

패턴수 = 0;

    for키 = 키(w)

patternCount = 패턴수 + 1;

patternCount_str = 스프린트f('%d', 패턴수);

pattern.addRecord(patternCount_str);

rec = PatternData.addRecord(key1, PatternCount_str);

rec.value = 바닥(r / w(키1));

끝


% 모델 인스턴스 생성for하위 직업

subCP = ws.addCheckpoint();

ws.addJobFromString(subModel).run(opt, subCP, cutstockData);


subMI = subCP.addModelInstance();


% 수정자 demdual 정의

demandDual = subMI.syncDB.addParameter('뎀듀얼', 1,'마스터의 이중 요구');

subMI.instantiate('mip를 사용하여 최소 z 가격 책정', opt, 피망 슬롯control.Modifier(demandDual));


% 발견신규패턴

패턴 추가됨 =참;

    그동안패턴추가됨

masterJob.run(opt, masterCP, cutstockData);


% gmssubMI.syncDB DB에 이중 복사

demandDual.clear();

        ford = masterJob.outDB.getEquation('수요').기록

rec =demandDual.addRecord(d1.key(1));

rec.value = d1.marginal;

끝


subMI.solve();


z_level = subMI.syncDB.getVariable('z').record.level;

        ifz_level < -0.00001

if PatternCount == maxpattern

fprintf('패턴이 벗어났습니다. 최대 패턴을 늘립니다(현재 [%d])\n', maxpattern);

패턴 추가됨 =거짓;

            그밖에

fprintf('새 패턴 값: %g\n', z_level);

패턴카운트 = 패턴카운트 + 1;

s = 패턴.addRecord(sprintf('%d', 패턴수));

                fory = subMI.syncDB.getVariable('y').기록

                    ify1.level > 0.5

rec = PatternData.addRecord(y1.key(1), s.key(1));

rec.value = round(y1.level);

끝

끝

끝

        그밖에

패턴 추가됨 =거짓;

끝

끝


% 해결최종MIP

opt.defines("solveMasterAs", "MIP");

masterJob.run(opt, cutstockData);


fprintf("최적의 해결책: %g\n", masterJob.outDB.getVariable("z").record.level);

    forxp = masterJob.outDB.getVariable("xp").기록

        ifxp1.level > 0.5

fprintf(" 패턴 [%s] [%g] 회: ", xp1.key(1), xp1.level);

            foraip = masterJob.outDB.getParameter("aip").기록

fprintf(" [%s] : [%g]\n", aip1.key(1), aip1.value);

끝

fprintf("\n");

끝

끝


관리되지 않는 리소스 % 정리

cutstockData.dispose();

subMI.dispose();

opt.dispose();

끝