Математические модели задач и их решение на ЭВМ

ЗАДАНИЕ № 1

Из пункта А в пункт Б ежедневно отправляются пассажирские и скорые поезда. Наличный парк вагонов разных типов, из которых ежедневно можно комплектовать данные поезда, и количество пассажиров вмещающихся в каждом вагоне приведены в таблице.







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

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







В данном случае неизвестными являются число скорых и пассажирских поездов Х1 и Х2

Составим математическую модель этой задачи.

Максимальное число пассажиров перевозимых данными поездами обозначим L. Тогда целевая функция будет иметь вид:



L= 0*(1*х1+1*х2)+58*(5*х1+8*х2)+40*(6*х1+4*х2)+32*(3*х1+1*х2) – max




Ограничение на искомое решение следующее:



1*х1+1*х2

5*х1+8*х2

6*х1+5*х2

3*х1+1*х2

Х1+х2<=10















ЗАДАНИЕ №2.



1. решить задачу геометрическим методом.

2. составить двойственную задачу для исходной.



2х1+5х2≥10

5х1+2х2≥10

3х1+4х2≤24

4х1+3х2≤24

Х1-2х2 ≤4

Z=3х1+х2→мах

Х1≥0;Х2≥ 0.

Х1+5x2>5

5x1+x2>5

X1+X2<7

3x1-4x2<12

-4x1+3x2<12

Z=4x1-3x2 – max

X1>0 X2>0



РЕШЕНИЕ



1. Поскольку рассматривается задача на максимум, то все ограничения следует привести к виду «≤». Для этого обе части первого и второго неравенств следует умножить на «-1». Получим: -  -2х1-5х2≤-10



-5х1-2х2≤-10

3х1+4х2≤24

4х1+3х2≤24

Х1≥0;Х2≥ 0.



2. Составим расширенную матрицу системы.



 -2 -5 -10

 -5 -2 -10

А1= 3 4 24

 4 3 24

 3 1 Z

3. Найти матрицу А1т, транспонированную к А1.



 -2 -5 3 4 3

А1т = -5 -2 4 3 1

 -10 -10 24 24 Z



4. Сформулируем двойственную задачу:



Z= -10у1 -10у2 +24у3 +24у4 → min.

-2 у1 - 5 у2 + 3 у3 + 4 у4≥3

-5у1 - 2 у2 + 4 у3 + 3 у4≥1

у1 ≥0; у2≥0; у3≥0; у4 ≥0.



ЗАДАНИЕ №3

Составить математическую модель задачи и решить ее на ЭВМ.

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

Данные для каждого варианта приведены

1.тарифы перевозок единицы груза от каждого поставщика каждому потребителю

2.запасы груза каждого поставщика

3.потребности в грузе каждого потребителя.











РЕШЕНИЕ



А1 + А 2 + А 3 + А 4 + А 5 = 30+20+10+27+30=117

В1 + В2 + В 3 + В 4 =30+40+50+10=130



Спрос превышает предложение и поэтому добавляем пятого фиктивного постивщика.130-117=13 Отсюда:



Х11+Х12+Х13+Х14+Х15 = 30

Х21+Х22+Х23+Х24+Х25 = 20

Х31+Х32+Х33+Х34+Х35 = 10

Х41+Х42+Х43+Х44+Х45 = 27

Х51+Х52+Х53+Х54+Х55 = 30

Х61+Х62+Х63+Х64+Х65=13

F = 7Х11+8Х12+5Х13+5Х14+5Х15+9Х16+1Х21+

+4Х22+2Х23+5Х24+9Х25+ 3Х31+5Х32+3Х33+8Х34+7Х35

+9Х36+2Х41+8Х42+7Х43+4Х44+5Х45+9Х46min.
























ЗАДАНИЕ №4

Представители одной фирмы могут принять по три стратегии. Матрица эффективности стратегий фирм представлена в таблице.

Определить верхнюю и нижнюю цену игры.

Найти седловую точку. В случае ее отсутствия составить двойственные задачи мат.програмирования.



К\С

С 1

С 2

С 3

К 1

1

7

2

К 2

5

4

8

К 3

4

6

3

K 4

1

3

2




РЕШЕНИЕ



Нижняя цена игры вычисляется α = maxi minj hij = maxi βj , где αi - наименьшее значение в i-той строке.







Верхняя цена игры вычисляется β = minj maxi hij = minj βj , где βj = =maxi hij - наибольшее значение в j-том столбце.



К\С

С 1

С 2

С 3



αi

К 1

3

7

3



3

К 2

8

1

5



1

К 3

2

6

4



2









α=

1

βj

8

7

5

β=

8



Седловая точка отсутствует, значит нужно составить двойственную задачу.



ЗАДАНИЕ №5

Имеются данные эффективности выпуска новой продукции при различных вариантах решений (стратегий) и различных состояниях среды (природы), таблица 1. Выбрать наилучшее решение, стратегию используя критерии:

Максимакса

Вальда

Сэвиджа

Гурвица (коэффициент пессимизма р=0,3)

