Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Мат.моделирование (методичка).doc
Скачиваний:
36
Добавлен:
14.11.2019
Размер:
4.09 Mб
Скачать

Лабораторная работа № 4.

Получение математической модели объекта по данным полного и дробного факторного эксперимента.

Цель работы: ознакомление с основами факторного эксперимента.

Задание:

  1. Ознакомится с теорией факторного эксперимента.

  2. Получить математическое описание зависимости выходного напряжения U2 трансформатора от входного напряжения U1 и сопротивления нагрузки R по результатам полного факторного эксперимента, проведенного на лабораторной установке «Однофазный трансформатор».

  3. Произвести статистический анализ полученной линейной модели. Сделать выводы об изучаемом явлении по модели.

  4. Получить математическое описание исследуемой зависимости U2=f(U1,R) по результатам планируемого расчета схемы замещения трансформатора. Применить схему ДФЭ

Основные понятия теории факторного эксперимента.

Математическая теория планирования эксперимента предлагает различные методы получения математического описание процесса в изученной системе [1,2]. Основой этих методов является возможность представления исследуемой (целевой) функции системы, зависящей от п независимых переменных параметров (факторов) системы x1,…xi….xn рядом Тейлора (1)

Ряд Тейлора, имеющий постоянные коэффициенты αо, αi, αik, αii, описывает исследуемую функцию в области базисной (исходной) точки xi* (i = 1,2,…n). Эта область ограничена предельными приращениями (шагами) переменных ± Δxi, которым соответствуют безразмерные нумерованные переменные равные ± 1. Нормированные переменные определяются по формуле (2)

(2)

В формуле (2) хi – текущее значение фактора

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

Определение коэффициентов ряда Тейлора по опытным данным у планируемого эксперимента (таблица 1) и определение числа членов ряда, достаточного для описания системы с заданной точностью, представляет задачу, решаемую статистическим методом факторного эксперимента.

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

Рис. 1

Полагают, что помеха е является случайной величиной, имеющей нормальный закон распределения с параметрами m{e}=0, σ2{e}=σ2{y}. В силу влияния помехи по экспериментальным данным у получают не истинное значение коэффициентов, а их статистические оценки. При этом математическая модель системы имеет вид

(3)

где ур – статистическая оценка предсказываемого значения функции у;

а0, аi, аik, аii – статистическая оценка коэффициентов α0, αi, αik, αii.

Для получения более точной математической модели системы эксперимент должен быть поставлен так, чтобы коэффициенты уравнения (3) являлись наилучшими оценками коэффициентов ряда Тейлора.

Наилучшие оценки коэффициентов ряда Тейлора дает метод наименьших квадратов [3], который, однако, не отвечает на вопрос об оптимальном объеме эксперимента. Теория планирования эксперимента, основанная на методе наименьших квадратов и применяющая методы дисперсионного анализа [4], предлагает такие программы (планы) факторных экспериментов, реализация которых дает наилучшие оценки коэффициентов ряда Тейлора при минимальном числе опытов.

Факторный эксперимент проводят на два этапа. На первом этап выдвигают гипотезу о виде уравнения (линейное, квадратичное, кубичное и д.т.), описывающем систему с определенной точностью. В соответствии с выдвинутой гипотезой о виде математической модели выбирают план эксперимента. По результатам эксперимента получают коэффициенты модели. Так как естественно желание получить наиболее простую модель, то сначала выдвигают гипотезу о линейной модели (4).

(4)

На втором этапе проводят статистический анализ значений коэффициентов a0, a1,…an с целью определения значимости влияния варьируемых факторов на функцию у. Затем приводят гипотезу об адекватности описания функции у принятым уравнением. В случае не состоятельности гипотезы о линейной аппроксимации функции выдвигают новую гипотезу о более сложной форме зависимости, например, квадратичной:

Для определения коэффициентов данного полинома реализуют факторный план второго порядка. Особенности планирования второго порядка излагаются в [1].

Получение математического описания процесса системы по данным эксперимента. Постановка «задачи регрессии»

