Главная страница
Навигация по странице:

  • , .

  • 1 вар записка. Корреляционный и регрессионный анализ экспериментальных данных


    Скачать 136.69 Kb.
    НазваниеКорреляционный и регрессионный анализ экспериментальных данных
    Дата19.11.2020
    Размер136.69 Kb.
    Формат файлаdocx
    Имя файла1 вар записка.docx
    ТипДокументы
    #151814

    МИНИСТЕРСТВО ОБРАЗОВАНИЯ И НАУИ РОССИЙСКОЙ ФЕДЕРАЦИИ

    Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования

    «ВОЛГОГРАДСКИЙ ГОСУДАРСТВЕННЫЙ ТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ»

    (ВолгГТУ)

    Кафедра «Технология машиностроения»

    Семестровая работа

    по дисциплине «Математическое моделирование

    технологических процессов

    тема: «Корреляционный и регрессионный анализ

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

    Вариант № 1


    Выполнил:

    студент гр. КТО-325

    Василенко Д.И.

    Проверил:

    Фролов Е. М.

    Волгоград 2018 г.

    Цель работы: провести корреляционный и регрессионный анализ результатов эксперимента.
    Задачи:

    1. Оценить степень тесноты взаимосвязи значений функции отклика со значениями независимой переменной. Определить достоверность данной оценки.

    2. Проанализировать условия применимости статистического анализа.

    3. Получить линейную регрессионную модель зависимость функции отклика от независимой переменной. Оценить значимость коэффициентов регрессии и адекватность уравнения регрессии.


    Задание:

    Таблица 1

    Вариант 1

    x

    y1

    y2

    y3

    1,12

    18,14

    18,29

    18,4

    2,24

    19,35

    19,5

    19,61

    3,36

    28,51

    28,66

    28,77

    4,48

    25,96

    26,11

    26,22

    5,6

    33,88

    34,03

    34,14

    6,72

    36,8

    36,95

    37,06

    7,84

    43,58

    43,73

    43,84

    8,96

    44,4

    44,55

    44,66

    10,08

    51,18

    51,33

    51,44

    11,2

    50,2

    50,35

    50,46

    12,32

    58,06

    58,21

    58,32

    13,44

    58,1

    58,25

    58,36

    14,56

    65,81

    65,96

    66,07

    15,68

    69,53

    69,41

    69,79

    16,8

    75,49

    75,37

    75,75

    17,92

    76,42

    76,3

    76,68

    19,04

    84,77

    84,65

    85,03

    20,16

    82,07

    81,95

    82,33

    21,28

    91,33

    91,21

    91,59

    22,4

    91,61

    91,49

    91,87

    23,52

    101,96

    101,84

    102,22

    24,64

    99,58

    99,46

    99,84

    25,76

    106,55

    106,43

    106,81

    26,88

    109,17

    109,05

    109,43

    28

    115,62

    115,5

    115,88

    29,12

    114,65

    114,53

    114,91

    30,24

    124,85

    124,73

    125,11

    Продолжение таблицы 1

    x

    y1

    y2

    y3

    49,28

    190,94

    190,89

    190,62

    50,4

    199,54

    199,49

    199,22

    51,52

    198,97

    198,92

    198,65

    52,64

    204,63

    204,58

    204,31

    53,76

    206,58

    206,53

    206,26

    54,88

    216,58

    216,53

    216,26

    56

    213,57

    213,52

    213,25

    57,12

    224,83

    224,78

    224,51

    58,24

    221,76

    221,71

    221,44

    59,36

    232,22

    232,17

    231,9

    60,48

    231,13

    231,08

    230,81

    61,6

    239,92

    239,87

    240,42

    62,72

    240,3

    240,25

    240,8

    63,84

    248,64

    248,72

    249,14

    64,96

    248,29

    248,37

    248,79

    66,08

    256,66

    256,74

    257,16

    67,2

    256,64

    256,72

    257,14

    68,32

    264,33

    264,41

    264,83

    69,44

    265,27

    265,35

    265,77

    70,56

    272,31

    272,39

    272,81

    71,68

    272,76

    272,84

    273,26

    72,8

    280,07

    280,15

    280,57

    73,92

    280,21

    280,29

    280,71

    75,04

    289,77

    289,85

    290,27

    76,16

    287,12

    287,2

    287,62

    77,28

    296,73

    296,81

    297,23

    78,4

    298,26

    298,34

    298,76


    1. Для оценки степени тесноты взаимосвязи значений функции отклика со значениями независимой переменной используется коэффициент парной корреляции, рассчитываемый по формуле ( – объем выборки; и – средние арифметические значения и в рассматриваемой выборке; , – стандартные отклонения и в рассматриваемой выборке).

    По результатам расчета коэффициент парной корреляции , что больше 0,5. Соответственно можно сделать вывод о наличии тесной взаимосвязи значений функции отклика со значениями независимой переменной.

    Так как объем рассматриваемой выборки N=69, то для оценки достоверности рассчитанного коэффициента корреляции проверяется нулевая гипотеза , согласно которой коэффициент корреляции генеральной совокупности равен нулю, посредством критерия Стьюдента. Результаты расчета приведен в таблице 2.
    Таблица 2

    Оценка достоверности

    к-та корреляции

    t расч

    292,1650045

    t критич

    1,667916115

    Вывод

    достоверен


    Поскольку критическое значение критерия Стьюдента при уровне значимости 0,05 и числе степеней свободы =69–2=67 меньше расчетного, то нулевую гипотезу об отсутствии линейной зависимости между значениями функции отклика и значениями независимой переменной следует отвергнуть и признать коэффициент парной корреляции достоверным.

    2. Применение корреляционного и регрессионного анализа правомерно и эффективно при соблюдении следующих условий:

    2.1. Параметр оптимизации – случайная величина с нормальным законом распределения. Для проверки используется критерий Пирсона (N – число интервалов, по которому строится эмпирический закон распределения; – вероятность попадания значения случайной величины в i-й интервал для теоретического закона распределения; – вероятность попадания значения случайной величины в i-й интервал для эмпирического закона распределения). Результаты расчетов приведены в таблице 3.
    Таблица 3

    Количество
    интервалов

    № интерв

    Частота

    Вероятность эксп

    Вероятность теор

    Хи практич

    9

    1

    8

    0,116

    0,071

    10,74056

     




    2

    8

    0,116

    0,065

    Хи критич

     




    3

    8

    0,116

    0,097

    48,30538

     




    4

    8

    0,116

    0,127

    Вывод

     




    5

    6

    0,087

    0,145

    соответствует

     




    6

    8

    0,116

    0,144

     

     

     

    7

    7

    0,101

    0,125

     







    8

    8

    0,116

    0,095










    9

    8

    0,116

    0,063





    Для ответа на вопрос о принятии нулевой гипотезы о том, что данная случайная величина подчиняется нормальному закону распределения, находится критическое значение критерия Пирсона в зависимости от уровня значимости 0,05 и числа степеней свободы k=s-1-r=9-1-2=6 (s – число интервалов, r – количество параметров предполагаемого распределения. В ходе сравнения наблюдаемого значения критерия Пирсона с теоретическим установили, что наблюдаемое значение критерия меньше теоретического, а соответственно рассматриваемая случайная величина распределена по нормальному закону.

    2.2. Дисперсии результатов параллельных опытов должны быть однородны.

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

    Оценка однородности дисперсий

    Критерий Кохрена практич

    0,162687512

    Критерий Кохрена теоретич

    0,189

    Вывод

    однородны


    Табличное значение критерия Кохрена определяется в зависимости от уровня значимости 0,05 и количества степеней свободы f=k-1=3-1=2 (где k ­ количество параллельных опытов). Поскольку расчетное значение критерия Кохрена меньше теоретического, дисперсии считаются однородными.

    2.3. Значения независимых переменных должны изменяться с пренебрежимо малыми ошибками по сравнению с ошибкой в определении функции отклика.

    Считаем, что эксперимент проводится на исправном оборудовании в соответствующих условиях и поэтому ошибки в задании независимых переменных пренебрежимо малы по сравнению с ошибками в определении функции отклика.

    2.4. Независимые переменные не должны коррелировать друг с другом.

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

    2.5. Процесс изменения зависимой переменной должен быть стационарным и случайным.

    Для проверки стационарности процесса строится случайная последовательность значений функции отклика, соответствующая порядку проведения измерений и разбивается на 7 равных отрезков. Для каждого из них рассчитывается дисперсию и с помощью критерия Кохрена устанавливается, являются ли значения дисперсий на каждом из отрезков оценками одной и той же генеральной дисперсии. Результат расчета приведен в таблице 5. Поскольку расчетное значение критерия Кохрена меньше теоретического для уровня значимости 0,05 и количества степеней свободы в знаменателе f=k-1=7-1=6 (где k ­ количество опытов в отрезке), дисперсии на каждом из отрезков являются оценками одной и той же генеральной дисперсии.


    Таблица 5

    Дисперсия по отрезку

    Критерий Кохрена практич

    101,657

    0,162687512

    104,679

    Критерий Кохрена теоретич

    99,173

    0,189

    93,488

    Вывод

    156,098

    однородны

    109,879

     

     

    131,992

    108,324

    54,210


    Затем проверяется гипотеза о равенстве всех генеральных средних. Чтобы она была справедлива, должно соблюдаться условие

    ( ; – табличное значение критерия Фишера; – общее среднее всех элементов отрезков при объединении их в одну выборку). Результаты расчета приведены в таблице 6.
    Таблица 6

    Среднеарифм по отрезку

    Критерий Фишера практич

    31,464

    9,228

    63,202

    Критерий Фишера теоретич

    95,927

    3,731589458

    128,311

    Вывод

    162,448

    не равны

    198,797

     

     

    231,531

    264,735

    290,611


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

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

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

    Данная гипотеза проверяется методом последовательных разностей (Критерий Аббе-Линника). Для чего по значениям выборки, расположенным в последовательности их наблюдения, образуется разностей между соседними членами: , . Затем рассчитывается критерий , в котором: , .

    Если ( – критическое значение при уровне значимости 0,95), то гипотеза о случайности выборки верна. Результаты расчета представим в таблице 7.
    Таблица 7

    тау расч

    0,002579371

    тау теор

    0,7759

    Вывод

    выборка случайна


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

    3. Коэффициенты линейной регрессии рассчитываются по формулам

    , . В рассматриваемом случае , и уравнение линейной регрессии выглядит следующим образом: . Построим график для сравнения измеренных значений функции отклика со значениями, полученными по уравнению регрессии (рис. 1).

    Значимость коэффициентов регрессии проверяется по критерию Стьюдента. Для этого рассчитывается ( – стандартное отклонение коэффициента регрессии) для каждого коэффициента и сравнивается с критическим значением критерия Стьюдента для уровня значимости 0,05 и числа степеней свободы . Коэффициент считается значимым, если . Так как во всех строках матрицы плана число параллельных опытов одинаково, то ( –дисперсия воспроизводимости). Поскольку дисперсии параллельных опытов однородны (установлено в пункте 2.2), то дисперсия воспроизводимости эксперимента равна ( – дисперсии результатов параллельных опытов).

    Результаты расчета представлены в таблице 8.
    Таблица 8

    Дисперсия воспроизводимости

    157,455

    Проверка b0

    t расч

    13,996

    t теор

    1,652

    Вывод

    к-т значим

    Проверка b1

    t расч

    4,196

    t теор

    1,656

    Вывод

    к-т значим


    На основе расчета можно сделать вывод, что оба коэффициента линейной регрессии значимы.

    Адекватность уравнения регрессии проверяется сравнением дисперсии выборки с дисперсией адекватности. Если дисперсия адекватности существенно отличается от дисперсии выборки, значит уравнение регрессии адекватно. Для оценки адекватности модели вычисляется опытное значение критерия Фишера: . Дисперсия выборки рассчитывается по формуле , а дисперсия адекватности –

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

    Результаты расчета представлены в таблице 9.
    Таблица 9

    Дисперсия адекватности

    5,621

    Критерия Фишера расч

    1256,285

    Критерий Фишера
    теоретич

    1,497

    Вывод

    уравнение адекватно


    На основе расчета можно сделать вывод, что уравнение регрессии дает адекватные результаты.


    Рисунок 1


    написать администратору сайта