При анализе переходных процессов

Задаём:

- диапазон времени 4m (4 миллисекунды); вывод расчётов Normal;

- максимальный шаг по времени 0.0001; начальные условия Zero (нулевые);

- количество точек 51; Ö Расчёт по постоянному току;

- температура List - 10, +50; Ö Автомасштаб;

- квадрат первый слева: масштаб по оси Х – линейный;

- квадрат второй слева: масштаб по оси Y – линейный;

- квадрат третий слева: цвет линии для графика функции – любой;

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

- в ячейке P – порядковый номер координатной плоскости для вывода графика функции – номер 1;

- в ячейке X Expression определено выражение для оси Х – время T;

- в ячейке Y Expression определено выражение для оси Y – значение напряжения в точке 5 схемы (выход усилительного каскада) – v(5);

- в ячейках X Range и Y Range определяются диапазоны представления на координатной плоскости значений аргумента и функции, соответственно. Однако, поскольку в нашем случае активирована опция Ö Автомасштаб, то значения ячеек X Range и Y Range системой MicroCAP 8 игнорируются.

Для запуска моделирования переходных процессов необходимо нажать кнопку Запуск. В результате на экране появится окно с двумя графиками зависимости выходного напряжения (в точке 5) усилительного каскададля значений температуры окружающей среды -10^оС и +50^оС при V 1 = 22В, R 4 = 40 Ом (рисунок 3.18).

температура окружающей среды -10^оС ()и +50^оС ();

V 1 = 22В (); R 4 = 40 Ом ()

Рисунок 3.18 – Графики зависимости выходного напряжения (точка 5) усилительного каскада

Выполним двойной щелчок на заголовке графиков УСКОЗЛ1-163.CIR Temperature=-10…50. В результате появляется диалоговое окно задания параметров оформления графиков Свойства (рисунок 3.19). В этом диалоговом окне можно сделать настройки, которые повышают удобство работы с графиками и помогают в последующем оформлении отчета. Выберем в диалоговом окне Свойства вкладку Цвета, шрифты, линии. На этой вкладке представлен список Объекты (рисунок 3.20). Редактирование типов, размеров и цветов элементов рисунка осуществляется при выборе соответствующего пункта из списка Объекты, например:

- редактирование типа, цвета и высоты шрифта текста;

- редактирование элементов координатной сетки;

- редактирование типа, цвета и высоты шрифта для меток на графике;

- редактирование толщины и цвета кривой и т.д.

Рисунок 3.19 – Верхняя часть диалогового окно Свойства

Рисунок 3.20 – Вкладка Цвета, шрифты и линии

диалогового окна Свойства

Схемы и графики, полученные в системе MicroCAP 8, можно передавать во внешние приложения с помощью команды Редактирование / Копировать в буфер обмена / Копировать видимую область в BMP формат. При этом графический образ копируется в буфер обмена Windows. Затем при помощи команды Правка / Вставить любого внешнего приложения, например текстового редактора Microsoft Word, скопированный графический образ размещается в требуемом месте.

Для снятия численных значений амплитуды с графиков зависимости выходного напряжения (точка 5) усилительного каскада следует воспользоваться режимом «электронного курсора». Для этого, находясь в окне отображения графиков зависимостей, следует нажать на пиктограмму (на рисунке 3.18 третья слева на панели инструментов). После этого на экране появляются изображения двух электронных курсоров, совмещенных с одним из графиков – выходным напряжением усилителя при температуре среды -10^оС(рисунок 3.21). Заметим, что электронные курсоры можно произвольно перемещать вдоль всего графика функции с помощью мыши или с помощью клавишей-стрелок и ®. Перемещение правого электронного курсора на экране происходит либо при нажатой правой клавиши мыши, либо при нажатой клавише Shift (в случае перемещения с помощью стрелок). Для нахождения локальных максимумов функции можно воспользоваться пиктограммой , расположенной на панели инструментов. При нажатии на эту пиктограмму левый электронный курсор перемещается на ближайший локальный максимум функции. Для того чтобы задействовать для поиска максимумов правый электронный курсор, необходимо его сначала активировать. Делается это правой кнопкой мыши; затем, также как и в предыдущем случае, нажатие на пиктограмму заставляет правый электронный курсор перемещаться на ближайший локальный максимум.