Байеса (вероятности для каждого состояния среды р1=0,2, р2=0,3, р3=0,3, р4=0,2)

Лапласа



ТАБЛИЦА 1.

ВАРИАНТЫ РЕШЕНИЙ

СОСТОЯНИЕ ПРИРОДЫ



П1

П2

П3

П4

А1

7

13

9

15

А2

15

8

11

12

А3

12

6

13

10

А4

11

10

15

14

А5

8

15,5

12

15







РЕШЕНИЕ



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



М = maximaxj hij = maxi Mi



Находим М=maxi hij, табл.2, т.е.максимальное значение в i-той строке.








ТАБЛИЦА 2.





М1= 15, М2= 15, М3=13, М4= 15, М5= 15,5.



Максимальное значение М = maxi Mi = 15,5, значит решение А5 оптимально.

Согласно критерию Вальда наиболее предпочтительным является решение, при котором W = maximinjhij = maxi Wi. Находим Wi = minjhij, т.е. минимальное значение W в i-той строке.












Максимальное значение W=10, следовательно решение А4 является наилучшим.

3. В соответствии с критерием Сэвиджа предпочтение отдается решению, для которого максимальные потери при различных вариантах обстановки окажутся минимальными, т.е. достигается значение:

S = minimaxj rij = miniSi.

Найдем матрицу потерь (табл.4 и 5): βj = maxi hij; rij = βj - hij.



ТАБЛИЦА 4. ВЫИГРЫШИ









ТАБЛИЦА 5. ПОТЕРИ.










Минимальное значение S = 7. Следовательно оптимальным решением является решение А5.

По критерию Гурвица предпочитается то решение, при котором G = maxi { minihij + (1- p) maxj hij } = maxi Gi.

Находим Gi = pWi + (1-p)Mi, р=0,3 по условию задачи.











Находим Gmax = 17,4 значит решение А2 является оптимальным.

Согласно критерию Байеса наилучшим является решение, при котором достигается максимум математического ожидаемого выигрыша (или минимум среднеожидаемого риска).

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

р1=0,2, р2=0,3, р3=0,3, р4=0,2. Определяем математическое ожидание выигрышей по каждому решению: МВ1 = ∑рihij.









Определяем максимум ожидаемого математического выигрыша. Он равен 12,85, что соответствует четвертому решению, которое, следовательно, и является оптимальным.

Определяем среднеожидаемый риск по каждому решению.



МРi = ∑pj rij











Определяем минимум среднеожидаемого риска. Он равен 2,3, что соответствует пятому решению, которое, следовательно, является оптимальным по данному критерию.

Определяем значения для каждого решения по критерию Лапласа.




ВЫИГРЫШИ:







Максимальный выигрыш составит 12,625 что соответствует 2-ому оптимальному решению.



ПРОИГРЫШИ:









Минимальный проигрыш составит 2,5, что соответствует 5-ому оптимальному решению.



ЗАДАНИЕ №6.

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

Построить линейную однофакторную модель зависимости расходов на питание от дохода семьи.

Определить коэффициент корреляции и оценить тесноту связи между доходами семьи и расходами на питание.

Определить коэффициент детерминации и коэффициент эластичности, объяснить их смысл.

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



Доходы семьи (х), тыс.грн.

2.2

3,6

4,2

5,8

6,7

7,9

8,6

10,6

Расходы на продукты (у)

1,2

2,0

2,6

2,9

3,1

3,9

4,5

5



РЕШЕНИЕ. Подготовим вспомогательную таблицу:



Табл 1





Табл 2







По формуле определим коэффициенты а0, и а1.



А0= ∑уi*∑xi^2-∑xiyi*∑xi / n*∑x^2-∑xi*∑xi

Ai=n*∑xiyi-∑xi*∑yi /n*∑x^2-∑xi*∑xi.









Тогда регрессионная модель, согласно формуле, запишется: 



Y^=А0+Аi*x



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







Для полученной модели определим:

А) коэффициент корреляции по формуле и оценим тесноту связи между доходами семьи и расходами на питание.



Xcp=∑xi/n Ycp=∑yi/n XYcp=∑xiyi/n 




Для этого вычислим средние значения доходов и расходов при помощи EXCEL. Расчеты приведены в табл 2



3. Хср= 49.6/8 = 6.2; Уср= 25.2/8 = 3.2 XcpУср=180,9/8 = 22,6.



Для вычисления среднеквадратических ошибок Sy, Sx имеем формулу:



Sy=√∑(yi-y^i)/n Sx=√∑(xi-x^)^2/n







Коэффициент корреляции вычислим по формуле:



rxy=xy^-x^*y^/sy*sx











3. Рассчитаем коэффициент детерминации: R2xy = 0,972111224. Значит, 97,2% величины расходов семьи на питание зависит от изменения доходов семьи, а остальные 2,8% связаны с изменением других, не включенных в модель факторов.









Вычислим коэффициент эластичности:



Эху=aix^/y^











С увеличением доходов семьи на 1% расходы на питание увеличатся в среднем на 0,8781%.

Найдем среднюю по модулю линейную относительную ошибку аппроксимации по формуле: d=1/n*∑(yi-y^) 











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



