О САЙТЕ
Добро пожаловать!

Теперь вы можете поделиться своей работой!

Просто нажмите на значок
O2 Design Template

ФЭА / АИТ / Лабораторная работа №3 по дисциплине МОДЕЛИРОВАНИЕ СИСТЕМ на тему: «Получение статической модели на основе пассивного элемента»

(автор - student, добавлено - 7-04-2014, 19:57)

СКАЧАТЬ:  lab3.zip [821,72 Kb] (cкачиваний: 32)  

СКАЧАТЬ:  3-laba-2.zip [226,91 Kb] (cкачиваний: 30)

 

 

Лабораторная работа №3

 

по дисциплине МОДЕЛИРОВАНИЕ СИСТЕМ

 

на тему:

 

«Получение статической модели на основе пассивного элемента»

 

 

Регрессионный и корреляционный анализ

1. Нахождение зависимости выходной температуры нефти ТВЫХ от давления нефти Р.

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

 

                                                                                                             Таблица 1

 


Твых, оС 

Р, МПа 

122,79

0,75

122,64

0,74

124,84

0,76

98,58

0,76

85,16

0,78

84,79

0,74

107,46

0,77

117,21

0,77

118,83

0,74

122,05

0,77

123,01

0,79

124,84

0,79

129,10

0,79

120,37

0,80

119,41

0,78

119,34

0,79

132,96

0,79

124,92

0,78

124,84

0,79

124,99

0,78

124,47

0,78

123,37

0,67

125,94

0,63

123,01

0,60

121,69

0,71

124,40

0,75

122,20

0,75

134,35

0,80

117,87

0,80

121,10

0,87

 

 

Весь диапазон изменения х на поле корреляции разобьем на равные интервалы=6. Найдем среднее значение ТВЫХ на каждом интервале по формуле:

,

где у - ТВЫХ, х – P, n – количество точек на интервале. Полученные 6 точек (по одной точке на каждый интервал) последовательно соединяем.

P

Т вых

0,63

98,58

0,68

11,04

0,73

118,34

0,78

118,34

0,83

121,81

0,88

124,84

 

 

 

 

По виду полученной линии регрессии определяем, что зависимость ТВЫХ от F  - линейная. Ниже приведена методика нахождения уравнения линии регрессии одной переменной  методом наименьших квадратов (МНК).

 

Методика нахождения коэффициентов в уравнении линейной регрессии от одного параметра.

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

 

по выборке объемом N.

Система нормальных уравнений для этого случая имеет вид:

 

или

 

 

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

 

 

 

 

Коэффициенты  проще найти по известному  из первого уравнения системы:                                     

,

где  - средние значения .

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

Для оценки линейной связи  вычисляется выборочный коэффициент корреляции r*:

,

где  - выборочные среднеквадратичные отклонения.

 

 

 

Из уравнения имеем:

 

 

Зависимость Твых от P (y = Твых, х = P).

=7129, =34.49, =45.39, =5420.094 N=60.

 

 

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

ТВЫХ=-12.32+173.36·P

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

 

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

Выборочный коэффициент корреляции:

(=854412.5211)

 

 

Параболическая регрессия

Если уравнение регрессии представляет собой полином некоторой степени, то при применении метода наименьших квадратов коэффициенты этого полинома находят решением системы линейных уравнений. Например, требуется определить по методу наименьших квадратов коэффициенты квадратичной функции – параболы второго порядка:

.

В этом случае

;               ;             

и система нормальных уравнений имеет вид:

 

 

 

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

 

Множественная корреляция 

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

.

В нашем случае к=2. То есть здесь мы имеем дело с гиперповерхностью регрессии, или иначе с поверхностью отклика.

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

;               ;                 ,

где yi0, x1i0, x2i0 – нормированные значения соответствующих факторов,

 - средние значения факторов,

sy, sx1, sx2 – среднеквадратичные отклонения.

;            ;       .

 

