Проведение вычислительного эксперимента

дипломная работа

2.3 ОПИСАНИЕ АЛГОРИТМА ПОИСКА ЭКСТРЕМУМА С ЗАПОМИНАНИЕМ ЭКСТРЕМУМА

Рассмотрим более подробно описанный выше алгоритм поиска экстремума с запоминанием экстремума.

Алгоритм запоминания экстремума заключается в использовании разности между текущим и экстремальным значением выходной величины для нахождения момента реверса системы и имеет вид:

= V1,V(0) = V1>0 (2.19)

Функция V(t) меняет знак на противоположный в моменты времени t , для которых выполняется соотношение:

с(tj) = -A0, с(tj)<0 , (2.20)

где

p(t) = q(t) - max q(ф),

0< ф < t,

A0 > 0 , с(t)=.

В силу этого алгоритма в вычислительном устройстве запоминается максимальное значение выходного сигнала q измерительного устройства, реализовавшееся до текущего момента времени t и непрерывно вычисляется функция:

с(t) = q(t) - max q(ф)

0 < ф ? t

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

В начальный момент t = 0,V(0) = V1 > 0 и координата х начинает увеличиваться. Если х(0) находится слева от точки х*, в которой достигается у* - максимум нелинейной характеристики, то точка х движется к х*, функция с(t) при тождественно равна 0, т.к. max q(t), 0 < ф ? t достигается в текущий момент времени t. При этом переключения функции V(t) не происходит пока точка х не перескочит положение х*, величина у, а значит и Z и q не начнут уменьшаться, функция с(t) станет отрицательной, т.к. max q(t) > q(t), 0 < ф ? t. В силу приведенного правила переключение произойдет при с(t) = -А0, когда точка х окажется на нисходящей ветви нелинейной характеристики на некотором расстоянии права от х* (рисунок 2.5).

Рисунок 2.5 - Поиск экстремума с запоминанием

После переключения V(t) начинается возвратное движение точки х, при котором она через некоторое время окажется слева от х и второе переключение произойдет, когда снова, уменьшаясь, функция с(t) не достигнет -А0. Далее возникает колебательное движение точки х в некоторой двусторонней окрестности точки х, а значит и выходной координаты Z относительно ее максимального значения (рисунок 2.6)

Если начальное положение х(0) точки х находится правее точки х*, то сначала точка х удаляется от х*, поскольку V(0) = V1>0, затем происходит переключение и дальнейший процесс аналогичен описанному выше.

Рисунок 2.6 - Колебательные движения в окрестности экстремума

Параметры V1 и А0 в законе управления выбираются так, чтобы обеспечить попадание величины Z(t) в заданную достаточно малую окрестность Zmax за приемлемое время и дальнейшее пребывание ее в этой окрестности, т.е. выполнении неравенства:

(2.21)

Чем больше величина V1, тем быстрее изменяется координата x(t), а значит и Z(t) и тем быстрее Z(t) может попасть в окрестность точки Zmax.

Однако при увеличении Vl возрастают динамические ошибки

воспроизведения входных сигналов линейной частью математической модели и измерительным устройством. Это может привести к “проскакиванию“ координатой Z(t) своего максимального значения Zmax и к дальнейшим колебаниям относительно Zmax со значительной амплитудой. При малых значениях V1 перемещение Z(t) в точку Zmax произойдет за значительное время.

2.4 АНАЛИЗ ТОЧНОСТИ РАБОТЫ СЭУ С ВЫДЕЛЯЕМОЙ НЕЛИНЕЙНОЙ ХАРАКТЕРИСТИКОЙ

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

Для начала рассмотрим метод оценивания точности функционирования СЭР со стационарной нелинейной характеристикой.

Пусть в алгоритме проведения экстремального эксперимента (законе управления) величина А0 определяется соотношением:

(2.22)

Предполагается, что возмущение ц0 и помеха ш0 в измерительном тракте равны нулю ц0 = ш0 = 0, тогда

(2.23)

