Адаптивный прогнозатор производственных показателей Советский патент 1991 года по МПК G06F15/46 

Фиг i

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

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

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

На фиг. 1 дано схема прогнозатора; на фиг. 2 - схема первого интегратора; на фиг. 3 - временная диаграмма его работы.

На фиг. 1 приняты следующие обозначения; x(t+ 0)- прогнозируемое нч интервале времени t, t + в значение пеовичной переменной; Q(t + в)- прогнозируемое значение расчетного интегрального показателя Q; Ki(t -f- G) - нормативны/, прогноз коэффициент пропорциональности; x(tj- фактическое значение первичной переменной в момент времени t; d - параметр идентификации: в (tK -1) - интервал прогнозиров - ния (разность между конечным t и текущим t моментами времени) t 1, t - и конечный моменты зрамзн - драницы интервала программирования рзбсгы производственной системы), а- параметр расчета интегрального показателя; Q(t) - заданное (плановое) значение интегрального показателя на момент t .

Адаптивный Прогнозатор содержит задатчик 1 программного изменения первичного параметра, с первого по седьмой фильтры 2-8 низкой частоты, с первого по седьмой блоки 9-15 вычитания, с первого по

шестой блоки 16-21 умножения, с первого по седьмой сумматоры 22-28,квадратор 29, первый 30 и второй 31 блоки деления, первый 32 и второй 33 интеграторы.

На фиг. 1 обозначены входы прогнозатора для измерения первичного параметра

34,для ввода коэффициента идентификации

35,интервала экстраполяции 36, конечного 37 и начального 38 моментов интервала прогнозирования, коэффициента расчета интегрального показателя 39, требующегося значения интегрального показателя 40 и выходы прогнозируемых значений первичного показателя 41, интегрального показателя 42 и расчетного коэффициента интенсивности

работы производственной системы 43. Функционально элементы, образующие прогно- зэтор. могут быть объединены в ряд блоков, котор1, обозначены на чертеже; идентификатор 44, окстраполятор 45, блок 46 пересчетного прогнозирования, блок 47 расчета ,гегрального показателя и блок 48 расчета коэффициента интенсивности работы производственной системы.

Интегратор 32 включает блоки 49 и 50

оперативной памяти, генератор 51 импульсов, блоки интегрирования 52, вычитания 53, реле 54 с контактами PI и Р2 и блок 55 умножения.

Прогнозатор работает следующим образом.

В момент пуска прогнозатора, совпадающего с начальным моментом f интервала планирования работы производственной системы, с ВЫХОДЕ- задатчика 1 на вход второго фильтра низкой частоты 3 и вход блока 10 начинает поступать сигнал x(t) заданной (плановой) траектории изменения первичной переменной. Одновременно на вход 34 поступает сигнал x(t) фактического значения

первичной переменной. Индентификатор 44, используя информацию x(t), x(t), формирует оценки коэффициентов модели задание - факт следующего вида; x°(t) K0(t) + k,x(t),

гдо xp(t) - расчетное значение первичной переменной; K0(t), ki(t) - оценки коэффициентов подели, пе-пумаемые на выходе идентификатора 44 в соответствии с алгоритмом