В новом масштабе имеем:

   ;       

Выборочный коэффициент корреляции при этом равен:

 

 

Уравнение регрессии между нормированными переменными не имеет свободного члена и принимает вид:

.  

Коэффициенты уравнения находятся из условия:

.

Условия минимума функции S определяются так же, как в случае зависимости от одной переменной:

;           …     

и система нормальных уравнений имеет вид:

 

 

Умножим левую и правую части уравнений на .

 

В результате при каждом коэффициенте  получается выборочный коэффициент корреляции r*. Принимая во внимание,  получаем систему нормальных уравнений:

 

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

Решив систему, рассчитывают коэффициент множественной корреляции R:

.

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

 

1. G от Pvh

Для определения вида уравнения регрессии построим эмпирическую линию регулятора по данным  таблицы:

 

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

 

Требуется определить по методу наименьших квадратов коэффициенты нелинейного уравнения регрессии по выборке объемом  N=60.

Находим коэффициенты a0, a1, a2:

 

 

В нашем случае:

 

 

Искомое уравнение регрессии:

 

Сравним графики полученного уравнения регрессии и эмпирической линии:

 

Для оценки нелинейной связи вычисляется коэффициент корреляции:

 

Вычислим :

 

В данном случае коэффициент корреляции =0,863. Т. к. коэффициент корреляции принадлежит интервалу [0,75; 1], то данный входной параметр влияет на выходной и его следует учитывать в уравнении множественной корреляции.

2. G от Tvh

Для определения вида уравнения регрессии построим эмпирическую линию регулятора по данным  таблицы:

 

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

 

 

Требуется определить по методу наименьших квадратов коэффициенты нелинейного уравнения регрессии по выборке объемом  N=60.

Находим коэффициенты b0, b1: 

 

В нашем случае:

 

Искомое уравнение регрессии:

 

Сравним графики полученного уравнения регрессии и эмпирической линии:

 

Для оценки линейной связи вычисляется выборочный коэффициент корреляции:

 

Или

 

Вычислим r:

 

В данном случае коэффициент корреляции =0,82. Т. к. коэффициент корреляции принадлежит интервалу [0,75; 1], то данный входной параметр влияет на выходной и его следует учитывать в уравнении множественной корреляции.

3. G от Gvh

Для определения вида уравнения регрессии построим эмпирическую линию регулятора по данным  таблицы:

 

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

 

 

Для определения коэффициентов в уравнении регрессии воспользуемся функцией MathCad – regress():

 

Искомое уравнение регрессии:

 

Сравним графики полученного уравнения регрессии и эмпирической линии:

 

Для оценки нелинейной связи вычисляется коэффициент корреляции:

 

В данном случае коэффициент корреляции =0,909. Т. к. коэффициент корреляции принадлежит интервалу [0,75; 1], то данный входной параметр влияет на выходной и его следует учитывать в уравнении множественной корреляции.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Вывод

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

 

Приложение 1

t-распределение Стьюдента. Значения t(ν3, q) в зависимости от числа степеней свободы ν3 и уровня значимости q=P[t > tν3,q]

Число степеней свободы k

Уровень значимости q

0,4

0,2

0,1

0,05

0,02

0,01

0,005

0,001

1

2

3

4

5

6

7

8

9

1

1,38

3,08

6,31

12,71

31,82

63,66

318,31

636,62

2

1,06

1,89

2,92

4,30

6,97

9,93

22,33

31,60

3

0,98

1,64

2,35

3,18

4,54

5,84

10,21

12,94

4

0,94

1,53

2,13

2,78

3,75

4,60

7,17

8,61

5

0,92

1,48

2,02

2,57

3,37

4,03

5,89

6,86

6

0,91

1,44

1,94

2,45

3,14

3,71

5,21

5,96

7

0,90

1,42

1,90

2,37

3,00

3,50

4,78

5,41

8

0,89

1,40

1,86

2,31

2,90

3,36

