apis/examples_matlab_control/tsp_8m_source.html

함수무료 슬롯(varargin)


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

    if나긴 > 0

wsInfo = 무료 슬롯control.WorkspaceInfo();

wsInfo.systemDirectory = varargin1;

ws = 무료 슬롯control.Workspace(wsInfo);

    그밖에

ws = 무료 슬롯control.Workspace();

끝


모델 =

    '$Title 여행하는 외판원 문제 '

    '$Ontext '

    '                                                                '

    'sub_tour 제거 제약 조건은 Python에 의해 생성됩니다.'

    '스크립트. MIP는 계속해서 해결되지만 무료 슬롯는 '

    'n개의 컷이 추가된 후에만 모델을 생성합니다.           '

    '                                                                '

    '$Offtext '

    '                                                                '

    '$tspdata를 설정하지 않은 경우 $abort ''tspdata가 설정되지 않음''                  '

    '                                                                '

    '2개 도시 설정'

    '    i(ii) 도시의 하위 집합 '

    '별칭 (ii,jj),(i,j,k);                                          '

    '                                                                '

    '매개변수 c(ii,jj) 거리 행렬;                             '

    '                                                                '

    '$gdxin %tspdata% '

    '$load ii c '

    '                                                                '

    '$nrCities를 설정하지 않은 경우 $set nrCities 20'

    'i(ii)$(ord(ii) < %nrCities%) = 예;                             '

    '                                                                '

    '변수 x(ii,jj) 결정 변수 - 여행 구간 '

    '          z 목적변수;                         '

    '이진변수 x; x.fx(ii,ii) = 0;                             '

    '                                                                '

    '방정식 목표 총 비용'

    '          rowsum(ii) 각 도시를 한 번만 떠나세요 '

    '          colsum(jj)는 각 도시에 한 번만 도착합니다.            '

    '                                                                '

    '* 할당 문제는 무료 슬롯의 완화입니다. '

    '목적.. z =e= sum((i,j), c(i,j)*x(i,j));                    '

    'rowsum(i).. sum(j, x(i,j)) =e= 1;                               '

    'colsum(j)..sum(i, x(i,j)) =e= 1;                               '

    '                                                                '

    '$cmax를 설정하지 않은 경우 $set cmax 2 '

    '컷 설정 /c0*c%cmax%/;                                           '

    '매개변수'

    '    acut(무료 슬롯ii,jj) 절단 제약 행렬 '

    '    rhscut(무료 슬롯) 제약 조건을 잘라냅니다. rhs;                            '

    '                                                                '

    'equation sscut(무료 슬롯) sub_tour 제거 컷;                   '

    'sscut(무료 슬롯).. sum((i,j), Acut(무료 슬롯i,j)*x(i,j)) =l= RHScut(무료 슬롯);  '

    '                                                                '

    'cc(컷) 이전 컷 설정; cc(컷) = 아니오;                        '

    '$if set cutdata 실행_로드 ''%cutdata%'', 참조, Acut, RHScut;   '

    '                                                                '

    'Acut(무료 슬롯i,j)$(cc(무료 슬롯) 아님) = eps;                              '

    'RHScut(컷)$(cc(컷) 아님) = 카드(ii);                           '

    '                                                                '

    '모델 할당 /all/;                                             '

    '                                                                '

    '옵션 optcr=0;                                                 '};

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


모델 인스턴스에 추가할 수 있는 컷 수 %

cutsPerRound = 10;

% 현재 삭감

curCut = 0;

% 삭감 한도for현재 모델 인스턴스(cmax = curCut + cutsPerRound)

c최대 = 0;


생성된 모든 컷을 수집하는 데 사용되는 데이터베이스 %

cutData = ws.addDatabase();

cc = cutData.addSet('cc', 1,'');

acut = cutData.addParameter('아컷', 3,'');

rhscut = cutData.addParameter('rhscut', 1,'');


% 도시 목록(i1, i2, i3, ...)

n = ;

하위 투어 = ;


    그동안 사실

% 생성신규컷 제한에 도달한 경우의 모델 인스턴스

        ifcurCut >= cMax

fprintf(',');

cMax = curCut + cutsPerRound;

cutData.export();


% 체크포인트 생성

tspJob = ws.addJobFromString(모델);

cp = ws.addCheckpoint();

opt = ws.addOptions();

opt.defines('nrcities', '20');

opt.defines('cmax', int2str(cMax - 1));

opt.defines('컷데이터', cutData.name);


%는에서 입력 데이터를 읽었습니다.'무료 슬롯gdx'

opt.defines('tspdata', ['"', 전체 파일(ws.systemDirectory,'API파일', '데이터', '무료 슬롯gdx'), '"']);

tspJob.run(opt, cp);


%는 n 목록을 한 번만 채웁니다.

            ifnumel(n) == 0

                fori = tspJob.outDB.getSet('나').기록

nend+1 = i1.key(1);

끝

끝


% 수정자 miAcut 및 miRhscut을 사용하여 모델 인스턴스를 인스턴스화

mi = cp.addModelInstance();

miAcut = mi.syncDB.addParameter('아컷', 3,'');

miRhscut = mi.syncDB.addParameter('rhscut', 1,'');

수정자 = 무료 슬롯control.Modifier(miAcut), 무료 슬롯control.Modifier(miRhscut);

mi.instantiate('mip min z 사용 할당', 수정자);

        그밖에

fprintf('.');

끝


% 모델 인스턴스 해결사용 중업데이트 유형을 누적하고 이후에 acut 및 rhscut 삭제

mi.solve(무료 슬롯control.globals.SymbolUpdateType.ACCUMULATE);

mi.syncDB.getParameter('아컷').clear();

mi.syncDB.getParameter('rhscut').clear();


%는 솔루션에서 그래프 정보를 수집합니다.

그래프 = 컨테이너.Map();

방문하지 않음 = n;

        fori = n

            forj = n

                ifmi.syncDB.getVariable('x').findRecord(i1,j1).level > 0.5

그래프(i1) = j1;

끝

끝

끝


% 모든 하위 투어를 찾고 acut 및 rhscut을 수정하여 필요한 컷을 추가합니다.

        그동안numel(방문하지 않음) ~= 0

ii = 방문하지 않음1;

하위 투어 = ii;

            그동안~strcmp(graph(ii), 방문하지 않음1)

ii = 그래프(ii);

하위 투어end+1 = ii;

끝


            fori = numel(notVisited):-1:1

for j = 1:numel(subTour)

if strcmp(notVisitedi, subTourj)

방문하지 않음(i) = [];

                        휴식

끝

끝

끝


% 두 데이터베이스에 컷 추가(cutData, mi.SyncDB)

for i = subTour

j = subTour의 경우

키 = ['c', int2str(curCut)], i1, j1;

rec = acut.addRecord(keys);

rec.값 = 1;

rec = miAcut.addRecord(keys);

rec.값 = 1;

끝

끝


키 = ['c', int2str(curCut)];

rec = rhscut.addRecord(키);

rec.value = numel(subTour) - 0.5;

rec = miRhscut.addRecord(키);

rec.value = numel(subTour) - 0.5;

cc.addRecord(키);

curCut = curCut + 1;

끝


        ifnumel(subTour) >= numel(n)

            휴식

끝

끝


fprintf('\nz=%g\n', mi.syncDB.getVariable('z').record.level);

fprintf('sub_tour: ');

    fori = 하위 투어

fprintf('%s -> ', 나1);

끝

fprintf('%s\n', 하위 둘러보기1);

끝