Файл: Золотарь, И. А. Экономико-математические методы в дорожном строительстве.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 19.10.2024
Просмотров: 141
Скачиваний: 0
|
|
|
|
|
Т а б л и ц а 11 |
|
Количество опытов п при |
|
Количество опытов п при |
||
к % |
£=0,95 |
е= 0 ,90 |
V, % |
£ = 0,95 |
£=0,90 |
|
|
||||
1 |
0,25 |
0,05 |
6 |
9,50 |
2,40 |
2 |
1,00 |
0,26 |
7 |
13,00 |
3,20 |
3 |
2,30 |
0,57 |
8 |
17,00 |
4,15 |
4 |
4,00 |
1,00 |
9 |
22,00 |
5,20 |
5 |
6,50 |
1,70 |
10 |
28,00 |
6,50 |
В математической статистике доказывается [32], что величина %о подчиняется следующему закону распределения:
г,п — 2 |
о |
Хп |
|
|
|
Р п (Хо)" Л-З |
|
(IV. 17) |
где Г (U) — гамма-функция.
Затем доказывается, что вероятность е относительной ошибки среднеквадратичного отклонения, не превышающей заданной ве
личины р, определяется следующим выражением: |
|
|
1 — ~ |
< Р ] = £(Р, л - 1 ) . |
(IV. 18) |
а |
|
|
Функция L (р, п — 1) табулирована (см. приложение 6).
Продолжим решение того же примера и установим с надежностью 0,9 дове
рительный интервал для полученного значения сг=0,81 ч. |
Следовательно, |
|
Из приложения 6 найдем, |
что при е=0,9 и n = 100 р=Ю,125. |
|
среднеквадратичное отклонение |
будет находиться в пределах от |
(1—0,125) сг= |
= (1—0,125)0,81 =0,74 ч до (1+0,125) а = (1 + 0,125)0,81 =0,91 ч. |
|
Следует иметь в виду возможность использования приближен ных формул для оценки доверительного интервала к статистиче
скому среднему I. |
30 и доверительной |
вероятности е = |
Так, при числе опытов |
||
= 0,95 можно использовать зависимость: |
|
|
t = t |
+ 2 , 0 ; . |
(IV.19) |
При е = 0,99 |
i n |
|
|
|
|
* = £ * ± 2 , 7 - ^ . |
(IV.20) |
|
|
i n |
|
С помощью распределения Стюдента может решаться и обрат ная задача — установление необходимого числа испытаний для по
лучения величин t (статистического среднего) и ю (статистическо го среднеквадратичного отклонения) с заданной надежностью, т. е.
62
отклоняющихся от достоверного в заданных наперед пределах. Так, для выше решенного приме ра ширина доверительного интер вала составляла 2а = 2-0,135 = = 0,27 ч. При этом доверительный уровень интервала был 0,90, а ко личество наблюдений за длитель ностью работы составляло 100.
Т а б л и ц а 12
№ |
* в р /’ |
№ |
ЛврГ |
|
образца |
образца |
кгс/см2 |
||
кгс/см 2 |
||||
1 |
280 |
4 |
270 |
|
2 |
310 |
5 |
285 |
|
3 |
315 |
6 |
330 |
Поставим теперь вопрос так: каково должно быть число наблюдений, что бы при том же уровне надежности 0,9 доверительный интервал был равен 0,1ч? На основе соотношения (IV.15) запишем
2£. Г-И— =0,1.
У п — 1
Можно полагать, что число наблюдений должно быть в этом случае значи
тельно больше 100. Тогда из таблицы приложения |
5 найдем |
£а = 1,64. |
|
Используя ранее вычисленное в примере |
значение |
сг=0,81 ч, получим |
|
0,81 |
|
|
|
2-1,64 ,----------=0,1, откуда |
700. |
|
|
У л — 1 |
|
|
|
Таким образом, при л=700 значение длительности работы ожидается в пре делах от 4,9—0,05=4,85 ч до 4,9+0,05=4,95 ч. При этом надежность этого ин тервала колебаний значений статистического среднего t составит 0,9.
Имеются и другие методы определения необходимого числа ис пытаний п для получения результатов с требуемой надежностью. Один из них основывается на предварительном (из рекогносциро вочной серии опытов) определении «меры изменчивости» V как от-
Т а б л и ц а 13
|
|
|
|
|
Значения е |
|
|
|
|
0,85 |
| |
0,90 |
| |
0,95 |
0,99 |
0,995 | |
0,999 |
|
|
|
|
|
Значения |
п |
|
|
0,05 |
207 |
|
270 |
|
384 |
663 |
787 |
1 082 |
0,04 |
323 |
|
422 |
|
600 |
1036 |
1231 |
1691 |
0,03 |
575 |
|
751 |
|
1067 |
1843 |
2 188 |
3 007 |
0,02 |
1295 |
|
1 691 |
|
2 400 |
4146 |
4 924 |
6767 |
0,01 |
5 180 |
|
6 764 |
|
9 603 |
16 587 |
19 699 |
27 069 |
П р и м е ч а н и е . Величина |
допустимой |
|
ошибки а дается долях от значения изучае- |
мой величины.
ношения, установленного в эксперименте среднеквадратичного от клонения к статистическому среднему, т. е.
К = 4 - 100%. |
(IV.21) |
t |
|
Затем по величине V и требуемой доверительной вероятности ре зультата е по табл. 11 устанавливается необходимое число опытов.
63
|
|
|
|
|
Таблица 14 |
№ |
Результат опыта |
{. |
(<г- о 2 |
||
опыта |
|||||
1 |
|
|
|
|
|
2 |
|
|
|
|
|
3 |
|
|
|
|
|
п |
|
|
|
|
|
|
|
п |
|
|
|
|
|
2 * / |
|
|
2 ( ^ г - 0 2 |
|
i = l |
|
Т = a f~ 2 |
;= 1 |
|
|
1 |
п |
- |
о 2 |
|
|
— 2 u = t |
||||
|
п |
/= 1 |
|
У |
|
|
|
|
|
||
|
|
а |
|
|
|
|
|
V 'n ^ i |
|
а |
_____ |
|
|
|
|
у п —1 |
Рассмотрим пример реализации этой методики. При испытании в 28-суточном возрасте серии бетонных кубиков в дорожной лаборатории получены результаты, приведенные в табл. 12. Определить, какое количество кубиков должно быть из готовлено и испытано для получения показателя прочности бетона на сжатие с
надежностью 0,95. |
_ |
|
|
1. |
Вычислим среднее значение R Bр: |
|
|
|
RBр = |
280 + 310 + 315 + 270 + 285 + 330 |
кгс/см2. |
|
----------------------- -т--------!--------------- = 299 |
||
2. |
Определим |
с помощью формулы (III.12) экспериментальное значение а: |
|
|
(280 — 299)2 + (310 — 299)2 + (315 — 299)2 + |
(270 — 299)2 + 43* |
|
|
|
6 |
|
|
|
(285—299)2 + (3 3 0 — 299)2 |
21,4 кгс/см2. |
|
|
------------------------------ ------------ = |
|
|
|
о |
|
3. Вычислим меру изменчивости
21,4
У= 2991О00/о=7-2°/о-
4.Из табл. 11 находим, что для получения Л вр с надежностью 0,95 нужно из готовить и испытать 14 кубиков.
Если, планируя эксперимент, мы не знаем меры изменчивости случайной величины, то количество необходимых опытов для суж дения об изучаемом явлении с необходимой надежностью резко возрастает. Это легко видеть из так называемой «таблицы доста точно больших чисел», где даются значения п при различных а и е (табл. 13). Поэтому всегда выгоднее провести серию предвари
тельных опытов, с тем чтобы определить а или V и воспользовать ся для установления необходимого числа опытов распределением Стюдента или табл. 11. В обоих этих случаях приходится вести обработку экспериментальных данных. При этом удобна форма за писи, приводимая в табл. 14.
64
Г па р а
V
л и н е й н о е п р о г р а м м и р о в а н и е
ПРИМЕНИТЕЛЬНО К ЭКОНОМИЧЕСКОМУ АНАЛИЗУ ДОРОЖНОГО СТРОИТЕЛЬСТВА
§ 10. Задачи, решаемые методом линейного программирования
Задачи линейного программирования являются разновидностью более широкого класса задач, называемых распределительными. Задачи распределения возникают в случаях, когда имеющихся ре сурсов недостаточно для выполнения каждой из работ с наиболь шим эффектом. При подобных обстоятельствах приходится забо титься о наилучшем использовании имеющихся ресурсов в целом. Рассмотрим в качестве примера приводимую ниже матрицу.
сц=0,4 |
с 1 2 = 0 ,1 |
« 1 3 = 0 >2 |
« 1 4 = 0 ,6 |
«15=0,9 |
<*!= 10 |
|
■*п |
*12 |
*13 |
*14 |
*15 |
||
|
||||||
С21= 0,6 |
С 2 2 = 0 ,4 |
с23=0,3 |
« 2 4 = 0 ,5 |
« -2 5 = 0 ,2 |
а 2= 12 (г) |
|
*21 |
*22 |
*23 |
*24 |
*25 |
||
|
||||||
с 31 = 0 ,5 |
с 32 = 0 ,2 |
« з з = 0 >6 |
с34 = 0 ,4 |
« 3 5 = 0 ,8 |
а3= 12 |
|
*31 |
*32 |
*33 |
*34 |
*35 |
|
|
й,=4 |
62=5 |
* 3 = 7 |
* 4 = 9 |
*5=9 |
|
|
|
|
U) |
|
|
|
Величины Cij выражают стоимости перевозки 1 м3 материала из
карьера i на участок работ j (руб.); а, дают запасы материала в карьерах, a bj — потребности в материалах для производства до рожных работ (тыс. м3) .
Рассматривая матрицу по столбцам, легко заметить, что на участки работ 1 , 2 и 3 выгоднее всего было бы вывозить материал
из карьера 1. Однако потребность в материале для упомянутых
участков составляет 4 + 5 + 7 = 16 |
тыс. м3, тогда |
как запасы мате |
риала в карьере 1 равны 10 тыс. |
м3. Очевидно, |
что часть материа |
ла на эти участки придется вывозить из других карьеров при более высокой стоимости перевозок. Аналогично можно сделать вывод о целесообразности вывозки необходимых для четвертого участка материалов из карьера 3, а для пятого участка — из карьера 2 .
Таким образом, минимальная общая стоимость перевозки со ставила бы при Х ц = 4; Xi2 = 5; X i 3 = 7; х 34 = 9 и *25 = 9:
3—1092 |
65 |
2 C/;.X (7 = ( c 11X n |
4 - C 12X i 2 + C13X 13 + C 34A :3 4+C 25^25) 1Q3 = |
= (0,4 -4 + 0 ,1 -5 + |
0,2-7 + 0,4 -9+ 0,7 -9) 103= 13400 руб. |
Подобный план перевозок был бы возможен при наличии сле дующих количеств материала в карьерах: ai=16; а^= 9 и а3 = 9. Однако наличие ресурсов в карьерах не соответствует этому ус ловию, и придется заботиться об отыскании решения, учитываю щего реальное наличие материалов в карьерах с возможной мини мальной стоимостью перевозок. Забегая вперед, отметим, что воз никающая при этом распределительная задача является транспорт ной задачей линейного программирования. В данном примере ее решение дается следующей матрицей:
1 |
5 |
4 |
0 |
0 |
ai=10 |
0 |
0 |
3 |
0 |
9 |
«2=12 |
3 |
0 |
0 |
9 |
0 |
«3=12 |
й ,=4 |
Ь2=5 |
Ь3=7 |
Ьа— 9 |
65= 9 |
|
Естественно, что общая стоимость перевозок будет больше вы численных выше 13 400 руб. и составит:
У,сих и = (0,4 -1 + 0 ,1 -5 + 0,2-4 + 0,3-3 + 0,7-9 +
+ 0,5 -3+ 0,4 -9) 103 = 14000 руб.
Для лучшего понимания математического смысла задач линей ного программирования рассмотрим следующую систему урав нений:
2xi + 3x2+ x3= 8 ; 1
х 1 + 2 х 2 + 2 х3 = 5. J |
( |
Система (V.1) является неопределенной, так как для отыскания трех неизвестных имеется всего два уравнения. Такая система в принципе имеет бесчисленное множество решений. Метод решения подобных систем состоит в приведении их к определенным за счет введения дополнительных условий. Так, если принять, что одно из неизвестных в системе (V.1) может равняться нулю, то получаем следующие три решения системы:
1 -е |
решение: х ^ О ; |
х 2= - ^ - |
|
; |
х 3= |
— |
i- ; |
||
'2 -е |
решение: |
х 2 = |
0 ; |
11 |
|
|
х 3= |
о |
; |
х х — -----; |
|
— |
|||||||
|
|
|
|
3 |
|
|
|
3 |
|
З-е решение: |
х3 = |
0 ; |
x ^ l ; |
х |
2 = |
2 . |
|
|
66