ЗАДАНИЕ №7.

По исходным данным из задачи 6 рассчитаем Se, Sa0, Sa1 по формулам. Для этого подготовим таблицу:














Se = √1/n-2*∑e^2

Sa0=Se*√ ∑x ^2/∑(xi-x^)^2 

Sa1 = Se*√ 1/∑(x-x^)^2



Согласно задаче имеем:

А0 = 0,3837079 А1 = 0,4461762. для вычисления фактических значений t-критерия воспользуемся формулами: ta0 = a0/ Sa0 = 1.84707; ta1 = 14,4617.

По таблице 1 приложения А найдем табличное значение t-критерия для степеней свободы df = 8-1-1 = 6 и уровня зависимости 6%,т.е. tтабл = 1,943.

При уровне значимости 6% имеет место неравенство: 

ta1 = 0,073525 ‹ tтабл = 1,943. Значит, с уверенностью 94% можно утверждать, что оценка А1 = 0,747263097 не является статистически значимой.

Аналогично проверим для другого параметра. ta0 = 1,743736 ‹ tтабл = 1,943, значит оценка А0 = 0,123251901 также не является статистически значимой.












Значимость уравнения регрессии в целом и коэффициента тесноты связи R2 определяется с помощью критерия Фишера. Значение оценки R2 получено в предыдущей задаче, R2 = 0,968583448. Фактическое значение Fфакт определяем по формуле: Fфакт = 184,9821.

Табличное значение Fтабл определяем по таблице: Fтабл = 5,99.

Поскольку Fфакт = 184,9821› Fтабл = 5,99, то с уверенностью 94% делается заключение о том, что уравнение регрессии в целом статистически значимо и статистически значим показатель степени связи R2, т.е. отвергается нулевая гипотеза о том, что R2 = 0.



ЗАДАНИЕ №8.

Имеются следующие исходные данные:

Годы

1997

1998

1999

2000

2001

2002

2003

2004

Объем реализации

10,84

11,12

10,6

11,31

11,62

12,0

12,73

11,12



Коэффициент достоверности аппроксимации для каждого типа линии тренда

Линейная у= 0,1795х – 347,71 R^2=0.4163





Логорифмическая у=359,19 Ln(x)-2718,8 R^2=0.1464







Степенная y=3E-102x^31.059 R^2=0.422







Экспонтенциальная у=4Е-13е^0.01558x R^2=0.4218

Как видно из рисунка в 2005г в сравнении с 2004г в среднем реализация продукции увеличилась на 0,42 млн. грн.


ЗАДАНИЕ №9.

Имеются данные испытаний нескольких величин по результатам обследования десяти статистически однородных филиалов фирмы, приведенные в таблице. х1- фондовооруженность, х2 – энерговооруженность, у – производительность труда.

Выполнить следующее:

Построить линейную регрессионную модель при помощи ПЭВМ.

Выполнить команду «Регрессия».

Определить по результатам команды «Регрессия» значение коэффициента множественной корреляции и детерминации.

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

Проверить статистическую значимость оценки степени достоверности взаимосвязи R2 и всей модели в целом.



РЕШЕНИЕ.

1. построить регрессионную модель.

2. выполнить команду «Регрессия», результаты которой показаны ниже.









Рис. Результаты команда «Регрессия»



Регрессионная модель принимает вид:



у^ = 0929087*2,9+ - 0,4502*4,5-3,246374



3. Согласно Рис коэффициенты множественной корреляции и детерминации,в данном случае R = 0,993689; R2 = 0,98742.

4. Статистическую значимость оценок параметров модели b,a1,а2 осуществим с помощью t-критерия. Для этого определим его табличное значение и его фактические значения для каждого из оцениваемых параметров. По таблице 1 приложения А при уровне значимости 1% найдем табличное значение t-критерия для степеней свободы df = 10-2-1 = 7 и уровня зависимости 7%,т.е. tтабл = 3,143.

Фактическое значение t-критерия для каждого из оцениваемых параметров смотрим на рисунке в столбце t-статистика в нашем случае:



t-a1= 15,73834 ta2= - 0,855361 tb=15,97697 



При уровне значимости 7% t-a1= 15,73834> tтабл имеет место равенство: Значит, с уверенностью 99% можно утверждать, что оценка А1 параметра модели является статистически значимой.

Условие ta2 = -0,855361< tтабл = 3,143 не выполняется, значит утверждаем, что этот критерий статистически не важен.

Условие tв = 15.97697> tтабл = 3,143 выполняется, значит и эта оценка статистически значима в модели.

5. Значимость уравнения регрессии в целом и коэффициента тесноты связи R2 определяем с помощью критерия Фишера. Фактическое Fфакт =274,684752

Табличное значение Fтабл определяем по таблице: Fтабл = 9,55. Условие Fфакт =274684752> Fтабл = 9.55 выполняется, поэтому с вероятностью 99% делается заключение о том, что R2 статистически значим, и уравнение регрессии в целом значимо, т.е. отвергается нулевая гипотеза R2 = 0.