Получение зависимости функции у от факторов системы в виде уравнения (регрессии) по результатам эксперимента представляет задачу, известную в математике как «задача регрессии» [5,6].

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

При выдвинутой гипотезе о линейной зависимости функции у от п факторов системы (4)

«задача регрессии» состоит в отыскании п + 1 неизвестных коэффициентов а0, а1,…аi,… аn по данным N измерений величины у. Это требует совместного решения системы N уравнений (5): (5)

где хij – значение i – фактора в j – опыте, x0j=1,

уj- результат измерения у в j – опыте.

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

ТАБЛИЦА 1

Номер опыта

Фактор

Исследуемая функция у

X1

X2

Xi

Xn

1

X11

X21

-

Xi1

-

Xn1

y1

-

-

-

-

-

-

-

-

J

X1j

X2j

-

Xij

-

Xnj

уj

-

-

-

-

-

-

-

-

N

X1N

X2N

-

XiN

-

XnN

yN

При числе опытов N, превосходящих число искомых коэффициентов п + 1, и невозможность рассматривать избыточные уравнения как противоречащие остальным в силу случайного характера измеряемой величины у существует многозначность в определении коэффициентов при решении системы (5).

Разрешение этой трудности дает применение принципа максимума правдоподобия, который может быть сформулирован следующим образом [7]:

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

Требование принципа максимума правдоподобии при решении «задачи регрессии» состоит в определении таких значений коэффициентов а0, а1,…аi,… аn, при которых вероятность получить наблюдаемые величины функции максимальна.

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

В следующем параграфе в рамках решения «задачи регрессии» излагается метод наименьших квадратов.

Применение метода наименьших квадратов для получения коэффициентов математической модели системы «задачи регрессии»

Для получения наилучших оценок коэффициентов математической модели системы в виде уравнения (3) применяется метод наименьших квадратов. Дальнейшее рассмотрение метода наименьших квадратов применительно к «задаче регрессии», а также статистические методы оценки его результатов ведутся в предположении:

1. Нормального закона распределения измеряемой величины у (2) при любой комбинации фиксированных уровней факторов xi (при номинальном распределении помехи е величина у имеет также нормальное распределение).

2. Пренебрежимо малой ошибки в установке уровней факторов xi по сравнению с ошибкой измерения величины у.

3. Равноточности измерения величины у во всех N опытах, то есть постоянство величины дисперсии σ2{y}=σ2{e} при любой комбинации уровней факторов x0, x1,…xi,… xn.

Требование метода наименьших квадратов при определении коэффициентов уравнения (3) по данным эксперимента состоит в минимизации суммы квадратов ошибок εj, обусловленных неточностью представления экспериментальных значений yj уравнением (3), дающим значения yp

,

где

(6)

Очевидно , что сумма

(7)

Экспериментальное значение S имеет место при выполнении следующих условий:

(8)

Систему уравнений (8) можно переписать в виде системы «нормальных» уравнений (9): (9)

При введении условных обозначений

,

система «нормальных» уравнений запишется сокращенно: ( (10)

Решение системы «нормальных» уравнений упрощается, если ввести матричную форму записи. Тогда система (5) запишется

(11)

Х – матрица факторов порядка

, (12)

А – матрица коэффициентов порядка

, (13) , (14)

Y – матрица , представляющая экспериментальные данные исследуемой функции у (14).

При этом система «нормальных» уравнений (10) примет вид:

, (15)

где ХТ – транспонированная матрица Х. , (16)

, (17)

При умножении обеих частей системы уравнений (15) на матрицу (ХТХ-1), обратную матрице ХТХ, получают формулу искомых коэффициентов:

(18)

При этом используют известное из матричной алгебры соотношение

Если обозначить элементы матрицы (ХТХ)-1 так, как показано в формуле (19)

(19)

и подставить (17) и (19) в формулу (18), то искомые коэффициенты найдутся из выражения (20) . (20)

Из полученной формулы (20) следует, что коэффициенты уравнения равны линейным комбинациям экспериментальных данных функции у, а именно

(21)

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

m{ai}=αiматаматическое ожидание коэффициента аi;

σ2{ai} – дисперсия коэффициента аi.

В [6] показано. что коэффициенты (21), полученные методом наименьших квадратов, обладают свойствами наилучших статистических оценок. А именно:

1) коэффициенты аi являются не смещенными оценками коэффициентов αi

