Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Шепелев, И. Г. Математические методы планирования и управления в строительстве конспект лекций

.pdf
Скачиваний:
12
Добавлен:
19.10.2023
Размер:
5.9 Mб
Скачать

мулам (2.4.5) и (2.4.6), определяем г| = 0,77 и е = 12,81. Кор­ реляционное отношение ц = 0,77, что больше г = 0,197, опре­ деленного для этих статистических данных: это означает, что степенная линия регрессии лучше аппроксимирует зависимость между выработкой и текучестью рабочей силы, чем прямая ли­ ния.

§ 2.5. Логарифмическая зависимость

Логарифмическая зависимость выражается формулой вида:

у = а + b lgx.

(2.5.1)

График логарифмической функции представлен на рис. 6. Для получения параметров логарифмической кривой нужно прологарифмировать наблюдения по х и, рассматривая их как

Рис. 6. График логарифмической функции

независимые переменные, определить параметры а и b по мето­ ду наименьших квадратов. При последующем использовании кривой для определения зависимой переменной необходимо ло­ гарифмировать независимую переменную и подставить ее зна­ чения в уравнение (2.5.1).

Рабочие формулы для определения параметров а и b имеют

вид:

 

 

 

где

А =

Е у-Е (Igx)2 — Е lgx-Ey-lgx;

 

 

B =

N E y l g x — Elgx-Ey;

(2.5.2)

 

Д =

N E (lg x)2 — (E lg x ) 2.

 

20

Корреляционные отношения р и е вычисляются но приведенным в предыдущем параграфе формулам.

Форма записи исходных данных для вычисления параметров линии регрессии и оценок этой линии приведены в табл. 5.

Таблица 5

 

■N-

 

'V»

У X lg X Уlg х lg

(У-У)2 (У-У)2

У-У

У-У

у-7

У

§ 2.6. Параболическая зависимость или многочлен n-й степени

Зависимость выражается формулой:

 

у = а + Ьх -ф-сх2.

(2.6.1)

Если степень независимого переменного равна двум, то это парабола второго порядка и т. д. Прямая зависимость является частным случаем многочлена.

Аппроксимация (подбор) параболической кривой осущест­

вляется также методом наименьших квадратов.

В целевую

функцию метода наименьших квадратов

£ (у — у)2-* min вмес­

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

уравнения

( 2. 6. 1)

-> min.

(2.6.2)

S = £(y — а Ь х сх-)

Возьмем частные производные от этого выражения по а, b и с

= — 2 £ (у — а — Ьх — сх2) = 0;

da

= — 2£ (у — а Ьх — сх2) х 0;

db

= — 22 (у — а — Ьх сх2) х 2 = 0.

dc

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

Na + 6£.x: + c£.x:2 — £у;

аЪ х ЬТ, х 2

ch х 3 = !>ух\

(2.6.3)

+

= £ у * 2.

 

21

Решая систему, находим любым известным методом параметры

параболы а, Ь и с.

 

(2.6.3) записать в матричной форме,

Если систему уравнений

a

b

c

 

N

2x

2x2

(2.6. 4)

l x

2д:2

2x3

2yx

2x2

2x3

2x*

2yx3

то параметры параболы можно определить из выражений:

_ Дь .

(2.6.5)

 

Д ’

дс

~д ’

где Д — определитель матрицы линейных уравнений; Да — определитель матрицы линейных уравнений, в котором

столбец, имеющий в своем составе «а», заменяют столбцом свободных членов;

Дь — определитель матрицы, в котором столбец, имеющий - в своем составе «Ь», заменен столбцом свободных чле­ нов;

Дс — определитель, в котором столбец, имеющий в своем со­ ставе «с», заменен столбцом свободных членов.

Определители матрицы (2.6.4) можно расписать в виде сле­

дующих выражений:

 

 

Д =

N £ х 2 Ел:4 +

ЕхЕх3 £х2 + £ х £ х2 £ х 3 -

 

-

£ х 2 £ x2 £ x 2 - (Ex)2 S x 4 — N (Ex3)2.

