Файл: Гинзбург Г.Л. Руководство по использованию метода наименьших квадратов в экономико-статистических расчетах учеб. пособие.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 27.07.2024

Просмотров: 102

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

а0= 0 ,2289062 •180—0,003906250 -5708= 18,906; +=0,003030303 •234=0,709;

as= —0,003906250-180+0,0001183712-5708=—0,027;

jTt= 18,9+0,71Н— 0.027*2 .

Произведем перенос осей координат по формуле

х' — t — 2{х к) — 1 = 2(х — 5) — 1 = — 11;

ух— 18,9 + 0,71(2* — 11) — 0,027 (4*3 — 44* +

121) =

 

= 7,8+2,61* — 0,11*2 .

 

Пример 2. В таблице 19 .приведены данные о производстве

изделия А по годам.

 

 

 

 

 

 

 

 

Таблица 19

 

Производство

 

 

 

 

Годы

изделия А

t

y t

y t2

yt

(в тыс. шт.)

 

 

 

 

 

 

У

 

 

 

 

1953

3

—4

—12

48

6,4

1954

10

—3

—30

90

7,9

1955

15

—2

—30

60

12,6

1956

21

—1

—21

21

20,3

1957

35

0

0

0

31,2

1958

42

1

42

42

45,3

1959

58

2

116

232

62,4

1960

81

3

243

729

82,7

1961

НО

4

440

1760

106,2

Итого

Ь„= 375

 

6, =748

*-.=2982

375,0

Здесь

Л—-I

п — 9, к = — — = 4 (нечетный случай).

Выровняем данный динамический ряд по параболе 2-го по­ рядка (т ~ 3):

а0=

0,2554112-375 — 0,0216450 •2882 = 31,23385;

+ =

0,0166667-748 = 12,46667;

+ =

— 0,0216450 •375 + 0,0032467 •2982 = 1,56493;

^t =

31,2+12.5f + 1,6/а

91


Произведем преобразование координат по формуле

x ' = t = x к — 1 = х — 5;

ух-= 31,2+12,5(jc—5)4-1 ,6(jc—5)e^=l .6jc®—3,5jc+8,7 .

§ 7. Нахождение параметров экспоненциальной формы связи

Экспоненциальному закону подчиняются такие важные для экономических исследований процессы, как рост народо­ населения, рост суммы 'накоплений по сделанному в сберега­ тельную кассу вкладу, рост числа научных работников и на­ учных публикаций, рост количества древесины на лесных массивах и т. д., т. е. такие процессы, в которых скорость или интенсивность изменения какой-либо величины пропорцио­ нальна самой этой величине.

Экономисту-исследователю часто приходится в своей ра­ боте анализировать реальные процессы, которые в опреде­ ленный период времени обнаруживают тенденцию к быстрому нарастанию и развиваются по закону, близкому к экспонен­ циальному. Это может быть поток телеграмм в предпразд­ ничные дни, опрос на железнодорожные и авиабилеты в пе­ риод летних отпусков, увеличение числа пассажиров на го­ родском транспорте в часы «пик», спрос на товары, входя­ щие в моду и т. д.

Для изучения подобных процессов экспоненциальный закон может оказаться наиболее подходящим. Идет ли речь о национальном доходе, о выплавке стали или о приплоде у кроликов, экономист должен уметь сделать научно обосно­ ванные выводы о тенденции развития этих явлений.

Получив статистические или экспериментальные данные и нанеся их на график, лучше всего на полулогарифмической бумаге, экономист, при определенном навыке, может устано­

вить, что данный процесс развивается по закону,

близкому

к экспоненциальному. После этого он может

приступить

к математической обработке имеющихся данных по предла­ гаемому ниже способу.

Рассмотрим суперпозицию показательной и полиномиаль­

ной функции

 

 

 

y { t ) = 10х» + М +

••• + хт-1 *т ~1t

(3,

20)

где t аргумент, принимающий дискретные значения

 

 