4,50

5,04

9

0,88

1,38

1,83

2,26

2,82

3,25

4,30

4,78

10

0,88

1,37

1,81

2,23

2,76

3,17

4,14

4,59

11

0,88

1,36

1,80

2,20

2,72

8,11

4,02

4,44

12

0,87

1,36

1,78

2,18

2,68

3,06

3,93

4,32

13

0,87

1,35

1,77

2,16

2,65

3,01

3,85

4,22

1

2

3

4

5

6

7

8

9

14

0,87

1,34

1,76

2,15

2,62

2,98

3,79

4,14

15

0,87

1,34

1,75

2,13

2,60

2,95

3,73

4,07

16

0,86

1,34

1,75

2,12

2,58

2,92

3,69

4,02

17

0,86

1,33

1,74

2,11

2,57

2,90

3,65

3,97

18

0,86

1 33

1,73

2,10

2,55

2,88

3,61

3,92

19

0,86

1,33

1,73

2,09

2,54

2,86

3,58

3,88

20

0,86

1,33

1,73

2,09

2,53

2,85

3,55

3,85

21

0,86

1,32

1,72

2,08

2,52

2,83

3,53

3,82

22

0,86

1,32

1,72

2,07

2,51

2,82

3,50

3,79

23

0,86

1,32

1,71

2,07

2,50

2,81

3,48

3,77

24

0,86

1,32

1,71

2,06

2,49

2,80

3,47

3,75

25

0,86

1,32

1,71

2,06

2,48

2,79

3,45

3,73

30

0,85

1,31

1,70

2,04

2,46

2,75

3,39

3,65

40

0,85

1,30

1,68

2,02

2,42

2,70

3,31

3,55

60

0,85

1,30

1,67

2,00

2,39

2,66

3,23

3,46

120

0,84

1,29

1,66

1,98

2,36

2,62

3,16

3,37

 

0,84

1,28

1,64

1,96

2,33

2,58

3,09

3,29

                               

 

 

 

 

 

 

 

 

 

Приложение 2

F -распределение Фише­ра. Значения F(ν4, ν3, q) в зависимости от числа степеней свободы ν4, ν3, и уровня значимости q=Р[F > F(ν4, ν3, q)]=0,05

 

ν3 \ ν4

1

2

3

4

5

6

7

8

9

 

1

2

3

4

5

6

7

8

9

10

 

1

161

200

216

225

230

234

237

239

241

 

2

18,5

19,0

19,2

19,2

19,3

19,3

19,4

19,4

19,4

 

3

10,1

9,55

9,28

9,12

9,01

8,94

8,89

8,85

8,81

 

4

7,71

6,94

6,59

6,39

6,26

6,16

6,09

6,04

6,00

 

5

6,61

5,79

5,41

5,19

5,05

4,95

4,88

4,82

4,77

6

5,99

5,14

4,76

4,53

4,39

4,28

4,21

4,15

4,10

7

5,59

4,74

4,35

4,12

3,97

3,87

3,79

3,73

3,68

8

5,32

4,46

4,07

3,84

3,69

3,58

3,50

3,44

3,39

9

5,12

4,26

3,86

3,63

3,48

3,37

3,29

3,23

3,18

10

4,96

4,10

3,71

3,48

3,33

3,22

3,14

3,07

3,02

11

4,84

3,98

3,59

3,36

3,20

3,09

3,01

2,95

2,90

12

4,75

3,88

3,49

3,26

3,11

3,00

2,91

2,85

2,80

13

4,67

3,80

3,41

3,18

3,03

2,92

2,83

2,77

2,71

14

4,60

3,74

3,34

3,11

2,96

2,85

2,76

2,70

2,65

15

4,54

3,68

3,29

3,06

2,90

2,79

2,71

2,64

2,59

16

4,49

3,63

3,24

3,01

2,85

2,74

2,66

2,59

2,54

17

4,45

3,59

3,20

2,96

2,81

