Математическое программирование и моделирование в экономике и управлении

Математическое программирование и моделирование в экономике и управлении

Министерство образования РФ

Санкт-Петербургская Лесотехническая академия им. С. М. Кирова

Кафедра: математических методов и моделирования в экономике и управлении

Курсовая работа по математическому программированию и моделирования в

экономике и управлении.

Выполнила: студентка ФЭУ, II курса, 4 группы

д/о, направление 521500

менеджмент

Гузеева Ольга

Зачётная книжка № 600033

Преподаватели: П. Н. Коробов, А. А. Моисеев

Санкт-Петербург

2002 год

Методология математического моделирования ассортиментной задачи (задачи

оптимизации программы выпуска продукции по ассортименту).

Этапы решения задач:

1. выбор проблемы решения;

2. постановка проблемы и разработка экономико-математической модели

(ЭММ);

3. выбор метода решения;

4. выполнение решения;

5. анализ результата и проведение эксперимента;

6. внедрение результата, полученного в результате опыта.

Задачи оптимизации:

1. обеспечение балансовой увязки между знаниями по выпуску продукции

разных видов и наличием производственных ресурсов (сырьё, материалы,

машинное время, трудовые ресурсы, энергия и т. п.);

2. обеспечение максимального экономического эффекта при использовании

производственных ресурсов;

3. проведение эксперимента (повторы решения при изменённых условиях,

чтобы выработать альтернативные варианты и выбрать из них наиболее

приемлемый).

Под оптимизацией программы выпуска продукции по ассортименту понимаются

такие объёмы выпуска различной продукции, которые обеспечивают получение

максимального экономического эффекта от реализации всей продукции.

Условия задачи: на предприятии имеются свободные ресурсы: сырьё, материалы,

машинное время, трудовые и т. п. В условии задачи известны фонды

производственных ресурсов на планируемый период, нормы их затрат на единицу

(десяток, сотню или комплект продукции), а также известны показатели

прибыли от реализации продукции. Найти программу выпуска продукции по

ассортименту, обеспечивающую максимальную суммарную прибыль от её

реализации.

|Виды |Фонды |Нормы затрат производственных |

|производственных |производственных |ресурсов на единицу продукции |

|ресурсов |ресурсов на | |

| |планируемый период | |

| | |Р1 …………….. Рj ……………… Рn |

|1 |bj | |

|. |. | |

|. |. | |

|. |. | |

|r |br |A=[arj]Rx n |

|. |. | |

|. |. | |

|. |. | |

|R |bR | |

|Критерий оптимальности |с1 ……………… сj ………………. cn |

j – индекс вида продукции;

Pj – виды продукции;

r – индекс вида производственных ресурсов (от 1 до R);

br – фонд r-производственного ресурса;

arj – норма затрат rj-производственного ресурса;

cj – критерий оптимальности; его сущность заключается в том, что это

экономический, технико-экономический показатель, который заложен в условии

задачи для суждения об оптимальности её решения;

xj –количество продукции Pj.

Х=(х1, х2…хj…xn) – оптимальная программа выпуска продукции по ассортименту.

Критерий оптимальности:

[pic]

Система ограничений:

[pic]

Суммарные затраты r-производственного ресурса на выполнение всех n видов

продукции не должен превышать фонды этого ресурса, которым предприятие

владеет на планируемый период.

Экономическое содержание и математическое моделирование распределительных

нетранспортных задач.

I. Известна программа выполнения продукции на период. Эта программа может

быть выполнена на разных станках, а также известны фонд эффективного

рабочего времени каждого исполнителя, часовая производительность каждого из

исполнителей при выработке каждого вида продукции. Известны затраты по

выполнению продукции у разных исполнителей.

i – индекс исполнителя (отдельной машины, рабочего, цеха, участка),

i=1,2…m;

j – индекс вида продукции (работы), j=1,2…n;

m – количество рабочих (станков);

n – число видов продукции (работ);

