Математические модели задач и их решение на ЭВМ
Значимость уравнения регрессии в целом и коэффициента тесноты связи R2 определяется с помощью критерия Фишера. Значение оценки R2 получено в предыдущей задаче, R2 = 0,968 583 448. Фактическое значение Fфакт определяем по формуле: Fфакт = 184,9821. Имеются данные испытаний нескольких величин по результатам обследования десяти статистически однородных филиалов фирмы, приведенные в таблице… Читать ещё >
Математические модели задач и их решение на ЭВМ (реферат, курсовая, диплом, контрольная)
ЗАДАНИЕ № 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+5×2>5
5x1+x2>5
X1+X2<7
3x1−4×2<12
— 4×1+3×2<12
Z=4×1−3×2 — 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. Найти седловую точку. В случае ее отсутствия составить двойственные задачи мат.програмирования.
КС | С 1 | С 2 | С 3 | |
К 1 | ||||
К 2 | ||||
К 3 | ||||
K 4 | ||||
РЕШЕНИЕ Нижняя цена игры вычисляется б = maxi minj hij = maxi вj, где бi — наименьшее значение в i-той строке.
Верхняя цена игры вычисляется в = minj maxi hij = minj вj, где вj = =maxi hij — наибольшее значение в j-том столбце.
КС | С 1 | С 2 | С 3 | бi | ||
К 1 | ||||||
К 2 | ||||||
К 3 | ||||||
б= | ||||||
вj | в= | |||||
Седловая точка отсутствует, значит нужно составить двойственную задачу.
ЗАДАНИЕ № 5
Имеются данные эффективности выпуска новой продукции при различных вариантах решений (стратегий) и различных состояниях среды (природы), таблица 1. Выбрать наилучшее решение, стратегию используя критерии:
1. Максимакса
2. Вальда
3. Сэвиджа
4. Гурвица (коэффициент пессимизма р=0,3)
5. Байеса (вероятности для каждого состояния среды р1=0,2, р2=0,3, р3=0,3, р4=0,2)
6. Лапласа ТАБЛИЦА 1.
ВАРИАНТЫ РЕШЕНИЙ | СОСТОЯНИЕ ПРИРОДЫ | ||||
П1 | П2 | П3 | П4 | ||
А1 | |||||
А2 | |||||
А3 | |||||
А4 | |||||
А5 | 15,5 | ||||
РЕШЕНИЕ
1. По критерию максимакса наилучшим признается решение, при котором достигается максимальный выигрыш равный М = 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 оптимально.
2. Согласно критерию Вальда наиболее предпочтительным является решение, при котором 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.
3. По критерию Гурвица предпочитается то решение, при котором G = maxi { minihij + (1- p) maxj hij } = maxi Gi.
Находим Gi = pWi + (1-p)Mi, р=0,3 по условию задачи.
Находим Gmax = 17,4 значит решение А2 является оптимальным.
4. Согласно критерию Байеса наилучшим является решение, при котором достигается максимум математического ожидаемого выигрыша (или минимум среднеожидаемого риска).
Вероятности для каждого состояния среды по условию задачи таковы:
р1=0,2, р2=0,3, р3=0,3, р4=0,2. Определяем математическое ожидание выигрышей по каждому решению: МВ1 = ?рihij.
Определяем максимум ожидаемого математического выигрыша. Он равен 12,85, что соответствует четвертому решению, которое, следовательно, и является оптимальным.
Определяем среднеожидаемый риск по каждому решению.
МРi = ?pj rij
Определяем минимум среднеожидаемого риска. Он равен 2,3, что соответствует пятому решению, которое, следовательно, является оптимальным по данному критерию.
5. Определяем значения для каждого решения по критерию Лапласа.
ВЫИГРЫШИ:
Максимальный выигрыш составит 12,625 что соответствует 2-ому оптимальному решению.
ПРОИГРЫШИ:
Минимальный проигрыш составит 2,5, что соответствует 5-ому оптимальному решению.
ЗАДАНИЕ № 6.
По экспериментальным данным опроса восьми групп семей о расходах на продукты питания, в зависимости от уровня дохода семьи, приведенным в таблице, требуется:
1. Построить линейную однофакторную модель зависимости расходов на питание от дохода семьи.
2. Определить коэффициент корреляции и оценить тесноту связи между доходами семьи и расходами на питание.
3. Определить коэффициент детерминации и коэффициент эластичности, объяснить их смысл.
4. Определить среднюю по модулю относительную ошибку аппроксимации и оценить точность построенной модели.
Доходы семьи (х), тыс.грн. | 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 | ||
РЕШЕНИЕ. Подготовим вспомогательную таблицу:
Табл 1
Табл 2
1. По формуле определим коэффициенты а0, и а1.
А0= ?уi*?xi2-?xiyi*?xi / n*?x2-?xi*?xi
Ai=n*?xiyi-?xi*?yi /n*?x2-?xi*?xi.
Тогда регрессионная модель, согласно формуле, запишется:
Y^=А0+Аi*x
Построим график зависимости и отметим экспериментальные точки.
2. Для полученной модели определим:
А) коэффициент корреляции по формуле и оценим тесноту связи между доходами семьи и расходами на питание.
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=v?(yi-y^i)/n Sx=v?(xi-x^)^2/n
Коэффициент корреляции вычислим по формуле:
rxy=xy^-x^*y^/sy*sx
3. Рассчитаем коэффициент детерминации: R2xy = 0,972 111 224. Значит, 97,2% величины расходов семьи на питание зависит от изменения доходов семьи, а остальные 2,8% связаны с изменением других, не включенных в модель факторов.
Вычислим коэффициент эластичности:
Эху=aix^/y^
С увеличением доходов семьи на 1% расходы на питание увеличатся в среднем на 0,8781%.
3. Найдем среднюю по модулю линейную относительную ошибку аппроксимации по формуле: d=1/n*?(yi-y^)
Коэффициент низкий что значит точность построения модели высока.
ЗАДАНИЕ № 7.
1. По исходным данным из задачи 6 рассчитаем Se, Sa0, Sa1 по формулам. Для этого подготовим таблицу:
Se = v1/n-2*?e2
Sa0=Se*v ?x ^2/?(xi-x^)^2
Sa1 = Se*v 1/?(x-x^)^2
Согласно задаче имеем:
А0 = 0,3 837 079 А1 = 0,4 461 762. для вычисления фактических значений t-критерия воспользуемся формулами: ta0 = a0/ Sa0 = 1.84 707; ta1 = 14,4617.
По таблице 1 приложения, А найдем табличное значение t-критерия для степеней свободы df = 8−1-1 = 6 и уровня зависимости 6%, т. е. tтабл = 1,943.
При уровне значимости 6% имеет место неравенство:
ta1 = 0,73 525 ‹ tтабл = 1,943. Значит, с уверенностью 94% можно утверждать, что оценка А1 = 0,747 263 097 не является статистически значимой.
Аналогично проверим для другого параметра. ta0 = 1,743 736 ‹ tтабл = 1,943, значит оценка А0 = 0,123 251 901 также не является статистически значимой.
2. Значимость уравнения регрессии в целом и коэффициента тесноты связи R2 определяется с помощью критерия Фишера. Значение оценки R2 получено в предыдущей задаче, R2 = 0,968 583 448. Фактическое значение Fфакт определяем по формуле: Fфакт = 184,9821.
Табличное значение Fтабл определяем по таблице: Fтабл = 5,99.
Поскольку Fфакт = 184,9821› Fтабл = 5,99, то с уверенностью 94% делается заключение о том, что уравнение регрессии в целом статистически значимо и статистически значим показатель степени связи R2, т. е. отвергается нулевая гипотеза о том, что R2 = 0.
ЗАДАНИЕ № 8.
Имеются следующие исходные данные:
Годы | |||||||||
Объем реализации | 10,84 | 11,12 | 10,6 | 11,31 | 11,62 | 12,0 | 12,73 | 11,12 | |
Коэффициент достоверности аппроксимации для каждого типа линии тренда
1) Линейная у= 0,1795х — 347,71 R2=0.4163
2) Логорифмическая у=359,19 Ln (x)-2718,8 R2=0.1464
3) Степенная y=3E-102x31.059 R2=0.422
4) Экспонтенциальная у=4Е-13е0.01558x R2=0.4218
Как видно из рисунка в 2005 г в сравнении с 2004 г в среднем реализация продукции увеличилась на 0,42 млн. грн.
ЗАДАНИЕ № 9.
Имеются данные испытаний нескольких величин по результатам обследования десяти статистически однородных филиалов фирмы, приведенные в таблице. х1— фондовооруженность, х2 — энерговооруженность, у — производительность труда.
Выполнить следующее:
1. Построить линейную регрессионную модель при помощи ПЭВМ.
2. Выполнить команду «Регрессия».
3. Определить по результатам команды «Регрессия» значение коэффициента множественной корреляции и детерминации.
4. Проверить статистическую значимость оценок параметров модели.
5. Проверить статистическую значимость оценки степени достоверности взаимосвязи R2 и всей модели в целом.
РЕШЕНИЕ.
1. построить регрессионную модель.
2. выполнить команду «Регрессия», результаты которой показаны ниже.
Рис. Результаты команда «Регрессия»
Регрессионная модель принимает вид:
у^ = 929 087*2,9+ - 0,4502*4,5−3,246 374
3. Согласно Рис коэффициенты множественной корреляции и детерминации, в данном случае R = 0,993 689; R2 = 0,98 742.
4. Статистическую значимость оценок параметров модели b, a1, а2 осуществим с помощью t-критерия. Для этого определим его табличное значение и его фактические значения для каждого из оцениваемых параметров. По таблице 1 приложения, А при уровне значимости 1% найдем табличное значение t-критерия для степеней свободы df = 10−2-1 = 7 и уровня зависимости 7%, т. е. tтабл = 3,143.
Фактическое значение t-критерия для каждого из оцениваемых параметров смотрим на рисунке в столбце t-статистика в нашем случае:
t-a1= 15,73 834 ta2= - 0,855 361 tb=15,97 697
При уровне значимости 7% t-a1= 15,73 834> tтабл имеет место равенство: Значит, с уверенностью 99% можно утверждать, что оценка А1 параметра модели является статистически значимой.
Условие ta2 = -0,855 361< tтабл = 3,143 не выполняется, значит утверждаем, что этот критерий статистически не важен.
Условие tв = 15.97 697> tтабл = 3,143 выполняется, значит и эта оценка статистически значима в модели.
5. Значимость уравнения регрессии в целом и коэффициента тесноты связи R2 определяем с помощью критерия Фишера. Фактическое Fфакт =274,684 752
Табличное значение Fтабл определяем по таблице: Fтабл = 9,55. Условие Fфакт =274 684 752> Fтабл = 9.55 выполняется, поэтому с вероятностью 99% делается заключение о том, что R2 статистически значим, и уравнение регрессии в целом значимо, т. е. отвергается нулевая гипотеза R2 = 0.