m{ai}=αi;

2) коэффициенты ai имеют минимальную дисперсию σ2{ai} среди класса оценок, равных линейным комбинациям экспериментальных данных.

Из формулы (21) следует однако и то, что значения коэффициентов ri зависит от числа членов уравнения (3), которое гипотетически принимается равным п + 1 перед проведением эксперимента (п – число существенно влияющих факторов на функцию у). При отказе от этой гипотезы необходим пересчет коэффициентов, что связано с дополнительными затратами труда и времени. Перерасчет коэффициентов устанавливается, если реализовать ортогональный план эксперимента [1].

Ортогональный план эксперимента.

План эксперимента называется ортогональным, если ортогональны между собой столбцы плана X1, X2,… Xn (табл.1), т.е.

при .

Тогда диагональные матрицы ХТХ и (ХТХ)-1

(23)

(24)

При подстановке формулы (24) в формулу (20) получают

(25)

откуда видно, что коэффициент аi(i=01,…n), равный

, (26)

не зависит от числа членов уравнения (3)

Дисперсия коэффициента σ2{ai} при этом находится согласно [6] из выражения (27)

где J – матрица из единиц порядка п×п [6].

Из (27) следует, что дисперсия любого коэффициента линейного уравнения (3) равна

(28)

то есть зависит от дисперсии опыта.

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

Полный факторный эксперимент.

В [1,2] показано, что ортогональный план эксперимента является оптимальным в широком смысле при выполнении следующих условий:

  1. Симметричное расположение варьируемых уровней факторов относительно центрального (среднего ) уровня , i=1,2,…n

  2. Равенство элементов диагональной матрицы

, i=1,2,…n

Реализация оптимального плана позволяет получить значения коэффициентов регрессии с минимальной дисперсией σ2{ai} при минимальном числе опытов.

Ортогональным планом является двухуровневый полный факторный эксперимент (ПФЭ) с числом опытов N=2п. Он содержит всевозможные неповторяющиеся комбинации уровней всех переменных факторов. В таблице 2 дана программа ПФЭ 22 с указанием нормированных значений факторов Xi max=+1, Ximin=-1, соответствующих максимальному ximax и минимальному ximin уровням факторов (рис.2.). ПФЭ 22 дан на рис. 3 в плоскости [x1, x2].

Рис. 2

Связь между величинами Хi и xi устанавливается формулой

,

Рис. 3.

где хi – значение фактора в именованных единицах;

хi* – базисное значение фактора;

Δхi – шаг изменения значения фактора

Из этой формулы следует, что при хi= хi Δхi нормированный уровень равен

Xi=±1.

По данным ПФЭ 22 получают коэффициенты линейного уравнения

(30)

Определяемые по формулам (26, 29):

(31)

Если матрицу X данного плана дополнить столбцом парного взаимодействия X1 X2

, (32)

то можно получить формулу коэффициента

В этом случае модель системы представляется неполным квадратичным уравнением

В таблице 3 дан ПФЭ 23

ТАБЛИЦА 2

j

X1

X2

Y

1

+1

+1

y1

2

-1

+1

y2

3

+1

-1

y3

4

-1

-1

y4

ТАБЛИЦА 3

j

X1

X2

X3

y

1

+1

+1

+1

y1

2

-1

+1

+1

y2

3

+1

-1

+1

y3

4

-1

-1

+1

y4

5

+1

+1

-1

y5

6

-1

+1

-1

y6

7

+1

-1

-1

y7

8

-1

-1

-1

y8

С помощью ПФЭ 23 получают математическую модель системы в виде уравнения

(34)

Коэффициенты модели определяют по формулам (26, 29):

