KR100310609B1 - Coefficient adjusting apparatus of proportional integral derivative(pid) controller - Google Patents

Coefficient adjusting apparatus of proportional integral derivative(pid) controller Download PDF

Info

Publication number
KR100310609B1
KR100310609B1 KR1019970020917A KR19970020917A KR100310609B1 KR 100310609 B1 KR100310609 B1 KR 100310609B1 KR 1019970020917 A KR1019970020917 A KR 1019970020917A KR 19970020917 A KR19970020917 A KR 19970020917A KR 100310609 B1 KR100310609 B1 KR 100310609B1
Authority
KR
South Korea
Prior art keywords
unit
value
control
coefficient
fuzzy
Prior art date
Application number
KR1019970020917A
Other languages
Korean (ko)
Other versions
KR19980084979A (en
Inventor
이봉국
김종환
김광춘
Original Assignee
이종수
엘지산전 주식회사
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 이종수, 엘지산전 주식회사 filed Critical 이종수
Priority to KR1019970020917A priority Critical patent/KR100310609B1/en
Publication of KR19980084979A publication Critical patent/KR19980084979A/en
Application granted granted Critical
Publication of KR100310609B1 publication Critical patent/KR100310609B1/en

Links

Images

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B11/00Automatic controllers
    • G05B11/01Automatic controllers electric
    • G05B11/36Automatic controllers electric with provision for obtaining particular characteristics, e.g. proportional, integral, differential
    • G05B11/42Automatic controllers electric with provision for obtaining particular characteristics, e.g. proportional, integral, differential for obtaining a characteristic which is both proportional and time-dependent, e.g. P.I., P.I.D.
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B13/00Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
    • G05B13/02Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
    • G05B13/0265Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric the criterion being a learning criterion
    • G05B13/027Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric the criterion being a learning criterion using neural networks only
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B13/00Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
    • G05B13/02Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
    • G05B13/04Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
    • G05B13/042Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance

Abstract

PURPOSE: A coefficient adjusting apparatus of proportional integral derivative(pid) controller is provided to automatically adjust a PID coefficient until rising time, overshot, and correction time as characteristic values of a system output indicating a system attribute arrives at the system output of a desired performance index. CONSTITUTION: A process monitoring station part(100) comprises an evolution calculation part(131), a management part(110), and a fuzzy professional device(120). The evolution calculation part(131) classifies a process to search an optimum regular base with regard to each classified process. The management part(110) interfaces with an outer system or a driver, and generates a control target set value of a control object. The management part(110) displays and manages the set value, a control value, and a process value. The fuzzy professional device(120) performs according to reasoning of a professional technique using a fuzzy rule when satisfying system response specifications of requiring change of a PID coefficient. A process control station part(200) directly collects data from the process via a field bus and controls the collected data. A process part(300) receives a control output value from a constant sampling period in the process control station part(200) and outputs a process value correspond to the received value.

Description

퍼지 전문가 장치를 이용한 피아이디 제어기의 계수조정 장치 및 그의 제어방법Coefficient Adjustment Device of PID Controller Using Fuzzy Expert Device and Its Control Method

본 발명은 퍼지 전문가 장치를 이용한 피아이디 제어기의 계수조정 장치 및 그의 제어방법에 관한 것으로, 특히 플랜트 제어 시스템에 적용되는 분산형 공정제어 시스템의 제어기법을 향상시키기 위해 전문가 기법을 도입하여 시스템의 운용을 보다 적응적으로 수행하는데 적당하도록 한 퍼지 전문가 장치를 이용한 피아이디 제어기의 계수조정 장치 및 그의 제어방법에 관한 것이다.BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to a coefficient adjusting device of a PID controller using a fuzzy expert device and a control method thereof. In particular, the system is operated by introducing expert techniques to improve the control method of a distributed process control system applied to a plant control system. The present invention relates to a coefficient adjusting device of a PID controller using a fuzzy expert device and to a method of controlling the same.

피아이디(PID : Proportional Integral Derivative) 제어기는 구조의 단순함에 비해 매우 안정하고, 견실한 특성을 가지고 있기 때문에 대부분의 산업 응용 부분에 많이 쓰이는데 일반적인 피아이디(PID) 제어기의 형태는 다음과 같다.Proportional Integral Derivative (PID) controller is used for most industrial applications because it is very stable and robust compared to the simplicity of structure. The form of general PID controller is as follows.

Figure kpo00002
-------- (Kc: 비례이득, Ti: 적분이득의 역수, Td: 미분이득)
Figure kpo00002
-------- (K c : proportional gain, T i : inverse of integral gain, T d : differential gain)

상기 피아이디(PID) 제어기는 분산제어 시스템에서 실질적인 제어를 담당하고 있는 제어 스테이션의 제어부에서 동작함으로써 제어치를 결정하며, 상기 피아이디(PID) 제어기의 여러 파라미터(비례, 적분, 미분계수) 설정은 분산제어 시스템의 감시부에 해당하는 감시 스테이션에서 수행한다.The PID controller determines a control value by operating in a control unit of a control station in charge of substantial control in a distributed control system, and setting of various parameters (proportional, integral, derivative) of the PID controller is performed. It is performed by the monitoring station corresponding to the monitoring part of the distributed control system.

도1은 일실시예로 종래 분산제어 시스템의 블록 구성도로서, 이에 도시된 바와 같이 분산제어 시스템의 상위 계층 스테이션으로 운전자의 휴먼인터페이스(Human Interface)를 담당하여 전체 공정의 동작상태를 모니터 화면을 통해 감시하는 기능과, 데이터의 저장 및 관리를 담당하는 제1∼N 피엠에스(PMS : Process Monitoring Station)(10A, 10B)부와; 공정으로 부터 필드 버스(Field Bus)를 통해 직접적인 데이터 수집과 제어를 수행하는 하위 계층 스테이션인 제1∼N 피씨에스(PCS : Process Control Staation)(30A ~ 30N)부와; 고급제어 알고리즘 및 특수한 기능을 갖는 알고리즘을 포함하고 있는 상위 계층의 에이씨에스(ACS : Advanced Control Station)(10C)로 구성된 것으로, 이와 같이 구성된 종래의 동작과정을 첨부한 도2를 참조하여 설명한다.FIG. 1 is a block diagram of a conventional distributed control system according to an embodiment. As shown in FIG. 1, the upper layer station of the distributed control system is in charge of the driver's human interface to monitor the operation state of the entire process. A first through NMS (PMS: Process Monitoring Station) (10A, 10B) which is in charge of monitoring and storing and managing data; First to NPC (PCS: Process Control Staation) units 30A to 30N which perform lower data collection and control directly through a field bus from a process; It is composed of an Advanced Control Station (ACS) 10C of an upper layer including an advanced control algorithm and an algorithm having a special function, which will be described with reference to FIG.

도2는 종래 피아이디 기법을 적용한 피엠에스/피씨에스부의 블록 구성도로서, 이에 도시된 바와 같이 제1∼N 피씨에스부(30A ~ 30N)에서 설정치(SV : Set Value)와, 공정치(PV : Process Value)가 가산기(31)에 공급되어 이들의 차이 즉, 오차(error)가 피아이디(PID) 제어기(32)에 공급되어 이로부터 일정한 샘플링 주기에서 제어 출력값(MV : Manipulated Value)이 공정부(40)와 공정관리 및 모니터부(12)로 출력되고, 제1∼N 피엠에스부(10A, 10B)는 상기 공정관리 및 모니터부(12)에 의해 제어 대상의 설정치와 제어치 및 공정치가 표시 및 관리되고, 설정치 발생부(11)에 의해 외부 시스템으로부터 간접적으로 발생한 신호가 선택되거나 운전자에 의해 데이터가 선택되게 함으로써, 운전자와의 인터페이스가 가능하게 된다.FIG. 2 is a block diagram of a PS / PCS unit to which the conventional PD technique is applied. As shown in FIG. 2, a set value (SV) and a process value (PV) of the first to NPCS units 30A to 30N are shown. : Process Value) is supplied to the adder 31, and the difference thereof, that is, error is supplied to the PID controller 32, from which the control output value (MV) is supplied in a constant sampling period. Outputted to the government unit 40 and the process management and monitor unit 12, and the first to N PMS units 10A and 10B are controlled by the process management and monitor unit 12, and control values and control values and process values. Is displayed and managed and the signal generated indirectly from the external system is selected by the set value generating unit 11 or data is selected by the driver, thereby enabling the interface with the driver.