(t—1,2,

a lQ, X,........Xm_,

(3,

21)

параметры, которые определим ниже.

92


После логарифмирования функции (3, 20) получим

ср(^)= lg (£)] = Х0 -)- \ t + ... -j- Xm_i tm~l

(3,21)

Параметры функции (3, 21) можно определить путем решения системы нормальных уравнений

5 0Х0 -Ь -j-. . . -f- Sm—i Xm_i — bQ

•SjXo

S aXj

+ S ^ m - l

— b x

 

 

 

(3,22)

5 mX0 +

^m+jXj -f- . . . -J- >S2m—2 Xm—1 — bm—i

где

 

 

 

5/ = 1^—)—2г—}—... -f- я 1—

0 ,1,2,. >2m—О

bi — \g у

+ ]g yx2‘+

•.. + lg Уп =

ig y} i‘

 

(i = 0, 1, 2, ... , m — 1) .

 

Для аппроксимации

экспериментальных

данных с по­

мощью функции (3, 20) предлагается алгоритм, с помощью

которого экономист-экспериментатор может получить пара­

метры уравнений динамики,

не только нерешая системы нор­

мальных уравнений (3, 22),

но даже и не выписывая самой

системы. Система (3, 22)

была решена в общем виде при ко­

личестве уравнений ш =2,

3, 4 и 5, после чего с помощью

ЭВМ М-22 были составлены специальные таблицы коэффи­ циентов, которые дают возможность сразу определять пара­ метры уравнений динамики по выведенным заранее форму­ лам (приложения 4 и 5). При решении системы (3, 22) было рассмотрено два случая: нечетный и четный.

Формулы для определения параметров экспоненциальной формы связи имеют такой же вид, как формулы, выведенные

в § 2 настоящей главы, только значения

 

Ь0 = Щ у , b ^ - m g y , b2 = I,t2\gy, bs =

T,t3lg y ,...

нужно находить с использованием таблиц

десятичных или

натуральных логарифмов '.

 

1 В данном случае таблицы десятичных логарифмов удобнее, так как

при работе с таблицами натуральных логарифмов нужно вычислять до­ полнительные поправки.

93


Теперь решим вопрос о том, как в каждом конкретном случае выбрать вид суперпозиции (3, 20), т. е. определить число т (количество уравнений). Ясно, что чем выше поря­ док аппроксимирующей суперпозиции, тем более сложную задачу придется решать для отыскания параметров аппрок­

симирующей формулы л0, . ,л,„_!, но и тем точнее бу­ дет произведено выравнивание по методу наименьших квад­

ратов отклонений логарифмов

'.

Ниже укажем каким значением числа т (для конкретной

задачи) следует ограничиться,

чтобы решение было более ра­

