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

Макарова Н.В. Статистика в Excel-1

.pdf
Скачиваний:
3
Добавлен:
11.04.2024
Размер:
11.91 Mб
Скачать

в столбце F с помощью инженерной функции МНИМ. ВЕЩ рассчитаны действительные части комплексных чисел (К^), а в столбце G с помощью инженерной функции МНИМ.ЧАСТЬ вычислены мнимые части комплексных чисел

(YJ.

Действительные и мнимые части рассчитанных в режиме «Анализ Фурье» комплексных чисел связаны с гармоническими коэффициентами следующими соотношениями:

Y Y Y

N N N

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

р^||||в|Я'^^'"^^^:-r-:fW:"^'

«0

0,00

 

щ

0,40

 

ai

0,09

 

^3

-0Л4

 

^4

-0,02

 

<^S

0,04

\

Ч

-0,08

 

^1

0,06

ч--0,27

а^

0,06

^10

-0,08

0,04

«12

-0,02

'^[З

-0,14

«М

0,09

^15

0,20

Таблица 17.7

•Ф^'^ШЩ. ^^^У

Ьу -0,12

ь.0,14

Ьъ

0,16

ЬА

-0,13

Ьь

0,44

Ьв

-0,42

bi

0,13

ь.0,00

Ь, -0,13

Ы0,42

Ьп

-0,44

Ьп

0,13

Ьп

-0,16

^14

-0,14

*15

0,12

340

Для нахождения теоретических значений У( необходимо от ре­ ального времени перейти к «радианному» времени по формуле

L=' N *

В табл. 17.8 приведены значения «радианного» времени (стол­ бец F), теоретические значения первых четырех гармоник (столб­ цы GJ) и их итоговая сумма (столбец К), соответствующая гармо­ нической модели:

у ^ = ^0 + aicost + bisint + Qjcoslt + bismlt + ^зСозЗ/ + b^sixvit +

-f a4cos4/ H- b^svcAt,

Таблица 17.8

isi

п

1

f/' - A ' ^ ; ;?.;•;?:

0

iWSI

2

щШ

3

pli

4

ршл

5

Wn2 6

f:::lrbL

7

\А72Л

8

шла

9

щтЩ

10

МШй

11

[1Ш1

12

0г:л

13

Г ''mi 14

15

 

 

МддЯИ^^Й^ It- V-,.J^^щш^н

 

 

t

^1

U2

 

^3

и.

 

 

0,000

0,402

0,087

 

-0,140

-0,022

0,328

1

0,393

0,324

0,161

 

0,093

-0,131

0,448

1

0,785

0,197

0,141

 

0,211

0,022

0,571

 

1,178

0,040

0,038

 

0,068

0,131

0,277

 

1,571

-0,123

-0,087

 

-0,159

-0,022

-0,391

 

1,963

-0,268

-0,161

 

-0,190

-0,131

-0,750

 

2,356

-0,372

-0,141

 

0,014

0,022

-0,477 1

2,749

-0,419

-0,038

 

0,200

0,131

-0,126

 

3,142

-0,402

0,087

 

0,140

-0,022

-0,197

 

3,534

-0,324

0,161

 

-0,093

-0,131

-0,387 1

3,927

-0,197

0,141

 

-0,211

0,022

-0,246

 

4,320

-0,040

0,038

 

-0,068

0,131

0,060

 

4,712

0,123

-0,087

 

0,159

-0,022

0,173

 

5,105

0,268

-0,161

I

0,190

-0,131

0,165

 

5,498

0,372

-0,141

 

-0,014

0,022

0,239

 

5,890

0,419

-0,038

1

-0,200

0,131

0,312

 

341

На рис. 17.4 представлены эмпирический фафик отклонений урожайности от основной тенденции (рад 1) и теоретические фафики первых четырех гармоник (ряды 2-5). На рис. 17,5 этот же эмпирический фафик показан вместе с итоговым теоретическим фафиком (ряд 6), полученным в результате суммирования первых четырех гармоник (рядов 2—5).

Рис. 17.4

Рис. 17.5

ПРИЛОЖЕНИЕ

Совместное использование режимов надстройки «Пакет анализа»

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

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

Пример П1. Предприятие «Импульс» за месяц произвело 2000 приборов, которым были присвоены заводские номера с 5001 по 7000 включительно. Все приборы изготавливаются по технической доку­ ментации, в соответствии с которой дисперсия чувствительности приборов не превышает 25 мкВ /м .

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

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

Необходимый объем выборки рассчитывается по формуле

343

гдеДд. - предельная ошибка выборки; а^ - дисперсия генеральной совокупности;

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

Для формирования контрольной выборки используем схему слу­ чайного повторного отбора при условии, что предельная ошибка вы­ борки не превышает 3 мкВ/м с уровнем надежности не менее 95 %.

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

r V

1,96^-25

.^^^ 1 W

« ч

п = —г— == г

= 10,68 «11

(приборов).

^\

3^

 

 

Примечание. В расчете необходимого объема выборки используется ко­ эффициент доверия /, для вычисления которого в Microsoft Excel предусмот­ рена функция СТЬЮДРАСПОБР {см. подразд. 6.3.8). Здесь коэффициент до­ верия / определяется по формуле = СТЬЮДРАСПОБР (0,05; 1999), где 0,05 =1—0,95 — требуемый уровень значимости, 1999 = 2000-1 — число степе­ ней свободы.

Таким образом, минимально допустимый объем выборки состав­ ляет 11 приборов. При меньшем объеме выборка не будет репрезен­ тативной.

После определения минимально допустимого объема выборки на рабочем листе Microsoft Excel в диапазон размером в 2000 ячеек (на­ пример, в В1:В2000) введем заводские номера приборов с 5001 по 7000 и сформируем контрольную выборку {см, главу 3).

Для быстрого ввода исходных данных (объем генеральной сово­ купности составляет все же 2000 ед.!) рекомендуем использовать та­ кой технический прием, как копирование ячеек с помощью правой клавиши мыши с последующей установкой через контекстное меню

арифметической прогрессии с шагом L

В результате будет сформирована выборка из 11 приборов с за­ водскими номерами: 5141, 5155, 5349, 5460, 5565, 5706, 5714, 5768, 6501,6771,6972.

344

Отобранные приборы прошли стендовые испытания, на которых были определены тактико-технические характеристики каждого прибора (табл. ПЛ).

 

 

 

 

Таблица ПЛ

 

 

 

 

О ^ : Ш

 

 

Тип

Чувствитель­

Частота

 

Заводской

распозна­

 

ферритовой

ность прибора,

 

номер прибора

ваемого

 

антенны

мкВ/м

 

 

сигнала, МГц

 

 

 

 

 

5141

ФА-77

90

10,07

 

5155

ФА-77М

96

9,73

 

5349

ФА-77

92

10,04

 

5460

ФА-77М

98

9,82

 

5565

ФА-77

86

10,57

Л0--

5706

ФА-77М

 

10,02

•\Ь

5714

ФА-77М

98

9,67

\Ш:

5768

ФА-77

90

9,98

ЛЗ:-

6501

ФА-77

86

10,51

14^

6771

ФА-77

92

9,92

шШЫ

6972

ФА-77М

90

9,93

На основании полученных из контрольной выборки значений характеристик приборов сделаем заключение о чувствительности приборов всей партии. Для этого с помощью режима «Описательная статистика» (см. главу 4) рассчитаем показатели, представленные в табл. П.2,

Во-первых, убедимся, что дисперсия чувствительности прибо­ ров не превышает 25 MKBVM^ (показатель Дисперсия выборки), а предельная ошибка выборки — 3 мкВ/м (показатель Уровень надеж­ ности).

Во-вторых, на основании рассчитанных по контрольной вы­ борке показателей (см, табл. П. 2) с уровнем надежности 95 % мож­ но предположить, что средняя чувствительность приборов всей партии будет находиться в пределах от 88,56 (=91,45-2,89) мкВ/м до 94,34 (=91,45+2,89) мкВ/м (пояснения к расчетам см, в подразд. 4.2).

345

 

 

 

 

 

Таблица п. 2

L -

 

.'• : • • • ' • • • . :

Я

•^i..---I.^::.....i

f .

 

4

/^3

Чувствительность

к - , : "

.•

;:;:•

 

 

[•'••'

.

ь

 

Среднее

91,45

Ь-'-

•'^•

i

Стандартная ошибка

1,30

Р -

 

Н-

--v^

Медиана

90,00

щ:

 

'<•)

'}

Мода

90,00

к

 

-^^

•:

4,30

 

 

 

 

Стандартное отклонение

[ v : .•

 

 

 

Дисперсия выборки

18,47

 

 

 

 

Эксцесс

-0,93

i|-

. .'U •,.'.:; Асимметричность

0,43

i

• .w: • ::!|Интервал

12,00

 

 

 

 

Минимум

86,00

fe-

.16

•••i:

Максимум

98,00

1

 

•17"-. '^,Сумма

1006,00

^:"

1 Я

•'^

Счет

11,00

Ё;;';,

; V

..•

Наибольший(1)

98,00

1 - 2 0

•}

Наименьший(1)

86,00

1-

-21.•'":•(:

Уровень надежности(95,0%)

2,89

В-третьих, коэффициент вариации

v= -^^-100% = --^^.100%«4,7% X 91,45

существенно меньше 40%, что свидетельствует о малой колеблемос­ ти признака в исследованной выборочной совокупности. Надеж­ ность средней подтверждается также и ее незначительным отклоне­ нием от медианы: 91,45 - 90,00 = 1,45.

В-четвертых, незначительное положительное значение коэффи­ циента асимметрии А^ позволяет говорить о том, что данное эмпири­ ческое распределение имеет несущественную правостороннюю асимметрию, а отрицательное значение эксцесса £)t - о его плоско-

346

вершинности, т.е. об отсутствии скопления членов ряда в центре распределения.

Следующим этапом проводимого исследования является уста­ новление степени и характера взаимосвязи предельной частоты рас­ познаваемого прибором сигнала К и его чувствительности X

Для оценки тесноты связи между двумя величинами чаще всего используются коэффициенты ковариации и корреляции (см, главу 13). Результаты расчетов этих коэффициентов приведены соответст­

венно в табл. П.З и П.4.

 

 

 

 

 

Таблица П.З

Ж;

^"^'"^^Щщ^IWT^Wc^W

 

 

 

Чувствительность

Частота

 

Чувствитель­

18,47

 

 

ность

 

 

Частота

-1,05

0,08

 

 

 

Таблица П. 4

pp..,,...,..,,..

••Л,^:,'';;,.;;.Е: Щ

Ц ^^^^щ-Ш?ЩгШМ'-^''''' -••"

'0г:-^:^^И

 

•73 .

 

Чувствительность

Частота

74; J Чувствитель­

1

 

 

ность

 

 

Частота

-0,86

1

Как видим, связь между предельной частотой распознаваемо­ го прибором сигнала У и его чувствительностью X является высо­ кой и обратной (г^у = —0,86), т. е. с повышением чувствительности прибора предельная частота распознаваемого им сигнала умень­ шается.

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

347

Спирмена. Для расчета этого коэффициента используем режим рабо­ ты «Ранг и персентиль» (см, главу 5), результаты выполнения которо­ го представлены в табл, П.5.

Таблица П. 5

 

D

• • • • / ' • : E . •

F

G

H

:1--У\-:^\

;; j " ' . ,.-•.^•.13

: ' 3 5 - ,

Точка Столбец]

Ранг

Процент

Точка

Столбец!

Ранг

процент

isM/::i

4

98

1

90,00%

5

10,57

1

100,00%

Шщ

7

98

1

90,00%

9

10,51

2

90,00%

шЩЩ

2

96

3

80,00%

1

10,07

3

80,00%

рШй

3

92

4

60,00%

3

10,04

4

70,00%

iiiii

10

92

4

60,00%

6

10,02

5

60,00%

•ill

1

90

6

30,00%

8

9,98

6

50,00%

iili

8

90

6

30,00%

11

9,93

7

40,00%

pifi

11

90

6

30,00%

10

9,92

8

30,00%

Iplljlll

6

88

9

20,00%

4

9,82

9

20,00%

iiiiii

5

86

10

,00%

2

9,73

10

10,00%

pill

9

86

10

,00%

7

9,67

11

,00%

По данным сгенерированной табл. П. 5 заполняем графы Ранг R^ и Ранг Ку(та6л. П. 6), на основании которых производим вычисления квадратов разности рангов (d]).

На заключительном этапе вычисляем коэффициент Спирмена по формуле

/2(/7^-1)

подставляя в которую исходные и рассчитанные данные, получим

61^/

а—^1|5_=л75.

р = 1 - /=1

п{п'-1)

11-(11-1)

348

.....:. ^

 

 

 

 

 

Таблица П.6

С

D

"'£'"'.

••••;.• р - :

-

': G^_ ..: .;|

 

 

 

 

'^"Ш-

Номер

Чувствитель­

Частота К,

Ранг

Ранг

 

Квадрат

прибора

ность X,

МГц

 

Ry

 

разности рангов

 

 

мкВ/м

 

 

 

 

dJ={R;,^Ryf

ттт

5141

90

10,07

6

3

 

9

•v^52:?::i

5155

96

9,73

3

10

 

49

,: Щ :-

5349

92

10,04

4

4

 

0

УШ1 5460

98

9,82

1

9

 

64

'/•0й

5565

86

10,57

10

1

 

81

: :S6: .i

5706

88

10,02

6

5

 

1

iill

5714

98

9,67

1

И

 

100

^^тщ

5768

90

9,98

6

6

 

0

йШ::*

6501

86

10,51

10

2

 

64

i i ll

6771

92

9,92

4

8

 

16

1Я§

6972

90

9,93

6

7

 

1

Шн

 

 

 

 

1 =

 

385

Пользуясь шкалой Чеддока (см. подразд. 13Л), можно констати­ ровать, что теснота связи между чувствительностью прибора X и зна­ чением предельной распознаваемой им частоты Уявляется высокой, что подтверждает сделанный ранее вывод.

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

Допустим, что в рассматриваемой ситуации не известен предпо­ лагаемый вид уравнения зависимости предельной частоты Уот чувст­ вительности приборах Для нахождения такого уравнения использу­ ем инструмент «Подбор линии тренда» из мастера диаграмм (см, под-

349