공정에 대하여 피아이디 제어기를 자동으로 튜닝하기 위하여 오프라인 튜닝과 온라인 튜닝 방법이 사용되는데, 오프라인 튜닝 방법은 초기에 결정된 피아이디 제어기 계수가 만족할만한 응답 특성(상승시간, 오버슈트, 정정시간 등)을 갖지않을 경우, 원하는 응답특성을 얻을때까지 단위 입력에 대한 응답결과를 반복, 수집하여 만족되는 응답특성에 이를 때의 피아이디 제어기의 계수를 최종 튜닝된 계수로써, 사용하는 방법으로서, 실제로 적용할 경우 운전자가 수많은 단위 입력에 대한 응답을 허용해야 하기 때문에 불편함을 초래할 수 있으며, 온라인 튜닝 방법은 초기에 피아이디 제어기 계수를 결정한 후, 이에대한 결과 응답시 샘플링 시간마다 피아이디 제어기 계수를 갱신하여 짧은 시간내에 상기 피아이디 제어기의 계수를 설정할 수 있어 시스템의 제어기 계수를 쉽고 용이하게 설정할 수 있다.In order to automatically tune the PID controller to the process, offline tuning and online tuning methods are used. The offline tuning method generates response characteristics (rising time, overshoot, settling time, etc.) that are satisfactory for the PID controller initially determined. If it is not available, it is a method of using the PID controller's coefficients as the final tuned coefficients when the response results of the unit input are collected and collected until the desired response characteristics are obtained. In this case, it may be inconvenient because the driver has to allow a response to a large number of unit inputs, and the online tuning method initially determines the PID controller coefficients, and then updates the PID controller coefficients at each sampling time in response to the results. The coefficient of the PID controller can be set within a short time. A system controller of the coefficients may easily be easily set.

상기와 같이 종래의 장치에 있어서는 첫 번째, 공정이 큰 시간 지연요소를 가질 때와; 두 번째, 실제 시스템은 고려하지 않은 모델특성(Unmodelled Dynamics)을 고려하지 못하고 선형화시키지 못한 비선형성이 크게 존재하는 경우와; 세 번째, 설정치 변화에 따라서 공정의 동작점이 변하는데, 이에 따라 공정의 동특성이 다르게 나타날 경우와; 네 번째, 제어 대상인 공정의 부하 변동이 심하거나 무시할 수 없는 외란(노이즈)이 제어 출력값이나 공정값에 존재하는 경우에 안정된 적용이 어려운 문제점이 있었다.In the conventional apparatus as described above, first, when the process has a large time delay element; Secondly, the actual system does not consider unmodeled dynamics and there is a large amount of nonlinearity that cannot be linearized; Third, the operating point of the process changes according to the change of the set value, and thus the dynamic characteristics of the process are different; Fourth, there is a problem that stable application is difficult in the case where the load variation of the process to be controlled is severe or a disturbance (noise) that cannot be ignored exists in the control output value or the process value.

따라서, 본 발명은 상기와 같은 종래의 문제점을 해결하기 위하여 창안한 것으로, 피아이디 제어기의 계수를 자동으로 조정하여 시스템의 속성을 나타내는 시스템 출력의 특성값인 상승시간, 오버슈트, 정정시간이 원하는 성능 지수의 시스템 출력에 이르기까지 자동으로 피아이디(PID) 계수를 적응시키는 퍼지 전문가 장치 및 제어방법을 제공함에 목적이 있다.Therefore, the present invention was devised to solve the above-mentioned conventional problems, and the rise time, overshoot, and settling time which are characteristic values of the system output indicating the property of the system by automatically adjusting the coefficients of the PID controller are desired. It is an object of the present invention to provide a fuzzy expert apparatus and control method that automatically adapts the PID coefficients up to the system output of the figure of merit.

도1은 종래의 분산제어 시스템의 구성도.1 is a block diagram of a conventional distributed control system.

도2는 종래 피아이디 기법을 적용한 피엠에스/피씨에스부의 블록 구성도.Figure 2 is a block diagram of a PS / PS unit applying a conventional PD technique.

도3은 본 발명 퍼지 전문가 장치를 이용한 피아이디 제어기의 계수조정 장치의 블록 구성도.Figure 3 is a block diagram of a coefficient adjusting device of the PID controller using the fuzzy expert apparatus of the present invention.

도4는 본 발명의 퍼지 전문가 장치의 흐름도.4 is a flow chart of a fuzzy expert apparatus of the present invention.

도5는 릴레이 궤환 시스템의 블록 구성도.5 is a block diagram of a relay feedback system.

도6은 릴레이의 디스크라이빙 함수.Fig. 6 is a disk driving function of a relay.

도7은 지글러-니콜스법에 의한 초기 피아이디 계수 결정표.7 is an initial PD coefficient determination table by the Ziegler-Nichols method.

도8은 본 발명에 따른 규칙 베이스표.8 is a rule base table according to the present invention.

도9는 도8의 퍼지집합 소속함수의 적합도.9 is a goodness of fit of the fuzzy set membership function of FIG. 8;

도10은 최적화된 규칙 베이스표.10 is an optimized rule base table.

도11은 본 발명을 이용한 모의실험 결과표.11 is a simulation result table using the present invention.

도12는 도3에서 진화연산부의 구성을 보인 블록 구성도.12 is a block diagram showing the configuration of the evolutionary calculation unit in FIG.

도13은 도12에서의 동작 흐름도.13 is an operational flowchart of FIG. 12;

*도면의 주요 부분에 대한 부호의 설명** Description of the symbols for the main parts of the drawings *

100 : 피엠에스부 110 : 관리부100: PMS 110: management unit

120 : 퍼지 전문가 장치 121 : 성능지수부120: fuzzy expert apparatus 121: performance index part

122 : 규칙 베이스부 123 : 패턴 추출부122: rule base section 123: pattern extraction section

124 : 모니터부 125 : 퍼지 추론부124: monitor unit 125: fuzzy inference unit

131 : 진화 연산부 200 : 피씨에스부131: Evolutionary Computing Unit 200: PS Part

202 : 가산기 201 : 피아이디 제어기202: adder 201: PID controller

300 : 공정부 301 : 규칙 정의부300: process unit 301: rule definition unit

302 : 개체 생성부 303 : 성능 평가부302: object generation unit 303: performance evaluation unit

304 : 비교 선택부 305 : 공정데이터 발생부304: comparison selection unit 305: process data generation unit

이와 같은 목적을 달성하기 위한 본 발명의 구성은, 도3에 도시한 바와 같이 공정을 분류하여 각 분류된 공정에 대해서 최적의 규칙베이스를 찾아내는 진화연산부(131)와, 외부 시스템이나 운전자와 인터페이스하여 제어 대상의 제어목표 설정치를 발생하고, 상기 설정치를 비롯하여 제어값과 공정값을 표시 및 관리하는 관리부(110)와, 피아이디 계수의 변경을 원하는 시스템 응답 사양(평가지수의 임계치, 속성치)을 만족할 때까지 퍼지규칙을 이용한 전문가 기법의 추론에 의해 수행하는 퍼지 전문가 장치(120)로 구성한 피엠에스부(100)와; 공정으로 부터 필드 버스(Field Bus)를 통해 직접적인 데이터 수집과 제어를 수행하는 하위 계층 스테이션인 피씨에스부(200)와; 상기 피씨에스부(200)에서 일정한 샘플링 주기에서 출력되는 제어 출력값을 입력받아 그에 해당하는 공정치를 출력하는 공정부(300)로 구성하여 된것을 특징으로 한다.The structure of the present invention for achieving the above object is, as shown in Figure 3, the evolutionary operation unit 131 for classifying the process to find the optimal rule base for each classified process, and interface with an external system or the driver The control unit 110 generates a control target set value of the control target and displays and manages the control value and the process value as well as the set value, and satisfies a system response specification (threshold value of the evaluation index, attribute value) for which the PID coefficient is changed. A PMS unit 100 configured as a fuzzy expert device 120 to perform by reasoning of expert techniques using fuzzy rules; A PC unit 200 which is a lower layer station performing direct data collection and control through a field bus from a process; The PC unit 200 may be configured to include a process unit 300 that receives a control output value output at a predetermined sampling period and outputs a process value corresponding thereto.