g.(i)-&(.))-(.V(

d + 5х (t )r

Ko(t) x(t) - Kix(t).

где и x(t) x(t) - x(t) - приращение фактического значения первичной переменной по отношению к ее сглаженному значению x(t);

6 x(t) - x(t) -x(t) - приращение заданного значения первичной переменной отно- сительно его сглаженного значения x(t);

(t) - сглаженное значение оценки коэффициента Ki(t) на момент времени t;

d - постоянный коэффициент алгоритма идентификации.

Сглаженные значения x(t), x(t) сигналов x(t), x(t) формируются на выходах первого и второго фильтров низкой частоты 2 и 3 соответственно. Приращения д x(t), d x(t) вычисляются блоками 9 и 10 соответственно. Сигнал 5x(t) поступает на первый вход блока 11, на второй вход которого подается сигнал fci(t) d x(t), формируемый на выходе блока 16 умножения. Сигнал разности д x(t) - fci(t)5 x(t) с блока 11 поступает на вход блока 30 деления. Сигнал д x(t) с выхода блока 10 подается на вход блока 16 умножения и на вход квадратора 29, который совместно с вторым сумматором 23 формирует сигнал {d + д x(t)}, поступающий на второй вход блока 30, причем константа d подается на вход сумматора 23 с входа устройства 35. На выходе блока 30 выдается сигнал отношения

(b(t)-ki(T)(b(t).

d+(5x(t)f

который складывается в первом сумматоре 22 с сигналом K(t), поступающим на вход сумматора 22 с выхода фильтра низкой час- тоты 4, на вход которого поступает результат суммирования, а именно сигнал оценки kt(t). Последний подается также на вход блока 17 умножения, на первый вход которого с выхода второго фильтра 3 низкой частоты поступает сигнал x(t).

Сигнал произведения Kix(t) с выхода фильтра 3 подается на вход блока 12. на другом входе которого имеется сигнал -x(t) с выхода первого фильтра 2 низкой частоты, На выходе блока 12 имеем сигнал Ј0(t) оценки коэффициента модели задание - факт. Сигналы k0(t), ki(t) подаются на соответствующие входы экстраполятора 45 (входы седьмого фильтра 8 низкой частоты и шестого блока 14 вычитания и входы пятого фильтра низкой частоты 6 и пятого блока 13 вычитания). Экстраполятор 45 реализует следую- щий лгоритм

K(t + в) а0 + Ь00;

Ki(t+ 6)-ai + bi0 , где Зо, ai, bo, bi - пэраметры экстраполятоРа 45;

K0(t+ в), Ki(t+ в)-экстраполированные

значения оценок коэффициентов Ko(t). Ki(t).

Параметры ai, bi экстраполятора 45 оцениваются с помощью двухслойной схемы фильтрации, включающей четвертый 5 и пятый 6 фильтры низкой частоты, соединенные между собой посредством блока 13. При этом на входы фильтра 6 и блока 13 подается сигнал Ki(t). С выхода фильтра 6 выдается сигнал а /, который подается на сумматор 24 и на вход блока 13, формирующего сигнал разности Ki(t) - ai, подаваемый на вход четвертого фильтра 5, на выходе которого выдается сигнал Ьь Аналогичным образом оцениваются параметры а0, bo экстраполятора. Соответствующая двухслойная схема фильтрации включает шестой 7 и седьмой 8 фильтры низкой частоты и блок 14. При этом на входы фильтра 8 и блока 14 подается сигнал Ko(t), а с выхода фильтра 8 выдается сигнал а , который подается на сумматор 25 и на вход блока 14, формирующего сигнал разности k0(t) - a0, подаваемый на вход шестого фильтра 7. на выходе которого выдается сигнал bo. На вторые входы третьего 18 и четвертого 19 блоков умножения с входа 36 подается сигнал интервала экстрапо- ляции в На выходе блоков 18 и 19 формируются сигналы произведений bi# , Ъ00 соответственно, которые подаются на входы сумматоров 24 и 25 соответственно. На выходах сумматоров 24 и 25 формируются сигналы экстраполированных значений коэффициентов K i(t + 0), K0(t + в) соответственно, которые подаются на входы блока 46 пересчетного прогнозирования блока 20умножения и вход сумматора 27. Блок 46 реализует следующий алгоритм:

х(т + в) K0(t + в) + Ki(t + 6). x(t + в);

x(t+ 6) x(t).

Сигнал x(t) подается на вход блока 20 с выхода второго фильтра 3, на выходе которого формируется сигнал K(t + x(t +0 ), который подается на вход сумматора 27, формирующего сигнал x(t + ф, выдаваемый на выход 41, а также на вход первого интегратора 32, являющийся одним из входов блока 47 расчета интегрального показателя. Сигнал момента времени t подается на другой вход интегратора 32 с входа 37. Блок 47 реализует алгоритм

к

Q(tK) x(t-r)dr + /x(t +0)6в

IH (t-tH), (tK-t),

где Q(t)-прогнозируемое на момент времени Значение интегрального показателя работы производственной системы; а- масштабный коэффициент;

x(t - т) - реализация входного сигнала x(t) устройства на отрезке времени t

x(t + в) - реализация прогнозируемого сигнала на интервале t, т.к .

На выходе интегратора 32 формируется

tkk сигнал / х (t + 0 ) d в который подается на

t

вход сумматора 26. На входы второго интегратора-33 подаются сигналы x(t) и t u с входов 34 и 38. На выходе интегратора 33

t формируется сигнал /x(t- r)dr котоt4

рый подается на вход сумматора 26. Выходной сигнал сумматора 26 подается на вход шестого блока 21 умножения, на другом входе которого установлен сигнал а подаваемый с входа 39. Результат умножения сигнал #(tK) подается на выход 42 устройства и на вход блока 48 расчета коэффициента интенсивности (вход блока 15). Блок 48 предназначен для определения таких значений коэффиента Ki(t) интенсивности на отрезке времени t, t , которые обеспечивали бы к моменту времени t равенство б(т Q(tK), где Q(tK) - заданное (плановое) значение интегрального показателя на конец интервала планирования. Блок 48 реализует следующий алгоритм

K-(t + 0)1.

Сигнал заданного значения Q(t ) подается с входа 40 на вход седьмого блока 15 вычитания и на блок 31 деления, С выхода блока 15 сигнал разности ДО Q(t ) - -Q(t ) подается на вход блока 31, на выходе которого формируется сигнал отношения AQ/Q(tк). Последний подается на вход сумматора 28, на другом входе которого имеется сигнал K I (t + в ), поданный с выхода сумматора 24. Результат суммирования с выхода сумматора 28 подается на выход 43.

Таким образом, на выходах устройства формируются три сигнала: сигнал прогнозируемого значения x(t+ в ) первичной произ- водственной переменной, сигнал прогнозируемого значения Q(t к) интегральной переменной на конец t интервала планирования, сигнал прогноза требуемого значения Ki (t + в ) коэффициента интенсивности работы производственной системы.

Интегратор 32 работает следующим образом.

Входной сигнал x(t + в ) поступает на вход блока 49 через контакт Pi, который замыкается с периодичностью, задаваемой генератором 51 импульсов. В моменты размыкания Р | на выходе блока 49 (в точке А) имеем постоянный сигнал U(A), который через замкнутый контакт Р2 поступает на первой вход блока умножения. На второй вход последнего с выхода блока вычитания по0 ступает сигнал разности tK - t, где tK - сигнал заданного момента окончания ин тервала прогнозирования, t - сигнал текущего момента времени, который формируется на выходе блока интегрирова5 ния прямоугольных импульсов. Генератор 51 вместе с блоком 52 играет роль датчика текущего времени, который включается в. момент г.уска устройства. Сигнал произведения U(A) - - 0 )tK - t, появляю0 на выходе блока умножения (в точке В), является оценкой интеграла

/x(t+$)d$. которая запоминается t

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

tK -t.

Временная диаграмма, поясняющая ра0 6oiy блоков интегратора 32 представлена на фиг 3. Точность оценивания интеграла в рассмотренной схеме определяется, главным образом, частотой сигнала прямоугольной волны, формируемого генератором 51.

5 Чем выше частота, тем точнее экстраполируется сигнал x(t + в ) и, как следствие, точнее оценивается искомый интеграл.

Таким образом, интегратор 32, содержащий описанные блоки, интегрирует прогно0 3HpveMbiu (экстраполируемый посредством простой сдвижки) сигнал x(t + в ). Формула изобретения Адаптивный прогнозатор производственных показателей, содержащий семь

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

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

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

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

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

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

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

производственной системы.

x(t+fl)

li-tl

Изобретение относится к автоматике и вычислительной технике и может быть использовано в системах управления производством. Цель изобретения - повышение точности прогнозирования и расширение функциональных возможностей прогноза- тора. Прогнозатор содержит задатчик 1 программного изменения первичного параметра, фильтры 2-8 низкой частоты, блоки 9-15 вычитания, блоки 16-21 умножения, сумматоры 22-28, квадратор 29. блоки 30 и 31 деления и интеграторы 32 и 33. За счет адаптации к ходу производственного процесса повышается точность прогнозирования его хода, включая и его обобщенные параметры.3 ил. 38 39 О Ј О ел 00 ел

Фиг.З