Бондаренко М.В.
Донецький національний технічний
університет
ОБГРУНТУВАННЯ ЗАКОНУ РОЗПОДІЛУ
ВИРОБНИЧОЇ СОБІВАРТОСТІ ДЛЯ ВУГЛЕДОБУВНИХ ПІДПРИЄМСТВ
Собівартість
вуглевидобутку є важливішим інтегральним критерієм, що характеризує не тільки
ефективність роботи вугільних підприємств, але і всього паливно-енергетичного
комплексу, включаючи й теплові
електростанції (ТЕС). Тому на даний момент головна ціль – встановлення основних
закономірностей формування цього економічного показника.
Фактична
виробнича собівартість складається з двох компонент: планованої
(детермінованої) та імовірнісної, величина якої обумовлена чинниками, вплив
яких не враховано в процесі планування через розмитість інформації або
випадкового їх прояву. Тому особливого значення набуває застосування
математичного апарату при розробці перспективних, оперативних і поточних
планів, оскільки висока інформаційна невизначеність у період їх формування
припускає оцінку імовірності прогнозованого рівню собівартості.
Згідно з
дослідженнями [2,3] в якості типового може бути
прийнято бета-розподіл, що описується системою рівнянь:


де α,β – статистичні параметри
розподілу,
хmin, xmax
– мінімальне
та максимальне значення випадкового параметру х.
Теоретичні
частоти розподілу можуть бути визначення за формулою:
f’=nР
де n – обсяг вибірки,
Р – оцінка імовірності влучення в заданий інтервал.
Імовірність
влучення в інтервал [xi,xi-1] можна визначити з умови:
Р=Р(хі)-Р(хі-1)
де 
Основою для
виконання аналізу гістограми фактичного розподілу виробничої собівартості як
випадкової величини служать помісячні дані про собівартість 1 т вугілля на
шахті «ім.О.О.Скочинського». (табл.1)
Таблиця 1 вихідні дані
|
Період |
Виробнича собівартість 1 т, грн |
Період |
Виробнича собівартість 1 т, грн |
||
|
Планова |
Фактична |
Плінова |
Фактична |
||
|
Січень 2007 |
287,06 |
350,73 |
Жовтень 2008 |
841,42 |
650,22 |
|
Лютий 2007 |
277,95 |
268,16 |
Листопад 2008 |
493,52 |
621,46 |
|
Березень 2007 |
260,33 |
260,39 |
Грудень 2008 |
470,67 |
757,42 |
|
Квітень 2007 |
278,38 |
298,92 |
Січень 2009 |
475,35 |
686,76 |
|
Травень 2007 |
259,57 |
277,72 |
Лютий 2009 |
702,94 |
739,17 |
|
Червень 2007 |
359,93 |
376,96 |
Березень 2009 |
437,47 |
627,63 |
|
Липень 2007 |
311,26 |
456,09 |
Квітень 2009 |
679,23 |
575,57 |
|
Серпень 2007 |
454,6 |
314,9 |
Травень 2009 |
682,37 |
635,57 |
|
Вересень 2007 |
280,66 |
340,55 |
Червень 2009 |
1481,26 |
1381,09 |
|
Жовтень 2007 |
673,78 |
773,5 |
Липень 2009 |
695,5 |
740,25 |
|
Листопад 2007 |
673,37 |
612,84 |
Серпень 2009 |
670,93 |
985,77 |
|
Грудень 2007 |
688,35 |
524,83 |
Вересень 2009 |
616,43 |
911,61 |
|
Січень 2008 |
856,81 |
552,59 |
Жовтень2009 |
732,31 |
806,81 |
|
Лютий 2008 |
737,22 |
577,21 |
Листопад 2009 |
674,73 |
750,78 |
|
Березень 2008 |
739,8 |
547,22 |
Грудень 2009 |
686,51 |
811,95 |
|
Квітень 2008 |
782,81 |
489,6 |
Січень 2010 |
662,81 |
729,53 |
|
Травень 2008 |
735,16 |
662,58 |
Лютий 2010 |
784,36 |
745,66 |
|
Червень 2008 |
792,86 |
572,69 |
Березень 2010 |
992,82 |
836,62 |
|
Липень 2008 |
743,11 |
858,79 |
Квітень 2010 |
1085,02 |
835,9 |
|
Серпень 2008 |
766,72 |
682,81 |
Травень 2010 |
1028,75 |
742,42 |
|
Вересень 2008 |
764,84 |
665,18 |
Червень 2010 |
1013,49 |
967,93 |
Для вивчення
форми розподілу вихідні дані представимо у вигляді інтервального ряду,
мінімальне значення якого дорівнює 200,00 грн , а максимальна – 1400 грн.
Довжина
інтервалу визначається за формулою Стерджеса:
![]()
де N -
число спостережень.
![]()
Емпіричні частоти інтервального ряду розподілу
представлені в табл.2.
Табл.2 Емпіричний ряд розподілу фактичного обсягу товарної продукції
|
Виробнича собівартість 1 т, грн |
Частота |
|
|
200 |
400 |
7 |
|
400 |
600 |
7 |
|
600 |
800 |
15 |
|
800 |
1000 |
6 |
|
1000 |
1200 |
0 |
|
1200 |
1400 |
1 |
Гістограма розподілу емпіричних частот і вид кривої, що описує
передбачувану функцію щільності розподілу розглянутого параметра, представлена
на рис. 1.
![]()