상기 퍼지 전문가 장치(120)는 공정의 과도응답이 끝난 후 과도응답을 분석하여 응답의 특징을 수치적인 값으로 출력하는 패턴 추출부(123)와; 제어기 계수를 바꾸는 규칙들을 가지고 있는 규칙 베이스부(122)와; 상기 패턴 추출부(123)의 출력과 각 패턴의 원하는 특징값과 상기 규칙 베이스부(122)의 규칙 베이스를 이용하여 피아이디 제어기 계수를 갱신하여 주는 퍼지추론부(125)와; 성능지수가 어느 임계치에 이를 때 상기 퍼지추론부(125)에서의 피아이디 계수의 갱신을 중지시키는 성능 지수부(121)와; 상기 성능 지수부(121)와 공정부(300)의 출력을 표시하는 모니터부(124)로 구성하여 된 것을 특징으로 한다.The fuzzy expert apparatus 120 includes a pattern extractor 123 for analyzing the transient response after the transient response of the process is finished and outputting the characteristic of the response as a numerical value; A rule base portion 122 having rules for changing controller coefficients; A fuzzy inference unit 125 for updating PID controller coefficients by using the output of the pattern extraction unit 123, desired feature values of each pattern, and a rule base of the rule base unit 122; A performance index unit 121 for stopping updating of the PID coefficients in the fuzzy inference unit 125 when the performance index reaches a certain threshold; Characterized in that it consists of a monitor 124 for displaying the output of the performance index unit 121 and the process unit 300.

본 발명의 제어방법은 주어진 분류의 공정에 대해 진화연산을 사용하여 규칙베이스를 최적화하는 제1 단계와; 대상 공정을 모델링하여 시스템을 분류하는 제2 단계와; 지글러-니콜스법에 따라 초기의 피아이디 계수를 결정하는 제3 단계와; 평가지수의 임계치 및 가중치를 결정하는 제4 단계와; 공정치를 이용하여 패턴을 추출하는 제5 단계와; 평가지수를 계산하고, 임계치의 초과 여부를 결정하여 퍼지추론에 의해 피아이디 계수를 갱신하는 제6 단계와; 갱신된 피아이디 계수를 이용하여 제어출력을 결정하고, 결정된 제어출력 및 공정치의 데이터를 저장하는 제7 단계로 이루어진 것을 특징으로 한다.The control method of the present invention comprises a first step of optimizing a rulebase using evolutionary computation for a given class of process; A second step of classifying a system by modeling a target process; Determining an initial PID coefficient according to the Ziegler-Nichols method; Determining a threshold value and a weight of the evaluation index; Extracting a pattern using a process value; Calculating an evaluation index, determining whether the threshold is exceeded, and updating the ID coefficient by fuzzy inference; And a seventh step of determining a control output using the updated PD coefficient and storing data of the determined control output and the process value.

그리고, 다른 제어방법은 주어진 분류의 공정에 대해 진화연산을 사용하여 규칙베이스를 최적화하는 제1 단계와; 대상 공정을 모델링하여 시스템을 분류하는 제2 단계와; 지글러-니콜스법에 따라 초기의 피아이디 계수를 결정하는 제3 단계와; 샘플링 주기마다 공정치를 이용하여 패턴을 추출하는 제4 단계와; 퍼지추론에 의해 피아이디 계수를 갱신하는 제5 단계와; 갱신된 피아이디 계수를 이용하여 제어출력을 결정하고, 결정된 제어출력 및 공정치의 데이터를 저장하는 제6 단계로 이루어진 것을 특징으로 한다.Another control method includes a first step of optimizing a rulebase using evolutionary operations for a given class of process; A second step of classifying a system by modeling a target process; Determining an initial PID coefficient according to the Ziegler-Nichols method; Extracting a pattern using a process value for each sampling period; Updating a PD coefficient by fuzzy inference; And a sixth step of determining a control output using the updated PD coefficient and storing data of the determined control output and the process value.

이하, 본 발명에 따른 실시예를 첨부한 도면을 참조하여 상세히 설명하면 다음과 같다.Hereinafter, exemplary embodiments of the present invention will be described in detail with reference to the accompanying drawings.

도4는 본 발명 퍼지 전문가 장치를 이용한 피아이디 제어기의 계수 조정장치의 동작 흐름도로서, 이에 도시한 바와 같이 계단응답이 단조증가 이거나 초기 부분만을 제외하고 단조증가라면 공정은 크게 두 가지 종류로 나눠질 수 있는데, 첫 번째 종류는 안정한(stable) 공정으로 동특성은 공정이득(Kp), 지연시간(L), 시상수(T)로 표현할 수 있고, 대표적인 안정한 공정의 전달함수는 다음식과 같다.4 is an operation flowchart of the coefficient adjusting device of the PID controller using the fuzzy expert apparatus of the present invention. As shown in FIG. 4, if the step response is monotonous increase or monotonic increase except for the initial portion, the process may be divided into two types. The first kind is a stable process, and its dynamic characteristics can be expressed as process gain (K p ), delay time (L), and time constant (T), and the transfer function of a typical stable process is as follows.

G(s) = Kp

Figure kpo00003
Figure kpo00004
------------------ (식1)G (s) = K p
Figure kpo00003
Figure kpo00004
------------------ (Equation 1)

두 번째 종류는 적분기를 갖는 불안정한 공정으로, 동특성은 공정이득(Kv), 지연시간(L), 시상수(Tv)로 표현할 수 있고, 불안정한 공정의 대표적인 전달함수는 다음식과 같다.The second type is an unstable process with an integrator. The dynamic characteristics can be expressed as process gain (K v ), delay time (L), and time constant (T v ), and the typical transfer function of an unstable process is as follows.

G(s) =

Figure kpo00005
---------------------- (식2)G (s) =
Figure kpo00005
---------------------- (Equation 2)

공정 분류의 기준이 되는 단위가 없는 정규화된 지연시간을 정의하면 정규화된 지연시간은 안정한 공정에 대해서 식3과 같이 지연시간과 시상수의 비로 정의되며, 정규화된 지연시간은 공정을 제어함에 있어서, 어려움의 척도로 사용되어질 수 있다. 즉, 작은

Figure kpo00006
을 갖는 공정은 제어하기 쉽고, 큰
Figure kpo00121
을 갖는 공정은 제어하기 힘들다.If you define the normalized delay time without unit as the basis of process classification, the normalized delay time is defined as the ratio of delay time and time constant for the stable process as shown in Equation 3, and the normalized delay time is difficult to control the process. Can be used as a measure of That is, small
Figure kpo00006
Process having a large, easy to control
Figure kpo00121
The process with is difficult to control.

θ1=

Figure kpo00007
--------------------------------- (식3)θ 1 =
Figure kpo00007
--------------------------------- (Eq. 3)

상기 식3과 유사하게 불안정한 공정에 대해서도 정규화된 지연시간을 정의하면 상기 식1의 전달함수[G(s)]가 불안정한 공정이라고 가정하면 상기 식2의 전달함수[s×G(s)]는 반대로 안정한 공정에 속하므로, 그 동특성을 지연시간(L), 시상수(Tv)로 표현하며, 적분기를 갖는 불안한 공정에 대하여 정규화된 지연시간을 정의하면 다음과 같다.Similarly to Equation 3, when the normalized delay time is defined for an unstable process, assuming that the transfer function [G (s)] of Equation 1 is unstable, the transfer function [s × G (s)] of Equation 2 is On the contrary, since it belongs to a stable process, its dynamic characteristics are expressed as delay time (L) and time constant (T v ), and the normalized delay time for an unstable process with an integrator is defined as follows.

θ2=

Figure kpo00008
-------------------------------- (식4)θ 2 =
Figure kpo00008
-------------------------------- (Equation 4)

정규화된 지연시간은 공정을 제어함에 있어서 어려움의 척도로 사용되므로 분류의 기준이 되며, 상기 정규화된 지연시간에 따라 다음과 같이 4가지의 안정한 공정과 2가지의 불안정한 공정으로 분류한다.Since the normalized delay time is used as a measure of difficulty in controlling the process, it is a standard of classification, and according to the normalized delay time, it is classified into four stable processes and two unstable processes as follows.