bi – фонд эффективного рабочего времени i-исполнителя в планируемом периоде

в часах;

?ij – часовая производительность j-продукции у i-исполнителя;

?=[ ?ij]mxn – известно;

sij – себестоимость производства единицы j-продукции у i-исполнителя;

S=[ sij] mxn – известно;

Pj – вектор показателей, которые характеризуют объёмы выпуска продукции

(выполнения работ) по всем видам – известно;

|Наименование |Фонд эффективного |P1 ………………… Pj …………………. Pn |

|исполнителя |рабочего времени | |

| | |производительность / себестоимость |

|1 |b1 | |

|. |. | |

|. |. | |

|. |. | |

|i |bi |?=[ ?ij]mxn / S=[ sij] mxn |

|. |. | |

|. |. | |

|. |. | |

|m |bm | |

Найти план распределения производственного задания по выпуску продукции

(выполнения работ) между исполнителями, при котором задание было бы

выполнено с минимальными суммарными затратами.

xij – затраты эффективного рабочего времени у i-исполнителя на произведение

j-продукции;

Х=[ xij]mxn – искомые величины.

Целевая функция:

[pic]

s’ij – себестоимость часового объёма выпуска продукции определённого вида

на определённом оборудовании.

Система ограничений:

[pic] – суммарные затраты эффективного рабочего времени на выполнение всех

видов работ не должен превышать фонда, которым располагает i-рабочий в

плановом периоде;

[pic] – суммарный объём выпущенной продукции j-вида у всех m исполнителей

должен быть равен производственному заданию;

[pic]

II. На предприятии известна программа выпуска продукции по видам, которая

может быть выполнена разными исполнителями (на разных участках). В условии

задачи известны: фонд эффективного рабочего времени каждого исполнителя в

плановом периоде, показатели норм затрат эффективного рабочего времени на

производство различных видов продукции на разном оборудовании, а также

прибыль от реализации единицы продукции, выработанной разными

исполнителями.

|Наименование |Фонд эффективного |P1 ………………… Pj …………………. Pn |

|исполнителя |рабочего времени | |

| | |нормы затрат / прибыль |

|1 |b1 | |

|. |. | |

|. |. | |

|. |. | |

|i |bi |A=[ aij]mxn / C=[ cij] mxn |

|. |. | |

|. |. | |

|. |. | |

|m |bm | |

i – индекс исполнителя (отдельной машины, рабочего, цеха, участка),

i=1,2…m;

j – индекс вида продукции (работы), j=1,2…n;

m – количество рабочих (станков);

n – число видов продукции (работ);

bi – фонд эффективного рабочего времени i-исполнителя в планируемом периоде

в часах;

aij – показатель нормы затрат на производство j-продукции у i-исполнителя;

A=[ аij]mxn – известно;

сij – показатель прибыли от единицы j-продукции у i-исполнителя;

С=[ сij] mxn – известно;

Pj – вектор показателей, которые характеризуют объёмы выпуска продукции

(выполнения работ) по всем видам – известно.

Требуется найти план распределения производственного задания между

исполнителями, при котором это задание было бы выполнено с максимальной

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

xij – объём (количество) j-продукции выработанной i-исполнителем;

Х=[ xij]mxn – искомые величины.

Целевая функция:

[pic]

Система ограничений:

[pic]

При решении этой системы линейных уравнений и неравенств, нужно найти такие

неотрицательные значения переменных, чтобы целевая функция принимала

максимальное значение.

Методология математического моделирования раскройной задачи (задачи

оптимизации программы раскроя материалов).

Пусть имеются ДСП стандартных размеров, из которых необходимо нарезать

m различных по размеру заготовок и деталей для производства мебели. ДСП

определённого размера может быть раскроена n способами (вариантами). По

каждому из возможных вариантов раскроя составляется соответствующая карта

раскроя, из которой видно, что при j (j=1,2…n) способе раскроя из одной