2,70

2,61

2,55

2,49

18

4,41

3,55

3,16

2,93

2,77

2,66

2,58

2,51

2,46

19

4,38

3,52

3,13

2,90

2,74

2,63

2,54

2,48

2,42

20

4,35

3,49

3,10

2,87

2,71

2,60

2,51

2,45

2,39

21

4,32

8,47

3,07

2,84

2,68

2,57

2,49

2,42

2,37

22

4,30

3,44

3,05

2,82

2,66

2,55

2,46

2,40

2,34

23

4,28

3,42

3,03

2,80

2,64

2,53

2,44

2,37

2,32

24

4,26

3,40

3,01

2,78

2,62

2,51

2,42

2,36

2,30

25

4,24

3,38

2,99

2,76

2,60

2,49

2,40

2,34

2,28

26

4,23

3,37

2,98

2,74

2,59

2,47

2,39

2,32

2,27

27

4,21

3,35

2,96

2,73

2,57

2,46

2,37

2,31

2,25

28

4,20

3,34

2,95

2,71

2,56

2,45

2,36

2,29

2,24

29

4,18

3,33

2,93

2,70

2,55

2,43

2,35

2,28

2,22

30

4,17

3,32

2,92

2,69

2,53

2,42

2,33

2,27

2,21

40

4,08

3,23

2,84

2,61

2,45

2,34

2,25

2,18

2,12

60

4,00

3,15

2,76

2,53

2,37

2,25

2,17

2,10

2,04

120

3,92

3,07

2,68

2,45

2,29

2,17

2,09

2,02

1,96

 

3,84

3,00

2,60

2,37

2,21

2,10

2,01

1,94

1,88

ν3 \ ν4

10

12

15

20

24

30

40

60

120

1

242

244

246

248

249

250

251

252

253

254

2

19,4

19,4

19,4

19,4

19,5

19,5

19,5

19,5

19,5

19,5

3

8,79

8,74

8,70

8,66

8,64

8,62

8,69

8,57

8,55

8,53

4

5,96

5,91

5,86

5,80

5,77

5,75

5,72

5,69

5,66

5,63

5

4,74

4,68

4,62

4,56

4,53

4,50

4,46

4,43

4,40

4,36

6

4,06

4,00

3,94

3,87

3,84

3,81

3,77

3,74

3,70

3,67

7

3,64

3,57

3,51

3,44

3,41

3,38

3,34

3,30

3,27

3,23

8

3,35

3,28

3,22

3,15

3,12

3,08

3,04

3,01

2,97

2,93

9

3,14

3,07

3,01

2,94

2,90

2,86

2,83

2,79

2,75

2,71

10

2,98

2,91

2,85

2,77

2,74

2,70

2,66

2,62

2,58

2,54

11

2,85

2,79

2,72

2,65

2,61

2,57

2,53

2,49

2,45

2,40

12

2,75

2,69

2,62

2,54

2,51

2,47

2,43

2,38

2,34

2,30

13

2,67

2,60

2,53

2,46

2,42

2,38

2,34

2,30

2,25

2,21

14

2,60

2,53

2,46

2,39

2,35

2,31

2,27

2,22

2,18

2,13

15

2,54

2,48

2,40

2,33

2,29

2,25

2,20

2,16

2,11

2,07

16

2,49

2,42

2,35

2,28

2,24

2,19

2,15

2,11

2,06

2,01

17

2,45

2,38

2,31

2,23

2,19

2,15

2,10

2,06

2,01

1,96

18

2,41

2,34

2,27

2,19

2,15

2,11

2,06

2,02

1 97

1,92

19

2,38

2,31

2,23

2,16

2,11

2,07

2,03

1,98

1,93

1,88

20

Ключевые слова -


ФНГ ФИМ ФЭА ФЭУ Яндекс.Метрика
Copyright 2021. Для правильного отображения сайта рекомендуем обновить Ваш браузер до последней версии!