1) 안정한 공정1) stable process

[분류 1]

Figure kpo00009
을 만족하는 공정[Category 1]
Figure kpo00009
Satisfying process

[분류 2]

Figure kpo00010
을 만족하는 공정[Category 2]
Figure kpo00010
Satisfying process

[분류 3]

Figure kpo00011
을 만족하는 공정[Category 3]
Figure kpo00011
Satisfying process

[분류 4]

Figure kpo00012
을 만족하는 공정[Category 4]
Figure kpo00012
Satisfying process

2) 불안정한 공정2) unstable process

[분류 5]

Figure kpo00013
을 만족하는 공정[Category 5]
Figure kpo00013
Satisfying process

[분류 6]

Figure kpo00014
을 만족하는 공정[Category 6]
Figure kpo00014
Satisfying process

초기의 피아이디 계수를 지글러-니콜스(Ziegler-Nichols)법에 따라 결정하고, 이때 피아이디 계수를 구하는 방법은 릴레이 응답에서 구한다.The initial PID coefficient is determined according to the Ziegler-Nichols method, and the method of obtaining the PI coefficient is obtained from the relay response.

도5는 릴레이 궤환시스템의 구성도이고, 도6은 릴레이의 디스크라이빙 함수이며, 도7은 지글러-니콜스법에 의한 초기 피아이디 계수 결정표로서, 이에 도시한 바와 같이 릴레이를 시스템에 직렬로 연결하면 다음식을 만족할 때 진동을 일으키케 된다.FIG. 5 is a block diagram of a relay feedback system, FIG. 6 is a disk driving function of a relay, and FIG. 7 is an initial PID coefficient determination table by the Ziegler-Nichols method. When connected, vibration occurs when the following equation is satisfied.

Figure kpo00015
또는
Figure kpo00016
----------- (식5)
Figure kpo00015
or
Figure kpo00016
----------- (Eq. 5)

상기 식5에서 N(A)는 릴레이의 디스크라이빙(describing) 함수이며, 이때의 진동주기(tc)와, 진폭(A)을 측정함으로써 대상 시스템에 대한 나이퀴스트(Nyquist) 곡선상의 어떤 특정한 점의 값을 알 수 있고, 이것을 이용하여 초기 이득을 결정한다.In Equation 5, N (A) is a describing function of the relay, and the oscillation period (t c ) and the amplitude (A) of the relay are measured on the Nyquist curve for the target system. The value of a particular point can be known and used to determine the initial gain.

높은 주파수 영역에서 적어도 180°의 위상뒤짐(Phase Lag)을 갖는 시스템에 릴레이를 되먹임 제어 루프를 넣으면 시스템의 출력은 진동을 하게되며, 릴레이의 진폭이 N≥0, M≥0 이라면 릴레이의 디스크라이빙 함수(N(A))는 다음과 같다.In the high frequency range, a relay loop is fed back into the system with at least 180 ° phase lag. The output of the system oscillates. If the relay amplitude is N≥0 and M≥0, the relay's disk The ice function N (A) is

N(A) =

Figure kpo00017
N (A) =
Figure kpo00017

여기에서,From here,

b1=

Figure kpo00018
b 1 =
Figure kpo00018

a1=

Figure kpo00019
a 1 =
Figure kpo00019

이므로, 릴레이의 디스크라이빙 함수는 다음과 같다.Therefore, the disk drive function of the relay is as follows.

Figure kpo00020
---------------------------- (식6)
Figure kpo00020
---------------------------- (Eq. 6)

릴레이의 되먹임을 가해 생긴 진동의 주기가 대략적으로 시스템의 임계주기(ultimate perriod)인 tc값이되고, 이때의 진폭을 A라할 때 상기 식5와 식6을 이용하여 다음식을 얻게된다.The period of vibration generated by the feedback of the relay is approximately t c , the critical period of the system (ultimate perriod), and when the amplitude is A, the following equations are obtained using Equations 5 and 6 above.

Figure kpo00021
Figure kpo00021

여기에서

Figure kpo00022
는 시스템의 임계이득(ultimate gain)이고, 릴레이 되먹임 실험을 통해 얻은
Figure kpo00023
,
Figure kpo00024
값을 이용하여 지글러-니콜스(Ziegler-Nichols)법에 따라 초기의 피아이디 계수를 결정하며, 이때 사용되는 주파수 영역의 지글러-니콜스(Ziegler - Nichols)법은 도7에 도시한 바와 같다.From here
Figure kpo00022
Is the critical gain of the system and is obtained from the relay feedback experiment.
Figure kpo00023
,
Figure kpo00024
The initial PID coefficient is determined according to the Ziegler-Nichols method using the value, and the Ziegler-Nichols method in the frequency domain used at this time is shown in FIG.

도3에서 퍼지 전문가 장치(120)는 패턴 추출부(123), 퍼지추론부(125), 규칙 베이스부(122) 등으로 구성되어 있는데, 제어방법이 온라인 튜닝 방법이냐 또는 오프라인 튜닝 방법이냐에 따라 수행하는 방법이 다르다.In FIG. 3, the fuzzy expert apparatus 120 includes a pattern extraction unit 123, a fuzzy inference unit 125, a rule base unit 122, and the like, depending on whether the control method is an online tuning method or an offline tuning method. How to do it is different.

오프라인 튜닝 방법이면 패턴 추출부(123)는 공정의 과도응답이 끝난 후 과도응답을 분석하여 응답의 특징을 수치적인 값으로 출력하며, 퍼지추론부(125)는 상기 패턴 추출부(123)의 출력과 각 패턴의 원하는 특정값과 진화연산에 의하여 최적화 되어있는 규칙 베이스(제어기 계수를 바꾸는 규칙들)를 이용하여 피아이디 제어기 계수를 갱신한다.In the offline tuning method, the pattern extractor 123 analyzes the transient response after the transient response is completed and outputs the characteristic of the response as a numerical value, and the fuzzy inference unit 125 outputs the pattern extractor 123. And the PID controller coefficients are updated by using a rule base (rules for changing controller coefficients) optimized by the specific value and the evolutionary operation desired for each pattern.

온라인 튜닝 방법이면 상기 패턴 추출부(123)는 공정의 샘플링 주기마다 응답의 특징을 수치적인 값으로 출력한다는 점이 다르다.In the online tuning method, the pattern extractor 123 differs in that it outputs the characteristic of the response as a numerical value for each sampling period of the process.

또한, 상기 퍼지 전문가 장치(120)는 오프라인 튜닝 방법이면 플랜트에 계단 입력을 가했을 때 플랜트 응답이 원하는 세 가지 속성치(상승시간, 오버슈트, 정정시간)가 최소화 되도록 피아이디 제어기(201)의 계수들을 동조하며, 패턴 추출부(123)는 과도응답이 끝난 후 세 가지의 응답특징을 추출한다.In addition, in the offline tuning method, the fuzzy expert apparatus 120 adjusts the coefficients of the PID controller 201 to minimize three property values (rising time, overshoot, and settling time) that the plant response desires when a step input is applied to the plant. Tuning, the pattern extractor 123 extracts three response features after the transient response is over.

이 세가지의 속성치는 다음 식과 같이 정의한다.These three attribute values are defined as follows.

Figure kpo00025
-------------------------- (식7-1)
Figure kpo00025
-------------------------- (Eq. 7-1)

상기 (식7-1)에서는

Figure kpo00026
은 상승시간,
Figure kpo00027
는 오버슈트,
Figure kpo00028
는 정정시간이고,
Figure kpo00029
는 패턴 추출부(123)의 출력이며,
Figure kpo00030
는 각 응답 속성의 원하는 값이고,
Figure kpo00031
는 각 응답 속성에 대한 가중치를 나타낸다.In Equation 7-1
Figure kpo00026
Is the rise time,
Figure kpo00027
Overshoot,
Figure kpo00028
Is the settling time,
Figure kpo00029
Is the output of the pattern extraction unit 123,
Figure kpo00030
Is the desired value for each response attribute,
Figure kpo00031
Denotes a weight for each response attribute.