температура окружающей среды -10^оС ()и +50^оС ();

V 1 = 22В (); R 4 = 40 Ом (); электронные курсоры совмещены

с графиком, соответствующим температуре среды -10^оС

Рисунок 3.21 – Реализация опыта №1 ПФЭ

Из рисунка 3.21 видно, что значения амплитуды выходного напряжения при температуре -10^оС составляют 460 мВ и 465 мВ. Более точные значения текущих координат электронного курсора отображаются под координатной плоскостью в столбцах Левый и Правый (первая строка): 459.990 мВ и 465.007 мВ.

Для перехода к измерению другой зависимости – выходного напряжения при температуре среды +50 ^оС – необходимо с помощью клавишей-стрелок ­ или ¯ переместить электронные курсоры на соседний график (рисунок 3.22). Из рисунка 3.22 видно, что значения амплитуды выходного напряжения при температуре +50^оС составляют 742.054 мВ и 748.991 мВ.

температура окружающей среды -10^оС ()и +50^оС ();

V 1 = 22В (); R 4 = 40 Ом (); электронные курсоры совмещены

с графиком, соответствующим температуре среды +50^оС

Рисунок 3.22 – Реализация опыта №2 ПФЭ

Подготавливаем в системе MicroCAP 8 схему, изображённую на рисунке 3.15, задав значения V 1 = 22В, R 4 = 80 Ом (, ). В этом случае при температуре окружающей среды -10^оС () значения амплитуды на выходе усилителя (точка 5) составляют 566.739 мВ и 565.848 мВ (рисунок 3.23), а при температуре окружающей среды +50^оС() значения амплитуды составляют 908.283 мВ и 911.786 мВ(рисунок 3.24).

температура окружающей среды -10^оСи +50^оС при V 1 = 22В () и

R 4 = 80 Ом (); электронные курсоры совмещены с графиком,

соответствующим температуре среды -10^оС ()

Рисунок 3.23 – Реализация опыта №5 ПФЭ

температура окружающей среды -10^оСи +50^оС при V 1 = 22В () и

R 4 = 80 Ом (); электронные курсоры совмещены с графиком,

соответствующим температуре среды +50^оС ()

Рисунок 3.24 – Реализация опыта №6 ПФЭ

Подготавливаем в системе MicroCAP 8 схему, изображённую на рисунке 3.15, задав значения V 1 = 24В, R 4 = 40 Ом (, ). В этом случае при температуре окружающей среды -10^оС () значения амплитуды составляют 691.959 мВ и 693.640 мВ (рисунок 3.25), а при температуре окружающей среды +50^оС() значения амплитуды составляют 791.565 мВ и 800.017 мВ (рисунок 3.26).

температура окружающей среды -10^оСи +50^оС при V 1 = 24В () и

R 4 = 40 Ом (); электронные курсоры совмещены с графиком,

соответствующим температуре среды -10^оС ()

Рисунок 3.25 – Реализация опыта №3 ПФЭ

Подготавливаем в системе MicroCAP 8 схему, изображённую на рисунке 3.15, задав значения V 1 = 24В, R 4 = 80 Ом (, ). В этом случае при температуре окружающей среды -10^оСзначения амплитуды составляют 847.627 мВ и 845.649 мВ (рисунок 3.27), а при температуре окружающей среды +50^оСзначения амплитуды составляют 961.79 мВ и 975.502 мВ(рисунок 3.28).

Полученные в результате анализа переходных процессов численные значения выходного напряжения y_{m g} = V 5_γ для различных сочетаний нормированных переменных (рисунки 3.21 – 3.28) заносим в матрицу планирования ПФЭ 2 ^п для п = 3 и числа измерений γ = 2 (таблица 3.5).

температура окружающей среды -10^оСи +50^оС при V 1 = 24В () и

R 4 = 40 Ом (); электронные курсоры совмещены с графиком,

соответствующим температуре среды +50^оС ()

Рисунок 3.26 – Реализация опыта №4 ПФЭ

