В. В. Куликов, соискатель
664074, Иркутск
Kulikov V. V., Postgraduate Student
Irkutsk, 664074
Н. Н. Куцый, д-р техн. наук, проф.
664074, Иркутск
Irkutsk, 664074
Е. А. Осипова, канд. техн. наук, доц.
664074, Иркутск
Irkutsk, 664074
В автоматических системах регулирования (АСР) промышленных процессов различного типа с запаздыванием при ограничении на значение перерегулирования широко применяется ПИД регулятор с реальным дифференцирующим звеном (далее ПИД регулятор). Как известно, признаком наличия большого запаздывания у объекта регулирования является отношение τob/Tob l 1, где τob, Tob — время запаздывания и максимальная постоянная времени объекта регулирования. При наличии в АСР большого запаздывания и ограничения на выходную координату параметрический синтез ПИД регулятора хорошо известными частотными методами становится затруднительным, что обусловливает интерес к разработке численных беспоисковых алгоритмов параметрической оптимизации, основанных на использовании функций чувствительности для определения градиента критерия оптимальности.
В данной работе предложен алгоритм численной оптимизации, вычисляющий исходя из минимума интегрального квадратичного критерия значения настраиваемых параметров ПИД регулятора в указанных АСР. Чтобы обеспечить отсутствие перерегулирования у результирующего переходного процесса, предлагается на этапе оптимизации вводить в АСР ограничение на регулирующее воздействие, которое, в свою очередь, учитывается путем добавления штрафной функции в интегральный критерий.
Предлагаемый алгоритм основан на методе сопряженных градиентов Полака—Поляка—Рибьера с его известными преимуществами. Составляющие градиента критерия оптимизации вычисляются с помощью таких функций чувствительности, которые позволяют получить все компоненты этого вектора без пробных поисковых вариаций настраиваемых параметров. Для вычисления шага оптимизации авторами реализован соответствующий алгоритм, в основе которого лежит градиентная процедура с применением функции чувствительности выходной координаты АСР к значению шага.
Сходимость предлагаемого алгоритма проверена с помощью вычислительной процедуры на основе зоны притяжения рекордных значений критерия оптимизации, которая определяется положительно определенной матрицей Гессе интегрального квадратичного критерия, построенного на разности усредненного и проверяемого переходных процессов.
In automatic control systems (ASR), industrial processes of various types with a delay with a limit on the amount of overregulation, a PID controller with a real differentiating link (hereinafter referred to as the PID controller) is widely used. As is known, a sign of the presence of a large delay in the object of regulation is the ratio τob/Tob ≥ 1, where τob is the value of the delay time and Tob is maximum time constant of the object of control. In the presence of a large delay and limitation in the ASR, the parametric synthesis of the PID controller by well-known frequency methods becomes difficult, which leads to interest in the development of numerical searchless algorithms for parametric optimization based on the use of sensitivity functions to determine the gradient of the optimality criterion. In this paper, an APO algorithm is formed that calculates, based on the minimum of the integral quadratic criterion, the values of the adjustable parameters of the PID controller in the specified ASR. In order to ensure that there is no re-regulation of the resulting transient process, the authors of this article propose to introduce a restriction on the regulatory effect into the automatic system at the optimization stage, which, in turn, is taken into account by adding a penalty function to the integral criterion. The proposed algorithm is based on the method of conjugate Polak—Polyak-Ribier gradients with its known advantages. The components of the gradient vector of the optimization criterion are calculated using such sensitivity functions that allow you to obtain all the components of this vector without trial search variations of the configurable parameters. To calculate the value of the optimization step, the authors implemented an appropriate algorithm based on a gradient procedure using the sensitivity function of the output coordinate of the ASR to the step value. The convergence of the generated APO algorithm was verified using a numerical procedure based on the zone of attraction of record values of the optimization criterion, which is determined by a positive-definite Hesse matrix of the integral quadratic criterion based on the difference between the averaged and the tested transients.
The authors declare that there are no conflicts of interest present.