плиты получается определённое количество (обозначим через aij) заготовок i

(i=1,2…m) вида (размера). По картам раскроя устанавливается также величина

отходов (площадь, вес, стоимость) при раскрое одной плиты j способом

(обозначим – сj). В задании на раскрой должно быть указано общее количество

заготовок каждого i вида (размера) – bi, которое необходимо нарезать из

плит, поступивших в раскрой (обозначим – R). В задаче требуется определить

оптимальный план раскроя ДСП, обеспечивающий минимальные отходы (или

минимальный расход раскраиваемых материалов), при условии выполнения

задания по выходу заготовок.

xj – количество ДСП, которое следует раскраивать с тем, чтобы нарезать

заданное число заготовок каждого вида, при этом суммарные отходы (или

суммарный расход плит) должны быть минимальными.

|Виды заготовок|Задание по раскрою |Способы раскроя |

| | |1 ……………………. j ………………….. n |

| | | |

|1 |b1 | |

|. |. | |

|. |. | |

|. |. | |

|i |bi |A=[ аij]mxn |

|. |. | |

|. |. | |

|. |. | |

|m |bm | |

|Отходы |C=[ cj] n |

Критерий оптимальности:

[pic]

Система ограничений:

[pic]

При решении этой системы линейных уравнений и неравенств, нужно найти такие

неотрицательные значения переменных, чтобы целевая функция принимала

минимальное значение.

Рассмотрим пример решения задачи оптимизации программы раскроя материалов

симплексным методом.

F=0.26x1+0.28x2+0.3x3+0.29x4=min

[pic]

F=0.26x1+0.28x2+0.3x3+0.29x4+0x5+0x6+0x7+0x8+0x9+M(y1+y2+y3+y4)=min

[pic]

| |B1 ………………………….. Bj ………………………………….. Bn |

| |b1 …………………………… bj ………………………………….. bn |

| |С=[ сij] mxn / Х=[ xij]mxn |

|A1 |a1 |c11|……………………. |c1j|…………………. |c1n|……………… |

| | | |x11………………… | |………x1j……… | |………….. x1n |

| | | | | | | | |

|. |. |. |. . |. |. . |. |. . |

|. |. |. |. |. |. |. |. |

|. |. |. |. . |. |. . |. |. . |

| | | |. | |. | |. |

| | | |. . | |. . | |. . |

| | | |. | |. | |. |

|Ai |ai |ci1|……………………. |cij|…………………. |cin|……………… |

| | | |xi1………………… | |………xij……… | |………….. xin |

| | | | | | | | |

|. |. |. | |. | |. | |

|. |. |. | |. | |. | |

|. |. |. | |. | |. | |

|Am |am |cm1|……………………. |c11|…………………. |c11|……………… |

| | | |xm1………………… | |………xmj……… | |…………..xmn |

| | | | | | | | |

Целевая функция:

[pic] (1)

Условие реализации продукции у каждого из поставщиков:

[pic] (2)

Условие обеспечения всех потребителей продукцией по их потребности:

[pic] (3)

Условие не отрицательности переменных:

[pic]

В решении системы линейных уравнений 2 и 3 необходимо найти такие не

отрицательные значения переменных, чтобы целевая функция принимала

минимальное значение.

m+n-1 – линейно независимых уравнений, ранг системы, r= m+n-1.

В каждом опорном плане должно быть m+n-1 базисных элементов (xij>0), если

таких переменных равно или больше, чем m+n-1, план называется

невырожденный; если одна или несколько базисных переменных равна нулю, то

такой план считается вырожденным.

Открытые транспортные задачи.

a) [pic]

[pic] (1)

[pic] (2)

[pic] (3)

[pic]

Bn+1: [pic] – потребность какого-то потребителя, находящегося за пределами

района (фиктивный потребитель).

[pic] (1)

[pic] (2)

[pic] (3)

[pic]

[pic]

сi, n+1=0 (i=1,2…m)

