Оборудование, технические и инструментальные средства
Для выполнения работы необходим персональный компьютер с установленной операционной системой MsWindows. Используется программа для ПК «ЛАБ1. Виртуальная лабораторная работа №1 по дисциплине МИОУ. Определение статической характеристики объекта». Описание хода выполнения работы. Исходные данные
1) Запустили программу ЛАБ1, зарегистрировались используя номер своей зачетной книжки. 2) Выбрали канал, по которому будет сниматься статическая характеристика (Y = f(T), Y = f(Q) или Y = f(F)). 3) Разбили интервал варьирования входной переменной на 10-15 равных частей. 4) Последовательно изменяя значение входной переменной от минимального до максимального, регистрировали значение выходной переменной Y. 5) Выполнили все необходимые с использование систем MsExcel и Mathcad.
Таблица 1.1 - Исходные данные
Результаты выполнения
Требуется аппроксимировать экспериментальные данные аналитическими зависимостями видов: 1)y = A+B·x 2)y = A·xB 3)y = A·10x·B 4)y = 1/(A+B·x) 5)y = x/(A+B·x) 6)y = A+B·lg(x) 7)y = A+B/x Для определения коэффициентов А и В воспользуемся (1.2), (1.3), а для нахождения критерия адекватности уравнением (1.4). Для этого предварительно (если необходимо) линеаризуем аппроксимирующие зависимости путем замены переменных, а затем вычислим нужные суммы: 1) y = A+B·x: y = y, x = x ΣXi= 2112; ΣYi= 536,02; ΣXi2= 416504; ΣXiYi = 109078,94 A = -58,845018B = 0,5602818 2) y = A·xB: y = lg(y), x = lg(x) ΣXi= 25,049892; ΣYi= 18,136324; ΣXi2= 57,103683; ΣXiYi = 41,454761 A = -4,3303918B = 2,6255856 3) y = A·10x·B: y = lg(y), x = x ΣXi= 2112; ΣYi= 18,136324; ΣXi2= 416504; ΣXiYi = 3547,229 A = 0,5132557B = 0,0059141 4) y = 1/(A+B·x): y = 1/y, x = x
ΣXi= 2112; ΣYi= 0,2787325; ΣXi2= 416504; ΣXiYi = 49,282514 A = 0,099244B = -0,0003849 5) y = x/(A+B·x): y = 1/y, x = 1/x ΣXi= 0,0589254; ΣYi= 0,2787325; ΣXi2= 0,0003248; ΣXiYi = 0,0016246 A = -0,0516398B = 14,370205 6) y = A+B·lg(x): y = y, x = lg(x) ΣXi= 25,049892; ΣYi= 536,02; ΣXi2= 57,103683; ΣXiYi = 1234,9894 A = -509,19191B = 244,99631 7) y = A+B/x: y = y, x = 1/x ΣXi= 0,0589254; ΣYi= 536,02; ΣXi2= 0,0003248; ΣXiYi = 2,6916826 A = 153,98126B = -19648,132
Таблица 1.2 – Результаты расчетов
На рисунке 1.3 показаны аналитические зависимости, описывающие эксперимент, а также показаны результаты аппроксимации данных, зависимость Yi от температуры, где видно, что график зависимости Y7 от температуры наиболее приближен к графику кривой, построенной по экспериментальным данным. Рисунок 1.3 - Аналитические зависимости, описывающие эксперимент, результаты аппроксимации данных Выводы
В ходе выполнения всех необходимых расчетов, сведения их в таблицу и построения нужных графических зависимостей стало понятно, что наиболее оптимальным образом экспериментальные данные описывает модель вида: y = A+B/x: y = y, x = 1/x
Лабораторная работа №2. «Однофакторный дисперсионный анализ объекта»
Цели и задачи работы
Целью работы является проверка: требуется проверить гипотезу о равенстве математических ожиданий для каждой конкретной серии опытов. А также научиться применять знания, полученные при изучении темы «Определение параметров регрессионной модели. Однофакторный дисперсионный анализ объекта управления».
Содержание работы. Теоретические основы
В любом эксперименте среднее значение наблюдаемых величин меняется в связи с изменением входных факторов, определяющих условия эксперимента, а также и случайных факторов (помех). Исследование влияния тех или иных факторов на изменчивость средних значений является задачей дисперсионного анализа. Дисперсионный анализ состоит в выделении и оценке отдельных факторов, вызывающих изменчивость изучаемой случайной величины. Для этого производится разложение суммарной выборочной дисперсии на составляющие, обусловленные независимыми факторами. Для того чтобы определить значимо ли влияние данного фактора необходимо оценить значимость соответствующей выборочной дисперсии в соответствии с дисперсией воспроизводимости, обусловленной случайными факторами. Предположим, что результат эксперимента зависит от некоторого одиночного фактора А, который принимает n различных значений (n-количество серий опытов). Для каждой серии опытов проводится m повторных наблюдений, результаты которых можно записать в следующем виде:
Y11 Y12 Y13... Y1m Y21 Y22 Y23... Y2m Y31 Y32 Y33... Y3m ............... Yn1 Yn2 Yn3... Ynm
На основе полученных статистических данных требуется проверить гипотезу о равенстве математических ожиданий для каждой конкретной серии. Если проверяемая гипотеза верна, то средние арифметические значения для всех серий практически не отличаются друг от друга, в противном случае предполагаемая гипотеза должна быть отвергнута. Обозначим через среднее значение i-й серии опытов, а через общее среднее значение для всех наблюдений:
Сущность дисперсионного анализа состоит в разложении суммы квадратов отклонений отдельных Yij от общего среднего на две суммы: Q - определяет общее отклонение значения каждого опыта (Yij) от среднего; QА - характеризует рассеяние, вызванное фактором А (выражение во вторых фигурных скобках); Qост - характеризует рассеяние вызванное случайными помехами (выражение в первых фигурных скобках).
Разделив суммы квадратов отклонений на соответствующие степени свободы получим следующие дисперсии:
σ2 = Q/f σА2= QA/f1 (2.2) σОСТ2 = Qост/f2
Число степеней свободы f = m·n -1 f1 = n - 1 f2 = n·(m-1) Проведение дисперсионного анализа состоит в сравнении оценок σА2и σОСТ2. Если гипотеза о том, что математические ожидания для каждой серии равны, верна, то σ А2 не должна существенно превышатьσОСТ2, что проверяется по критерию Фишера:
F = σА2/σОСТ2 (2.3)
Если F < Fкр, то различие между σА2 и σОСТ2 можно считать несущественным, т.е. влияние фактора А сравнимо с влиянием случайных помех. Если F > Fкр, то различие между σА2 и σОСТ2 существенно, т.е. фактор А оказывает влияние на выходную величину. Значение Fкр определяют по квантилям распределения Фишера, при уровне значимости α ("альфа") и степеням свободы f1 и f2: Fкр = f(α, f1, f2)
Воспользуйтесь поиском по сайту: ©2015 - 2024 megalektsii.ru Все авторские права принадлежат авторам лекционных материалов. Обратная связь с нами...
|