온라인 튜닝 방법이면 플랜트에 계단 입력을 가했을 때 플랜트 응답이 원하는 세 가지 속성치(오차, 오차율, 오차합)가 최소화 되도록 피아이디 제어기(201)의 계수들을 동조하며, 패턴 추출부(123)는 샘플링 주기마다 세 가지의 응답특징을 추출한다. 이 세가지의 속성치는 다음 식과 같이 정의한다.In the online tuning method, when the stair input is applied to the plant, the plant response tunes the coefficients of the PID controller 201 so that three desired attribute values (error, error rate, and error sum) are minimized, and the pattern extraction unit 123 performs a sampling period. Each of the three response features is extracted. These three attribute values are defined as follows.

e(K) = (yr(K)-y(K))e (K) = (y r (K) -y (K))

Δe(K) = (e(K)-e(K-1)) ------------------------ (식7-2)Δe (K) = (e (K) -e (K-1)) ------------------------ (Equation 7-2)

Figure kpo00032
=
Figure kpo00033
Figure kpo00032
=
Figure kpo00033

상기 (식7-2)에서는 e(K)는 오차, Δe(K)는 오차율,

Figure kpo00034
는 오차합, yr(K)는 설정치, y(K)는 현재의 공정 출력치를 나타낸다.In Equation 7-2, e (K) is an error, Δe (K) is an error rate,
Figure kpo00034
Is the sum of errors, y r (K) is the set point, and y (K) is the current process output.

도8은 본 발명에 따른 규칙 베이스표이고, 도9는 도8의 퍼지집합 소속함수의 적합도로서, 이에 도시한 바와 같이 퍼지규칙은 3개의 출력 퍼지변수를 가지고 있는데, 오프라인 튜닝 방법 또는 온라인 튜닝 방법에 따라 출력 퍼지변수가 다르게 적용된다.FIG. 8 is a rule base table according to the present invention, and FIG. 9 is a goodness of fit of the fuzzy set membership function of FIG. 8. As shown in FIG. 8, the fuzzy rule has three output fuzzy variables. Depending on the output fuzzy variable is applied differently.

도8a에서와 같이 오프라인 튜닝 방법이면 전건부(antecedent)에 해당하는 입력 퍼지변수는

Figure kpo00035
,
Figure kpo00036
,
Figure kpo00037
이고, 후건부(consequent)에 해당하는 출력변수는 ΔKc, ΔTi,
Figure kpo00038
ΔTd이며, 이때 퍼지집합 변수는 GOOD, BAD이고, 다음과 같은 의미를 갖는다.In the offline tuning method as shown in FIG. 8A, the input fuzzy variable corresponding to the antecedent is
Figure kpo00035
,
Figure kpo00036
,
Figure kpo00037
The output variables corresponding to the consequent parts are ΔKc, ΔT i,
Figure kpo00038
ΔT d , where the fuzzy set variables are GOOD and BAD and have the following meanings.

GOOD : 원하는 속성을 만족함GOOD: Satisfy the desired property

BAD : 원하는 속성을 만족하지 못함BAD: does not satisfy desired attributes

즉, 원하는 속성보다 좋으면 GOOD이고, 좋지 않으면 BAD를 의미하며, 각 집합의 소속함수는 도9a에 도시한 바와 같은 형태이고, 도8a에 도시한 바와 같이 9개의 규칙(

Figure kpo00039
,
Figure kpo00040
,∼,
Figure kpo00041
)을 사용하였다.That is, if the property is better than the desired property, GOOD, and if not, it means BAD, and the membership function of each set is as shown in FIG. 9A, and as shown in FIG.
Figure kpo00039
,
Figure kpo00040
, ~,
Figure kpo00041
) Was used.

도8b에서와 같이 온라인 튜닝 방법이면 전건부(antecedent)에 해당하는 입력 퍼지변수는 오차, 오차율, 오차합이고, 후건부(consequent)에 해당하는 출력변수는 피아이디 계수의 변화량인ΔKc, ΔTi,

Figure kpo00042
ΔTd이며, 이때 퍼지집합 변수는 포지티브, 제로, 네거티브이고, 다음과 같은 의미를 갖는다.A rear-line tuning method as in 8b around the conditional input fuzzy variables for the (antecedent) is error, error rates, error sum, and then conditional output variable corresponding to the (consequent) is blood ID coefficient change amount of ΔKc, ΔT i of ,
Figure kpo00042
ΔT d , where the fuzzy set variable is positive, zero, negative and has the following meaning.

포지티브(Positive) : 정(+) 방향Positive: Positive Direction

제로(Zero) : 제로Zero: Zero

네거티브(Negative) : 부(-) 방향Negative: negative direction

즉, 입력에 대하여 퍼지집합 변수인 포지티브(Positive)는 입력이 양인 값인 집합, 제로(Zero)는 입력이 대략 0인 집합, 네거티브(Negative)는 입력이 음의 값인 집합을 의미하며, 각 집합의 소속함수는 도9b에 도시한 바와 같은 형태이고, 도8a에 도시한 바와 같이 27개의 규칙(

Figure kpo00043
,
Figure kpo00044
, ∼,
Figure kpo00045
)을 사용하였는데,
Figure kpo00046
는 실수값인 싱글톤(Singleton)이다.That is, a positive fuzzy set variable with respect to the input means a set of positive inputs, zero means a set of approximately zero inputs, and a negative means a set of negative inputs. The membership function has the form shown in FIG. 9B, and as shown in FIG. 8A, the 27 rules (
Figure kpo00043
,
Figure kpo00044
, ~,
Figure kpo00045
),
Figure kpo00046
Is the real value Singleton.

피아이디 제어기의 계수들(오프라인 튜닝 방법, 온라인 튜닝 방법)은 다음의 식에 의해서 구해진다.The coefficients of the PID controller (offline tuning method, online tuning method) are obtained by the following equation.

Figure kpo00047
Figure kpo00047

Figure kpo00048
Figure kpo00048

Figure kpo00049
Figure kpo00049

여기서,

Figure kpo00050
,
Figure kpo00051
,
Figure kpo00052
는 각각 피아이디 제어기의
Figure kpo00053
번째 계수이고,
Figure kpo00054
,
Figure kpo00055
,
Figure kpo00056
는 각각 피아이디 제어기 계수들을 갱신하는 도3의 퍼지추론부(125)의 출력이다.here,
Figure kpo00050
,
Figure kpo00051
,
Figure kpo00052
Each of the PID controllers
Figure kpo00053
Th coefficient,
Figure kpo00054
,
Figure kpo00055
,
Figure kpo00056
Is the output of the fuzzy inference unit 125 of FIG. 3, which updates the PID controller coefficients, respectively.

전체 시스템의 성능지수는 도3의 성능지수부(121)에서 관리하는데, 이는 오프라인 튜닝 방법과 온라인 튜닝 방법에 의하여 구분되며, 오프라인 튜닝 방법에 의한 성능지수는 다음 식과 같다.The performance index of the entire system is managed by the performance index unit 121 of FIG. 3, which is classified by the offline tuning method and the online tuning method, and the performance index by the offline tuning method is as follows.

성능지수 =

Figure kpo00057
----------------------- (식9-1)Figure of merit =
Figure kpo00057
----------------------- (Eq. 9-1)

상기 (식9-1)에서

Figure kpo00058
는 평가지수의 가중치를 나타내며, 성능지수가 어느 임계치에 이를 때 피아이디 계수의 갱신을 중지하도록 하고,
Figure kpo00059
는 패턴 추출부(123)의 출력이다.In the above formula (9-1)
Figure kpo00058
Denotes the weight of the evaluation index, stops updating the PID coefficients when the performance index reaches a certain threshold,
Figure kpo00059
Is the output of the pattern extraction unit 123.

온라인 튜닝 방법에 의한 성능지수는 다음 식과 같다.The figure of merit by online tuning method is as follows.

성능지수 =

Figure kpo00060
---------------------- (식9-2)Figure of merit =
Figure kpo00060
---------------------- (Eq. 9-2)

상기 (식9-2)에서