циональным (т = 1 интереса

не представляет, так как

#=10*» =const — прямая, параллельная оси t).

Рассмотрим возможность применения аппроксимирующей формулы (3, 20) и поясним наши соображения на примерах.

§ 8. Применение экспоненциальной формы связи при анализе конкретных экономических процессов

Пример 1. В таблице 21 приведены данные о производ­ стве изделия А по годам.

У

5 t -4 -3 - i -i о V 2 3 к

Рис. 17.

1 Здесь сумма квадратов отклонений для исходных переменных! не

минимизируется, что в значительной мере ухудшает смысл аппроксимации.

94

Здесь

т= 2, п = 9. & = 4

Решение. Нанесем экспериментальные точки (у) на гра­ фик, используя полулогарифмическую бумагу (см. рис. 17 и

У

18). Очевидно, что для аппроксимаци подходит суперпози­ ция первого и второго порядка (прямая на рис. 17 и квад­ ратичная парабола на рис. 18). В таблице 20 приведены дан­ ные для обоих случаев.

95

Годы

Производство изделийА шт.тыс( .)

t

 

 

1959

 

3

—4

1960

 

10

—3

1961

 

15

—2

1962

 

21

—1

1963

 

35

0

1964

 

42

I

1965

 

58

2

1966

 

81

3

.1967

 

ПО

4

Итого

375

0

Т а б л и ц а 20

е

lg у

tig У

У1

t*lgy

yt

16

0,47712

—1,90848

5,49

7,63392

3,9

9

1,00000

—3,00000

8,13

9,00000

7,4

4

1,17609

—2,35218

12,02

4,70436

13,2

1

1,32222

— 1,32222

17,78

1,32222

21,7

0

1,54407

0,00000

26,30

0,00000

33,6

1

1,62325

1,62325

38,90

1,62325

48,5

4

1,76343

3,52686

57,54

7,05372

65,5

9

1,90849

5,72547

85,11

17,17641

82,6

16

2,0413°

8,16556

125,90

32,66224

97,3

60

ь,=

Ьг =

81,17612

S у/=

 

= 12,85606

= 10,45826

 

 

=373,7

По формулам (§ 2) имеем

Хо = 0,111111-12,85606= 1,42845 X, = 0,016667-10,45826 = 0,17429

Аналитическое уравнение экспоненциальной суперпозиции по­ лучит вид

y ( i ) = Ю1-434-0'17* или lg [t/(f)] = 1,43 + 0,Ш

По формулам (§ 2) имеем:

Х0 = 0,25541 •12,85606 — 0,021645 -81,17612= 1,5265;

Xt = 0,01667 •10,45826 = 0,17430;

Х ,=—0,02164-12,85606+0,0032467-81,17612=—0,01472.

Уравнение суперпозиции будет:

1Q 1,526+ 0,174/-0,0147/’

или

lg [y(t)] = 1,526 + 0,174* — 0.0147<*

Пример 2. Для определения размера, динамики и струк­ туры национального дохода и анализа факторов его роста произведем выравнивание данных о произведенном нацио­ нальном доходе СССР за период 1960—'1969 >гг-

96


 

 

 

 

 

Т а б л и ц а 21

 

 

Произведенный национальный доход

 

 

 

В фактически

 

 

1п у,

Годы

t

действовав­

\пу

t\пу

ших ценах

=5,25+033/

 

 

(млрд, руб)

 

 

 

1960

—9

145,0

4,9768

—44,7912

4,9561

1961

—7

152,9

5,0297

—35,2079

5,0219

1962

—5

164,6

5,1035

—25,5175

5,0877

1963

—3

168,8

5,1288

—15,3864

5,1535

1964

—1

181,3

5,2001

— 5,2001

5,2193

1965

1

193,5

5,2653

5,2653

5,2851

1966

3

207,4

5,3347

16,0041

5,3509

1967

5

225,5

5,4183

27,0915

5,4161

1968

7

244,1

5,4976

38,4832

5,4825

1969

9

261,7

5,5672

50,1048

5,5483

 

1

 

 

 

 

м

И

О

Ьо =52,5220 6, = 10,8458

Ьг =52,5220

Здесь т= 2 (количество уравнений), п=10 (количество экс­

периментальных точек), k = — =

Если нанести экспериментальные точки (у) на график, используя полулогарифмическую бумагу (в натуральных ло­ гарифмах), то эти точки будут расположены вблизи прямой линии.

Это значит, что для аппроксимации нужно избрать супер­ позицию первого порядка

у =

(3, 23)

Если бы на графике точки оказались расположенными вблизи квадратичной параболы, то для аппроксимации сле­ довало бы избрать суперпозицию второго порядка

y — eK+ht+i.j* .

,(3, 24) и т. д.

По формулам (§ 2) имеем:

>ч, = 0,1 52,5220 = 5,2522;

л,=0,00(30) •10,8458=0/32863.

Подставим значения найденных параметров в уравнение

у = ех°+x>f и получим эмпирическое уравнение экспонен­ ты, выравнивающей уровни национального дохода.

7—984

97