Для составления планов двухуровневого полного факторного эксперимента не нужен каталог. План для п факторов получают в виде таблицы на п столбцов и N=2n строк. Частота изменения уровней +1 и –1 в каждом столбце уменьшается в два раза по сравнению с предыдущим столбцом. Несмотря на то, что планы ПФЭ 2п являются оптимальными, объем их слишком велик при достаточно большом числе факторов

(п = 10). Значительную экономию в числе опытов дают планы типа дробного факторного эксперимента.

Дробный факторный эксперимент

План, составляющий 1/2т часть объема полного факторного эксперимента 2п, называется дробным факторным экспериментом 2п-т. При т=1, 2, 3… объем ДФЭ составляет 1/2, 1/4, 1/8, … часть объема ПЭФ. Примечание ДФЭ возможно в случаях, когда имеется предварительная информация о слабых т зависимостях между факторами. Тогда столбцы взаимосвязей используются для задания изменений уровней других т факторов.

Например, полный двухфакторный эксперимент с числом опытов N=22=4 (табл. 5) можно применить для трехфазной задачи (n=3), если известно априори, что слаба взаимосвязь X1 X2. Программа изменения третьего фактора X3 задается столбцом X1 X2 , т.е.

Х31 Х2

При применении ПФЭ в случае трех факторов приводится N=23=8 опытов. Таким образом, ДФЭ сокращает в этом случае объем эксперимента в 2 раза.

ТАБЛИЦА 5.

Опыт

j

X1

X2

X3=X1 X2

y

1

+1

+1

+1

y1

2

-1

+1

-1

y2

3

+1

-1

-1

y3

4

-1

-1

+1

y4

ТАБЛИЦА 6.

Опыт

j

X1

X2

X3

X4= X1X2

X5= X1X3

X6= X3X2

y

1

+1

+1

+1

+1

+1

+1

y1

2

-1

+1

+1

-1

-1

+1

y2

3

+1

-1

+1

-1

+1

-1

y3

4

-1

-1

+1

+1

-1

-1

y4

5

+1

+1

-1

+1

-1

-1

y5

6

-1

+1

-1

-1

+1

-1

y6

7

+1

-1

-1

-1

-1

+1

y7

8

-1

-1

-1

+1

+1

+1

y8

Таблица ПФЭ 23 применяется для получения линейной модели для семи факторов, если слабы взаимосвязи X1X2, X1X3X2X3, X1X2X3. Столбцы взаимосвязи в этом случае используются как программы изменения факторов

X4=X1X2, X5 =X1X3, X6 =X2 X3, X7=X1X2X3

Полученный таким образов план для семи факторов является ДФЭ 27-4. Объем плана составляет часть ПФЭ 27. Подробно рассматриваются данные планы в методических указаниях [10].

Статистический анализ результатов полного и дробного факторного экспериментов.

Статистический анализ результатов факторного эксперимента состоит

  1. В проверке значимости коэффициентов модели.

  2. В проверке адекватности уравнения f(X1, Xn)

Проверка коэффициента ai на значимость состоит в проверке гипотезы о равенстве нулю истинного значения коэффициента αi. Принимая эту гипотезу, находят доверительный интервал для коэффициента αi, т.е. αi=0. Принимая эту гипотезу, находят доверительный интервал для коэффициента ε{αi},внутри которого в силу влияния помехи, находится истинное значение αi с вероятностью Р. обычно принимают Р = 0,95 – 0,98. Доверительный интервал для коэффициента ai, имеющего нормальное распределение, равен

, (36)

где σ{ai} - оценка среднеквадратичного отклонения величины ai, определяемая по ошибке эксперимента σ2{y}

, (37)

tT – значение определяемое по таблице t-распределения Стьюдента [5] в зависимости от степени свободы и вероятности Р. Здесь γ – число повторений N опытов в плане эксперимента.

Ошибка эксперимента σ2{y} находится по значениям у в результате повторения всех N опытов плана γ раз (табл. 7). Обычно берут γ=2-4 [1,2]. Тогда

(38)

В формуле (38)