б) [pic]

[pic] (1)

[pic] (2)

[pic] (3)

[pic]

Аn+1: [pic] – фиктивный поставщик.

[pic] (1)

[pic] (2)

[pic] (3)

[pic]

[pic]

Ограничение транспортных возможностей.

а) xij=0 => cij=М, где М»0;

б) 0 ? хij ? dij

dij – характеризует транспортные возможности между i-поставщиком и j-

потребителем.

Тогда поставщик Аi условно делится на Аi` и Аi``, при этом ai`=dij и ai``=

ai`-dij, cij`=cij и cij``=М, где М»0.

Рассмотрим пример решения транспортной задачи методом потенциалов.

| |В1 |В2 |В3 |В4 |В5 |Ui |

| |200 |250 |275 |255 |120 | |

?11=-1

?12=0

?13=M-11

?21=6

?24=7

?25=5

?31=6

?33=M-6

?34=7

?35=5

?41=0

?42=-4

?44=4

?52=13

?54=0

?55=2

| |В1 |В2 |В3 |В4 |В5 |Ui |

| |200 |250 |275 |255 |120 | |

?11=-1

?12=4

?13=M-11

?21=6

?22=4

?24=7

?25=5

?31=2

?33=M-10

?34=3

?35=1

?41=0

?44=4

?52=7

?54=3

?55=2

| |В1 |В2 |В3 |В4 |В5 |Ui |

| |200 |250 |275 |255 |120 | |

?12=5

?13=M-10

?15=1

?21=6

?22=4

?24=6

?25=5

?31=2

?33=M-10

?34=2

?35=1

?41=0

?44=3

?52=7

?54=2

?55=2

F=7x1+10x2+Mx3+6x4+7x1+10x2+Mx3+6x4+9x5+5x6+6x7+8x8+8x9+6x10+11x11+

+10x12+6x13+11x14+9x15+7x16=min

при ограничениях:

[pic]

F=7*45+6*155+5*125+6*125+6*125+11*25+9*125+6*255=6300

Оптимальный план поставок для деревообрабатывающих предприятий,

обеспечивающий минимальные транспортные затраты в сумме 6300000 руб.,

заключается в следующем:

1-ое лесозаготовительное предприятие поставляет 45 т. м3 1-ому

деревообрабатывающему предприятию;

1-ое – 4-ому: 255 т. м3;

2-ое – 2-ому: 125 т. м3;

2-ое – 3-ему: 125 т. м3;

3-е – 2-ому: 125 т. м3;

3-е – 3-ему: 25 т. м3;

у 3-го предприятия остаётся запас в 120 т. м3;

4-е – 1-ому: 155 т. м3;

4-е – 3-ему: 125 т. м3;

имеется альтернативный приведённому план поставок при тех же транспортных

издержках:

1-ое – 4-ому: 255 т. м3;

2-ое – 2-ому: 125 т. м3;

2-ое – 3-ему: 125 т. м3;

3-е – 1-ому: 25 т. м3;

3-е – 2-ому: 125 т. м3;

у 3-го предприятия остаётся запас в 120 т. м3;

4-е – 1-ому: 130 т. м3;

4-е – 3-ему: 150 т. м3.

Оптимизация замены оборудования. Динамическое программирование в

планировании производством и управлении им.

Под динамическим программированием понимается вычислительный метод,

опирающийся на аппарат рекуррентных соотношений.

Динамическое программирование – планирование многошагового процесса,

при котором на каждом шаге решения, оптимизируется только этот шаг. Идея

динамического программирования заключается в том, что отыскание множества

переменных, что имело место в линейном программировании, заменяется на

многократное отыскание одной или очень небольшого числа исходных

переменных.

Весь процесс динамического программирования планируется в виде

составления функциональных уравнений, которые решаются на каждом шаге.

Под функциональными уравнениями понимаются такие уравнения, в которых

выражается функциональная зависимость между множеством функций – это