температура окружающей среды -10^оСи +50^оС при V 1 = 24В () и

R 4 = 80 Ом (); электронные курсоры совмещены с графиком,

соответствующим температуре среды -10^оС ()

Рисунок 3.27 – Реализация опыта №7 ПФЭ

температура окружающей среды -10^оСи +50^оС при V 1 = 24В () и

R 4 = 80 Ом (); электронные курсоры совмещены с графиком,

соответствующим температуре среды +50^оС ()

Рисунок 3.28 – Реализация опыта №8 ПФЭ

Таблица 3.5 – Матрица планирования ПФЭ 2 ^п для п = 3 и γ = 2 (ПФЭ 2³)

m			-1	…………	ym 1 = V 5 m 1	ym 2 = V 5 m 2	т (ym) = (уm 1 + уm 2)/2 = т (V 5 т)
	-1	-1	-1	…	459.99	465.007	т (у 1) = 462.499
	+1	-1	-1	…………	742.054	748.991	т (у 2) = 745.523
	-1	+1	-1	…	691.959	693.64	т (у 3) = 692.799
	+1	+1	-1	…………	791.565	800.017	т (у 4) = 795.791
	-1	-1	+1	…	565.848	566.739	т (у 5) = 566.293
	+1	-1	+1	…………	908.283	911.786	т (у 6) = 910.034
	-1	+1	+1	…	846.627	845.649	т (у 7) = 846.138
	+1	+1	+1	…	961.79	975.502	т (у 8) = 968.646

Среднее значение т (у_т) в каждой т -ой строке в таблице 3.5 вычислено по формуле (3.22):

.

Для вычисления коэффициентов полинома (3.12):

используются выражения (3.16), в которых следует принять N = 8, а величины у_т – равными средним измеренным значениям т (у_т), взятым из таблицы 3.5 (у_т = т (у_т)). Тогда:

Если пренебречь квадратичными членами, то полином (3.12) для нашего примера примет вид:

. (3.12б)

Перейдём от коэффициентов В, которые получаются из опыта, к коэффициентам β в натуральном масштабе в полиноме (3.8). При использовании полиномов первого порядка расчёт проводим по формулам (3.15):

;

; ;

.

Полином первого порядка с коэффициентами β в натуральном масштабе для нашего примера примет вид:

V 5 = (-1325.193 + 3.551· Т + 77.378· V 1 + 3.716· R 4) мВ.

Значения строчной дисперсии вычислены по формуле (3.23), которая при γ = 2 примет вид:

.

D (y ₁) = (459.99 – 462.499)² +(465.007 – 462.499)² = 12.585мВ²;

D (y ₂) = (742.054 – 745.523)² + (748.991 – 745.523)² = 24.061 мВ²;

D (y ₃) = (691.959 – 692.799)² + (693.64 – 692.799)² = 1.413 мВ²;

D (y ₄) = (791.565 – 795.791)² + (800.017 – 795.791)² = 35.718 мВ²;

D (y ₅) = (565.848 – 566.293)² + (566.739 – 566.293)² = 0.397 мВ²;

D (y ₆) = (908.283 – 910.034)² + (911.786 – 910.034)² = 6.136мВ²;

D (y ₇) = (846.627 – 846.138)² + (845.649 – 846.138)² = 0.478 мВ²;

D (y ₈) = (961.79 – 968.646)² + (975.502 – 968.646)² = 94.009 мВ².

Исходя из строчной дисперсии, по формуле (3.24) находим дисперсию воспроизводимости:

Зная дисперсию воспроизводимости, находим дисперсию полиномиальных коэффициентов B_i:

.

По формуле (3.30) вычисляем среднее квадратическое отклонение коэффициентов:

σ(B_i) = [ D (B_i)]^0.5 = 2.731^0.5= 1.653мВ.

Значения коэффициентов уравнения (3.12б) существенно больше. Следовательно, можно полагать, что все оставшиеся коэффициенты полинома значимы.

Дополнительно производим проверку значимости с использованием критерия Стьюдента. Для уровня значимости, равного a = 0.05 и степеней свободы М = (N ·γ – N) = 8 · 2 – 8 = 8, из таблицы 3.3 находим t _ТАБЛ(УЗ) = 2.306.