Figure kpo00061
는 평가지수의 가중치를 나타내며, 성능지수가 어느 임계치에 이를 때 피아이디 계수의 갱신을 중지하도록 하고, 여기서
Figure kpo00062
은 상승시간(Rise Time),
Figure kpo00063
는 오버슈트(Overshoot),
Figure kpo00064
은 정정시간(Settling Time)이다.In the above formula (9-2)
Figure kpo00061
Denotes the weight of the evaluation index, and when the performance index reaches a certain threshold to stop updating the PID coefficients, where
Figure kpo00062
Is the rise time,
Figure kpo00063
Means overshoot,
Figure kpo00064
Is the settling time.

본 발명의 시스템 운용을 보다 적응적으로 수행하는데 적당하도록 도8의 규칙 베이스가 최적화되어야 하고, 최적화의 대상은 도8a에서의 9개의 규칙(

Figure kpo00065
,
Figure kpo00066
,∼,
Figure kpo00067
)과 도8b에서의 27개의 규칙(
Figure kpo00068
,
Figure kpo00069
, ∼,
Figure kpo00070
)이며, 이를 최적화하기 위해 진화연산(EC : Evolutionary Computation)을 사용하는데, 상기 진화연산은 자연 진화에서의 서로 다른 측면들을 강조하는 유전 알고리즘(GA : Genetic Algorithm), 진화기법(ES : Evolutionary Strategy), 진화 프로그래밍(EP : Evolutionary Programming) 등으로 구분되며, 상기 유전 알고리즘은 염색체 단위에서의 진화를 중요시하고, 상기 진화기법은 개체단위에서의 행동의 변화를, 상기 진화 프로그래밍은 종 단위에서의 행동변화를 중요한 현상으로 여기고 있으며, 이들의 공통적인 개념은 최적의 개체만이 살아남는다는 자연계의 적자 생존의 원리를 최적화기법으로 사용한다.The rule base of FIG. 8 should be optimized to be more adaptive to perform the system operation of the present invention, and the object of optimization is the nine rules (Fig. 8A).
Figure kpo00065
,
Figure kpo00066
, ~,
Figure kpo00067
) And the 27 rules in FIG.
Figure kpo00068
,
Figure kpo00069
, ~,
Figure kpo00070
Evolutionary Computation (EC) is used to optimize this, which is a genetic algorithm (GA) that emphasizes different aspects of natural evolution, and an evolutionary strategy (ES). , Evolutionary Programming (EP), etc., wherein the genetic algorithm emphasizes evolution at the chromosome level, and the evolutionary technique is a change in behavior at the individual level, and the evolutionary programming is a change in behavior at the species level. Are considered to be important phenomena, and their common concept uses the principle of survival of the fittest in nature that only optimal individuals survive.

일반적으로 퍼지논리는 전문가의 지식을 효과적으로 표현할 수 있어서 피아이디 제어 계수를 제어시스템이 원하는 특성을 만족하도록 하는 전문가 장치를 구성할 때 적합하지만 전문가로부터 체계적으로 정보를 얻기가 어렵기 때문에 진화연산을 사용하여 전문가 장치의 지식을 최적화하는 기능을 제공한다.In general, fuzzy logic is suitable for constructing an expert device that can effectively express the knowledge of experts so that the PID control coefficients satisfy the desired characteristics of the control system, but it is difficult to obtain information systematically from the expert. It provides the ability to optimize the knowledge of expert devices.

도12는 도3에서 진화연산부의 구성을 보인 블록 구성도이고, 도13은 도12에서의 동작 흐름도로서, 이에 도시한 바와 같이 규칙의 최적화를 위하여 먼저 각 규칙에 대한 개체들의 초기화(P개의 부모개체를 랜덤하게 생성하고, 변수의 개수가 N개이면 1개체는 N개의 실수배열[(

Figure kpo00071
Figure kpo00072
,
Figure kpo00073
),
Figure kpo00074
∈{1,∼,P}]로 구성)를 하고, 각각의 부모개체에 대해서 부모개체 평가함수(J)를 계산하는데, 평가함수는 다음 식과 같다.FIG. 12 is a block diagram showing the construction of the evolutionary operation unit in FIG. 3, and FIG. 13 is an operation flowchart of FIG. 12. As shown in FIG. If objects are generated randomly and the number of variables is N, one object is N real arrays [(
Figure kpo00071
Figure kpo00072
,
Figure kpo00073
),
Figure kpo00074
∈ {1, ~, P}], and the parent object evaluation function (J) is calculated for each parent object. The evaluation function is as follows.

J =

Figure kpo00075
----------------- (식10)J =
Figure kpo00075
----------------- (Eq. 10)

여기에서

Figure kpo00076
는 상수이고,
Figure kpo00077
은 상승시간(Rise Time),
Figure kpo00078
는 오버슈트(Overshoot),
Figure kpo00079
은 정정시간(Settling Time) 이다.From here
Figure kpo00076
Is a constant,
Figure kpo00077
Is the rise time,
Figure kpo00078
Means overshoot,
Figure kpo00079
Is the settling time.

상기 (식10)에 의한 평가함수(J)를 계산한 다음 자손개체를 생성하는데, 자손 생성에는 다음 식을 이용하여 생성한다.After calculating the evaluation function (J) according to Equation 10, the offspring are generated, and the offspring is generated using the following equation.

Figure kpo00080
=
Figure kpo00081
Figure kpo00082
+
Figure kpo00083
Figure kpo00080
=
Figure kpo00081
Figure kpo00082
+
Figure kpo00083

Figure kpo00084
=
Figure kpo00085
------------ (식11)
Figure kpo00084
=
Figure kpo00085
------------ (Eq. 11)

상기 (식11)에서

Figure kpo00086
,
Figure kpo00087
이다.In the above formula (11)
Figure kpo00086
,
Figure kpo00087
to be.

자손개채를 생성하고 나서, 자손개체에 대한 평가함수

Figure kpo00088
를 계산하고, 부모/자손개체 모두 포함하여 2×P개의 개체에 대해 확률적으로 선택된 비교(Q개의 무작위로 선택된 상대방과 비교)를 하여 자신이 상대방보다 낮은 목적함수를 가지게되면 승(win)수를 한점씩 올리며, 또한 2×P개의 개체 중에서 승수가 높은 P개의 개체의 다음 세대의 부모로서 선택하고, 종료조건을 만족하는지를 판단하여 만족하면 종료하고, 만족하지 못하면 다시 자손개체 생성부터 다시 시작한다.Evaluate function of offspring after creating offspring
Figure kpo00088
If we have a lower objective function than the other party, we compare the randomly selected targets (compared with Q randomly selected opponents) against 2 × P individuals, including both parents and grandchildren. Is increased by one point, and is selected as the parent of the next generation of P entities with a high multiplier among 2 × P entities, and if it satisfies the termination condition, it is terminated if it is satisfied. .

일예로 다음과 같은 6가지 공정에 대해서 상기 진화 프로그래밍을 사용하여 상기 도8의 규칙( r1∼ r9)( r1∼ r27)을 최적화하였는데, 안정한 공정의 특성을 포괄하는 4가지의 공정모델과 불안정한 공정의 특성을 포괄하는 2가지의 공정모델을 사용하여 나타내면,For example, the evolutionary programming was used to optimize the rules (r 1 to r 9 ) (r 1 to r 27 ) of FIG. 8 for the following six processes. Using two process models covering the model and the characteristics of the unstable process,

가) 안정한 공정의 경우A) for stable processes

Figure kpo00089
Figure kpo00090
Figure kpo00089
Figure kpo00090

Figure kpo00091
Figure kpo00092
Figure kpo00091
Figure kpo00092

Figure kpo00093
Figure kpo00094
Figure kpo00093
Figure kpo00094

Figure kpo00095
Figure kpo00096
Figure kpo00095
Figure kpo00096

나) 불안정한 공정의 경우B) unstable processes

Figure kpo00097
Figure kpo00098
Figure kpo00097
Figure kpo00098

Figure kpo00099
Figure kpo00100
Figure kpo00099
Figure kpo00100

로 나타내는데, 상기 식들은 공정에서 발생하는 대부분의 경우를 포함하고 있어 다양한 공정에 적용할 수 있다.The equations include most of the cases occurring in the process and can be applied to various processes.

도11은 본 발명을 이용한 모의실험 결과표로서, 이에 도시한 바와 같이 공정 모델들에 대해서 다음과 같은 진화 프로그래밍을 사용하여 최적화 하였다.FIG. 11 is a simulation result table using the present invention. As shown in FIG. 11, the process models were optimized using the following evolutionary programming.