Рис.1 Гістограма розподілу емпіричних частот фактичної собівартості 1 т
вугілля
Візуальний аналіз гістограми дозволив припустити вид теоретичної функції
розподілу та розрахувати її параметри:
а) вид кривої, що описує функцію щільності розподілу розглянутого
показника, дозволяє припустити, що
,
. Допустимо, що параметри бета-розподілу дорівнюють
,
, тоді функція щільності має вигляд:
(*)
б) математичне очікування дорівнює:
![]()
Графік даної функції на досліджуваному інтервалі [200;1400] має вигляд
(див.рис.2), відповідний гістограмі розподілу емпіричних частот, представлений
на рис.1.

Рис. 2 Щільність бета-розподілу на інтервалі [200;1400]
при
, ![]()
Скориставшись можливостями програми Microsoft Excel, можна визначити значення
інтегральної функції щільності бета-імовірності Р(х). Для оцінки істотності відхилень емпіричних і теоретичних
частот розраховується критерій згоді Пірсонаа і визначається його критичне
значення для відповідної довірчої імовірності. Результати відповідних
розрахунків представлено в табл.3.
Критичне значення критерію Пірсона для довірчої імовірності 0,99 виявилося
нижче розрахункового, що підтверджує можливість використання досліджуваного
бета-розподілу для прогнозування зміни техніко-економічного показника (рис.3).
![]()
![]()
![]()
![]()
![]()
![]()
![]()
![]()
![]()
![]()

Рис.3 Гістограма частот та теоретична функція щільності імовірностей
Гістограма, що приведена до виду щільності – це гістограма, в якій вектор (Z) числа потраплянь
квазістатистики у відповідний інтервал гістограми розмірністю N замінено вектором fі:
![]()
де
- шаг гістограми.
Наприклад, ![]()
Таким чином, розподіл фактичної собівартості 1 т вугілля щонайкраще
описується за допомогою функції (*).
Таблиця 3. Перевірка відповідності емпіричного розподілу бета-розподілові
за умови
, ![]()
|
Фактична виробна собівартість 1 т вугілля, грн |
Р(хі-1) |
Р(хі) |
Імовірність |
Теоретична частота, f’ |
Емпі-рична частота, f |
(f- f’)2/ f’ |
|
|
200,00 |
400,00 |
0 |
0,131944 |
0,131944 |
5,541666 |
7 |
0,383772 |
|
400,00 |
600,00 |
0,131944 |
0,407407 |
0,275462 |
11,569444 |
7 |
1,804739 |
|
600,00 |
800,00 |
0,407407 |
0,6875 |
0,280092 |
11,763888 |
15 |
0,890217 |
|
800,00 |
1000,00 |
0,6875 |
0,888888 |
0,201388 |
8,458333 |
6 |
0,714491 |
|
1000,00 |
1200,00 |
0,888888 |
0,983796 |
0,094907 |
3,986111 |
0 |
3,986111 |
|
1200,00 |
1400,00 |
0,983796 |
1 |
0,016203 |
0,680555 |
1 |
0,149943 |
|
|
|
Обсяг вибірки |
42 |
Розрах.критерій(Σ) |
7,929273 |
||
|
|
|
|
|
Критич. знач. для імовірності 0,99 |
16,81 |
||
2. Голенко Д.И. Статестические методы в экономических
системах.-М.:»Статистика», 1970.-320
3.Евдокимов Ф.И. Механизм управления затратами на основе оптимизации
параметров процес сов угледобычи/ Евдокимов Ф.И., Кучер В.А //Наукові праці
Донецького національного технічного університету. Серія:економічна. Випуск 22.
– Донецк, ДонДТУ,2000.С.152.-160