PID Tuning via Desired Step Response Curve Fitting

기술 요약: 원하는 단계 응답 곡선 피팅을 통한 PID 조정 (PID-SRCF)

문제 제기
비례 - 적분 - 미분 (PID) 제어기의 조정은 제어 공학의 핵심 과제로 남아 있으며, 정착 시간, 오버슈트, 강인성 등 특정 시간 영역 성능 기준을 충족시키기 위해 계수 ($K_p, K_i, K_d$) 를 계산해야 합니다. 지그글러 - 니콜스, 람다 조정, 극점 배치, MATLAB 의 독점 도구 등 수많은 분석적 및 경험적 방법이 존재하지만, 이러한 접근 방식들은 종종 특정 플랜트 모델 (예: 시간 지연이 있는 1 차 시스템), 휴리스틱 규칙, 또는 사용자 정의 과도 사양과 정확히 일치하지 않을 수 있는 암시적 목표 응답에 의존합니다. 또한 극점 배치와 같은 분석적 방법은 폐루프 영점의 영향을 간과하여 예기치 않은 오버슈트를 초래할 수 있으며, 독점 도구는 최적화 목표에 대한 투명성이 부족할 수 있습니다.

방법론
본 논문은 실제 폐루프 단계 응답과 사용자 정의된 원하는 단계 응답 사이의 $L_2$-노름 오차를 직접 최소화하는 수치 최적화 접근법인 PID 단계 응답 곡선 피팅 (PID-SRCF) 방법을 제안합니다.

• 목적 함수: 이 방법의 핵심은 시뮬레이션 시간 $T_{sim}$ 동안 원하는 궤적 $y_{desired}(t)$과 실제 시스템 응답 $y_{PID}(t)$ 사이의 제곱 평균 제곱근 오차 (RMSE) 를 최소화하는 제약 조건이 있는 비선형 최적화 문제입니다:
  $$\min_{K_p, K_i, K_d} \sqrt{\int_{0}^{T_{sim}} |y_{desired}(t) - y_{PID}(t)|^2 dt}$$
• 원하는 응답 정의: 사용자는 비진동 거동을 위한 1 차 시스템 (예: 시간 지연이 있는 1 차 시스템) 과 같은 전달 함수로 목표 응답을 정의하거나, 2 차 시스템에 대한 성능 매개변수 (정착 시간 $T_s$ 및 백분율 오버슈트 $PO$) 를 지정할 수 있습니다. 알고리즘은 이에 해당하는 감쇠비 ($\zeta$) 와 고유 진동수 ($\omega_n$) 를 자동으로 계산합니다.
• 최적화 알고리즘: 이 문제는 순차 2 차 계획법 (SQP) 알고리즘을 사용하는 MATLAB 의 fmincon 솔버를 통해 해결됩니다. 최적화는 양수 PID 게인 ($K_p, K_i, K_d \geq 0$) 을 보장하도록 제약되며, 편향을 피하기 위해 0 에서 초기화됩니다.
• 검증: 최적화 후, 극점 위치를 검사하여 폐루프 안정성을 확인하고 최대 감도 ($M_s$) 지표를 통해 강인성을 평가합니다.

주요 기여