Если в начальный момент времени t=0 выполняются соответствующие установившемуся режиму при отсутствии помех начальные условия:

(2.24)

то через заведомо попадает и в дальнейшем при t> уже не выйдет из окрестности максимума нелинейной характеристики у* = f(x*) = Zmax, заданной неравенством:

?А1 (t?),

где - момент второго (после начала работы системы) переключения управления V(t).

Здесь

А1= 2А0 + f1V1д1 +ц0д2.

В нашем случае:

А1 = 2A0+f1V1д1, (2.25)

Используемые в соотношениях величины f1, V1 определены в формулах (2.6)

и (2.19). Параметры д1, , д2, выражается через характеристики корней характеристических уравнений, соответствующих дифференциальным уравнениям линейной части и измерительного устройства.

Обозначим через Рj (j =) корни характеристического уравнения

N(p) = an+1pn + ... + а2р + а1 =0 модели линейной части производственного процесса.

Тогда д1 определяется как:

для объекта 1-го порядка:

(2.26)

Для объекта 2-го порядка:

(2.27)

Для объекта 3-го порядка:

(2.28)

Обозначим через ( j =) - корни характеристического уравнения для измерительного устройства, тогда определяется как:

(2.29)

Передаточная функция измерительного устройства имеет вид:

(2.30)

Характеристики большинства объектов не остаются с течением времени неизменными, а изменяются иногда со значительными скоростями. Кроме того, иногда характеристики объектов не определяются аналитически, а затраты на экспериментальное определение могут быть очень велики. В то же время любая характеристика аппроксимируется отрезками достаточно простых кривых - парабол не выше второго порядка, причем вид аппроксимирующей параболы (коэффициенты управления) зависит от значения входной координаты. Пределы измерения коэффициентов в большинстве случаев определяются без значительных затрат. Учитывая, что при поиске экстремума входная координата изменяется во времени по определенному или случайному закону, характеристику любого экстремального объекта можно представить квадратичной параболой, дрейфующей во времени.

Выделяют два основных случая изменения характеристик во времени: на объект воздействует помеха; происходит регулярный дрейф характеристики. Обоим случаям нетрудно найти аналоги в теории обычных систем автоматического регулирования.

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

Уравнение дрейфующей параболы имеет вид:

у = С2 -(х- sin(w * t))2 + С1, (2.31)

где с1= --1,

с2 = 2,

щ = 0.01.

На рисунке 2.7 представлены характеристики объекта, соответствующие времени t1 = 0, t2 = 25,t3 = 50.

Их экстремумы равны:

х1* = 0 х2* =0.3 х3*=0.5

y1 = 2 y2* =2 y3*=2

Рисунок 2.7 - Горизонтальный дрейф нелинейной характеристики

Исследуем точность функционирования СЭР с дрейфующей нелинейной характеристикой.

Из-за нестационарности НХ естественным показателем точности функционирования представляется величина модуля разности текущего значения выхода управляемого объекта и текущего значения максимума нелинейной характеристики:

е(t) = z(t) - m(t). (2.32)

Конечно, в силу ряда технологических причин, величина е (0) при t=0 может быть достаточно большой, и анализировать точность СЭР целесообразно в “квазистационарном ” режиме “захвата” экстремума нелинейной характеристики, наступающем после первого переключения u(t) в момент t1. Поскольку сигнал управления имеет релейный характер, то в описываемых условиях “квазистационарный ” режим будет колебательным.

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

Пусть такая система описывается дифференциальным уравнением

R(w) = S(v)

или передаточной функцией

(2.33)

где v(t) - внешнее воздействие, w(t) - выходной сигнал;

V(p),W(p) - их изображения по Лапласу при нулевых начальных условиях для самих функций и их производных.

Назовем максимальным нормированным отклонением выходного сигнала w(t) (k-ой производной w(k)(t) выходного сигнала) системы соответственно величины

максимальной нормированной ошибкой воспроизведения входного сигнала системой величину

В этих определениях предполагается, что система до момента включения t = 0 находится в состоянии покоя, т.е. w(t) = v(t) = 0 при t ? 0.

