CC BY
117
Моделирование производства культивируемых грибов с использованием аппарата производственных функций
А.А. Набоких
Простейшую модель производства можно представить как некоторую систему, перерабатывающую различные виды ресурсов в готовую продукцию (рис. 1).
Рис. 1. - Простейшая модель производства
В качестве ресурсов для производства культивируемых грибов служат:
1) сырье (мицелий, солома, торф, пр.);
2) трудовые ресурсы;
3) энергозатраты;
4) научно-исследовательские ресурсы;
5) технологические ресурсы;
6) транспортные ресурсы и др.
Производственная функция (ПФ) описывает соотношение, существующее между количеством используемых в производстве различных факторов (ресурсов) и максимальным количеством продукта, который может быть произведен при помощи данных факторов. ПФ является математической моделью процесса производства продукции в данной экономической системе и выражает устойчивую закономерную количественную зависимость между объемными показателями ресурсов и оборота. Производственные функции лежат в основе моделирования деятельности самых разнообразных производственных структур и систем от отдельных предприятий и организаций до регионов, отраслей и экономики страны в целом.
Поэтому к выбору типа функции, привлекаемой информации и методам ее обработки должны быть предъявлены высокие требования в отношении объективности отражения реальных закономерностей общественного производства. Связь между затратами ресурсов и выпуском продукции является характерной и относительно устойчивой для экономических систем любого уровня, технология производства в которых не претерпевает существенных скачков. Отражение и количественное определение этой связи и составляет экономическое содержание понятия производственной функции.
В общем виде производственная функция имеет вид: у = /(хь*2,...,хп), где у -выпуск продукции, *1,*2,...,хп - факторы производства. Выпуск продукции у является
кумулятивным результатом действия всех факторов производства.
К настоящему времени накоплен большой опыт построения производственных функций для различных объектов [1,2]. В зависимости от масштабов объекта исследования, используемых технологий производства, взаимозаменяемости производственных факторов, могут быть использованы разнообразные функции.
В начале работы рассмотрим двухфакторную производственную функцию с предельно агрегированными характеристиками. Для исследуемого производства в качестве таковой возьмем ПФ с одним выходом - выпуском продукции (У) и с двумя входными воздействиями - затратами капитальных (К) и трудовых (Ь) ресурсов. Структурно такая агрегированная модель грибоводческого предприятия представлена на
рис. 2.
Рис. 2. - Структурное представление модели грибоводческого предприятия
Математическое моделирование функционирования грибоводства в Кировской области выполнено для периода 2002-2011 гг. на основе реальных исходных данных. В качестве показателя производства была принята величина годового объема произведенных культивируемых грибов в Кировской области (в тоннах), в качестве показателя капитальных ресурсов выбрана величина основных фондов по полной учетной стоимости (в фактических ценах на конец года) (в тыс. руб.), под показателем «трудовые ресурсы» понимается среднегодовая численность промышленно-производственного персонала. Математическая обработка эмпирических данных производилась с помощью инструмента анализа «Регрессия» программы Microsoft Excel.
Формализовано двухфакторная ПФ в общем виде определяется зависимостью
Y = F(K, L) .
Первоначально в качестве модели была выбрана двухфакторная однородная производственная функция Кобба-Дугласа с предельно агрегированными
характеристиками:
Y(t) = A • Ka(t) • L1-a (t) (1)
В (1) А - коэффициент, характеризующий масштабную эффективность системы. Большей величине параметра А при фиксированных значениях факторов производства соответствует большее значение Y(t). Следовательно, и производственный процесс, описываемый такой функцией, более эффективен. Кроме того, А - коэффициент, учитывающий влияние факторов, не вошедших в уравнение ПФ. Величины а и в=1-а -функции чувствительности (эластичности), характеризующие относительный вклад капитала и труда в грибное производство. Коэффициенты а и в показывают, на сколько процентов изменится в среднем результативный признак с изменением значения соответствующего фактора на 1% при неизменности воздействия других факторов.
Функция Кобба-Дугласа чаще всего используется для описания среднемасштабных хозяйственных субъектов, характеризующихся устойчивым, стабильным
функционированием, когда вовлечение дополнительной единицы ресурса приносит эффект, пропорциональный средней производительности имеющегося ресурса.
Общепринятым подходом при идентификации параметров ПФ является использование исходных статистических данных в натуральных единицах измерения. Однако практика расчетов показывает, что при этом разброс значений масштабных коэффициентов А может составлять несколько порядков. Поэтому в работе наряду с натуральными единицами использованы исходные данные в безразмерной форме. Сведение к последним осуществлялось путем приведения временного ряда к первому году интервала идентификации. Кроме того, идентификация моделей и расчеты проведены для несглаженных и сглаженных исходных данных. Сглаживание производилось методом скользящего среднего по трем точкам с целью устранения влияния случайных выбросов.
При оценке множественной регрессии для обеспечения статистической надежности требуется, чтобы число наблюдений, по крайней мере, в три раза превосходило число оцениваемых параметров [3]. Идентификация ПФ осуществлялась методом минимизации среднего квадратического отклонения модельных зависимостей от эмпирических данных.
У к
ПФ (1) сводится к линейной 1п— = 1п А + а1п—, параметры которой оценивались
Ь X
методом наименьших квадратов. Качество моделей и их параметров определялось значениями множественного коэффициента корреляции Я, коэффициента детерминации Я , ^-критерия Стьюдента, К-критерия Фишера.
Результаты оценки параметров и качества моделей однородной производственной функции Кобба-Дугласа при уровне значимости 5% представлены в таблице 1.
Таблица №1
Параметры модели и коэффициенты Данные в натуральных единицах Несглаженные данные (безразмерная форма) Сглаженные данные (безразмерная форма)
А 4,0986Е-10 2,07283 1,38748969
а 4,27682682 4,276827 5,667408
Я 0,77950062 0,779501 0,91934019
Я2 0,60762121 0,607621 0,845186
Р-значение А 0,009064 0,269498 0,5215
Р-значение а 0,007849 0,007849 0,001234
К 12,38846 12,38846 32,75628
Значимость К 0,007849 0,007849 0,001234
Высокие значения множественного коэффициента корреляции говорят о достаточно тесной связи между функцией и всеми включенными в модель факторами-аргументами. Коэффициент детерминации (особенно в случае сглаженных данных) свидетельствует о достаточно полном охвате полученной моделью всех существенно влияющих на функцию факторов. На другие факторы приходится только 8,1%.
Если вероятность (значимость К) меньше заданного уровня значимости, то гипотеза о незначимости регрессии (т.е. гипотеза о том, что все коэффициенты функции регрессии равны нулю) отвергается и считается, что регрессия значима. В нашем случае регрессия значима во всех моделях.
В столбце Р-Значение вычисляются уровни значимости, соответствующие значениям критериальных характеристик. Если вычисленный уровень значимости меньше заданного уровня, то принимается гипотеза о значимом отличии коэффициента от нуля; в противном случае принимается гипотеза о незначимом отличии коэффициента от нуля. В нашем случае коэффициент А незначимо отличается от нуля во второй и третьей моделях.
Полученные коэффициенты свидетельствуют о том, что однородная
производственная функция У(^) = А • Ка (I) • Ь1-а(^) неудовлетворительно описывает анализируемые процессы грибоводства, модели не обладают требуемыми аппроксимативными свойствами, и, следовательно, не пригодны для описания
функционирования грибоводческого предприятия в рассматриваемый период.
Далее были изучены возможности моделирования функционирования системы на основе двухфакторной неоднородной ПФ Кобба-Дугласа:
У(0 = А • Ка(Г) • Ьв(1), где а + РФ 1. (2)
ПФ (2) сводится к линейной 1п У = 1п А + а1п К + Д1п Ь, параметры последней оценивались методом наименьших квадратов.
Приведенная зависимость обладает более гибкими аппроксимативными
свойствами, так как имеет три адаптивных настроечных параметра А, а, Д
Идентификация параметров модели (2) была проведена аналогично однородной ПФ Кобба-Дугласа.
Полученные численные характеристики моделей для неоднородной
производственной функции Кобба-Дугласа приведены в таблице 2.
Таблица №2
Параметры модели и коэффициенты Данные в натуральных единицах Несглаженные данные (безразмерная форма) Сглаженные данные (безразмерная форма)
А 6,8586Е-11 1,22398053 0,79413566
а 2,25810493 2,25810493 2,69214467
в 2,86804963 2,86804963 3,60832425
X 5,1261545 5,12615456 6,30046892
Я 0,88181452 0,88181452 0,97168122
я2 0,77759685 0,777597 0,944164
Р-значение А 0,004499 0,0453873 0,0499991
Р-значение а 0,0191807 0,0191807 0,0107424
Р-значение в 0,0453251 0,0453251 0,0329155
Р 12,23719 12,23719 42,2743
Значимость Р 0,005188 0,005188 0,000737
Статистические показатели качества аппроксимации, приведенные в таблице 2, показывают, что описательные свойства модели (2) улучшились по сравнению с моделью (1): коэффициент детерминации увеличился, повысилась адекватность моделей.
Для грибоводческого сектора Кировской области как в таблице 1, так и в таблице 2 величины а и в достаточно высокие, что указывает на эффективное использование как капитальных, так и трудовых ресурсов. Изменение выпуска грибной продукции более чувствительно к изменениям объема использования фактора производства труд.
В таблице 2 помимо факторных эластичностей а и в приведены значения полной эластичности Х=а+в, характеризующие интегральную, масштабную эффективность грибного производства. Поскольку Х>1, имеем возрастающую отдачу от масштаба.
Графики реальных и модельных траекторий изображены на рис. 3, 4.
-У эмпирическое
У натуральные единицы измерения
■У неслаженные данные, безразмерная величина
-У сглаженные данные, безразмерная величина
Рис. 3. - Графики реальных и модельных данных для однородной ПФ Кобба-Дугласа
123456789 10
У эмпирическое
У натуральные единицы измерения
У несглаженные данные, безразмерная величина
У сглаженные данные,безразмерная величина
Рис. 4. - Графики реальных и модельных данных для неоднородной ПФ Кобба-Дугласа По данным модели (2) (сглаженные данные, безразмерная форма) рассчитаны
средняя капиталоотдача У^) = 0,794 • К1,692 ^) • Ь3,608^) и средняя производительность
К (V )
У (t)
труда = ^’^^4 • К2,в92^) • Ь2,608^), которые показывают, сколько тонн
культивируемых грибов приходится на 1 тыс. руб. основных фондов и 1 чел./час
ду (t)
затрачиваемого труда. Предельная капиталоотдача ^) = 2,137 • К 1,692(^ • Ь3,608^) и
предельная производительность труда дУ^) = 2,864 • К2,692 (^ • Ь2,608^) показывают, на
дЬ^ )
сколько дополнительных тонн продукции увеличится производство культивируемых грибов, если объем соответствующего фактора производства вырастет на одну единицу при неизменности других.
По данным наилучшей модели У^) = А • Ка ^) • Ьв ^), а + в ^ 1 (сглаженные данные, безразмерная форма) были построены изокванты, отражающие зависимость между факторами, входящими в производственную функцию при постоянном значении самой функции. Общее уравнение изокванты )) имеет вид:
У °,277^ ) т (t) =__ _____^_____
() 0,938 • К 0,746(0.
Изоквантами являются гиперболы, асимптоты которых представляют оси координат.
К изокванте для среднего значения количества произведенной продукции за исследуемый период была построена изокоста, представляющая собой касательную к изокванте, в точке, соответствующей минимуму суммы затрат на капитал и трудовые ресурсы:
Ь(К($)) = 4,148 - 0,9996К(t) .
Угловой коэффициент изокосты равен предельной норме замещения и характеризует возможность замены затрат на трудовые ресурсы капитальными при сохранении производства в регионе. Расчеты показывают, что абсолютное значение коэффициента предельной нормы замещения приблизительно равно единице. То есть при фиксированном значении объема производства культивируемых грибов рост капитальных вложений способен обеспечить пропорциональное снижение количества используемых
трудовых ресурсов, и наоборот, рост используемых трудовых ресурсов пропорционален снижению капитальных затрат.
Научно-технический прогресс (НТП) при построении производственной функции
может быть учтен с помощью введения множителя научно-технического прогресса вп, где параметр у характеризует темп прироста под влиянием НТП:
У (г) = А • ка(г) • Ьв(г) • в71 (3)
Такая динамическая производственная функция включает нейтральный, т.е. не материализованный в одном из факторов, технический прогресс.
Логарифмируя зависимость (3), получим линейную зависимость 1п У = 1п А + а1п К + р\п Ь + ^ , параметры которой определяются с помощью метода наименьших квадратов. Полученные численные характеристики модели приведены в таблице 3.
___________________________________________________________________________Таблица №3
Параметры модели и коэффициенты Данные в натуральных единицах
А 0,2870943
а 1,33760997
в 4,38723714
У 0,48817298
Я 0,92840367
я2 0,86193337
Р-значение А 0,0492399
Р-значение а 0,0482289
Р-значение в 0,0216244
Р-значение у 0,0104067
Р 12,48576
Значимость Р 0,00545
Как видно из таблицы 3, уравнение (3) значимо на высоком уровне по F-критерию, как и его параметры по t-критерию. Высоким является коэффициент детерминации. Следовательно, изменение выпуска готовой продукции грибоводства в зависимости от затрат ресурсов и влияния НТП на 86% объясняется построенной моделью (3).
Как и в моделях (1) и (2), эффективность использования трудовых и капитальных ресурсов велика. Положительное значение параметра у указывает на современные технологии в этой отрасли.
Графики моделируемой производственной функции (3) и реального выпуска продукции грибоводства представлены на рис. 5.
450 400 350 300 250 200 150 100 50 0
123456789 10
Рис. 5. - Графики реальных и модельных данных для ПФ Кобба-Дугласа с НТП
-Y эмпирическое Y расчетное
Для описания функционирования грибоводческого предприятия может быть использована и производственная функция Леонтьева. Указанная функция предназначена для моделирования строго детерминированных технологий, не допускающих отклонения от технологических норм использования ресурсов на единицу продукции. Обычно используются для описания мелкомасштабных или полностью автоматизированных производственных объектов.
Функция Леонтьева имеет вид:
У = шт
(4)
где коэффициенты К0 и Ь0 выражают количество соответствующего ресурса, необходимого для производства единицы продукта. Расчеты показали, что К0=233,178 тыс.руб., Ь0=1,644 чел./час.
При использовании производственной функции (4) имеет постоянная отдача от масштаба. Эластичность замены факторов по любому ресурсу в такой модели равна нулю.
Литература:
1.Клейнер, Г.Б. Производственные функции: теория, методы, применение [Текст] / Г.Б. Клейнер. - М.: Финансы и статистика, 1986. - 239 с.
2.Лотов, А.В. Введение в экономико-математическое моделирование [Текст] / А.В. Лотов. - М.: Наука, 1986. - 392 с.
3.Замков, О.О., Толстопятенко, А.В., Черемных, Ю.В. Математические методы в экономике: Учебник. 2-е изд. [Текст] / О.О. Замков, А.В. Толстопятенко, Ю.В. Черемных. - М.: Дело и сервис, 1999. - 368 с.
