- •ВВЕДЕНИЕ
- •1. ЭКСПЕРИМЕНТ КАК ПРЕДМЕТ ИССЛЕДОВАНИЯ
- •1.1. Понятие эксперимента
- •1.2. Виды экспериментальных исследований
- •1.3. Математическая модель объекта исследования
- •1.4. Факторы
- •1.5. Параметр оптимизации
- •Вопросы и задания для самоподготовки
- •2. ПРОСТЫЕ СРАВНИВАЮЩИЕ ЭКСПЕРИМЕНТЫ
- •2.1. Предварительная обработка экспериментальных данных
- •2.2. Статистические гипотезы
- •2.5. Сравнение двух рядов наблюдений
- •2.6. Сравнение двух дисперсий
- •2.7. Проверка гипотезы о законе распределения
- •2.7.1. Общие сведения
- •2.7.2. Пример проверки гипотезы о нормальном законе распределения экспериментальных данных
- •Вопросы и задания для самоподготовки
- •3. АНАЛИЗ РЕЗУЛЬТАТОВ ПАССИВНОГО ЭКСПЕРИМЕНТА
- •3.1. Элементы дисперсионного анализа
- •3.1.1. Общие сведения
- •3.1.2. Пример применения однофакторного дисперсионного анализа
- •3.2. Характеристика видов связей между рядами наблюдений
- •3.4. Определение тесноты связи между случайными величинами
- •3.5. Парная линейная корреляция
- •3.6. Статистическое изучение корреляционной связи
- •3.6.2. Исключение из массива первичной информации промахов
- •3.6.4. Измерение степени тесноты связи, оценка ее существенности
- •3.8. Линейная множественная регрессия
- •Вопросы и задания для самоподготовки
- •4. МНОГОФАКТОРНЫЕ ЭКСПЕРИМЕНТЫ
- •4.1. Общие сведения
- •4.2. Полный факторный эксперимент
- •4.2.1. Кодирование факторов
- •4.2.2. Матрицы планирования эксперимента
- •4.2.3. Рандомизация опытов
- •4.2.4. Проведение эксперимента
- •4.2.5. Проверка однородности дисперсии параллельных опытов, воспроизводимости эксперимента
- •4.2.6. Расчет коэффициентов регрессии, проверка их значимости
- •4.2.7. Проверка адекватности модели
- •4.2.8. Пример применения планов первого порядка полного факторного эксперимента
- •4.3. Дробный факторный эксперимент
- •4.4. Применение плана первого порядка дробного факторного эксперимента для исследования технических устройств
- •4.4.1. Теоретическая модель усилителя
- •4.5. Планы второго порядка
- •4.5.1. Ортогональные планы
- •4.5.2. Ротатабельные планы второго порядка
- •5.1. Метод покоординатной оптимизации
- •5.3. Симплексный метод планирования
- •ПРИЛОЖЕНИЯ
4.4. Применение плана первого порядка дробного факторного эксперимента для исследования технических устройств
Исследование технических объектов обычно связано с их математическим описанием, которое позволяет выбрать структуру и параметры устройства на стадии разработки, а также организовать наиболее эффективные режимы функционирования в процессе эксплуатации.
Для ряда устройств построению эмпирической модели предшествует получение теоретической модели. Сравнение результатов эксперимента с теоретической зависимостью позволяет объяснить смысл изучаемого явления и показать преимущества эксперимента.
В качестве примера рассмотрим усилитель с положительной и отрицательной обратными связями. Схема усилителя собрана на ин-
Д |
|
|
|
|
тегральной микросхеме серии К140УД6 (рис. 4.3). Входной сигнал |
||||
UВХ подается на инвертирующий вход. Отрицательная обратная связь |
||||
создается с помощью сопротивлений R и R . Сопротивления R , |
R |
4 |
||
1 |
И2 |
3 |
|
образуют положительную обратную связь, которая при определенных
KU UВЫХ UВХ . |
|
(4.36)А |
|
|
|
|
R1 1 |
|
R2 |
|||||||||
Задача исследован |
я состоит в |
UВХ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
получении и анал зе коэфф ц ента |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
б |
|
|
|
|
I1 |
|
|
|
|
|
|
|
|
|
|||
усиления KU как функц |
сопро- |
|
|
|
|
|
|
|
|
2 |
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
тивлений отрицательнойии положи- |
|
|
|
|
|
|
|
|
|
|
|
|
UВЫХ |
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
тельной обратных связей. Необхо- |
|
|
|
|
R3 |
|
|
R4 |
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
||||||||
димо также сравнить теоретическую |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
С |
|
|
|
|
|
3 |
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
I3 |
|
|
|
|
|
|
|
|||||
и эмпирическую зависимости KU от |
|
|
|
|
|
|
|
|
|
|
|
|||||||
сопротивлений R1, R2, R3 , R4. |
|
|
|
|
Рис. 4.3. Усилитель |
|||||||||||||
При исследовании |
подобного |
с отрицательной и положительной |
||||||||||||||||
устройства следует учитывать, что |
|
|
|
|
обратными связями |
соотношениях между сопротивлениями повышает коэффициент усиления по напряжению схемы
в зависимости от параметров отрицательной и положительной обратных связей коэффициент усиления
KU изменяется в широких пределах, а при некоторых значениях сопротивлений схема может оказаться неработоспособной.
Для решения поставленной задачи исследования получим теоретическую и эмпирическую модели усилителя, а затем сравним их.
77
4.4.1. Теоретическая модель усилителя
Для схемы, собранной на операционном усилителе, при выводе зависимости KU f R1,R2,R3,R4 принимаются следующие допуще-
ния: коэффициент усиления в разомкнутом состоянии и его входное сопротивление бесконечно велики. Следовательно, потенциалы в точках 1, 3 (см. рис. 4.3) одинаковы:
|
|
|
|
1 |
3, |
|
|
|
|
|
|
|
|
(4.37) |
|||
а входной ток усилителя равен нулю. |
|
|
|
|
|
|
|
|
|
||||||||
На основании закона Ома можем записать уравнения |
|
||||||||||||||||
|
|
|
|
1 |
UВХ I1R1; |
|
|
|
|
|
|
(4.38) |
|||||
|
|
|
|
3 I3R3. |
|
|
|
|
|
|
|
|
(4.39) |
||||
Токи I1, I3 также определяются на основе закона Ома: |
|
||||||||||||||||
|
|
|
|
I1 UВХ UВЫХ |
R1 R2 ; |
(4.40) |
|||||||||||
|
|
|
|
I3 UВЫХ |
R3 R4 . |
|
|
|
|
(4.41) |
|||||||
Подставляя в уравнения |
(4.38), (4.39) выражения для токов |
||||||||||||||||
|
|
|
|
|
|
А |
Иполучаем уравнение, свя- |
||||||||||
(4.40), (4.41) с учетом соотношения (4.37), |
|||||||||||||||||
зывающее выходное и входное напряжения: |
|
|
|
|
|
||||||||||||
|
|
|
R1 |
|
|
|
|
ДR3 |
|
|
|
||||||
UВЫХ |
|
|
|
UВХ U |
ВЫХ |
|
|
UВЫХ . |
(4.42) |
||||||||
R R |
|
R R |
|||||||||||||||
|
|
1 |
|
2 |
|
|
|
|
|
|
|
3 |
|
4 |
|
|
|
и |
|
|
|
|
|
|
|
|
|
|
|
|
|||||
После преобразования соотношения (4.42) получим формулу |
|||||||||||||||||
для коэффициента ус лен я |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
С |
|
бU |
1 R |
R |
|
|
|
|
|
||||||||
KU |
ВЫХ |
|
|
|
3 |
4 |
|
. |
|
|
(4.43) |
||||||
|
|
|
R |
|
R R |
R |
|
|
|||||||||
U |
ВХ |
|
|
|
|||||||||||||
|
|
|
|
|
3 |
4 |
1 |
2 |
|
|
|
|
Выражение (4.43) является общим, так как из него могут быть получены известные соотношения коэффициентов усиления по напряжению для инвертирующего и неинвертирующего усилителей. При R4 (разомкнута положительная обратная связь) находим коэффициент усиления для инвертирующего усилителя:
K R |
R ; |
(4.44) |
|
U |
2 |
1 |
|
при R2 (разомкнута отрицательная обратная связь) получаем ко- |
|||
эффициент усиления для неинвертирующего усилителя: |
|
||
K 1 R |
R . |
(4.45) |
|
U |
4 |
3 |
|
Усилитель устойчиво работает при отрицательных значениях KU , т.е. условие устойчивой работы определяется неравенством
R1 R2 R3 R4 . |
(4.46) |
78
KU
R3/R4<R1/R2
1
R3/R4
0
R3/R4=R1/R2
-R2/R1
R3/R4>R1/R2
ний R3 R4 при R1 R2 |
const. |
Рис.4.4. Зависимость коэффициента |
|
Коэффициент |
усиления |
усиления по напряжению |
|
от параметров схемы |
|||
при |
|
||
|
|
||
R1 R2 R3 R4 |
(4.47) |
|
принимает неопределенное значение , т.е. имеет место разрыв |
||||
функции. В области, где |
|
|
И |
R4 |
KU 0, с уменьшением отношения R3 |
||||
абсолютная величина KU |
уменьшается и при R3 R4 0 модуль дос- |
|||
|
|
Д |
|
|
|
А |
|
|
|
б |
|
|
|
тигает наименьшего значения, определяемого выражением (4.44). При
R3R4 R1R2 усилительинера отоспособен. Поэтому область исследования задается неравенством (4.46).
4.4.2. Эмп р ческая модель усилителя
В качестве выходной величины (отклика) примем коэффициент
усиления по напряжению, в качестве факторов – сопротивления отри- |
|
цательной и положительнойС |
обратных связей R1, R2, R3 , R4. Из дру- |
гих возможных факторов на коэффициент усиления KU может оказывать влияние напряжение питания E 15 В. Предварительный однофакторный эксперимент показал, что при отклонениях напряжения питания на 10% от номинального значения коэффициент усиления практически не изменяется.
Целью эксперимента является построение экспериментальной зависимости KU f R1,R2,R3,R4 . Получение такой зависимости диктуется тем, что при построении теоретической зависимости (4.43)
79
принималось допущение об идеальности операционного усилителя К140УД6, на котором собран усилитель напряжения.
По результатам экспериментального исследования можно оценить, насколько реальный коэффициент усиления KU отличается от теоретического и в каких пределах он будет изменяться при отклонении факторов R1, R2, R3 , R4 от их номинальных значений.
В процессе эксперимента коэффициент усиления определяется как отношение выходного и входного напряжений по формуле (4.37). Оба напряжения измеряются комбинированным цифровым прибором Щ4310. Согласно паспортным данным основная относительная погрешность 0 измерения постоянного напряжения определяется в процентах выражением
обусловленная отклонениями температурыИи напряжения питания от
|
|
|
|
|
|
0 |
|
|
UN |
|
(4.48) |
|
|
||||
0,5 0,5 |
1 , |
||||
|
|
UX |
|
|
|
где UN – предел измерения; UX |
|
Д |
|
||
– измеренное значение. |
|
||||
Кроме основной, имеет место дополнительная погрешность, |
|||||
|
А |
|
|
номинальных значений. Дополнительная погрешность используемого
|
б |
|
прибора согласно паспортным данным может составлять |
|
|
|
Д 1,5 0 . |
(4.49) |
и |
|
|
Таким образом, о щая погрешность средства измерения |
|
|
С |
0 Д . |
(4.50) |
Входное сопрот влен е RВХ прибора при измеряемых напряжениях 2 мВ – 20 В RВХ 100 МОм. Входное напряжение постоянного тока UВХ 13 мВ снимается с источника ВСП-50. Величины варьируемых сопротивлений устанавливаются с помощью измерительных магазинов сопротивления Р33 с погрешностью (%)
|
6 |
|
R |
N |
|
|
|
R 0,2 6 10 |
|
|
|
, |
(4.51) |
||
|
|
|
|||||
|
|
RX |
1 |
||||
|
|
|
|
|
|
где RN – наибольшее значение сопротивления магазина (100 кОм); RX – номинальное значение включенного сопротивления.
Выбор интервалов варьирования и кодирование факторов. Ин-
тервалы варьирования факторов определяются в предварительном эксперименте из условия существенного (на 30 – 40 %) изменения отклика – коэффициента усиления. Заданное изменение отклика имеет
80
место при изменении сопротивлений примерно на 10 % от номинального значения. В качестве номиналов сопротивлений приняты их значения, соответствующие KU 92,5. Номинальные значения (основные уровни) факторов и интервалы их варьирования представлены в табл. 4.17.
Таблица 4.17
Условия эксперимента
Уровень варьируемых |
Кодовое |
R1, кОм |
R2, кОм |
R3, кОм |
R4, кОм |
||
факторов |
обозначение |
X1 |
X2 |
X3 |
X4 |
||
Нижний уровень |
–1 |
9,5 |
470 |
9 |
1000 |
||
Верхний уровень |
+1 |
10,5 |
530 |
11 |
1200 |
||
Основной уровень |
0 |
|
10 |
500 |
10 |
1100 |
|
Интервал варьирова- |
X |
u |
0,5 |
30 |
1 |
100 |
|
ния |
|||||||
|
|
|
|
|
План и результаты эксперимента. В предварительном экспе-
рименте выяснено, что взаимодействие факторов X1X3X4 , т.е. коэф-
фициент b |
, является статистически незначимым. Поэтому в качест- |
134 |
И |
ве плана эксперимента можно принять план ФЭ 24 1 с генерирую-
щим соотношением X |
2 |
X |
1 |
X |
3 |
X |
4 |
. При этом определяющий кон- |
|
|
|
|
|
|
Д |
||||
траст имеет вид 1 X1X2X3X4 . |
Такой контраст обеспечивает сме- |
||||||||
шивание линейных эффектов сАэффектами тройных взаимодействий и |
парных взаимодейств й между со ой. Так как тройные взаимодейст-
б |
4 1 |
вия обычно стат ст чески незначимы, то принимаемый ДФЭ 2 |
|
и |
|
следует считать допуст мым. План и результаты измерения отклика
Ku y UВЫХ |
UВХ представлены в табл. 4.18. |
|
Обработка результатов ДФЭ 24 1. Результаты эксперимента, |
||
т.е. значения коэффициентаС |
усиления, получены делением выходного |
|
напряжения |
(показания |
вольтметра) на входное напряжение |
UВХ 13мВ. Повторные измерения дают тот же результат, поэтому в каждой строке плана проводится одно измерение отклика. Дисперсию воспроизводимости в таком случае следует определять по метрологическим характеристикам средства измерения. Основную и дополнительную погрешности измерения вольтметра Щ4310 определяем по формулам (4.48) – (4.50). В эксперименте выходное напряжение измерялось на пределе UN 2 В и составляло UВЫХ 1 1,5 В. Следова-
81
тельно, основная относительная погрешность прибора 0 1%, дополнительная Д 1,5 0 1,5%, общая погрешность 2,5%.
Таблица 4.18
План и результаты ДФЭ 24 1
Номер |
X0 |
X1 |
|
X2 |
X3 |
|
X4 |
|
X1X3 |
|
X2X3 |
X3X4 |
|
y |
~ |
|
~ |
|||
опыта |
|
|
|
|
|
|
y |
|
y y |
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
1 |
|
|
+ |
– |
|
+ |
– |
|
– |
|
+ |
|
– |
+ |
|
110,0 |
110,7 |
|
-0,7 |
|
2 |
|
|
+ |
+ |
|
– |
– |
|
– |
|
– |
|
+ |
+ |
|
73,8 |
74,1 |
|
-0,3 |
|
3 |
|
|
+ |
– |
|
– |
+ |
|
– |
|
– |
|
– |
– |
|
107,6 |
107,7 |
|
-0,1 |
|
4 |
|
|
+ |
+ |
|
+ |
+ |
|
– |
|
+ |
|
+ |
– |
|
112,3 |
11,5 |
|
0,8 |
|
5 |
|
|
+ |
– |
|
– |
– |
|
+ |
|
+ |
|
+ |
– |
|
77,7 |
76,5 |
|
1,2 |
|
6 |
|
|
+ |
+ |
|
+ |
– |
|
+ |
|
– |
|
– |
– |
|
80,0 |
80,3 |
|
-0,3 |
|
7 |
|
|
+ |
– |
|
+ |
+ |
|
+ |
|
– |
|
И |
|
113,8 |
113,9 |
|
-0,1 |
||
|
|
|
|
|
|
+ |
+ |
|
|
|||||||||||
8 |
|
|
+ |
+ |
|
– |
+ |
|
+ |
|
+ |
|
– |
+ |
|
76,5 |
77,3 |
|
-0,8 |
|
9 |
|
|
– |
0 |
|
0 |
0 |
|
0 |
|
|
Д |
0 |
|
91,5 |
94,0 |
|
-2,5 |
||
|
|
|
|
|
|
0 |
0 |
|
|
|||||||||||
UX |
|
Абсолютная |
погрешность |
измерения |
U UX /100%, |
где |
||||||||||||||
– максимальное значение выходного напряжения. В нашем слу- |
||||||||||||||||||||
чае |
|
|
|
|
|
|
б |
|
|
|
|
|
|
|
|
|
||||
U 0,04 В. Для доверительной вероятности PД 0,95 величина |
||||||||||||||||||||
доверительного интервала |
U |
2S , |
где SU |
– ошибка воспроизво- |
||||||||||||||||
|
|
|
|
|
|
и |
|
|
|
U |
|
|
|
|
|
|
|
|
||
|
|
|
А |
|
||||||||||||||||
димости результата |
змерения выходного |
напряжения. Отсюда |
||||||||||||||||||
SU 0,02 В. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
Коэффициент ус лен я определяется косвенным методом, |
по- |
этому среднее квадрат ческое отклонение результата измерения –
ошибка эксперимента S y |
– вычисляется по формуле |
|
|||||||
|
S y K0 |
SU |
|
|
SU |
|
|
(4.52) |
|
|
|
|
|
|
|||||
С |
U UВЫХ0 |
UВХ |
|
|
|||||
и составляет S y 1,5. |
|
|
|
|
|
|
|
|
|
Значения коэффициентов уравнения регрессии |
|
||||||||
|
~ |
b1X1 b2X2 b3X3 b4X4 |
|
||||||
|
y b0 |
||||||||
|
b13X1X3 |
b23X2X3 b34X3X4 |
(4.53) |
||||||
вычисляются по формуле (4.13). |
|
|
|
|
|
|
|||
Расчеты дают следующие результаты: |
|
|
|||||||
b0 94; |
b1 8,2; |
b2 10,1; |
b3 8,6; |
b4 7,0; |
|||||
b13 0,15; |
b23 0,44; b34 |
0,44. |
|
|
82
Ошибка определения коэффициентов bu вычисляется по формуле (4.18) и составляет
Sb S y |
N 1,5 |
8 0,5. |
Для выявления значимости коэффициентов уравнения регрессии |
||
строим доверительный интервал шириной |
2 b 2tT Sb 2 1,86 0,5 1,86.
Табличное значение t-критерия Стьюдента определяем для уровня значимости 0,05 и числа степеней свободы d.f N 8.
Коэффициенты, характеризующие взаимодействие факторов, оказались статистически незначимыми, поэтому уравнение регрессии имеет вид
~ |
|
|
|
|
|
|
|
|
|
|
8,6X3 7X4 . |
(4.54) |
y 94 8,2X1 10,1X2 |
||||||||||||
Из уравнения (4.54) видно, что коэффициент усиления убывает с |
||||||||||||
увеличением сопротивлений R1, R4 и возрастает с увеличением R2, |
||||||||||||
R3 , что согласуется с теоретической формулой (4.43). |
|
|||||||||||
Проверим адекватность |
|
полученной математической |
модели |
|||||||||
|
|
|
|
А |
~ |
|
||||||
(4.54) экспериментальным данным. РасчетныеИзначения y и разности |
||||||||||||
~ |
|
|
|
|
|
|
|
|
|
|
|
|
y y приведены в табл. 4.16. Дисперсия адекватности (4.22) |
|
|||||||||||
2 |
1 |
|
N |
|
|
Д2 |
|
|||||
|
|
|
|
|||||||||
|
|
|
|
|
|
~ |
|
|
||||
Sад |
|
|
i 1 yi yi |
0,85; |
|
|||||||
8 4 |
|
|||||||||||
дисперсия воспроизводимости |
|
|
|
|
|
|
|
|
|
|||
|
S2 y 1,52 2,25; |
|
||||||||||
т.к. числовое значен е д |
бсперсии воспроизводимости больше число- |
|||||||||||
вого значения дисперс |
адекватности, расчетное значение критерия |
|||||||||||
и |
|
|
|
|
|
|
|
|
|
|||
Фишера определяется по формуле |
|
|
|
|
|
|
||||||
F |
S2 y |
2,25 |
|
|
||||||||
|
|
|
|
|
|
|
2,647. |
|
||||
|
|
|
|
|
|
|
||||||
С |
|
|
Sад2 |
|
|
|
0,85 |
|
|
|
||
Для d.f1 8 4 4; |
d.f2 8 |
|
|
и 0,05 критическое значение |
FT 3,838. Так как F FT , то уравнение (4.54) адекватно описывает функцию отклика.
По полученной модели могут быть вычислены коэффициенты влияния (чувствительности) факторов на коэффициент усиления и определены числовые характеристики KU как случайной величины в зависимости от разброса сопротивлений резисторов. Для этого в уравнении (4.54) необходимо перейти от кодированных факторов к физическим переменным. Тогда
83