yiq – значение у в j-м опыте q-й серии (в таблице 7 уj1 и уj2)

ˉyj – среднее значение у в j-м опыте (по данным γ серии опытов).

(39)

При проверке значимости коэффициентов, значение коэффициента ai сравнивают с его доверительным интервалом ε{ai}. Если модуль коэффициента |ai| превышает доверительный интервал ε{ai}, т.е.

|ai|> ε{ai} (40)

ТАБЛИЦА 7

Опыт

j

Х1

Х2

Опытные данные двух серий (γ=2)

Среднее значение

Расчетное значение

yj1

yj2

ypj

1

2

3

4

+1

-1

+1

-1

+1

+1

-1

-1

y11

y21

y31

y41

y12

y22

y23

y24

yp1

yp2

yp3

yp4

то гипотеза о равенстве нулю коэффициента αi отвергается. Коэффициент ai признается значимым, и слагаемое aixi остается в уравнении системы в силу существенного влияния фактора xi на функцию y.

Если выполняется неравенство

|ai|< ε{ai} (41)

то это означает, что коэффициент ai незначим. Слагаемое aixi исключается из уравнения.

После проверки коэффициентов модели на значимость проверяется гипотеза об адекватности представления в области исследования функции y(x1, xn) принятым полиномом.

Проверка этой гипотезы основана на сравнении ошибки от неадекватности уравнения S2ад{y} с ошибкой эксперимента S2{y}. Для чего находится оценка дисперсии неадекватности

(42)

где п3 – число значимых коэффициентов уравнения,

уpj – значение функции у, полученное по уравнению в j-м опыте,

- экспериментальное значение функции у в j-м опыте.

При проверке адекватности уравнения применяется критерий Фишера [5]. Для этого находят величину F, имеющую F распределений Фишера

(43)

Величину F сравнивают с табличным значением FT, определяемым при степенях свободы ν1=N-n3 и ν2=N-(γ-1) при вероятности Р = 0,95 – 0,98.

В случае

F ≤ FT

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

F >FT

уравнение неадекватно. Тогда отказываются от выдвинутой гипотезы о моделировании системы полученным уравнением [10].

Порядок выполнения работы.

1. Получить зависимость выходного напряжения U2 трансформатора от входного напряжения U1 и сопротивления нагрузки R в виде линейного уравнения

где функция у соответствует U2, фактором X1=U; X2=R.

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

U1* ± ∆U1 , R* ±∆R

здесь U1*, R* - базисные значения факторов (согласовать с преподавателем);

U1, ∆R - шаги изменения значений факторов; принять что ∆U1 =20% U1* ∆R=20% R*

2.Составить таблицу 1, где представлен полный факторный эксперимент 22 (столбцы 2,3)

ТАБЛИЦА 1.

Опыт

j

Фактор

Исследуемая функция у

Среднее значение

Значение уpj , полученное по уравнению (1)

Х1

Х2

yj1

yj2

1

2

3

4

+1

-1

+1

-1

+1

+1

-1

-1

y11

y21

y31

y41

y12

y22

y23

y24

yp1

yp2

yp3

yp4

Нормированным уравнениям факторов, равным +1 и – 1, соответствуют следующие значения факторов в именованных единицах

+1= U1*+∆U1 ; R*+ ∆R

-1= U1* - ∆U1 ; R* - ∆R

3.Проведение эксперимента:

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

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

в) после первой серии измерения провести вторую серию, результаты записать в пятом столбце таблицы 1.

г) найти среднее значение выходного напряжения в каждом опыте (столбец 6 таблицы 1).

4.По полученным опытным данным рассчитать коэффициенты а0, а1, а2 линейного уравнения (1)

(1)

Формулы для коэффициентов таковы: (2)

5.Проверить зависимость коэффициентов ai. Для этого следует найти оценку дисперсии исследуемой величины S2( ) по результатам повторных опытов:

(3)

Тогда доверительный интервал для коэффициентов ai равен:

, (4)

где tT –величина, взятая по таблице t-распределения (распределение Стьюдента) [1] при уровне значимости α=0,05 и степени свободы ν=N(γ-1) .