Тогда:

.

Полученное число 3.811 меньше коэффициентов B_i полинома (3.12б):

B ₀ = 748.465 мВ; B ₁ = 106.533 мВ; B ₂ = 77.378 мВ; B ₃ = 74.313 мВ.

Это согласно (3.28) является подтверждением того, что все коэффициенты полинома значимы.

Производим проверку адекватности путём сравнения результатов эксперимента т (у_{т g}) для разных сочетаний переменных с теми результатами, которые при том же сочетании переменных получены для параметра после расчета по полиному (3.12б). Результаты расчета помещаем в таблицу 3.6.

Таблица 3.6 – Сравнение результатов эксперимента т (у_{т g}) с результатами, которые получены после расчета по полиному

m	т (ут)	ут РАСЧ	[ т (ут) – ут РАСЧ]^2
	462.499	748.465 + 106.533·(-1)+ 77.378·(-1)+ 74.313·(-1) = 490.242	769.681
	745.523	748.465 + 106.533·(+1)+ 77.378·(-1)+ 74.313·(-1) = 703.308	1782.09
	692.799	748.465 + 106.533·(-1)+ 77.378·(+1)+ 74.313·(-1) = 644.998	2284.98
	795.791	748.465 +106.533·(+1)+ 77.378·(+1)+ 74.313·(-1) = 858.064	3877.94
	566.293	748.465 + 106.533·(-1)+ 77.378·(-1)+ 74.313·(+1) = 638.867	5266.88
	910.034	748.465 + 106.533·(+1)+ 77.378·(-1)+ 74.313·(+1) = 851.933	3375.80
	846.138	748.465 + 106.533·(-1)+ 77.378·(+1)+ 74.313·(+1) = 792.623	2757.81
	968.646	748.465+ 106.533·(+1)+ 77.378·(+1)+ 74.313·(+1) = 1006.69	1447.29

Приведём результаты расчета выходного параметра у_{т РАСЧ} по этому полиному с нормированными переменными:

. (3.12б)

Затем по формуле (3.31) вычисляем дисперсию адекватности:

,

где т (у_{т g}) – среднее значение результата эксперимента, взятое из т -ой строки таблицы 3.5; у_{т РАСЧ} – результаты расчета по полиному (3.12б) для того же сочетания переменных; d – количество членов, оставленных в полиноме.

Для D _АД(y)принимается (N – d) степеней свободы, а для D _B(у)берется N ·(γ – 1) степеней свободы.

Для проверки адекватности по формуле (3.32) вычисляем случайную величину F _АД распределения Фишера:

.

Адекватность полинома установим с помощью таблицы 3.4 распределения Фишера при сравнении величины F _АДс F _ТАБЛ. Если F _АД< F _ТАБЛ, то принимается решение об адекватности полинома результатам опыта. В противном случае принимается решение о его неадекватности.

Для рассматриваемого примера:

- количество членов, оставленных в полиноме d = 4;

- количество степеней свободы для дисперсии адекватности:

(N – d) = 8 – 4 = 4, т.е. М ′ = 4;

- количество степеней свободы для дисперсии воспроизводимости:

N ·(γ – 1) = 8·(2 – 1) = 8, т.е. М = 8;

- из таблицы 3.4 при уровне значимости a = 0.05 с учетом того, что М = 8, a М ′ = 4, получаем F _ТАБЛ = 6.04, а при уровне значимости a = 0.01 F _ТАБЛ = 14.8.

Как видно, F _АД = 246.714 > F _ТАБЛ = 6.04 при уровне значимости a = 0.05 и F _АД = 246.714 > F _ТАБЛ = 14.8 при уровне значимости a = 0.01. Следовательно, в отличие от предыдущей задачи полученный линейный полином не адекватен результатам опыта с вероятностью ошибки в этом решении не меньше, чем 1 –a = 0.99.

Для обеспечения адекватности в полином следует ввести коэффициенты взаимодействия. В рассматриваемом примере можно вычислить коэффициенты взаимодействия B_ij 2-го порядка [20]:

(3.18)

и получить неполный полином 2-го порядка. Как упоминалось ранее, получение коэффициентов при из результатов ПФЭ в нашем случае невозможно.

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

С помощью вспомогательной таблицы 3.8 вычислим относительные отклонения измеренных значений от расчётных (δ%) при использовании линейного полинома и от расчётных значений при использовании неполного полинома 2-го порядка (δ% _УТОЧ). Эти отклонения рассчитаем по формулам:

(3.33)

и

, (3.34)

где

. (3.35)

По данным таблицы 3.8 видно уменьшение относительных отклонений измеренных значений от расчётных при использовании полинома 2-го порядка (δ%_УТОЧ < δ%).

Таблица 3.7 – Вспомогательная таблица для вычисления коэффициентов взаимодействия 2-го порядка

m			-1	…………
	-1	-1	-1	…	462.499	462.499
	+1	-1	-1	…………	-745.523	745.523
	-1	+1	-1	…	-692.799	-692.799
	+1	+1	-1	…………	795.791	-795.791
	-1	-1	+1	…	566.293	-566.293
	+1	-1	+1	…………	-910.034	-910.034
	-1	+1	+1	…	-846.138	846.138
	+1	+1	+1	…	968.646	968.646
	=

Таблица 3.8 – Вспомогательная таблица для вычисления относительных отклонений измеренных значений от расчётных

m				т (ут)	ут РАСЧ	δ%		δ%УТОЧ
	-1	-1	-1	462.499	490.242	6.00	490.242 - 50.158 + 7.236 = 447.319	3.282
	+1	-1	-1	745.523	703.308	5.66	703.308 + 50.158 + 7.236 = 760.702	2.036
	-1	+1	-1	692.799	644.998	6.90	644.998 + 50.158 - 7.236 = 687.92	0.704
	+1	+1	-1	795.791	858.064	7.83	858.064 - 50.158 - 7.236 = 800.67	0.613
	-1	-1	+1	566.293	638.867	12.8	638.867 - 50.158 - 7.236 = 581.473	2.68
	+1	-1	+1	910.034	851.933	6.39	851.933 + 50.158 - 7.236 = 894.855	1.668
	-1	+1	+1	846.138	793.623	6.21	793.623 + 50.158 + 7.236 = 851.017	0.577
	+1	+1	+1	968.646	1006.69	3.93	1006.69 - 50.158 + 7.236 = 963.767	0.504

Вычислим коэффициенты взаимодействия 2-го порядка b _ij в натуральном масштабе:

; (3.36)

;

.

В двух последних формулах шаг изменения температуры окружающей среды Δ х ₁ = Δ Т = ±30^оС, шаг изменения напряжения питания Δ х ₂ = Δ V 1 = ±1В, шаг изменения сопротивления нагрузки Δ х ₃ = Δ R 4 = ±20 Ом.

Неполный полином 2-го порядка с коэффициентами βв натуральном масштабе для нашего примера примет вид:

Для анализа точности и стабильности выходного параметра вычислим по формуле (3.21) коэффициенты k_{i ОТН} в уравнении отклонений:

и коэффициенты взаимодействия в уравнении отклонений

, (3.37)

гденоминальная величина:

у ₀ = V5 ₀ = (829.269 + 828.989)/2 = 829.129 мВ

это значение выходного параметра V5, определённое из графика (рисунок 3.29) для номинальных значений исходных параметров: температуры окружающей среды 20^оС (х _1,0), напряжения питания V 1₀ = 23В(х _2,0) и R 4₀ = 60 Ом(х _3,0).

;

;

;

;

.

Получим явный вид уравнения отклонений неполного полинома 2-го порядка:

. (3.38)

Подставим численные значения величин коэффициентов k_{i ОТН} и k_{ij ОТН}:

температура окружающей среды 20^оС(х _1,0),

напряжение питания V 1₀ = 23В(х _2,0) и R 4₀ = 60 Ом (х _3,0)

Рисунок 3.29 – График зависимости выходного напряжения

усилительного каскададля номинальных значений

Дата добавления: 2014-12-27