сущность и отличие динамического программирования от линейного.

Содержание проблемы и сущность алгоритма решения.

Процесс решения задачи осуществляется следующим способом. Берётся

период в N лет. К этому времени оборудование отработало некое количество

лет и пришло t0 возраста.

Решение задачи начинается с последнего N-го года, составляется пара

функциональных уравнений в предположении, что пришло старое оборудование

без замены:

1) рассчитывается доход от эксплуатации оборудования при замене;

2) рассчитывается доход от эксплуатации оборудования в течение

года при условии его старения.

Вторая гипотеза: к N-ому году оборудование могло прийти замененным в каком-

то году, тогда составляется пара уравнений, в которых определяется доход за

год от эксплуатации единицы оборудования при условии замены или сохранения

оборудования.

Шаг второй: рассматриваем (N-1) год.

Рассматриваются две гипотезы:

. пришло старое оборудование без замены;

. пришло оборудование, которое было заменено.

Шаг третий: рассматривается (N-2) год при двух гипотезах, составляются

уравнения, рассчитывается доход.

Решение продолжается по всем шагам. На первом году будет одна гипотеза,

что пришло старое оборудование, используемое t0 лет.

Составление функциональных уравнений.

Под критерием оптимальности может быть принят любой экономический

показатель, если он хорошо подготовлен, т.е. он должен быть отчищен от

факторов, не зависящих от работы оборудования.

r(t) – стоимость продукции, созданной единицей оборудования возраста t лет

за год.

U(t) – затраты на содержание в течение года единицы оборудования возраста t

лет.

С(t) – затраты на замену единицы оборудования возраста t лет (затраты на

приобретение, отладку за вычетом остаточной стоимости старого

оборудования).

i – год установки нового оборудования.

Доход замены оборудования рассчитывается:

f’=r(t)-U(t)-C(t)

Доход от сохранения оборудования:

f’’=r(t)-U(t)

Если f’>f’’, то оборудование необходимо заменить, если f’?f’’ – оставить.

Шаг 1-ый: N-ый год.

Гипотеза 1: пришло старое оборудование возраста N+t0 лет.

Тогда доход за N-ый год при условии замены или сохранения оборудования:

[pic]

Гипотеза 2: пришло новое оборудование.

[pic]

Возьмём N-t-й год:

[pic]

Шаг 2-ой: (N-1)-ый год.

Рассчитывается суммарный условный доход, при условии замены или сохранения.

Гипотеза 1: пришло старое оборудование.

[pic]

Гипотеза 2: пришло новое оборудование.

[pic]

Рассмотрим пример решения задачи о замене оборудования.

Исходная информация по старому оборудованию (t0=7):

|Показатель |Значение показателей на единицу оборудования возраста (лет) в |

| |тыс. руб. |

| |8 |9 |10 |11 |12 |

|r(t) |100 |87 |74 |67 |60 |

|U(t) |27 |32 |39 |42 |40 |

|C(t) |135 |148 |150 |165 |172 |

Исходная информация по новому оборудованию:

|Показатель |Значение показателей на единицу оборудования возраста (лет) в |

| |тыс. руб. |

| |0 |1 |2 |3 |4 |

|r1(t) |135 |105 |85 |80 |75 |

|U1(t) |12 |15 |20 |22 |25 |

|C1(t) |- |152 |160 |170 |180 |

|r2(t) |125 |100 |90 |84 | |

|U2(t) |13 |15 |17 |20 | |

|C2(t) |- |132 |142 |152 | |

|r3(t) |136 |120 |116 | | |

|U3(t) |15 |16 |19 | | |

|C3(t) |- |156 |162 | | |

|r4(t) |145 |135 | | | |

|U4(t) |20 |17 | | | |

|C4(t) |- |180 | | | |

|r5(t) |162 | | | | |

|U5(t) |35 | | | | |

I этап (5 год):