(2.6.6)

J X = Ey E x2 Ex4 + Ex3 Ex £x2y + E xy Ex2 Ex3 -

 

-

£ x 2y (Ex2)2 - Exy ExEx4 — E у (E x 3)2.

(2.6.7)

R b = N E xy Ex4 +

Ey Ex3 Ex2 + Ex Ex2 E x 2y -

 

- N Ex3 E x 2y -

Ex4 Ex Ey - (Ex2)2 E xy.

(2.6.8)

Д с = N Ex2 Ex2y + Sx Ex3 Ey -f Ex Ex2 Exy —

 

-

N Ex3 Exy -

E x 2y Ex2 Ey — (S x2)2 Ey.

(2.6.9)

22

 

 

Т а б л и ц а 6

У

X X2 X 3 Л-‘

ху

Х2У

Форма записи и обработки исходных данных для параболи­ ческой зависимости представлена в табл. 6.

Оценка точности аппроксимации параболы также произво­ дится по корреляционному отношению rj и ошибке аппроксима­

ции е.

§ 2.7. Корреляционные зависимости периодического вида

Кривые периодического вида могут найти широкое распрост­ ранение при аппроксимации зависимостей многих экономичес­ ких явлений во времени. Наблюдения времени могут быть пред­ ставлены в виде равноотстоящих переменных х, выраженных & тригонометрической форме.

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

■2ic;