Величины Д, д называются точностными показателями качества (ТПК) системы.

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

Если

(2.34)

Заметим еще, что

Д(w, G(p)) = Д(w, pG(p)).

Максимальная нормированная ошибка д(G(p)) и максимальные нормированные отклонения Д(w(k), G(p)) также определяются и при условиях

w(t) = v(t) = с при t ? 0.

В этом случае под Д(w, G(p)) понимается максимальное возможное

отклонение от величины с.

Опишем дополнительные условия, накладываемые на поведение нестационарной нелинейной характеристики. Требования унимодальности f(t,x) при любом фиксированном t недостаточно для получения содержательных утверждений о точности функционирования СЭР. Положим

r = x-xm(t). (2.35)

y = f(t,х) = f(t,r + хm(t)) ? g(t, r) (2.36)

Функция g(t,r) при фиксированном t описывает движение по нелинейной характеристике, а при фиксированном r - её деформацию во времени. Функция хm (t) отображает перемещение абсциссы точки максимума нелинейной характеристики; при этом g(t,0) = m(t).

Используя введенные выше обозначения (2.34) и (2.36), положим:

2.37

Априорная информация о поведении нелинейной характеристики исчерпывается соотношениями:

у = с2х2 + с1 - уравнение стационарной нелинейной характеристики,

r0 = -0.5 - начальное положение х,

r = x-xm(t),

Отсюда

y = c2(x-xm(t))2 +с1,

хm = б sin(щt).

Таким образом

у = С2 * (х - sin(w * t))2 + С1, a g = С2 * r2 + С1.

Для объекта 1-го порядка:

Где б = (б- положительная часть вещественного корня уравнения ф1p+1=0). Тогда

Для объекта 2-го порядка:

Соответствующее Gy (р) характеристическое уравнение имеет корни

Для объекта 3-го порядка:

Представим ПФ Gy(р) в виде:

Тогда соответствующее Gv2(p) характеристическое уравнение имеет корни:

р12 =-б(1 ± iм),

Если описанная выше СЭР до начала работы в момент /==0 находилась в состоянии покоя:

q(t) = z(t) = y(t) = f(0, х(0)) = g(0, r(0)), x(t) = x(0) при t < 0,

r(0) = x(0) - xm (0) < 0, r(0) ? [-r0,0); (2.38)

величины А, ф0 заданы и выполнены условия

k > О, В < r0, (2.39)

то после первого переключения сигнала управления u(t) в момент t1 > 0 отклонение текущего значения абсциссы НХ от текущего значения абсциссы её максимума и отклонение текущего значения выхода СЭР z(t) от текущего значения m(t) максимума НХ заведомо удовлетворяют неравенствам

(2.40)

Таблица 1 - параметры системы для объекта 1-го порядка

Таблица 2 - параметры системы для объекта 2-го порядка

Более сложен вопрос о выборе желаемого значения величины х1 - скорости перемещения абсциссы нелинейной характеристики x(t). При малом х1, абсцисса xm(t) текущего экстремума нелинейной характеристики будет удаляться от x(t); большое х1, увеличивает величину А0, связанную с динамическими погрешностями линейной части управляемого объекта, а значит и величину В. Существует, следовательно, оптимальное значение скорости х1.

Выполняя расчет точностных показателей качества для разных порядков объекта управления (Приложение Б), мы получаем данные для минимизации величины В по параметру х1 которые приведены в таблицах 1, 2, 3.

Таблица 3 - параметры системы для объекта 3-го порядка

Минимизация величины В по параметру х приводит при х оценкам:

- для объекта 1-го порядка - |г(/)| < В = 0.282, z(t) - m{t) < 0.106;

- для объекта 2-го порядка - |г(/)| < В = 0.364, z(t) - m(t) < 0.176;

- для объекта 3-го порядка - |г(/)| < В = 0.455, z(t) - m(t) < 0.275 .

j = 0.055 к

Делись добром ;)