|Оборудова-н|Возраст|Условие максимального дохода за 5 год |

|ие | | |

| | |Формула |Расчёты |Полити-к|

| | | | |а |

|Ст |12 |[pic] |165-35-172= -42| |

| | | | |Сохран. |

| | | |60-40=20 | |

|Н |1 |[pic] |165-35-180= -50| |

| | | | |Сохран. |

| | | |135-17=118 | |

|Н |2 |[pic] |130-162= -32 | |

| | | |116-19=97 |Сохран. |

|Н |3 |[pic] |130-152= -22 | |

| | | |84-20=64 |Сохран. |

|Н |4 |[pic] |130-180= -50 | |

| | | |75-25=50 |Сохран. |

II этап (4, 5 год):

|Оборудова-н|Возраст|Условие максимального дохода за 5 год |

|ие | | |

| | |Формула |Расчёты |Полити-к|

| | | | |а |

|Ст |11 |[pic] |=145-20-165| |

| | | |+ |Замены |

| | | |+118=78 | |

| | | | | |

| | | |=67-42+ | |

| | | |+20=45 | |

|Н |1 |[pic] |=125-156+ | |

| | | |+118=87 | |

| | | |=120-16+ |Сохран. |

| | | |+97=201 | |

|Н |2 |[pic] |=243-142= | |

| | | |=101 | |

| | | |=90-17+64= |Сохран. |

| | | |=137 | |

|Н |3 |[pic] |=243-170= | |

| | | |=73 | |

| | | |=80-22+50= |Сохран. |

| | | |=108 | |

III этап (3, 4, 5 год):

|Оборудова-н|Возраст|Условие максимального дохода за 5 год |

|ие | | |

| | |Формула |Расчёты |Полити-к|

| | | | |а |

|Ст |10 |[pic] |=136-15-150| |

| | | |+ |Замены |

| | | |+201=172 | |

| | | | | |

| | | |=74-39+ | |

| | | |+78=113 | |

|Н |1 | [pic] |=322-132= | |

| | | |=190 | |

| | | |=100-15+ |Сохран. |

| | | |+137=222 | |

|Н |2 |[pic] |=322-160= | |

| | | |=162 | |

| | | |=85-20+108=|Сохран. |

| | | | | |

| | | |=173 | |

IV этап (2, 3, 4, 5 год):

|Оборудова-н|Возраст|Условие максимального дохода за 5 год |

|ие | | |

| | |Формула |Расчёты |Полити-к|

| | | | |а |

|Ст |9 |[pic] |=125-13-148| |

| | | |+ | |

| | | |+222=186 | |

| | | | |Сохран. |

| | | |=87-32+ | |

| | | |+179=227 | |

|Н |1 | [pic] |=334-152= | |

| | | |=182 | |

| | | |=105-15+ |Сохран. |

| | | |+173=263 | |

V этап (1, 2, 3, 4, 5 год):

|Оборудова-н|Возраст|Условие максимального дохода за 5 год |

|ие | | |

| | |Формула |Расчёты |Полити-к|

| | | | |а |

|Ст |8 |[pic] |=135-12-135| |

| | | |+ | |

| | | |+263=251 | |

| | | | |Сохран. |

| | | |=100-27+ | |

| | | |+227=300 | |

Оптимальная политика отношения к оборудованию, обеспечивающая максимальную

прибыль 300 тыс. руб., заключается в следующем: в 1 год сохранить

оборудование, при этом доход составит (300-263)=37 тыс. руб.; во 2 год –

сохранить, при доходе (263-172)=91 тыс. руб.; в 3 год – заменить, при

убытке (172-201)=55 тыс. руб.; в 4 год – сохранить, при доходе (201-97)=104

тыс. руб.; в 5 год – сохраняем, при доходе 97 тыс. руб.

-----------------------

125

0

0

125

0

125

195

70

75

120

45

0

200

155

0

45

70

25

80

125

0

25

25

0

155

130

125

150