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

«Ростовский государственный строительный университет»

Кафедра «Высшей математики»

РАСЧЕТНАЯ РАБОТА

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

Вариант №6

Выполнил: ст. гр. МПГС-

ХХХХ

Руководитель: доцент кафедры

ХХХХ

Ростов-на-Дону

2014г.

Планирование эксперимента

По мере роста сложности исследуемых процессов и явлений возрастают затраты на аппаратуру и проведение эксперимента. B ходе испытаний собирается большое количество данных, требующих обработки и анализа. Широкое применение экспериментальных методов привело к созданию теории эксперимента, которая призвана дать ответы на следующие вопросы:

1. Как нужно организовать эксперимент, чтобы наилучшим образом решить поставленную задачу (в смысле затраты времени и средств или точности результата);

2. Как следует обрабатывать эксперимент, чтобы получить максимальное количество информации об исследуемом объекте или явлении;

3. Какие обоснованные выводы можно сделать об исследуемом объекте по результатам эксперимента.

Рассматривается общая модель для исследования качества технической системы. Ha такой объект действуют четыре группы факторов:

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

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

Ha практических занятиях рассматриваются модели полного факторного эксперимента (ПФЭ) в так называемых линейных моделей для случая одного

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

фактора, двух факторов, трех факторов, четырех факторов и так далее. Показывается, что число экспериментов должно быть N = 2*, где к - число

в различных пределах и иметь различную размерность, что усложняет обработку эксперимента и не обеспечивает нужную степень точности. Поэтому от размерных X_j переходят к безразмерным X_i с помощью соотношения используемых факторов. При числе факторов три и выше имеет место избыточность числа экспериментов по сравнению с числом определяемых коэффициентов многочлена «линейной аппроксимации», что приводит к необходимости использования метода наименьших квадратов (MHK) для их нахождения.

Для простоты рассмотрим случай трех факторов и ограничимся «линейной моделью»:

у = а + fr,x, + b₂x₂ + b₃x₃ + C₁X₁X₂ + C₂X₁X₃ + C₃X₂X₃ .

Предположим, что нами проведен ПФЭ с числом повторения m = 3 экспериментов в каждой его вершине (табл. 4):

Таблица 4

Составляем Excel – таблицу (табл.5)

Видно, что векторы X₁, X₂, X₃, X₁ • X₂, X₁ • X₃, X₂ • X₃ ортогональны.

Y-F

Поэтому все коэффициенты определяются соотношениями

Заполняем вначале заголовки таблицы - ячейки A1:Q1 ,A12 и

A17:A23. Затем вводим номера экспериментов в ячейках A3:A10 . Вводим исходные данные эксперимента в ячейки B3:D10 и H3:J10 . B ячейки E3 , F3 и G3 записываем соответственно формулы: =B3*C3 , =B3*D3 и =C3*D3 . Выделяем ячейки E3:G3 и производим протяжку формул до ячейки GlO . B ячейках КЗ - Q3 записываем соответственно формулы:

=(H3+I3+J3)/3, =K3*B3 , =K3*C3 , =K3*D3 , =K3*E3 , =K3*F3 , =K3*G3 . Выделяем диапазон K3:Q3 и производим протяжку формул до ячейки QlO . B ячейках C17 - C23 записываем формулы, описанные в ячейках A17 - A23 . B ячейке B12 находим сумму диапазона B3 - BlO и производим протяжку формулы до ячейки Q12 . B результате этих действий в ячейках C17 - C23 получаем коэффициенты интересующего нас многочлена:

у = 12,564 - 0,794x₁ + 0,792х₂ + 0,495х₃ + 1,725 х_хх₂ - 0,973*,¾ - 2,416х₂х₃

Отметим, что при записи аппроксимирующего многочлена мы в коэффициентах оставили три знака после запятой.

Целью получения аппроксимирующего многочлена является возможность нахождения оптимальных значений определяющих факторов. При большом числе факторов это является достаточно сложной задачей. При трех же факторах и, если нас устраивает невысокая точность их определения, нахождение minj) и maxj) можно определить методом простого перебора.

Программ на VBA , которая определяет оптимальные значения у,

представлена для нашего аппроксимирующего многочлена ниже.

В ячейках КЗ-Q3 записываем соответственно формулы: =(H3+I3+J3),

=K3*B3 , =K3*C3 , =K3*D3 , =K3*E3 , =K3*F3 , =K3*G3 . Выделяем диапазон K3:Q3 и производим протяжку формул до ячейки QlO . B ячейках C17 - C23 записываем формулы, описанные в ячейках A17 - A23 . B ячейке B12 находим сумму диапазона B3 - BlO и производим протяжку формулы до ячейки Q12 . B результате этих действий в ячейках C17 - C23 получаем коэффициенты интересующего нас многочлена:

у = 12,564 - 0,794x₁ + 0,792х₂ + 0,495х₃ + 1,725 х_хх₂ - 0,973*,¾ - 2,416х₂х₃

Отметим, что при записи аппроксимирующего многочлена мы в коэффициентах оставили три знака после запятой.

Целью получения аппроксимирующего многочлена является возможность нахождения оптимальных значений определяющих факторов. При большом числе факторов это является достаточно сложной задачей. При трех же факторах и, если нас устраивает невысокая точность их определения, нахождение minj) и maxj) можно определить методом простого перебора.

Программ на VBA , которая определяет оптимальные значения у,

представлена для нашего аппроксимирующего многочлена ниже:

Private Sub workbook_open()

Const h = 0.01

Const a = 17.875, b1 = 2.293, b2 = 1.762, b3 = -1.091

Const c1 = 1.656, c2 = -2.183, c3 = -1.214

Dim x1, x2, x3, y As Single

Dim x1min, x2min, x3min, ymin As Single

Dim x1max, x2max, x3max, ymax As Single

x1max = 10: x2max = 10: x3max = 10: ymax = -100000

x1min = 10: x2min = 10: x3min = 10: ymin = 100000

For x1 = -1 To 1 + h / 2 Step h

For x2 = -1 To 1 + h / 2 Step h

For x3 = -1 To 1 + h / 2 Step h

y = a + b1 * x1 + b2 * x2 + b3 * x3 + c1 * x1 * x2 + c2 * x1 * x3 + c3 * x2 * x3

If y < ymin Then

ymin = y

x1min = x1: x2min = x2: x3min = x3

End If

If y > ymax Then

ymax = y

x1max = x1: x2max = x2: x3max = x3

End If

Next x3

Next x2

Next x1

Range("A1").Value = "Ymin ="

Range("B1").Value = ymin

Range("A2").Value = "X1min ="

Range("B2").Value = x1min

Range("A3").Value = "X2min ="

Range("B3").Value = x2min

Range("A4").Value = "X3min ="

Range("B4").Value = x3min

Range("A6").Value = "Ymax ="

Range("B6").Value = ymax

Range("A7").Value = "X1max ="

Range("B7").Value = x1max

Range("A8").Value = "X2max ="

Range("B8").Value = x2max

Range("A9").Value = "X3max ="

Range("B9").Value = x3max

End Sub

Для нашего примера оптимальные значения будут (табл. 6):

Таблица 6