Файл: Золотарь, И. А. Экономико-математические методы в дорожном строительстве.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