Если абсолютная величина коэффициента ai превышает доверительный интервал ε(ai)

|ai|> ε{ai}

то коэффициент ai значим и слагаемое aixi остается в уравнении.

Если

|ai|< ε{ai}

То коэффициент ai незначим и слагаемое aixi исключается из уавнения.

6. Проверить адекватность уравнения (1). Для этого нужно найти величину S2ад{y}

(7)

где п3 – число значимых коэффициентов уравнения,

уpj – значение функции у, полученное по уравнению в j-м опыте,

- экспериментальное значение функции у в j-м опыте.

После чего найти F-отношение

(8)

Сравнить F-отклонение с табличным значением FT, полученным по таблице F-распределения (распределение Фишера) при α=0,05 ν1=N-n3 и ν2=N(γ-1)γ

Если F<FT , то увеличение адекватно и может быть использовано для анализа изучаемой зависимости U2=f(U1,R).

7. По полученной математической модели сделать выводы:

а) какова степень влияния факторов U1 и R на выходное напряжение U2 ;

б) как следует изменить значения U1* ,R*, чтобы влияние в области с большим значением U2

Указание к п. 7 б): перемещение к оптимальной области с U2max осуществить в направлении градиента функции U2=f(U1,R). Для этого использовать полученные коэффициенты a1 и a2, являющиеся оценками составляющих градиента функции, а именно:

,

Новые базисные уровни факторов U1 и R найти в этом случае по формулам:

U1**=∆U1* +n∆U1; (10)

R**= R* +m∆R (10)

Где ∆U1 и ∆R - шаги изменения факторов;

n и m - числа шагов, которые нужно сделать, чтобы движение точки [U1*, R*] происходило в направлении градиента функции U2.

Числа шагов n и m находят так: задают произвольно, например число шагов n (3…5) и получают m из соотношения

(11)

После определения новых уровней факторов U1**, R** провести ПЭФ 22 с центром в этой точке. Если линейная модель в окрестности точки [U1**, R**] неадекватна, то реализуют в этой области эксперимент по плану второго порядка для получения модели в виде квадратичного уравнения: (12)

Для этого проводят опыты по схеме центрального или ротатабельного композиционного планирования эксперимента [U1**, R**]. если линейная модель в области с центром [U1*, R*] адекватна, то продолжают поиск оптимальной области.

8. Получить математическое описание зависимости из расчета схемы замещения трансформатора по плану дробного факторного эксперимента U2=f(U1,R). Переменными величинами считать:

напряжение на входе U1=U120%U1* (фактор х1)

сопротивление нагрузки R=R*±20%R* (фактор х2)

суммарное индуктивное сопротивление рассеяния храсрас-*±10%храс* (фактор х31х2), где храсрас1+х'рас2.

Получить линейное уравнение

(13)

После этого проверить значимость коэффициентов и адекватности уравнения. Получить формулы для коэффициентов ai, которые смешаны с коэффициентами взаимодействий a13, a24, a123.

Если уравнение (9) адекватно, то это означает, что коэффициенты a13, a24, a123 несущественны.

Контрольные вопросы.

  1. Какова идея наименьших квадратов?

  2. Какой план эксперимента является ортогональным и каковы его достоинства?

  3. Покажите таблицу полного факторного эксперимента 23. Какое уравнение получают по данным ПФЭ 23?

  4. Что представляет собой дробный факторный эксперимент? Покажите таблицу ДФЭ , где х41х2.

  5. Как осуществляются смешанные оценки коэффициентов?

  6. В чем состоит оценка значимости коэффициентов и адекватности матаматической модели?

Литература:

  1. Круг Г.К. Статистические методы в инженерных исследованиях: Лабораторный практикум. В.П. Бородюк и др. – М.: В.Ш., 1983.

  2. Веников В.А. Теория подобия и моделирования. – М.:В.Ш., 1971.

Артемева Е.М., Круг Н.К. Факторный эксперимент: Учебное посоьие.- И., ИМИ, 1982.