CC BY
151
УДК 681.516.5 Дата подачи статьи: 1.11.16
Б01: 10.15827/0236-235Х.030.1.045-050 2017. Т. 30. № 1. С. 45-50
СИСТЕМЫ АВТОМАТИЧЕСКОГО УПРАВЛЕНИЯ ОБЪЕКТАМИ С ЗАПАЗДЫВАНИЕМ: РОБАСТНОСТЬ, БЫСТРОДЕЙСТВИЕ, СИНТЕЗ
В.З. Тхан, аспирант, [email protected]; Д.Ю. Берчук, аспирант, [email protected] (Национальный исследовательский Томский политехнический университет, просп. Ленина, 30, г. Томск, 634050, Россия)
Рассматривается задача синтеза систем автоматического управления объектами с запаздыванием. Особенность ее постановки связана с формированием и решением уравнения синтеза численным методом без аппроксимации передаточной функции звена запаздывания. Этим создаются потенциальные возможности для повышения точности синтеза регуляторов. В работе они раскрываются за счет привлечения численного подхода к синтезу регуляторов - вещественного интерполяционного метода. Метод позволяет найти стабилизирующий регулятор по модели желаемой системы автоматического управления, известному неаппроксимированному математическому описанию объекта управления и принятой структуре регулятора. Вторая особенность задачи также связана с принятым методом и звеном запаздывания и реализуется на этапе получения желаемой модели.
Для усиления демонстрационных возможностей и наглядности подхода в работе использован минимальный набор показателей - перерегулирование и время установления переходного процесса. Модель желаемой системы совместно с неаппроксимированным описанием объекта управления и принятой структурой регулятора позволили сформировать более точное уравнение синтеза по сравнению с традиционным подходом, основанным на приближении модели звена запаздывания дробно-рациональным выражением.
Рассмотрены численный способ решения уравнения синтеза, а также исследование процедуры синтеза и свойств робастности систем с запаздыванием. На этой основе предложен алгоритм достижения максимального быстродействия систем автоматического управления с запаздыванием в условиях поддержания перерегулирования в заданных пределах и показан механизм ухудшения робастности синтезированной системы при увеличении запаздывания сигнала в объекте управления даже в условиях подстройки параметров регулятора под изменяющееся время запаздывания.
Результаты проиллюстрированы расчетными примерами.
Ключевые слова: синтез систем управления с запаздыванием, робастность, быстродействие.
Системы автоматического управления (САУ) объектами с запаздыванием составляют обширный и важный для практики класс динамических систем. Расчет таких систем имеет специфические особенности в силу двух обстоятельств. Во-первых, это физические особенности, когда реакция объекта задерживается на время т относительно момента подачи сигнала на вход объекта. Во-вторых, это же явление вносит свои особенности в расчетные процедуры САУ из-за наличия в их математическом описании трансцендентной составляющей e -рт, отражающей существование запаздывания т. С позиций синтеза систем с запаздыванием обе особенности учитываются при построении регуляторов на основе популярного предиктора Смита [1]. Подход эффективен, но имеет свои ограничения [2, 3]. В частности, в условиях переменного запаздывания система управления может значительно изменять свои свойства вплоть до потери устойчивости. К этому можно добавить, что параметры дробно-рациональной части объекта практически всегда также изменяются с течением времени, что ухудшает качество управления. В настоящей работе делается попытка поиска средств, позволяющих раскрывать робастные свойства САУ с запаздыванием на этапе синтеза регулятора и представлять их в наглядной форме. Аналитическое решение такой задачи, очевидно, не-
возможно из-за трудностей действий с функций e -рт, входящей в модель САУ. По этой причине в работе используется один из численных методов, позволяющих оперировать описаниями САУ с запаздыванием.
Выбор численного метода расчета систем с запаздыванием
При решении задач автоматического управления основным математическим аппаратом является преобразование Лапласа, которое позволяет заменить операции в области времени более простыми действиями в области изображений. Получение таких изображений осуществляется по формуле прямого преобразования
Р(Р) = Р = 5+У® • (1)
Она удобна для аналитических преобразований, однако переход к численным действиям в общем случае затруднителен из-за вычислительных сложностей. Поэтому привлекают частные случаи преобразования (1), вытекающие из двух возможных трансформаций комплексной переменной р=5+/ю. Первый вариант определяется условием 5=0, при котором комплексная переменная вырождается в чисто мнимую р=/®. Получаемые изображения F(j®) соответствуют частотному методу на основе
преобразования Фурье, они легко преобразуются в численные модели вида F{jЮi), /=1, 2, ... [4]. Второй вариант вытекает из преобразования (1) при условии ю=0, когда комплексная переменная вырождается в вещественную p=5 [5, 6]. Получаемые функции-изображения
F(5) = J7 f (t)e-"a'dt, 5 e [С, да], С > 0:
(2)
являются вещественными функциями вещественной переменной, что важно с позиций объема вычислений и практических задач управления. Функции F(5) позволяют достаточно просто перейти к дискретным представлениям F(5i), /=1, 2, ..., привлечь для действий над ними хорошо разработанные численные методы. В отношении САУ с запаздыванием, которые являются предметом рассмотрения в настоящей работе, отметим: модели вида F(5i), /=1, 2, ..., встречают существенно меньшие трудности при их преобразовании по сравнению с аналитическом подходом, так как действия осуществляются в вещественной арифметике - над числами e -5и, /=1, 2, ...,
В формуле вещественного преобразования (2) предусмотрено ограничение переменной 5 снизу величиной 00, что обеспечивает сходимость интеграла для широкого класса функций-оригиналов ^), охватывающих практически все важные виды сигналов динамических систем.
Синтез регуляторов САУ с запаздыванием численным методом
Выбранный численный метод позволяет находить регуляторы автоматических систем, управляющих так называемыми особыми объектами, в том числе с запаздыванием, с распределенными параметрами [5, 6]. Воспользуемся такой возможностью для рассмотрения быстродействия САУ с запаздыванием и их робастных свойств. С этой целью рассмотрим конкретные задачи, в рамках которых можно выполнять численные расчеты, анализировать результаты и делать обобщающие выводы.
В качестве примера управления использовались экструдеры, предназначенные для изготовления пластиковой нити для 3Б-принтеров [7], а также близкие к ним объекты [8], управление которыми сводится к типовой схеме, представленной на рисунке 1.
Система содержит регулятор, объект управления и цепь обратной связи. Для объекта [8] элементы САУ представлены описаниями
Жн (Р) = ■ 0,5
p
4 p 2 + 8 p +1
р (p> p(a,p +1)
K, = 1.
По технологии требуется поддерживать перерегулирование ст в заданных пределах:
стз - Дст < CT < СТз + Дст.
(3)
Обозначения стз, Дст, ст в (3) соответствуют заданному значению перерегулирования Стз =5 %, достигнутой величине ст и допустимому отклонению Дст = [стз - ст|, Дст = +1 %. Еще одно требование -время установления переходного процесса должно быть минимальным:
i ^ min. (4)
В этих условиях требуется синтезировать регулятор, не прибегая к аппроксимации звена запаздывания, рассмотреть робастные свойства системы, а также выяснить зависимость этих свойств от величины запаздывания.
Для решения задачи составим уравнение син-
Wp (pW„ (р) . Входящая в левую
теза
w (p) s
1 + W. ( pWo ( p)kK
часть уравнения желаемая передаточная функция Ж* (р) найдется по исходным данным. Неопределенность возникает лишь при задании желаемого времени установления ^ , приближенным ориентиром для которого являются коэффициенты знаменателя модели объекта. В условиях численного метода и итерационного приближения к искомому решению неточность такого рода не имеет существенного значения. Более того, как будет ясно из дальнейшего, эта неоднозначность и итерационная процедура позволяют получить множество решений, из которых можно выбрать наиболее подходящее, привлекая дополнительные критерии. Желаемое время переходного процесса ^ будем рассматривать в качестве переменного параметра, обеспечивая его минимизацию при выполнении других условий. Для начальных расчетов оно было принято равным 115 с, что позволило получить методом [5] желаемую передаточную функцию:
Ж(р)= +1 .
ж 794р2 + 42 р +1)
В правой части уравнения синтеза содержится модель объекта управления, которая известна, и передаточная функция регулятора, вид которой выбран на основе априорных данных:
W р ( p ) =
b1 p + b0
Р(а1 Р +1)
Для упрощения уравнения синтеза, сопоставляющего желаемую и синтезируемую системы в за-
мкнутом состоянии, переидем к равносильному уравнению = Жр(р)-Ж„(р). описывающему
эти же системы в разомкнутом состоянии. В левой части уравнение содержит передаточную функцию
р, , р + 0,048 разомкнутой системы Жр (р) =-, кор (37,81р +1)
торая находится однозначно по уже известной функции Жжз (р).
В результате получаем расчетное соотношение изображений по Лапласу: р + 0,048
p(37,81p +1)
-, biP + bo
0,5
(5)
,-3 p
р(а1 р +1) 4р2 + 8р +1
Дальнейшие шаги могут быть сделаны в рамках известных инженерных методов. Они связаны с заменой передаточной функции звена запаздывания на дробно-рациональное выражение, адаптированное к классическим методам. Другой вариант -использование принятого численного метода. В последнем случае появляется возможность исключения аппроксимации звена запаздывания, что позволяет надеяться на повышение точности решения уравнения. Для его реализации перейдем в уравнении (5) к вещественным изображениям. Такой переход для устойчивых звеньев и систем можно сделать формально заменой переменных: р^-5 [6]. Очевидно, что все полюсы функций в соотношении (5) расположены в левой полуплоскости, так что можно воспользоваться формальным приемом, что приводит к форме 5 + 0,048
8(37,818 +1)
_ Ь 8 + b
0,5
(6)
5(а5 +1) 452 + 85 +1
Для численного решения уравнения (6) заменим входящие в него функции их численными представлениями вида F(5I■), /=1, 2, 3, в соответствии с рекомендациями [6]. В результате получим систему линейных алгебраических уравнений (СЛАУ):
5 + 0,048
8, (37,818, +1)
~ Щ + bo
0,5
(7)
-• e~3S', i = 1, 2,3.
5, Ц5, +1) 45,2 + 85,.+1
При формировании СЛАУ для упрощения задачи узлы интерполирования приняты равномерными: 5/ = ■51, /=1, 2, 3. В этом случае первый узел 51 можно использовать в качестве инструментальной переменной, что позволяет из множества решений выбрать такое, которое в наибольшей степени отвечает предъявляемым требованиям, в частности по критерию (3).
По итогам принятых уточнений констатируем: имеются заданное значение перерегулирования стз, условие (4) и два переменных параметра - ^ и 51,
с помощью которых необходимо достичь заданных показателей качества САУ. Для выстраивания приоритетов будем считать перерегулирование ведущим параметром, учитывая его интегральный смысл, так как он отражает не только качество работы САУ, но и косвенно запас устойчивости. Такой подход известен и успешно используется [9]. Для объединения трех выделенных параметров воспользуемся формой ст=/(51) для различных ^ •
Это позволит оптимизировать быстродействие синтезируемых систем, оценить их робастность [10], а также влияние величины запаздывания объекта управления на эти свойства САУ.
Быстродействие САУ с запаздыванием
Быстродействие создаваемой системы при прочих равных условиях будет зависеть от желаемого времени установления переходного процесса, задаваемого при формировании эталонной передаточной функции. На рисунке 2 приведены результаты синтеза регуляторов для нескольких значений ^ .
Итоги расчетов даны в виде графиков ст=/(51) для нескольких значений желаемого времени установления: = 100 ^ /уж2 = 115 ^ /уж3 = 130 c. Видно, что положение кривой ст=/(51) зависит от задаваемой величины ^ • Изменяя ее, можно найти такое
положение графика ст=/(51), при котором перерегулирование ст соответствует заданному значению стз с заранее определенной точностью. Таким образом, на этапе синтеза регулятора можно обеспечить выполнение условий (1) и (2). Для принятого расчетного варианта, как показывает рисунок 2, приемлемой оказывается итерация, соответствующая времени = 115 а
-38
Обратим внимание на важную для практики возможность синтезировать регуляторы, которые ограничивают перерегулирование не только сверху, как это обычно бывает, но и снизу. Подчеркнем, что речь идет не о простейших САУ, а о системах с запаздыванием.
Робастность САУ с запаздыванием
Робастное свойство САУ многогранно, поэтому обычно ограничиваются какими-либо конкретными показателями системы. В данном случае за такой показатель вновь примем перерегулирование в силу его обобщающих свойств по отношению к показателям устойчивости и качества САУ. В [10] показано, что при определенных условиях изменение величины 51 приближенно соответствует изменению значений физических параметров САУ и, следовательно, коэффициентов ее передаточной функции. Это позволяет оценивать робастность системы по функции ст=Д51) - величине той ее части, которая удовлетворяет условию (3). Для расчетов и сравнительного анализа вариантов используют количественную оценку робастности, которая определяется максимальной величиной интервала [51тт, 51тах], внутри которого выполняется условие (3). Обозначим такую оценку как
Гст= 52тт - 51тах, (8)
указывая индексом ст на отношение показателя ро-бастности к перерегулированию.
Воспользуемся (8) для исследования влияния величины запаздывания на робастные свойства САУ. С этой целью сохраним неизменной дробно-рациональную часть модели объекта управления и структуру регулятора, а величину запаздывания будем изменять, оценивая робастность синтезируемых САУ по величине Гст. Для оценивания итогов вновь обратимся к наглядной зависимости ст=Д51). Результаты расчетов для нескольких величин запаздывания представлены на рисунке 3.
Рис. 3. График зависимости показателя робастности Га от величины запаздывания
Fig. 3. A graph of dependence of a robustness indicator on a delay value
График показывает, что увеличение запаздывания существенно ухудшает робастные свойства САУ по перерегулированию. При этом подчеркнем, что график получен в условиях настройки регулятора САУ на каждое новое изменение времени запаздывания. В тех случаях, когда такая подстройка не выполняется, робастность снижается намного сильнее.
Есть основания предполагать, что показатели качества и устойчивости имеют примерно такую же зависимость от величины запаздывания. Представленные сведения свидетельствуют о том, что в рамках рассматриваемого численного подхода имеется возможность оценивать влияние изменения запаздывания на свойства САУ, в частности на перерегулирование. С этих позиций можно дать практическую рекомендацию по выбору величины переменной 51: она должна находиться в середине интервала [51тт, 51тах]. Возможности расширения робастных свойств САУ увеличиваются, если перерегулирование ограничивается только сверху, что имеет место во многих практических задачах. При этом достаточно сместить рабочую точку 51 в сторону начала интервала [51тт, 51тах].
Повышение точности безаппроксимационного решения уравнения синтеза САУ с запаздыванием
Использование неаппроксимированного уравнения синтеза САУ с запаздыванием для повышения точности таких систем привлекательно и в каких-то случаях может оказаться эффективным. Проверим эту возможность, используя прежний объект управления. Для сравнения с традиционным подходом, основанным на приближении звена запаздывания рациональной дробью, выберем два достаточно сильно различающихся по точности приближения звена запаздывания е хр:
К
Жр, (Р) = —0^ , (9)
ap+1
wnp2 ( p) = ■
Kp + bo
2 , (10)
а2р + аР + 1
Вычислим неизвестные коэффициенты численным методом [6]. Для структур (9) и (10) получены две модели объекта управления:
0,5 1
WH1 (5) =
452 + 85 +1 4,865p + 1
Wh2 (5) = .
0,5
-0.752p +1
452 + 85 +1 2,816р2 + 2,177р +1
Для получения представления о точности приближения воспользуемся наглядным представлением ошибки в области вещественных изображений как функции А №(5) = е -рх - №пр(5), где №пр(5) -аппроксимированная вещественная передаточная
функция звена запаздывания. Результаты аппроксимации представлены на рисунке 4.
Рис. 4. Графики ошибок аппроксимации AW1(5) = e -35 - Wnp1(S), AW2(5) = e -35 - W^)
Fig. 4. The approximation errors graphs AW1(8) = e -35 - Wnp1(8), AW2(S) = e -35 - W^)
Сопоставляя графики, можно констатировать, что погрешность аппроксимации в случае функции второго порядка W2(p) значительно меньше по сравнению с максимально простой функцией первого порядка Wi(p). Об этом можно судить не только по графикам, но и по количественным оценкам, например, вида AW = max | AW(5) |.
Таким образом, имеются три модели объекта управления: точная и две приближенные, у одной из которых сравнительно высокая точность. Наличие трех различных по точности моделей позволит установить влияние погрешности представления модели на результаты синтеза САУ. Для оценивания результатов были получены функции ст=/(51), соответствующие принятым моделям объектов. Они отвечают требованию по перерегулированию (3), которое является обязательным требованием синтеза САУ. Выяснилось, что вопреки ожиданиям робастность по перерегулированию для всех трех случаев примерно одинакова. Заметно худшие результаты получены для простейшей модели (8) в том смысле, что робастность оказалась ошибочно завышенной за счет увеличения верхней границы 5max.
В то же время быстродействие САУ оказалось существенно зависимым от точности модели объекта. С понижением точности модели время установления ty синтезированной САУ увеличивается. Для принятых моделей оно составляет соответственно 149, 151 и 180 с. Это время отсчитывалось для среднего значения инструментальной переменной 51 = (51max - S1mm)/2, когда робастность каждой системы максимальна. Отметим, что при снижении требований по робастности можно уменьшить время установления ty, опираясь на функцию ст=Д51).
Заключение
В работе рассмотрены вопросы синтеза регуляторов САУ объектами, имеющими запаздывание, направленные на повышение робастных свойств этих систем. Для проведения исследований был выбран численный путь, позволяющий снизить трудности, связанные с трансцендентной формой модели звена запаздывания. Этот выбор позволил
- составить и решить уравнение синтеза, в котором модель звена запаздывания не требует аппроксимации, что потенциально обеспечивает повышенную точность синтеза САУ;
- использовать возможность синтеза САУ, в которых перерегулирование находится в заданных пределах, а быстродействие близко к максимальному;
- установить наглядную зависимость, отражающую ухудшение робастности САУ при увеличении запаздывания.
Подход может быть реализован не только в рамках синтеза САУ с запаздыванием, но и при построении самонастраивающихся САУ, в которых изменяется как запаздывание, так и другие параметры объекта управления. Для систем этого класса имеются положительные предпосылки [12].
Литература
1. Смит О.Дж.М. Автоматическое регулирование; [пер. с англ. под ред. Е.П. Попова]. М.: Физматгиз, 1963. 847 с.
2. Филимонов А.С. Спектральная декомпозиция динамических систем с запаздываниями: Теория и применения: дисс. д. технич. наук. М.: РУДН, 2003. 426 с.
3. Громов Ю.Ю., Земской Н.А., Лагутин А.В., Иванова О.Г., Тютюнник В.М. Системы автоматического управления с запаздыванием. Тамбов: Изд-во ТГТУ, 2007. 76 с.
4. Козлов О.С., Скворцов Л.М. Решение задач теории автоматического управления в программном комплексе МВТУ // Мехатроника, автоматика, управление. 2016. Т. 17. № 1. С. 3-10.
5. Анализ и оптимальный синтез на ЭВМ систем управления; [под ред. А.А. Воронова и И.А. Орурка]. М.: Наука, 1984. 344 с.
6. Гончаров В.И. Вещественный интерполяционный метод синтеза систем автоматического управления. Томск: Изд-во ТПУ, 1995. 100 с.
7. Тхан Вьет Зунг, Берчук Д.Ю. Экспериментальное определение передаточной функции нагревателя экструдера // Молодежь и современные информационные технологии: сб. тр. XIII Междунар. науч.-практич. конф. Томск: Изд-во ТПУ, 2016. С. 261-262.
8. Михалевич С.С., Байдали С.А., Чурсин Ю.А. Графический метод расчета параметров настройки регуляторов, приводящих систему к заданным показателям устойчивости // Промышленные АСУ и контроллеры. 2013. N° 10. С. 35-39.
9. Гончаров В.И., Лиепиньш А.В., Рудницкий В.А. Синтез робастных регуляторов низкого порядка // Изв. РАН. Теория и системы управления. 2001. № 4. С. 36-43.
10. Крутова И.Н. Формирование алгоритма управления итерационным процессом настройки параметров в системе с упрощенной эталонной моделью // Автоматика и телемеханика. 1998. № 2. С. 72-84.
11. Денисенко В.В. Компьютерное управление технологическим процессом, экспериментом, оборудованием. М.: Горячая линия-Телеком, 2009. 608 с.
12. Goncharov V. Rudnicki V. Real interpolation method in automatic control systems self-adjustment problem // Systems Science, vol. 36, no. 3, 2010, pp. 35-37.
Software & Systems
DOI: 10.15827/0236-235X.030.1.045-050
Received 1.11.16 2017, vol. 30, no. 1, pp. 45-50
TIME DELAY AUTOMATIC CONTROL SYSTEMS: ROBUSTNESS, RESPONSE TIME, SYNTHESIS
V.D. Than1, Postgraduate Student, [email protected] D.Yu. Berchuk1, Postgraduate Student, [email protected]
1 National Research Tomsk Polytechnic University, Lenin Ave. 30, Tomsk, 634050, Russian Federation
Abstract. The article considers synthesis of time delay automatic control systems (ACS). A synthesis statement is connected to equation solving by a numerical technique without approximation of a time delay element. This issue has a potential to improve synthesis accuracy of regulators. A real interpolation method (RIM) has number of advantages in synthesis of time delay ACS. RIM allows finding a regulator that compensates time delay using a targeted model, a defined non-approximated mathematical model of object and a regulator structure. Moreover, RIM allows building transfer function with direct quality factors on the stage of developing a targeted model. It's simpler and more useful for a developer.
The article demonstrates potentials of the method using minimal quantity of factors: overshoot and settings time. The targeted model with a defined non-approximated mathematical model and a regulator structure allows building accurate synthesis equation, which is useful for numerical research of synthesis and robust properties of time delay systems. The article proposes an algorithm of a maximal control rate for a time delay system in terms of specified overshoot. It also shows a negative effect of increasing time delay to robustness of the system. A result is graphically illustrated and proved by a number of examples.
Keywords: synthesis of time delay control systems, robustness, performance.
1. Smith O.J. Automatic control (Russ.ed.: Popov E.P. Moscow, Fizmatgiz Publ., 1963, 847 p.).
2. Filimonov A.S. Spektralnaya dekompozitsiya dinamicheskikh sistem s zapazdyvaniyami: teoriya i primeneniya [Spectral Decomposition of dynamical systems with delays: Theory and applications]. Dr.Sc. thesis. Moscow, PFUR Publ., 2003.
3. Gromov Yu.Yu., Zemskoy N.A., Lagutin A.V., Ivanova O.G., Tyutyunnik V.M. Sistemy avtomaticheskogo up-ravleniya s zapazdyvaniem [Automatic control systems with delay]. Tambov, TSTU Publ., 2007, 76 p.
4. Kozlov O.S., Skvortsov L.M. Solving problems in the theory of automatic control software package MVTU.Me-khatronika, avtomatika, upravlenie [Mechatronics, Automation, Control]. 2016, vol. 17, no. 1, pp. 3-10 (in Russ.).
5. Analiz i optimalny sintez na EVM sistem upravleniya [Analysis and Optimum Synthesis of Computer Control Systems]. A.A. Voronov, I.A. Orurk (Eds.). Moscow, Science, 1984, 344 p.
6. Goncharov V.I. Veshchestvenny interpolyatsionny metodsinteza sistem avtomaticheskogo upravleniya [Synthesis of electromechanical actuators for industrial robots]. Tomsk, TPU Publ., 2002, 100 p.
7. Than V.Z., Berchuk D.Yu. Experimental determination of the transfer function of the extruder heater. Youth and modern IT: Proc. 13 Int. Science and Practice Conf. of Students, Postgraduates and Young Scientists. Tomsk, 2015, pp. 261-262 (in Russ.).
8. Mikhalevich S.S., Baydali S.A., Chursin Yu.A. Graphical method for calculating regulatory settings that specifies the system to the parameters of stability. Promyshlennye ASUi kontrollery [Industrial ACS and controllers]. 2013, no. 10, pp. 35-39 (in Russ.).
9. Goncharov V.I., Liepinsh A.V., Rudnitsky V.A. The synthesis of low order robust controllers. Izvestiya Ross. akademii nauk. Teoriya i sistemy upravleniya [Jour. of Computer and Systems Sciences International]. 2001, no. 4, pp. 36-43.
10. Krutova I.N. Formation of a control algorithm for setting up iterative process in a system with a simplified reference model. Avtomatika i telemekhanika [Automation and Remote Control]. 1998, no. 2, pp. 72-84 (in Russ.).
11. Denisenko V.V. Kompyuternoe upravlenie tekhnologicheskim protsessom, eksperimentom, oborudovaniem [Computer control of a production process, experiment equipment]. Moscow, Goryachaya Liniya Telecom Publ., 2009, 608 p.
12. Goncharov V., Rudnicki V. Real interpolation method in automatic control systems self-adjustment problem. Systems Science. 2010, vol. 36, no. 3, pp. 35-37.
References
