Отделите модель оптимизации от данных

Чтобы получить масштабируемую, переиспользуемую задачу оптимизации, создайте задачу таким образом, чтобы отделить данные задачи от структуры модели.

Предположим, что у вас есть проблема многопериодного планирования с несколькими продуктами. Периоды времени указаны в векторе, periodsи продукты в векторе строки, products.

periods = 1:10;
products = ["strawberry","cherry","red grape",...
    "green grape","nectarine","apricot"];

Чтобы создать переменные, которые представляют количество продуктов, используемых в каждом периоде, используйте операторы, которые берут размеры из данных. Для примера:

usage = optimvar('usage',length(periods),products,...
    'Type','integer','LowerBound',0);

Чтобы позже изменить периоды времени или продукты, необходимо изменить данные только в periods и products. Затем можно запустить тот же код для создания usage.

Другими словами, чтобы поддерживать гибкость и допускать повторное использование, не используйте оператор, который имеет жестко закодированные размеры данных. Для примера:

usage = optimvar('usage',10,6,... % DO NOT DO THIS
    'Type','Integer','LowerBound',0);

Такое же внимание уделяется выражениям, а также переменным. Предположим, что затраты на продукты указаны в матрице данных costs, размера length(periods)-by- length(products). Чтобы симулировать допустимые данные, создайте случайную целочисленную матрицу соответствующего размера.

rng default % For reproducibility
costs = randi(8,length(periods),length(products));

Лучшая практика состоит в том, чтобы создавать выражения затрат, которые берут размеры из данных.

costPerYear = sum(costs.*usage,2);
totalCost = sum(costPerYear);

Таким образом, если вы когда-либо измените размеры данных, операторы, которые создают costPerYear и totalCost не изменяйте. Другими словами, чтобы поддерживать гибкость и допускать повторное использование, не используйте оператор, который имеет жестко закодированные размеры данных. Для примера:

costPerYear = optimexpr(10,1); % DO NOT DO THIS
totalcost = 0;
for yr = 1:10 % DO NOT DO THIS
    costPerYear(i) = sum(costs(i,:).*usage(i,:));
    totalcost = totalcost + costPerYear(i);
end

Похожие темы