(2.7.1>

12

В течение года можно получить 12 наблюдений экономического

показателя у и у2, уз, .... у 12. Тогда зависимость

величины у от

времени х можно выразить уравнением:

 

m

 

у — а0-ф- 2 (ак cos kx -f- Ьк sin kx),

(2.7.2)'

k=l

 

здесь Go, Gk и &k — коэффициенты линии регрессии, число их рав­

но

2m + 1. Если N > 2 т + 1, коэффициенты Gk и Ьк находятся

по

методу наименьших квадратов. Целевая функция имеет вид:

 

N

m

(2.7.3)

 

S = 2 [у — ао— 2 (ak cos + sin &х )2 “*■min-

 

N=1

к—1

 

 

Для

вычисления

неизвестных параметров уравнения а0,

Gk,

и т.

д. необходимо продифференцировать выражение (2.7.3).

в

частных производных относительно ао, Gk, Ьк и т. д., приравнять, полученные производные нулю и составить систему ортогонадь-

2$

ных уравнений. При этом необходимо принять во внимание, что функции 1, cos kx, sin kx, ( k — 1, m),

. N

где m — составляют систему ортогональных по отношению к

ряду (2.7.1) функций, а сумма

 

 

N

 

2N-*tc

О, при k Ф /;

 

\

 

2N-1-*

. / , п , N

(2.7.4)

,cos—- — cos-

N

N=1

N

—, при

к = 1 ф 0 ф

 

 

 

 

 

Благодаря этому свойству решение нормальных уравнений ока­ зывается сравнительно простым [3]. Для вычисления парамет­ ров уравнения (2.7.2) имеем:

N

аО-— У-

N

у cos kx,

(2.7.5)

N

*k ^

»

sin^ .

N- l

Вкачестве примера составим корреляционное уравнение за­ висимости поставок леса от времени года. Данные подекадных поставок в процентах к плану представлены в табл. 7.

Если зависимость аппроксимировать тригонометрической кривой с k =■4, то расчетные данные удобно расположить в ви­ де табл. 8. Подставим данные табл. 7 в форму табл. 8, просум­

мируем столбцы 2, 3, 4, 5, 6, 7, 8, 9, 10. Получим суммы:

Ss =

+ 456,8

= + 388,698 S 4 = — 628,962 Sr, = — 704,669

Se =

- 495,376

= — 829,764 Se =

- 208,947 So =

- 1039,841

S io -= + 531,813.

 

 

 

 

В соответствии с формулами (2.7.5)

определим параметры урав­

нения периодического типа

 

 

ап

456,8

=

12,68

 

 

36

 

 

 

 

 

 

 

 

а I — 388,698

=

21,59

628,962

- 34,94

 

 

18

 

 

18

 

24

 

 

 

 

 

 

 

 

 

 

 

Т а б л и ц а 7

Наблю­

x ,

У

Наблю­

X ,

У

Наблю­

х ,

У

Наблю­

X ,

У

дения

ДН И

дения

дни

дения

дни

дения

дни

i

10

- 1 7 ,8

2

20

- 4 5 ,0

3

30

+60,0

4

40

+ 16,2

5

50

+47,5

6

60

—45,8

7

70

+ 13,2

8

80

+26,2

9

90

+81,0

X

У

X

ч

COS

с

 

<*3

 

 

у

 

1

2

3

4

10

100

—51,5

19

190

 

- 4 7 ,7

28

280

+96,0

и

п о

- 3 8 ,5

20

200

 

—44,7

29

290

+ 115,0

12

120

—14,5

21

210

 

+ 136,0

30

300

+330,0

13

130

+47,0

22

220

 

- 9 ,8

31

310

+89,5

14

140

—24,6

23

230

 

+ 16.0

32

320

+260,0

15

150

- 9 ,7

 

24

240

 

—25,0

33

330

+62,6

16

160

- 5 ,1

 

25

250

 

+ 10,1

34

340

—100,0

17

170

—46,5

26

260

 

+43,8

35

350

-100,0

18

180

- 8 2 ,0

27

270

 

—93,0

36

360

-100,0

 

 

 

 

 

 

 

 

 

 

Т а б л и ц а 8

ч

Ч

Ч

Ч

Ч

ч

 

 

 

 

 

 

СМ

со

тГ

 

 

 

 

 

1У—У

 

со

Tf*

 

 

 

 

 

С<3

с

*5

С

«О

«О

У

(у-у) (У-У)2

(У-У)

(У-У)3

г г

с>

 

О

 

О

Ss

 

 

Ss

 

 

 

 

 

 

5

6

7

8

9

10

и

12

13

14

15

16

t

I.

2,5

S7

2S.

X9

"io

2 »

2 13

S l4

^15

216

а 2 — -

704,669

39,14

b2 =

-

495,376

=

— 27,52

 

 

18

 

 

 

18

 

 

а я =

-

829,764 -

- 46,09

b3=

-

208,947

=

- 11,60

 

 

18

 

 

 

18

 

 

a k =

-1039,841

— — 57,76

b* =

531,813 _

29,54,

 

 

18

 

 

 

18

 

 

Уравнение имеет вид:

 

 

у =

12,68 + 21,59 cos л: — 34,94 sin л: — 39,14 cos 2л: —

— 27,52

sin 2л: — 46,09 соэЗл:— 11,60

эшЗл:— 57,76

cos 4л; +

 

+ 29,54 sin 4л:.

 

(2.7.6)

Подставив значения cos л:, sin л: в уравнение (2.7.6), получим

расчетные значения зависимого переменного у. Заполнив ими столбец 11, табл. 8 и рассчитав значения для столбцов 12—16, получим необходимые данные для вычисления корреляционно­

го отношения г| и ошибки аппроксимации е.

Вычислим

ri

S (у - 3 + _

. / j _

13597,3160

= 0,97

S ( y _ + 2 ~

К

301230,1945

 

 

-•5,69645-100% = 15,84%.

Г Л А В А I I I

МЕТОДЫ ОЦЕНКИ

§ 3.1. Генеральная и выборочная совокупности

Группа предметов или явлений, объединенных каким-либо общим признаком или свойством качественного или количест­ венного характера, называется статистической совокупностью. Например, отчетные данные о выработке на одного трудящего­ ся, о себестоимости, уровне механизации и т. д.

Генеральная совокупность — это такая случайная совокуп­ ность, когда можно предположить, что группа объединяет очень большое, неограниченное число наблюдений. Применительно к экономическим явлениям генеральная совокупность обозначает наблюдения по всем возможным предприятиям с родственной технологией за длительный период времени, в том числе и не­ который период времени в будущем.

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

Чтобы правильно судить о всей совокупности, необходимо иметь представительную выборочную совокупность. Представи­ тельная выборка называется репрезентативной. Например, если при исследовании выработки взять только передовые предприя­ тия, то результаты исследования не будут характерны для всех предприятий и по средним значениям, полученным из такой вы­ борки, невозможно судить о всех предприятиях.

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

Вэкономике наиболее представительной является такая вы­

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

ченные по выборке, истинным характеристикам и взаимосвязям, которые могли быть получены при изучении генеральной сово­ купности.

В практике часто приходится решать задачу: достаточно ли число наблюдений в выборке для того, чтобы судить по ней о генеральной совокупности.

Если обозначить:

N — число наблюдений в выборке,

S — среднеквадратичное или стандартное отклонение, полу­ ченное по результатам выборки,

а — вероятность приближенного равенства о0 ~ S, где Оо — стандарт генеральной совокупности,

е — точность наблюдений,

k = N — 1— число степеней свободы,

qs — аргумент функции распределения вероятности. Аргумент qs определяется из выражения е = qs • S. Из функции L (qs&) = а вероятность а зависит только от q и k. Можно с за­ данной вероятностью а построить доверительный интервал, в котором будет находиться стандартное отклонение а. Этот ин­ тервал можно записать в виде следующего соотношения:

 

 

 

Р =

( S - q sS < 3o< S + qs S}=«.

 

(3.1.1)

 

 

 

L(qs6)

 

 

 

 

 

Пример определения достаточного числа наблюдений. Опре­

делить

Идост при

заданной

вероятности а =

0,95,

S = 0,12,

е =

0,06.

Определяем qs —

 

= 0,5. По таблице функ­

ции

L(qs&),

(приложение

1)

находим k = 14,

следовательно

'N =

k +

1 =

15.

 

 

 

 

 

Такой метод оценки репрезентативности выборки пригоден

для

любого

N. Если N >

20,

то можно пользоваться

упрощен­

ным методом, основанным на гипотезе нормального распределе­ ния случайной величины. При этом

 

N >

t‘2 S2

(3.1.2)

 

- ^ ,

 

 

£2 X2

 

где ta — аргумент,

характеризующий вероятность нормального

распределения в

интегральной

функции

распределения Ф (t);

ta определяется по таблицам

функции

Ф (t) (приложение II)

взависимости от заданной вероятности а.

х— среднее значение случайной величины в выборке. Если выполняется соотношение (3.1.2), то число наблюдений доста­ точно.

28

§ 3. 2. Оценка представительности коэффициента корреляции

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

Ошибка выборочного коэффициента корреляции определя­ ется по формуле

<3 -2 Л )

У N - 1.

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

■ ■ г-

.

(3.2.2)

] / n -

1

 

Полагая коэффициент корреляции распределенным нормально,

можно утверждать, что

 

 

 

Р {г — К ar < р < г + t, зг) = 2 Ф (t),

(3.2.3)

при Р =

0,9

t =

1,643,

 

Р =

0,95

t =

1,96,

 

Р =

0,99

t =

2,576.

 

Особенно интересна проверка так называемой нулевой гипо­ тезы. Известно (см. § 2.2), что если коэффициент корреляции по модулю больше нуля, то между двумя случайными величинами имеется связь. Однако г, определенный по частичной выборке, отличается от истинного коэффициента корреляции р. Может быть, что при | г j > 0 , р — 0 тогда связь, установленная по час-, тичной выборке, в генеральной совокупности отсутствует. Для ответа на вопрос: есть ли связь в генеральной совокупности, за­ даются гипотезой, что такой связи нет и проверяют эту гипотезу.

Будем считать, что р = 0, тогда

ор = —===г.

 

 

 

I / n -

г .

 

Р jг -

~ 7 = = r < Р <

г + ,

) = а.

(3.2.4)

I

>AN - 1

J/N —1 J

 

При проверке нулевой гипотезы интересен только нижний пре­ дел

г -------^ == < р;

(3.2.5)

/ N - 1

 

29

Соседние файлы в папке книги из ГПНТБ