개체군을 초기화하고 초기 개체군 즉, 부모 세대를 [-1, 1]에서 임으로 20개를 생성하는 제1 단계와; 개체군을 다음과 같은 평가함수에 의해 평가하는 제2 단계와;A first step of initializing the population and generating twenty initial populations, ie, parent generations from [-1, 1]; A second step of evaluating the population by the following evaluation function;

J = a

Figure kpo00101
(
Figure kpo00102
은 패턴 추출부의 출력, a는 상수)J = a
Figure kpo00101
(
Figure kpo00102
Is the output of the pattern extractor, a is a constant)

새로운 개체군, 즉 자손 세대를 부모 세대로부터 생성할 때는 다음 식과 같이 부모세대에 외란을 주어 생성하는 제3 단계와;A third step of generating a new population, ie, a generation of offspring from the parent generation, by disturbing the parent generation as in the following equation;

자손세대 = 부모세대 + N(0,J)Offspring = parent generation + N (0, J)

새로운 개체군을 평가하는 제4 단계와; 10개의 평가비교 대상을 확률적인 선정을 통해서 승의 계수(평가치(J)가 작은 것이 승이 된다)로 하는 제5 단계와; 모든 개체군을 승수가 많은 순서대로 서열화하고, 그 중에서 20개를 선택하는 제6 단계와; 만족한 해이면 종료하고, 그렇지 않으면 상기 제3 단계로 돌아가는 제7 단계로 이루어져 있다.A fourth step of evaluating the new population; A fifth step of making ten coefficients of evaluation (probably the smallest value of evaluation J) by means of probabilistic selection; A sixth step of sequencing all populations in order of multiplier, and selecting twenty among them; If it is a satisfactory solution, it ends, and if not, it returns to the third step.

상기 진화 프로그램을 1000세대 동안 적용하여 도10과 같은 규칙 베이스를 얻어 이 규칙 베이스를 사용하여 모의실험을 하여 도11의 결과를 얻었다.The evolution program was applied for 1000 generations to obtain a rule base as shown in FIG. 10, and simulated using the rule base to obtain the result of FIG.

이상에서 설명한 바와 같이, 본 발명 퍼지 전문가 장치를 이용한 피아이디 제어기의 계수 조정장치 및 방법은 지글러-니콜스(Z-N)법을 적용했을 때의 상승시간, 정정시간보다 작은 값을 갖도록 하고, 오버슈트의 경우 '0'이 되도록하며, [분류 2]의 경우 오버슈트가 작으므로 오버슈트를 '0'이 되도록하는 것과 동시에 지글러-니콜스(Z-N)법을 적용했을 때의 상승시간, 정정시간보다 향상되게 하는 효과가 있다.As described above, the apparatus and method for adjusting the coefficients of the PID controller using the fuzzy expert apparatus of the present invention have a value smaller than the rise time and the settling time when the Ziegler-Nichols (ZN) method is applied, In the case of [Category 2], the overshoot is small so that the overshoot is '0' and the rise time and correction time when Ziegler-Nichols (ZN) method is applied. It is effective.

Claims (7)

공정을 분류하여 각 분류된 공정에 대해서 최적의 규칙베이스를 찾아내는 진화연산부와, 외부 시스템이나 운전자와 인터페이스하여 제어 대상의 제어목표 설정치를 발생하고, 상기 설정치를 비롯하여 제어값과 공정값을 표시 및 관리하는 관리부와, 피아이디 계수의 변경을 원하는 시스템 응답 사양(평가지수의 임계치, 속성치)을 만족할 때까지 퍼지규칙을 이용한 전문가 기법의 추론에 의해 수행하는 퍼지 전문가 장치로 구성한 피엠에스부와; 공정으로 부터 필드 버스(Field Bus)를 통해 직접적인 데이터 수집과 제어를 수행하는 하위 계층 스테이션인 피씨에스부와; 상기 피씨에스부에서 일정한 샘플링 주기에서 출력되는 제어 출력값을 입력받아 그에 해당하는 공정치를 출력하는 공정부로 구성하여 된 것을 특징으로 하는 퍼지 전문가 장치를 이용한 피아이디 제어기의 계수조정 장치.The evolutionary operation unit that classifies the process and finds the optimal rule base for each classified process, and interfaces with an external system or an operator to generate the control target setpoint for the control object, and displays and manages the control value and the process value as well as the setpoint. And a fuzzy expert device configured to perform by inference of an expert technique using fuzzy rules until a system response specification (a threshold value and an attribute value of an evaluation index) desired to change the PD coefficient is satisfied. A PCS unit, which is a lower layer station that performs data collection and control directly through a field bus from a process; And a control unit configured to receive a control output value output at a predetermined sampling period from the PC unit and output a process value corresponding to the control unit. 주어진 분류의 공정에 대해 진화연산을 사용하여 규칙베이스를 최적화하는 제1 단계와; 대상 공정을 모델링하여 시스템을 분류하는 제2 단계와; 지글러-니콜스법에 따라 초기의 피아이디 계수를 결정하는 제3 단계와; 평가지수의 임계치 및 가중치를 결정하는 제4 단계와; 공정치를 이용하여 패턴을 추출하는 제5 단계와; 평가지수를 계산하고, 임계치의 초과 여부를 결정하여 퍼지추론에 의해 피아이디 계수를 갱신하는 제6 단계와; 갱신된 피아이디 계수를 이용하여 제어출력을 결정하고, 결정된 제어출력 및 공정치의 데이터를 저장하는 제7 단계로 이루어진 것을 특징으로 하는 퍼지 전문가 장치를 이용한 피아이디 제어기의 계수조정 제어방법.A first step of optimizing the rulebase using evolutionary operations for a given class of process; A second step of classifying a system by modeling a target process; Determining an initial PID coefficient according to the Ziegler-Nichols method; Determining a threshold value and a weight of the evaluation index; Extracting a pattern using a process value; Calculating an evaluation index, determining whether the threshold is exceeded, and updating the ID coefficient by fuzzy inference; And a seventh step of determining a control output using the updated PD coefficient and storing data of the determined control output and the process value. 주어진 분류의 공정에 대해 진화연산을 사용하여 규칙베이스를 최적화하는 제1 단계와; 대상 공정을 모델링하여 시스템을 분류하는 제2 단계와; 지글러-니콜스법에 따라 초기의 피아이디 계수를 결정하는 제3 단계와; 샘플링 주기마다 공정치를 이용하여 패턴을 추출하는 제4 단계와; 퍼지추론에 의해 피아이디 계수를 갱신하는 제5 단계와; 갱신된 피아이디 계수를 이용하여 제어출력을 결정하고, 결정된 제어출력 및 공정치의 데이터를 저장하는 제6 단계로 이루어진 것을 특징으로 하는 퍼지 전문가 장치를 이용한 피아이디 제어기의 계수조정 제어방법.A first step of optimizing the rulebase using evolutionary operations for a given class of process; A second step of classifying a system by modeling a target process; Determining an initial PID coefficient according to the Ziegler-Nichols method; Extracting a pattern using a process value for each sampling period; Updating a PD coefficient by fuzzy inference; And a sixth step of determining a control output using the updated PD coefficient and storing data of the determined control output and the process value. 제1항에 있어서, 상기 진화연산부는 전문가 장치의 규칙 구조를 정의하는 규칙 정의부와; 각 규칙에 대한 개체군을 생성하는 개체 생성부와; 각 개체에 대하여 지정된 성능에 대한 평가를 하는 성능 평가부와; 상기 성능 평가부에서 평가된 개체에 대하여 각각을 비교하여 정의된 선택기준에 의하여 부모개체로 선택되도록 하는 비교 선택부와; 성능 평가를 위한 공정의 데이터를 발생하는 공정데이터 발생부로 구성함을 특징으로 하는 퍼지 전문가 장치를 이용한 피아이디 제어기의 계수조정 장치.The apparatus of claim 1, wherein the evolution calculation unit comprises: a rule defining unit defining a rule structure of the expert apparatus; An entity generating unit generating an entity group for each rule; A performance evaluation unit for evaluating a specified performance for each individual; A comparison selection unit which compares each of the individuals evaluated by the performance evaluation unit and selects the parent object based on a defined selection criterion; A coefficient adjusting device of a PID controller using a fuzzy expert device, characterized in that the process data generator for generating data of the process for performance evaluation. 제1항에 있어서, 상기 퍼지 전문가 장치는 공정의 과도응답이 끝난 후 과도응답을 분석하여 응답의 특징을 수치적인 값으로 출력하는 패턴 추출부와; 제어기 계수를 바꾸는 규칙들을 가지고 있는 규칙 베이스부와; 상기 패턴 추출부의 출력과 각 패턴의 원하는 특징값과 상기 규칙 베이스부의 규칙 베이스를 이용하여 피아이디 제어기 계수를 갱신하여 주는 퍼지추론부와; 성능지수가 어느 임계치에 이를 때 상기 퍼지추론부에서의 피아이디 계수의 갱신을 중지시키는 성능 지수부와; 상기 성능 지수부와 공정부의 출력을 표시하는 모니터부로 구성함을 특징으로 하는 퍼지 전문가 장치를 이용한 피아이디 제어기의 계수조정 장치.The fuzzy expert apparatus of claim 1, further comprising: a pattern extractor configured to analyze the transient response after the transient response of the process is completed and output a characteristic of the response as a numerical value; A rule base portion having rules for changing controller coefficients; A fuzzy inference unit for updating a PID controller coefficient by using the output of the pattern extracting unit, a desired feature value of each pattern, and a rule base of the rule base unit; A performance index unit for stopping updating of the PID coefficients in the fuzzy inference unit when the performance index reaches a certain threshold; And a monitor unit for displaying outputs of the performance index unit and the process unit. 제2항에 있어서, 상기 제6 단계는 과도응답이 끝난 후 세 가지의 응답특징(상승시간, 오버슈트, 정정시간)을 추출하는 단계와; 상기 세 가지의 응답특징으로부터 성능지수를 계산하여 성능지수가 원하는 사양을 만족하는지를 검사하는 단계와; 상기 단계의 검사에 의해 만족하면 수행을 끝내고, 만족하지 않으면 선택된 규칙을 사용하여 퍼지추론을 통해 피아이디 계수를 갱신하는 단계를 포함하여 이루어진 것을 특징으로 하는 퍼지 전문가 장치를 이용한 피아이디 제어기의 계수조정 제어방법.3. The method of claim 2, wherein the sixth step comprises: extracting three response features (rising time, overshoot, correction time) after the transient response is finished; Calculating a figure of merit from the three response features and checking whether the figure of merit satisfies a desired specification; If satisfied by the check of the above step is finished, if not satisfied by using a selected rule updating the PID coefficients through the fuzzy inference coefficient adjustment of the PID controller using a fuzzy expert device, characterized in that Control method. 제3항에 있어서, 상기 제5 단계는 샘플링 주기마다 세 가지의 응답특징(오차, 오차율, 오차합)을 추출하는 단계와; 상기 세 가지의 응답특징으로부터 성능지수를 계산하여 성능지수가 원하는 사양을 만족하는지를 검사하는 단계와; 상기 단계의 검사에 의해 만족하면 수행을 끝내고, 만족하지 않으면 선택된 규칙을 사용하여 퍼지추론을 통해 피아이디 계수를 갱신하는 단계를 포함하여 이루어진 것을 특징으로 하는 퍼지 전문가 장치를 이용한 피아이디 제어기의 계수조정 제어방법.4. The method of claim 3, wherein the fifth step comprises: extracting three response features (error, error rate, error sum) for each sampling period; Calculating a figure of merit from the three response features and checking whether the figure of merit satisfies a desired specification; If satisfied by the check of the above step is finished, if not satisfied by using a selected rule updating the PID coefficients through the fuzzy inference coefficient adjustment of the PID controller using a fuzzy expert device, characterized in that Control method.
KR1019970020917A 1997-05-27 1997-05-27 Coefficient adjusting apparatus of proportional integral derivative(pid) controller KR100310609B1 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
KR1019970020917A KR100310609B1 (en) 1997-05-27 1997-05-27 Coefficient adjusting apparatus of proportional integral derivative(pid) controller

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
KR1019970020917A KR100310609B1 (en) 1997-05-27 1997-05-27 Coefficient adjusting apparatus of proportional integral derivative(pid) controller