1. 통합 프레임워크: 본 논문은 PID-SRCF 가 고전적 조정 기법을 포괄하는 "메타 방법"으로 작용한다고 주장합니다. 특정 원하는 전달 함수를 선택함으로써, 이 방법은 지그글러 - 니콜스 (4 분의 1 감쇠 응답을 목표로 함), 람다 조정 (FOTD 응답을 목표로 함), 그리고 극점 배치 (임계 감쇠 응답을 목표로 함) 의 결과를 재현할 수 있습니다.
2. 명시적 사양: 전통적인 방법에서 목표 응답이 조정 규칙에 암시적으로 내재되어 있는 것과 달리, PID-SRCF 는 엔지니어가 원하는 폐루프 동역학을 명시적으로 정의할 수 있게 하여 설계 사양과 해법 방법을 분리합니다.
3. 오픈 소스 구현: 저자들은 독점 블랙박스 도구와 대조적으로 투명한 오픈 소스 MATLAB 구현을 제공하여 최적화 로직의 검토 및 수정을 가능하게 합니다.
4. 강인성 분석: 본 연구는 비오버슈트 응답 (1 차 또는 임계 감쇠 2 차) 을 지정하는 것이 본질적으로 최대 감도 ($M_s$) 를 견고한 산업 임계값 ($1 < M_s < 2$) 내에서 제한함을 분석적으로 증명합니다. 구체적으로, 임계 감쇠 2 차 목표는 $M_s$를 $2/\sqrt{3} \approx 1.155$로 산출합니다.

결과
다양한 시스템 차수와 조정 시나리오에 걸쳐 비교 평가가 수행되었습니다:

• 지그글러 - 니콜스 사례: 3 차 프로세스의 경우, PID-SRCF 는 지그글러 - 니콜스 반응 곡선 방법과 거의 동일한 PID 계수 및 성능 지표 (정착 시간, 오버슈트, IAE) 를 달성하여 지그글러 - 니콜스가 암시적으로 특정 2 차 응답 형태를 목표로 한다는 것을 확인했습니다.
• 람다 조정 사례: 1 차 시간 지연 (FOTD) 프로세스의 경우, PID-SRCF 로 유도된 PI 계수는 극점 - 영점 상쇄가 포함된 람다 조정 방법의 계수와 밀접하게 일치하여 표준 산업 관행을 재현하는 능력을 검증했습니다.
• 극점 배치 사례: 1 차 프로세스에 대해 비오버슈트 응답을 목표로 할 때, PID-SRCF 는 분석적 극점 배치보다 낮은 게인 값을 산출했습니다. 폐루프 영점을 무시하는 분석적 방법은 9.2% 의 상당한 오버슈트를 초래한 반면, PID-SRCF 는 2.5% 의 오버슈트만으로 목표를 달성했습니다.
• 고차 시스템: 3 차 프로세스의 경우, PID-SRCF 는 MATLAB 의 독점 "제어 시스템 조정기"와 거의 동일한 PI 계수를 산출했습니다 (네 번째 유효 숫자만 다름). 이러한 유사성은 두 도구 모두 폐루프 응답과 참조 모델 사이의 편차를 최소화하기 때문이며, PID-SRCF 는 특허받은 솔버 (systune) 가 아닌 오픈 소스 솔버 (fmincon) 를 통해 이를 달성한다는 점이 기인합니다.

의의 및 주장
본 논문은 PID-SRCF 가 휴리스틱 분석 방법과 독점 소프트웨어 모두에 대한 유연하고 계산 효율적인 대안을 제공한다고 주장합니다. 그 주요 의의는 다음과 같은 능력에 있습니다:

• 원하는 응답 모델을 단순히 변경함으로써 여러 개의 서로 다른 조정 방법을 단일 통합 알고리즘으로 대체할 수 있습니다.
• 고차 시스템에서 실패할 수 있는 근사에 의존하지 않고 과도 응답 특성 (오버슈트 및 정착 시간) 에 대한 정밀한 제어를 제공합니다.
• 비오버슈트 목표를 선택할 때 감도를 자연스럽게 제한함으로써 강인성을 보장합니다.
• 오픈 소스 구현을 통해 투명성과 접근성을 제공하여 실무자가 라이선스 제약 없이 조정 로직을 이해하고 수정할 수 있게 합니다.

저자들은 근본적인 최적화 원칙 (참조 모델에 대한 오차 최소화) 이 새로운 것은 아니지만, PID-SRCF 의 명시적 공식화와 오픈 소스 가용성이 기존 방법을 효과적으로 대체할 수 있는 실용적이고 통합된 산업용 PID 조정 도구를 제공한다고 결론지었습니다.