Publications (2)

Publication Number Publication Date
KR19980084979A KR19980084979A (en) 1998-12-05
KR100310609B1 true KR100310609B1 (en) 2001-12-17

Family

ID=37530936

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1019970020917A KR100310609B1 (en) 1997-05-27 1997-05-27 Coefficient adjusting apparatus of proportional integral derivative(pid) controller

Country Status (1)

Country Link
KR (1) KR100310609B1 (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2011131078A1 (en) * 2010-04-20 2011-10-27 杭州和利时自动化有限公司 Data management method and system for field bus instrument management system
KR101242678B1 (en) * 2011-02-28 2013-03-12 고려대학교 산학협력단 PID control method of changing parameters adaptively and apparatus thereof
KR101640931B1 (en) * 2016-02-17 2016-07-20 한국기계연구원 automatic PID gain tuning method in articulated robot

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117555224A (en) * 2024-01-09 2024-02-13 冠县新瑞实业有限公司 Wheat gluten production control method based on decoupling control system

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2011131078A1 (en) * 2010-04-20 2011-10-27 杭州和利时自动化有限公司 Data management method and system for field bus instrument management system
KR101242678B1 (en) * 2011-02-28 2013-03-12 고려대학교 산학협력단 PID control method of changing parameters adaptively and apparatus thereof
KR101640931B1 (en) * 2016-02-17 2016-07-20 한국기계연구원 automatic PID gain tuning method in articulated robot

Also Published As

Publication number Publication date
KR19980084979A (en) 1998-12-05

Similar Documents

Publication Publication Date Title
Li et al. A new methodology for designing a fuzzy logic controller
Das et al. Performance comparison of optimal fractional order hybrid fuzzy PID controllers for handling oscillatory fractional order processes with dead time
US7219087B2 (en) Soft computing optimizer of intelligent control system structures
CN111045326A (en) Tobacco shred drying process moisture prediction control method and system based on recurrent neural network
Sharma Designing and modeling fuzzy control Systems
Khan et al. Neufuz: Neural network based fuzzy logic design algorithms
CN116914751B (en) Intelligent power distribution control system
CN113325721A (en) Model-free adaptive control method and system for industrial system
Kim et al. Design of an optimal fuzzy logic controller using response surface methodology
KR100310609B1 (en) Coefficient adjusting apparatus of proportional integral derivative(pid) controller
Babuška et al. Promising fuzzy modeling and control methodologies for industrial applications
Gholamian et al. A hybrid intelligent system for multiobjective decision making problems
CN114202063A (en) Fuzzy neural network greenhouse temperature prediction method based on genetic algorithm optimization
Motter et al. Predictive multiple model switching control with the self‐organizing map
Pan et al. Enhancement of fuzzy PID controller with fractional calculus
Gebreel An adaptive interactive multi-objective optimization approach based on decision neural network
CN112733372B (en) Fuzzy logic strong tracking method for load modeling
Abraham Beyond integrated neuro-fuzzy systems: Reviews, prospects, perspectives and directions
JP2781886B2 (en) Control device and control method
Yahiaoui et al. Constrained genetic algorithm-based computer-aided control system design-fixed versus variable size population
Melin et al. Two-rule-based linguistic fuzzy controllers
KR0121103B1 (en) Control method and apparatus for multi-reference expert type fuzzy logic
Navarro et al. Fuzzy logic implementation of industrial controllers
Medel et al. Real-time neuro-fuzzy digital filtering: a technical scheme
Ibeas et al. Representations of Multi‐Model Based Controllers by Using Artificial Intelligence Tools

Legal Events

Date Code Title Description
A201 Request for examination
E701 Decision to grant or registration of patent right
GRNT Written decision to grant
FPAY Annual fee payment

Payment date: 20130515

Year of fee payment: 13

FPAY Annual fee payment

Payment date: 20140620

Year of fee payment: 14

FPAY Annual fee payment

Payment date: 20150702

Year of fee payment: 15

LAPS Lapse due to unpaid annual fee