KR102605904B1 - 임의 교란을 가진 강성 우주선을 위한 자세 추적 제어 방법 및 장치 - Google Patents
임의 교란을 가진 강성 우주선을 위한 자세 추적 제어 방법 및 장치 Download PDFInfo
- Publication number
- KR102605904B1 KR102605904B1 KR1020230045514A KR20230045514A KR102605904B1 KR 102605904 B1 KR102605904 B1 KR 102605904B1 KR 1020230045514 A KR1020230045514 A KR 1020230045514A KR 20230045514 A KR20230045514 A KR 20230045514A KR 102605904 B1 KR102605904 B1 KR 102605904B1
- Authority
- KR
- South Korea
- Prior art keywords
- value
- attitude
- control
- spacecraft
- disturbance
- Prior art date
Links
- 238000000034 method Methods 0.000 title claims abstract description 36
- 230000007423 decrease Effects 0.000 claims description 9
- 238000004364 calculation method Methods 0.000 claims description 7
- 238000004590 computer program Methods 0.000 claims description 3
- 229910052709 silver Inorganic materials 0.000 claims description 2
- 239000004332 silver Substances 0.000 claims description 2
- 230000007613 environmental effect Effects 0.000 abstract description 7
- 238000010586 diagram Methods 0.000 description 35
- 230000006870 function Effects 0.000 description 27
- 239000011159 matrix material Substances 0.000 description 10
- 238000005516 engineering process Methods 0.000 description 8
- 238000004088 simulation Methods 0.000 description 5
- 230000001052 transient effect Effects 0.000 description 3
- 230000004888 barrier function Effects 0.000 description 2
- 230000006641 stabilisation Effects 0.000 description 2
- 238000011105 stabilization Methods 0.000 description 2
- 230000000087 stabilizing effect Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000003190 augmentative effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 230000000593 degrading effect Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 239000000428 dust Substances 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 239000000446 fuel Substances 0.000 description 1
- 230000014509 gene expression Effects 0.000 description 1
- 229910052760 oxygen Inorganic materials 0.000 description 1
- 230000001144 postural effect Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000012887 quadratic function Methods 0.000 description 1
- 230000002441 reversible effect Effects 0.000 description 1
- 238000004904 shortening Methods 0.000 description 1
- 229910052717 sulfur Inorganic materials 0.000 description 1
- 230000001960 triggered effect Effects 0.000 description 1
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64G—COSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
- B64G1/00—Cosmonautic vehicles
- B64G1/22—Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
- B64G1/24—Guiding or controlling apparatus, e.g. for attitude control
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64G—COSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
- B64G1/00—Cosmonautic vehicles
- B64G1/22—Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
- B64G1/24—Guiding or controlling apparatus, e.g. for attitude control
- B64G1/244—Spacecraft control systems
- B64G1/245—Attitude control algorithms for spacecraft attitude control
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Chemical & Material Sciences (AREA)
- Combustion & Propulsion (AREA)
- Radar, Positioning & Navigation (AREA)
- Aviation & Aerospace Engineering (AREA)
- Automation & Control Theory (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
본 발명은 임의 교란을 가진 강성 우주선을 위한 자세 추적 제어 방법 및 장치에 관한 것으로서, 더욱 상세하게는 시스템의 불확실성과 환경 교란을 가진 우주선에 대하여 사전 정의된 시간 내에 기준 자세를 추적하여 제어하는 방법 및 그러한 방법을 수행하는 장치에 관한 것이다.
본 발명에 의하면, 시스템의 불확실성, 외부 장애 및 가상 제어의 시간 파생을 포함한 전체 장애에 대한 초기 추정 오차가 무한대인 경향이 있더라도, 사전 정의된 시간 동안 정확하게 우주선에 대한 기준 자세를 추적하여 제어할 수 있는 방법 및 장치를 제공한다.
본 발명에 의하면, 시스템의 불확실성, 외부 장애 및 가상 제어의 시간 파생을 포함한 전체 장애에 대한 초기 추정 오차가 무한대인 경향이 있더라도, 사전 정의된 시간 동안 정확하게 우주선에 대한 기준 자세를 추적하여 제어할 수 있는 방법 및 장치를 제공한다.
Description
본 발명은 임의 교란을 가진 강성 우주선을 위한 자세 추적 제어 방법 및 장치에 관한 것으로서, 더욱 상세하게는 시스템의 불확실성과 환경 교란을 가진 우주선에 대하여 사전 정의된 시간 내에 기준 자세를 추적하여 제어하는 방법 및 그러한 방법을 수행하는 장치에 관한 것이다.
비선형 항, 관성 행렬 불확실성, 마찰 토크 및 외부 교란은 우주선 자체 시스템의 불확실한 부분으로 인식된다. 불확실한 부분이 존재하는 상황에서 원하는 성능을 달성하기 위해 자세 추적 컨트롤러를 개발하는 것은 어려운 작업이다.
지금까지 슬라이딩 모드 제어, 백스텝 제어, 최적 제어, 이벤트 트리거 제어, 모델 예측 제어과 같은 수많은 관찰자 기반 자세 제어가 만들어 졌지만, 이러한 제어기들은 교란이 일정하거나 무시할 만한 변화를 가지고 있다는 가정에 기반하며 추정 오차가 무한한 시간 범위 내에서 원점으로 수렴되도록 보장한다.
실용적인 관점에서 볼 때, 어떤 우주선이든 초기 조건이 아무리 크더라도 제한된 시간 내에 제어를 완수해야 하기 때문에 고정 시간 제어가 필수적이다. 수렴 시간은 제어 이득의 선택에 크게 좌우됨에 따라 과도적이고 안정적인 상태에서 성능을 저하시키지 않으면서 특정한 안착 시간을 달성하는 것은 필수적이나, 그와 같이 불확실한 부분이 존재하는 상황에서 고정된 시간 내에 안정된 자세 제어를 수행하는 것은 매우 어려운 실정이다.
본 발명은 이와 같은 문제점을 해결하기 위해 창안된 것으로서, 시스템의 불확실성, 외부 장애 및 가상 제어의 시간 파생을 포함한 전체 장애에 대한 초기 추정 오차가 무한대인 경향이 있더라도, 사전 정의된 시간 동안 정확하게 우주선에 대한 기준 자세를 추적하여 제어할 수 있는 방법 및 장치를 제공하는데 그 목적이 있다.
이와 같은 목적을 달성하기 위하여 본 발명에 따른 임의 교란을 가진 강성 우주선을 위한 자세 추적 제어 방법은, (a) 우주선에 대하여 측정된 자세값(이하 '측정 자세값'이라 한다)과, 목표 자세값과 차이값(이하 '자세 오차'라 한다)를 산출하는 단계; (b) 상기 자세 오차 및, 성능 함수로부터 가상 자세 제어값을 산출하는 단계; (c) 상기 측정 자세값 및, 우주선의 자세를 제어하기 위한 제어 토크의 피드백 값으로부터 우주선에 대하여 작용하는 외란(disturbance)의 추정값을 산출하는 단계; 및, (d) 상기 측정 자세값, 상기 가상 자세 제어값 및 상기 외란 추정값으로부터 새로운 제어 토크 값을 산출하는 단계를 포함하고, 상기 성능함수는, 시간 t=0으로부터 감소하여, 유한 시간 내에 특정 값으로 수렴하고, 이후 일정한 값으로 유지된다.
상기 단계(d) 이후, (e) 상기 단계(d)에서 산출된 제어 토크 값을, 우주선의 자세를 제어하는 제어부로 전달하는 단계를 더 포함할 수 있다.
상기 성능함수는, 시간 t=0으로부터 감소하여, 유한 시간 내에 특정 값으로 수렴하고, 이후 일정한 값으로 유지될 수 있다.
상기 가상 자세 제어값을 구하는 파라미터에는, 상기 제어 토크값 산출시, 해당 제어 토크값이 특정 수렴시간(이하, '제어 수렴시간'이라 한다) 내에 특정값으로 수렴하여 이후 일정한 값으로 유지되는 경우, 그 제어 수렴시간 값을 더 포함할 수 있다.
상기 가상 자세 제어값은, 로 산출될 수 있고, 상기 ρi는 성능함수 값, , , β는 제어 파라미터이고, 이며, 는 제어 파라미터, Tc는 상기 제어 수렴시간이다.
상기 외란 추정값의 산출은, (c1) 미지 입력 선형 시스템(unknown input linear system, UILS)을 정의하는 단계; (c2) 상기 UILS의 상태변수의 추정값의 미분 식과, 상기 UILS 식으로부터 상기 UILS에 포함된 미지 입력값(unknown input)을 구하는 단계; 및, (c3) 상기 미지 입력값으로부터 외란 추정값을 산출하는 단계를 포함할 수 있다.
상기 UILS 식은, 일 수 있고, 여기서 는 상기 UILS의 상태변수, k1, k2는 상수 파라미터, 은 미지 입력값이고, 상기 상기 UILS의 상태변수의 추정값의 미분 식은,
이며, 여기서 로서 추정 오차이며, 상기 외란 추정값 은,
로 산출되며, 여기서 이고,
이며, 이때 로서, 목표 자세값()과 측정 자세값()의 오차를 나타내며, 이고, 상기 추정 오차는 수렴시간(이하 '외란 관찰 수렴시간'이라 한다) 내에 특정값으로 수렴함으로써 상기 외란 관찰 수렴시간 내에 상기 외란 추정값이 산출될 수 있다.
상기 제어 토크값은,
로 산출될 수 있고, 여기서 는 상기 단계(c)에서 산출된 외란 추정값이다.
상기 제어 토크는 상기 제어 수렴시간 내에 특정값으로 수렴할 수 있다.
상기 제어 수렴시간 는, 상기 외란 추정값 산출시의 추정 오차의 수렴시간(이하 '외란 관찰 수렴시간'이라 한다) 보다 큰 값일 수 있다.
본 발명의 다른 측면에 따르면, 임의 교란을 가진 강성 우주선의 자세를 추적 제어하기 위한 장치(이하 '자세 추적 제어 장치'라 한다)는, 적어도 하나의 프로세서; 및 컴퓨터로 실행가능한 명령을 저장하는 적어도 하나의 메모리를 포함하되, 상기 적어도 하나의 메모리에 저장된 상기 컴퓨터로 실행가능한 명령은, 상기 적어도 하나의 프로세서에 의하여, (a) 우주선에 대하여 측정된 자세값(이하 '측정 자세값'이라 한다)과, 목표 자세값과 차이값(이하 '자세 오차'라 한다)를 산출하는 단계; (b) 상기 자세 오차 및, 성능 함수로부터 가상 자세 제어값을 산출하는 단계; (c) 상기 측정 자세값 및, 우주선의 자세를 제어하기 위한 제어 토크의 피드백 값으로부터 우주선에 대하여 작용하는 외란(disturbance)의 추정값을 산출하는 단계; 및, (d) 상기 측정 자세값, 상기 가상 자세 제어값 및 상기 외란 추정값으로부터 새로운 제어 토크 값을 산출하는 단계가 실행되도록 하고, 상기 성능함수는, 시간 t=0으로부터 감소하여, 유한 시간 내에 특정 값으로 수렴하고, 이후 일정한 값으로 유지된다.
본 발명의 또 다른 측면에 따르면, 임의 교란을 가진 강성 우주선의 자세를 추적 제어하기 위한, 컴퓨터로 판독 가능한 비일시적 저장 매체에 저장된 컴퓨터 프로그램은, 비일시적 저장 매체에 저장되며, 프로세서에 의하여, (a) 우주선에 대하여 측정된 자세값(이하 '측정 자세값'이라 한다)과, 목표 자세값과 차이값(이하 '자세 오차'라 한다)를 산출하는 단계; (b) 상기 자세 오차 및, 성능 함수로부터 가상 자세 제어값을 산출하는 단계; (c) 상기 측정 자세값 및, 우주선의 자세를 제어하기 위한 제어 토크의 피드백 값으로부터 우주선에 대하여 작용하는 외란(disturbance)의 추정값을 산출하는 단계; 및, (d) 상기 측정 자세값, 상기 가상 자세 제어값 및 상기 외란 추정값으로부터 새로운 제어 토크 값을 산출하는 단계가 실행되도록 하는 명령을 포함하고, 상기 성능함수는, 시간 t=0으로부터 감소하여, 유한 시간 내에 특정 값으로 수렴하고, 이후 일정한 값으로 유지된다.
본 발명의 또 다른 측면에 따르면, 임의 교란을 가진 강성 우주선의 자세를 추적 제어하기 위한 장치(이하 '자세 추적 제어 장치'라 한다)는, 우주선에 대하여 측정된 자세값(이하 '측정 자세값'이라 한다)과, 목표 자세값과 차이값(이하 '자세 오차'라 한다)를 산출하는 자세 오차 산출부; 상기 자세 오차 및, 성능 함수로부터 가상 자세 제어값을 산출하는 가상 제어부; 상기 측정 자세값 및, 우주선의 자세를 제어하기 위한 제어 토크의 피드백 값으로부터 우주선에 대하여 작용하는 외란(disturbance)의 추정값을 산출하는 외란 관찰부; 및, 상기 측정 자세값, 상기 가상 자세 제어값 및 상기 외란 추정값으로부터 새로운 제어 토크 값을 산출하는 자세 추적 제어부를 포함하고, 상기 성능함수는, 시간 t=0으로부터 감소하여, 유한 시간 내에 특정 값으로 수렴하고, 이후 일정한 값으로 유지된다.
본 발명에 의하면, 시스템의 불확실성, 외부 장애 및 가상 제어의 시간 파생을 포함한 전체 장애에 대한 초기 추정 오차가 무한대인 경향이 있더라도, 사전 정의된 시간 동안 정확하게 우주선에 대한 기준 자세를 추적하여 제어할 수 있는 방법 및 장치를 제공하는 효과가 있다.
도 1은 강성 우주선 자세 시스템에 사용되는 좌표 프레임을 나타내는 도면.
도 2는 성능 함수 ρ(t)의 일 실시예를 나타내는 도면.
도 3은 본 발명의 자세 추적 제어 장치의 구성을 나타내는 블럭 다이어그램.
도 4는 도 3의 자세 추적 제어 장치가 수행하는 자세 추적 제어 방법의 순서도.
도 5는 본 발명에 의해 수행된 결과로서, 다양한 초기조건 설정시 자세 추적 오차를 나타내는 도면.
도 6은 본 발명에 의해 수행된 결과로서, 다양한 초기조건 설정시 각속도 추적 오차를 나타내는 도면.
도 7은 본 발명에 의해 수행된 결과로서, 다양한 초기조건 설정시 시변 이득을 나타내는 도면.
도 8은 본 발명에 의해 수행된 결과로서, 다양한 초기조건 설정시 제어 입력을 나타내는 도면.
도 9는 본 발명에 의해 수행된 결과로서, 다양한 수렴시간 설정시 자세 추적 오차를 나타내는 도면.
도 10은 본 발명에 의해 수행된 결과로서, 다양한 수렴시간 설정시 각속도 추적 오차를 나타내는 도면.
도 11은 본 발명에 의해 수행된 결과로서, 다양한 수렴시간 설정시 시변 이득 λ를 나타내는 도면.
도 12는 본 발명에 의해 수행된 결과로서, 다양한 수렴시간 설정시 제어 입력 τ를 나타내는 도면.
도 13은 본 발명과 기존 기술의 간의 성능 비교를 나타내는 도면.
도 14는 본 발명과 기존 기술의 간의 자세 추적 오차를 비교하여 나타내는 도면.
도 15는 본 발명과 기존 기술의 간의 각속도 추적 오차를 비교하여 나타내는 도면.
도 16은 본 발명과 기존 기술의 간의 제어 입력 τ를 비교하여 나타내는 도면.
도 17은 추정 오차를 나타내는 도면.
도 18은 총 교란과 그 추정값을 나타내는 도면.
도 2는 성능 함수 ρ(t)의 일 실시예를 나타내는 도면.
도 3은 본 발명의 자세 추적 제어 장치의 구성을 나타내는 블럭 다이어그램.
도 4는 도 3의 자세 추적 제어 장치가 수행하는 자세 추적 제어 방법의 순서도.
도 5는 본 발명에 의해 수행된 결과로서, 다양한 초기조건 설정시 자세 추적 오차를 나타내는 도면.
도 6은 본 발명에 의해 수행된 결과로서, 다양한 초기조건 설정시 각속도 추적 오차를 나타내는 도면.
도 7은 본 발명에 의해 수행된 결과로서, 다양한 초기조건 설정시 시변 이득을 나타내는 도면.
도 8은 본 발명에 의해 수행된 결과로서, 다양한 초기조건 설정시 제어 입력을 나타내는 도면.
도 9는 본 발명에 의해 수행된 결과로서, 다양한 수렴시간 설정시 자세 추적 오차를 나타내는 도면.
도 10은 본 발명에 의해 수행된 결과로서, 다양한 수렴시간 설정시 각속도 추적 오차를 나타내는 도면.
도 11은 본 발명에 의해 수행된 결과로서, 다양한 수렴시간 설정시 시변 이득 λ를 나타내는 도면.
도 12는 본 발명에 의해 수행된 결과로서, 다양한 수렴시간 설정시 제어 입력 τ를 나타내는 도면.
도 13은 본 발명과 기존 기술의 간의 성능 비교를 나타내는 도면.
도 14는 본 발명과 기존 기술의 간의 자세 추적 오차를 비교하여 나타내는 도면.
도 15는 본 발명과 기존 기술의 간의 각속도 추적 오차를 비교하여 나타내는 도면.
도 16은 본 발명과 기존 기술의 간의 제어 입력 τ를 비교하여 나타내는 도면.
도 17은 추정 오차를 나타내는 도면.
도 18은 총 교란과 그 추정값을 나타내는 도면.
이하 첨부된 도면을 참조로 본 발명의 바람직한 실시예를 상세히 설명하기로 한다. 이에 앞서, 본 명세서 및 청구범위에 사용된 용어나 단어는 통상적이거나 사전적인 의미로 한정해서 해석되어서는 아니되며, 발명자는 그 자신의 발명을 가장 최선의 방법으로 설명하기 위해 용어의 개념을 적절하게 정의할 수 있다는 원칙에 입각하여 본 발명의 기술적 사상에 부합하는 의미와 개념으로 해석되어야만 한다. 따라서, 본 명세서에 기재된 실시예와 도면에 도시된 구성은 본 발명의 가장 바람직한 일 실시예에 불과할 뿐이고 본 발명의 기술적 사상을 모두 대변하는 것은 아니므로, 본 출원시점에 있어서 이들을 대체할 수 있는 다양한 균등물과 변형예들이 있을 수 있음을 이해하여야 한다.
먼저, 본 발명의 특징을 개략적으로 설명하면 다음과 같다.
1. 제약된 사전 정의된 시간(predefined-time) 백스텝 제어를 수행하도록 설계되었다. 시스템 궤적(system trajectory)이 초기 조건과 제어 이득에 관계없이 유한한 시간 내에 0으로 수렴하는 것을 보장한다.
2. 커맨드 필터나 조각별(piece-wise) 연속 함수를 사용하지 않고도 특이점 문제를 해결하기 위하여, 가상 자세 제어에 이차 함수가 사용된다.
3. 자세 추적 오차에 대한 제약을 만족시키기 위해 가상 자세 제어에 탄형 BLF(tan-type barrier Lyapunov function)가 포함된다. 따라서 제약 조건은 간단한 부등식 조건을 만족시키는 방식으로 처리한다.
4. 사전 정의된 시간(predefined-time) 외란 관측부를 사용하여 가상 제어의 미분과 전체 외란(disturbance)의 미분을 정확하게 재구성한다. 외란 관측부의 안정화 시간은 초기 추정 오차가 무한대인 경향이 있더라도 항상 사전 정의된 시간보다 작다. 또한 이 시간은 외란 관측부의 조정 가능한 이득으로 결정된다.
이하에서는 각 도면을 참조하여 본 발명에 대하여 설명해 나가기로 한다.
도 1은 강성 우주선 자세 시스템에 사용되는 좌표 프레임을 나타내는 도면이다.
원형 궤도에서 움직이는 강체 우주선의 자세 운동 방정식은 자세 운동학과 우주선 역학으로 표현할 수 있다
1-1. 자세 운동학
강성 우주선 자세 시스템의 경우, 그림 1에 도시된 3개의 주요 좌표 프레임이 사용된다. 이 그림에서 궤도 기준, 관성 및 차체 고정 프레임은 각각 FO, FI 및 FB로 표시된다. 강체 우주선의 궤도 위치를 결정하기 위해 지구 중심 관성 프레임 FI가 사용된다. 오일러 각도는 이해하기 쉽기 때문에 자세 표현과 운동학 방정식을 제시하는 데 사용된다.
를 롤-피치-요() 회전 순서를 통해 얻은 궤도 기준 프레임 FO에 대한 고려된 우주선의 오일러 각도로 하고, 는 관성 기준 프레임 FI에 대한 관성 회전 속도를 나타내며 차체 고정 프레임 FB로 표현된다. 이때, 운동 방정식은 다음과 같이 주어진다.
여기서 벡터 와 행렬 은 각각 다음과 같이 정의된다.
여기서 C와 S는 각각 코사인과 사인을 의미하고, 는 궤도 속도를 나타낸다.
이기 때문에 행렬 가 가역적이기 위해서는 t≥0에 대한 조건 를 만족해야 하며, 여기서 와 은 자연수의 집합이다. 이를 위해 초기 피치 각도θ(0)는 가 되도록 제한되어야 한다. 자세 기동 중에 컨트롤러는 을 얻을 수 있어야 한다.
1-2. 강체 우주선 역학
오일러 모멘트 방정식에 기초하여, 강체 우주선의 역학은 다음과 같이 구할 수 있다.
여기서 는 로 표현된 우주선의 관성 행렬이며, 여기서 와 는 각각 관성 행렬의 공칭 부분(nominal part)과 불확실한 부분을 나타낸다. 게다가, 와 는 각각 제어 토크 및 환경 교란 토크를 나타냅니다. 는 스큐 대칭 행렬(skew-symmetric matrix)을 의미한다.
이때, 불확실한 관성 행렬 는 로 경계가 있으며, 여기서 a1은 양수이지만 알 수 없는 상수이다. 또한 환경 교란(environmental disturbance)은 경계가 있다. 즉, 여기서 a2는 양수이지만 알려지지 않은 상수이다. 또한 오일러 각도와 회전 속도는 측정 가능하며 피드백에 사용할 수 있다고 가정한다.
식(1)의 미분은 다음과 같은 결과가 된다.
우주선 자세 시스템 식(1) 및 식(3)을 다음과 같은 동적 시스템으로 변환하는 것은 간단하다.
여기서,
이다.
자세 추적 오차를 로 정의하자. 여기서 는를 가진 두번 미분 가능한 시변 기준 궤적(twice differentiable time-varying reference trajectory)이다. 이때 그리고 이며, 우주선 자세 시스템 식(5)는 다음과 같이 표현될 수 있다.
여기서 이다.
본 발명에서, 강성 우주선 자세 시스템 식(6)을 위한 사전 정의된 시간(predefined-time) 수렴을 갖는 탄형 BLF 기반 백스텝 제어(tan-type BLF-based backstepping control)를 이용하여 다음과 같은 제어 목표를 달성한다.
1. 우주선 자세는, 시스템의 불확실성 및 환경 교란에도 불구하고, 미리 정의된 시간 내에 기준 자세 궤적 를 추적한다,
2. 과도 상태(예: 최대 오버슈트) 및 정상 상태(예: 최종 추적 오차) 설정에서 자세 추적 오차 x1에 대한 사전 정의된 성능이 미리 지정되어 있다.
다음은, 본 발명의 자세 추적 제어에 대한 본격적인 설명에 앞서, 비선형 시스템의 안정성을 입증하기 위한 양의 확정 랴푸노프 함수(Positive-definite Lyapunov Function)에 대한 몇가지 명제를 기술한다.
명제1) 다음과 같은 비선형 시스템을 고려한다.
다음을 만족하는 양의-확정 랴푸노프 함수가 존재한다고 가정한다.
여기서 및 이다. 이때, 시스템 (7)의 원점은 특정 안착 시간 Tc로 규정된, 시간 안정적인 특징을 가진다.
랴푸노프 안정성의 개념에 기초하여, 양의 확정 랴푸노프 함수(Positive-definite Lyapunov Function) 로 표시된 비선형 시스템(7)의 에너지가 증가하지 않는 경우, 즉, 이면 시스템 상태 x(t)가 발산하지 않고 시스템이 안정적이다. 게다가, 조건이 만족된다면, 시스템은 점근적으로 안정적이고 시간이 무한대인 경향이 있을 때 시스템 상태 x(t)는 0으로 수렴한다. 전통적인 랴푸노프 안정성 이론은 무한한 시간 지평선에 걸쳐 시스템의 안정성을 보장한다는 점에 주목한다. 그러나, 명제1은 시스템 상태 가 유한한 시간 Tc 내에 0으로 수렴할 수 있는 충분한 조건을 제시한다.
명제2) 비선형 시스템(7)을 고려한다. 랴푸노프 함수 가 다음과 같은 부등식을 만족한다면
,
여기서 , 시스템 (7)은 실질적으로 규정된 시간 안정적(prescribed-time stable)이고, 여기서 잔차 집합은 다음과 같이 표현된다
여기서 파라미터 는 인 수렴 시간을 나타낸다. 이기 때문에, 와 이다. 와 로 정의하면, 의 값은 a<b 인지 또는 b<a 인지를 결정한다. 더 정확하게 말하면, 만약 이면 a>b이고 아닐 경우 그 반대이다.
명제3) 모든 와 에 대해, 다음과 같은 부등식이 성립한다:
명제4) 이고 인 경우,
이다.
도 2는 성능 함수 ρ(t)의 일 실시예를 나타내는 도면이다.
먼저 다음과 같은 유한 시간 규정 성능 함수(FTPPF, finite-time prescribed performance function)를 정의한다. 즉,
정의1) 양의 및 비증가 함수 ρ(t)는 이면 유한 시간 규정 성능 함수(FTPPF, finite-time prescribed performance function)이고, 임의의 에 대하여 이며, 여기서 Tf는 시스템 궤적이 정상 상태 영역에 도달하는 데 필요한 시간이고 ρT는 최대 정상 상태 경계이다.
이때, 자세 추적 오차에 대한 성능 요구 사항을 충족하기 위해 다음을 보장해야 한다.
실제로 성능 함수 는 설계자가 일시적(예: 최대 허용 오버슈트) 및 정상 상태(예: 최대 정상 상태 경계)에서 추적 오차에 대해 유리한 성능 특성을 지정할 수 있도록 한다. 함수 (11)는 (정의1)에 주어진 요건을 충족하며, FTPPF로 사용된다.
여기서 양의 상수 ρ0, c1, c2는 Tf, ρT, ρ0의 값에 따라 적절하게 선택된다. 더 정확히 말하면, 와 는 각각 b1, b2가 양의 홀수와 짝수 정수가 되도록 한다.
수식(11)의 예가 도 2에 나와 있다. ρ(0)=25, ρT=1, Tf=8, b1=5, b2=10으로 선택하면 가 나온다. 함수 는 에서 시작하여 t=8에서 궁극값(ultimate value) 로 수렴하여 그 이후에 계속 유지되는 것을 알 수 있다.
설계 절차를 간소화하기 위해 추적 오차와 관련된 변환된 오차 변수를 로 정의한다. 여기서 이다. 따라서, 출력 제약 조건(10)은 모든 t > 0에 대해 를 만족시키는 것과 같다. 새로운 변환 오차에 기초하여, 동적 시스템(6)은 다음과 같이 다시 쓸 수 있다.
도 3은 본 발명의 자세 추적 제어 장치(100)의 구성을 나타내는 블럭 다이어그램이고, 도 4는 도 3의 자세 추적 제어 장치(100)가 수행하는 자세 추적 제어 방법의 순서도이다.
먼저 도 3의 각 구성과 도 4의 순서도를 전체적으로 설명한 이후에, 그 세부적 내용에 대하여 설명하기로 한다.
자세 오차 산출부(110)는 이전 단계에서 측정된 우주선의 자세값(이하, '측정 자세값'이라 한다)과, 목표 자세값의 차이(이하, '자세 오차'라 한다)를 산출한다(S410).
가상 제어부(120)는 이러한 자세 오차 및, 성능 함수(performance function)로부터 가상 자세 제어값(virtual attitude control)을 산출한다(S420).
외란 관찰부(130)는 측정 자세값과, 후술하는 자세 추적 제어부(140)에서 산출하는 제어 토크(control torque) 값으로부터 우주선에 작용할 외란(disturbance)의 추정값(이하, '외란 추정값'이라 한다)을 산출한다(S430).
자세 추적 제어부(140)는 측정 자세값, 가상 자세 제어값 및 외란 추정값 등으로부터, 사전 정의된 시간 내에 우주선의 자세를 안정화 시키기 위한 제어 토크(control torque) 값을 산출한다(S440).
우주선의 제어부(150)는, 입력 받은 제어 토크 값을 이용하여 우주선의 자세 제어를 수행한다(S450). 이때 우주선에는 외부의 복합적인 교란(composite interference)이 가해지며, 이러한 교란 입력 및 제어 토크 값에 의한 자세 제어 후의 우주선에 대하여 직접 측정(S460)된 자세값은, 전술한 바와 같이 다시 자세 오차 산출부(110), 자세 추적 제어부(140), 외란 관찰부(130) 등으로 자세 추적 제어를 위하여 제공된다(S470).
이하에서는, 지금까지 설명한 도 3의 각 구성과 도 4의 순서도에 대한 전체적인 내용에 포함되어 있는 구체적인 내용을 상세하게 설명하기로 한다.
앞서 명제1, 명제2를 통하여 랴푸노프 함수 V에 대하여 기술한 바 있다. 이를 다시한번 간략히 정리하면, 랴푸노프 함수의 경우 다음과 같은 3가지의 조건을 만족하면 asymptotically stable하게 된다. 예를 들어 다음의 V함수는 x에 대한 함수이고 다음과 같은 3가지 조건을 만족하게 되면 x는 0으로 수렴하게 된다.
특히, 본 발명에서는 명제1, 명제2와 같은 조건을 만족시키는 경우 x가 0으로 수렴할 뿐만 아니라, 특정 시간, 예를 들어 Tc 내에 수렴하도록 한다는 것을 보여준다.
이하에서는 랴푸노프 함수를 특정한 함수로 정하고, 그러한 랴푸노프 함수에 의하여 식 (12)의 λ, 즉, 전술한 바와 같이 인 λ가 특정 시간 내에 0으로 수렴함울 보임으로써, 이에 따라 결국 자세 오차인 가 그 특정 시간 내에 0으로 수렴하게 된다는 것을 보이게 된다.
- Tan-type BLF(barrier Lyapunov function) 자세 제어 설계
이하에서는 백스텝 제어 접근법에 기초한 제한된 자세 제어를 설명한다. 일반적으로 백스텝 전략은 랴푸노프 함수의 선택과 피드백 제어의 설계를 혼합하는 재귀적 절차이다. 그것은 전체 시스템에 대한 설계 문제를 (심지어 스칼라) 하위 시스템에 대한 일련의 설계 문제로 분해한다. 또한, 백스텝은 슬라이딩 모드 제어, 피드백 선형화 및 이득 스케줄링과 같은 다른 방법에서 발생하는 것보다 덜 제한적인 조건에서 안정화 및 궤적 추적 제어 문제를 모두 해결할 수 있다.
1단계: 랴푸노프 함수는 임의로 정할 수 있는데, 본 발명에서는 이러한 랴푸노프 함수를 tan-type으로 식(13)과 같이 구성한다.
여기서 tan(·)은 탄젠트 함수를 나타낸다. 랴푸노프 함수(13)는 집합 에서 미분 가능하다. 의 도함수는 다음과 같다.
여기서 와 는 i번째 대각선 요소가 0이 아닌 유일한 요소가 되도록 일관된 차원(consistent dimension)을 가진 대각선 행렬을 보여준다. 또한, 는 중간 오차(intermediate error)를 나타내고, 는 다음과 같이 표현되는 가상 자세 제어(virtual attitude control) 값을 나타낸다.
여기서 이고, 는,
로서, 이고, Tc는 수렴시간을 나타낸다.
이러한 가상 자세 제어값(15)를 (14)로 대입함으로써, 다음을 획득하게 된다.
전술한 도 3의 가상 제어부(120)는, 자세 오차 산출부(110)에서 산출한 자세 오차 및 성능함수 ρi(t)에 의하여 오차 변수 λi(t)를 로 구하고, 이로부터 를 구하여, 위 조건을 만족하는 제어 파라미터인 , 및 수렴시간 Tc로부터 식(15)에 의하여 가상 자세 제어값 α를 구한다(S420).
본 발명의 가상 제어값의 차별점은 다음과 같다.
가상 제어값 (15)의 도함수는 다음과 같다.
여기서,
이고,
이다. 따라서
이다.
종래의 백스테핑(backstepping) 방법에 의하면, 가상 자세 제어는,
와 같이 설계되었다. 그러나 이 경우의 에는 일 때
이므로 특이점(singularity)의 문제가 발생한다. 그러나 본 발명에서는 식 (15)의 시간 도함수에 항을 포함하고, 이로써 특이점(singularity)의 문제는 제거되는 장점이 있다.
2단계: 다음과 같이 증강 랴푸노프 함수를 구성한다.
여기서 이다. V2의 도함수를 구하면,
전술한 바와 같이, 이므로,
이 성립한다. 여기서 는 전체 외란값인데, 이는 후술하는 바와 같이 외란 관찰부(130)에 의해 추정된다.
실제 제어 입력, 즉 제어 토크값은
에 의해 구해지는데, 이러한 제어 토크값 (21)은, 도 3 및 도 4에 대한 개략적 설명 부분에서 전술한 바와 같이, 자세 추적 제어부(140)에 의해 산출된다(S440). 여기서 이고, 는 외란 추정값이다.
즉, 자세 추적 제어부(140)는 측정 자세값(), 가상 자세 제어값(α) 및 외란 추정값() 등으로부터, 사전 정의된 시간 내에 우주선의 자세를 안정화 시키기 위한 제어 토크(control torque) 값을 산출한다(S440).
즉, z 는 가상 자세 제어값(α)을 이용하여 으로부터 구해지고, 외란 추정값()은 후술하는 바와 같이 외란 관찰부(130)에 의해 산출된다.
또한 , 등은 가상 제어부(120)의 가상 자세 제어값(α) 산출시 설명한 바와 같이, 자세 오차 와 관계된 값으로서, 로부터 직접 구할 수도 있고, 또는 가상 제어부(120)로부터 계산된 값을 넘겨받아 식(21)에서 사용할 수도 있다. 식(21)에 의해 구해진 제어 토크(τ)는 우주선의 제어부(150)로 전달된다.
식(21)을 식 (19)에 대입하면,
이고, 여기서 이다.
이하에서는 외란 관찰부(130)가 외란 추정값 를 산출하는 과정(S430)을 설명한다.
이 부분에서는 식(20)의 총 교란 Td를 정확하게 추정하기 위해 사전 정의된 시간(predefined-time) 내에 외란을 관찰하는 역할을 외란 관찰부(130)가 수행한다. 외란 관찰부(130)는 초기 추정 오차와 관계없이 미리 지정된 시간 동안 추정 오차가 원점으로 수렴되도록 보장한다. 또한, 이 추정 시간은 관찰자에서 조정 가능한 이득으로 명시적으로 결정된다.
먼저, 식(20)을 다음과 같이 다시 쓰면 다음과 같다.
여기서 이고, k1은 양의 상수이다.
다음의 보조 시스템(auxiliary system)은 미지 입력 선형 시스템(UILS, unknown input linear system)을 구성하고 미지 입력 관찰자의 이론적 개념을 사용하기 위해 정의된다.
에 의해 z와 za 사이의 불일치(discrepancy)를 정의한다. 그러면, UILS는 다음과 같이 주어진다.
여기서 는 UILS의 상태, y는 시스템 출력, k2 > 0이다. 구성된 UILS(25)에서 은 미지 입력(unknown input)의 역할을 하며, 다음과 같은 미리 정의된 시간(predefined-time) 외란 관찰자에 의해 추정된다.
여기서 는 추정 오차(estimation error)로서 로 표현되고, 이며, Te는 미리 정의된 시간(predefined-time) 외란 관찰자의 안착 시간(settling time)을 나타낸다.
추정 오차의 역학은 다음과 같이 표현할 수 있다.
연속적으로 미분 가능한 양의 확정 랴푸노프 함수(positive definite Lyapunov function)를 로 정의한다. 이 경우 의 시간 도함수는 다음과 같이 주어진다.
부등식 (28)은 다음과 같이 다시 쓸 수 있다.
여기서 ε는 음이 아닌 상수이다.
미분방정식 (29)를 풀면
와 같이 나타내어진다. 따라서 다음 부등식을 얻는다.
여기서 이고 임을 주목해야 한다. 이것은 초기 추정 오차가 무한대에 가깝더라도, 수렴시간은 Te를 넘지 않는다는 것을 보장한다.
따라서 미지 입력 은 미리 정해진 시간 Te 내에 다음 식과 같이 에 의해 추정된다.
을 기반으로, 총 외란 는,
로 재구성된다.
를 정의함으로써, 식 (33)을 이용하여,
를 얻는다.
- 안정성 분석(Stability Analysis)
지금까지 설명한 바와 같은 자세 추적 제어부(140)와 외란 관찰부(130)를 가지는 본 발명의 폐루프(closed loop) 자세 추적 제어 장치(100)의 미리 정해진 시간 안정성(predefined-time stability)은, 다음의 정리로 확립된다.
정리1) 불확실한 우주선 자세 시스템(6)을 고려한다. 외란 관찰부(130)와 자세 추적 제어부(140)가 각각 방정식(26) 및 (21)로 설계된 경우, 자세 추적 기동은 사전 정의된 시간 Tc 내에 수행되며, 자세 추적 오류 x1에 대한 제약 조건이 충족된다.
(증명)
최종 랴푸노프 함수를 다음과 같이 정의한다.
을 시간에 대하여 미분하고 식 (22) 및 식 (28)을 이용하면,
를 얻는다.
식 (34)로부터, 임을 알 수 있다. 따라서, Young 부등식에 따르면,
를 얻는다. 따라서, 식 (36)은 다음과 같이 간소화된다.
로 선택하자. 외란 관찰부(130)의 안착 시간(settling time, 즉, Te)은 자세 추적 제어부(140)의 안착 시간(즉, Tc)보다 적절하게 작아야 하므로, 이 된다. 따라서, 식 (38)은, 다음과 같이 다시 쓸 수 있다.
명제 2에 근거하여, 본 발명의 폐루프 형태의 자세 추적 제어 장치(100)는, 사전 정의된 시간 안정적(predefined-time stable)이며, 오차 궤적에 대한 제약이 충족되는 중에,고정된 시간 동안 기준 궤적을 추적한다고 결론 내릴 수 있다. 이것으로 증명이 완료되었다.
이하에서는 도 5 내지 도 18을 참조하여, 본 발명에 대한 시뮬레이션 결과를 설명하기로 한다.
명시되어 있는 사전 정의 시간 제약 제어 프레임워크 (Predefined-Time Constrained Contorl Framework)를 검증하기 위하여 강체 우주선에서 수치 시물레이션이 수행되었다.
공칭관성행렬(Nominal Inertia Matrix)은 이다. 관성 행렬의 불확실성은 일 때 로 정의된다.
상수항 는 연료 누출로 간주할 수 있다. 관성 불확실성과 더불어 이 단위 스텝일 때 다음과 같은 환경적 외란이 우주선에 작용할 수 있다:
두 번째로 고려된 외란은 우주먼지와의 접촉으로 인해 발생 할 수 있다.
관심 영역의 고해상도 이미지를 촬영하는 것은 우주선의 중요한 임무중 하나이다. 이를 수행하기 위하여 우주선은 각과 같은 기준 자세 궤적(reference attitude trajectory)으로 빠르고 정확한 자세 추적기동을 수행해야 한다. 명시한 제어기와 observer의 매개변수(parameter)는,
그리고 로 선택되었다.
먼저 도 5 내지 도 8에서는, 본 발명에 의해 수행된 결과로서, 다양한 초기조건을 설정한 경우의 시뮬레이션 결과를 보여준다. 특히 도 5 내지 도 8에서는, 4가지의 초기조건을 고려한다.
주요 목표는 폐쇄 루프 시스템(Closed- Loop System)의 수렴시간이 초기조건에 의존하지 않고 시스템 궤적이 내에서 안정화 된다는 것을 보여주는 것이다. 네 가지 초기조건은 다음과 같다.
도 5는 다양한 초기조건 설정시 자세 추적 오차를 나타내는 도면이고, 도 6은 다양한 초기조건 설정시 각속도 추적 오차를 나타내는 도면이며, 도 7은 다양한 초기조건 설정시 시변 이득을 나타내는 도면이고, 도 8은 다양한 초기조건 설정시 제어 입력을 나타내는 도면이다.
도 5 및 도 6에서 자세 및 각도 추적 오차 궤적은 수렴시간의 한계, 즉 보다 작은 8초 전에 안정화 된다. 이는 모든 초기조건에 해당되며 주어진 초기조건에 대해 이전에 궤적이 안정화 된다는 주장을 뒷받침한다. 또한 초기조건에 상관 없이 자세 추적 오차에 대한 목표 성능을 달성하였다.
도 7에 따르면, 자세 추적 오차가 성능함수(즉, )에 접근 하면 함수 가 1로 향하는 경향을 보이며, 결과적으로 tan 함수가 증가함을 알 수 있다. 도 8에서 자세추적오차궤적(attitude tracking error trajectory)이 허용된 영역의 경계와 접촉할 경우, 제어입력이 증가하여 제약조건을 넘지 않도록 한다.
다음으로, 도 9 내지 도 12에서는, 본 발명에 의해 수행된 결과로서, 다양한 수렴시간을 설정한 경우의 시뮬레이션 결과를 보여준다. 특히 도 9 내지 도 12에서는, 4가지의 수렴시간을 고려한다.
주요 제어 목표는 수렴 시간이 다를 경우 컨트롤러가 시스템의 궤적을 0으로 수렴하도록 강제하는 방법을 평가하는 것이었다. 이를 위해 초기 조건은 deg 및 deg/s로 취하였다. 관찰자와 제어기의 파라미터는 Tc와 Tf를 제외하고는 이전 부분과 동일하였으며, 4가지 경우를 고려하였다.
(1),
(2) ,
(3) ,
(4) .
도 9는 본 발명에 의해 수행된 결과로서, 다양한 수렴시간 설정시 자세 추적 오차를 나타내는 도면이고, 도 10은 본 발명에 의해 수행된 결과로서, 다양한 수렴시간 설정시 각속도 추적 오차를 나타내는 도면이며, 도 11은 본 발명에 의해 수행된 결과로서, 다양한 수렴시간 설정시 시변 이득 λ를 나타내는 도면이고, 도 12는 본 발명에 의해 수행된 결과로서, 다양한 수렴시간 설정시 제어 입력 τ를 나타내는 도면이다.
도 9와 도 10을 기준으로 Tc와 Tf의 값이 작으면 수렴이 빨라진다. 분명히, 수렴 시간이 더 작아지는 경우에도 자세 추적 오차에 대한 규정된 성능은 여전히 충족된다. 이는 제안된 제어 체계가 추적 오차에 대한 제약 조건을 위반하지 않고 사전 정의된 시간 전에 시스템을 안정화할 수 있음을 확인한다. 그러나, 도 11에 나타낸 것과 같이, 수렴 시간이 작아지면, 시변 이득 가 상승하여, 1에 가까워질 수 있다. 결과적으로, 도 12에 나타낸 것과 같이, 필요한 제어 노력이 증가한다. 실제 적용에서 제어 에너지의 중요성 때문에 수렴 속도와 필요한 제어 노력 사이에 절충이 이루어져야 한다.
다음으로, 도 14 내지 도 16에서는, 본 발명의 제어와 기존 기술을 이용한 제어 간의 성능 비교 결과가 도시되어 있다.
본 발명의 제어 체계는 이하 논문 1에서 교란 관찰자 기반 사전 정의 시간 자세 추적 제어(DOBPTATC, disturbance-observer-based predefined-time attitude tracking control) 및 강력한 고정 시간 자세 제어(RFTAC, robust fixed-time attitude control)와 비교된다. DOBPTATC와 RFTAC는 견고한 우주선을 위한 고정된 시간 수렴을 가진 교란 관찰자 기반 자세 추적 제어이다. 따라서 비교에 적합할 수 있다. 초기 조건, 관성 불확실성, 환경 교란 및 기준 궤적은 이하 논문 2에 제시된 것으로 간주된다.
논문 1. Cao, L.; Xiao, B.; Golestani, M. Robust fixed-time attitude stabilization control of flexible spacecraft with actuator uncertainty. Nonlinear Dyn. 2020, 100, 2505-2519. [CrossRef]
논문 2. Xiao, B.; Wu, X.; Cao, L.; Hu, X. Prescribed Time Attitude Tracking Control of Spacecraft with Arbitrary Disturbance. IEEE Trans. Aerosp. Electron. Syst. 2022, 58, 2531-2540. [CrossRef]
RFTAC 및 DOBPTATC와 비교하여 본 발명의 제어 방법의 성능을 추가로 분석하기 위해 3가지 제어 법칙에 따른 수렴 시간과 정상 상태 추적 오차를 도 13에 정리하였다. 도 13에서, 안착 시간은 자세 추적 오차가 영역 에 진입하고 그 이후에 머무르는 시간을 보여준다. 제안된 제어 방식은 정착 시간을 단축하여 수렴률을 향상시키고, 를 감소시켜 점 정확도(point accuracy)를 향상시키며, 를 감소시켜 점 안정성(point stability)을 동시에 향상시키는 것을 관찰할 수 있다. 따라서 시뮬레이션 결과는 우주선의 자세 추적 기동을 달성하는 데 있어 본 발명의 제어 프레임워크의 우수한 성능을 확인한다.
도 14는 본 발명과 기존 기술의 간의 자세 추적 오차를 비교하여 나타내는 도면이며, 도 15는 본 발명과 기존 기술의 간의 각속도 추적 오차를 비교하여 나타내는 도면이고, 도 16은 본 발명과 기존 기술의 간의 제어 입력 τ를 비교하여 나타내는 도면이다.
본 발명의 제어 방법은 정착 시간을 줄이면서 훨씬 더 빠른 수렴을 초래하는 것으로 나타났다. 더욱 정확하게는 본 발명의 제어 장치, RFTAC 및 DOBPTATC에서의 오차 궤적은 각각 대략 8초, 14초, 19초 후에 안정화되었다. 이는 본 발명의 제어 방법이 RFTAC 및 DOBPTATC에 대해 수렴 시간을 42.8%, 57.9% 단축한다는 것을 의미한다. 또한 RFTAC 및 DOBPTATC와 달리 제안된 제어 방법은 자세 추적 오차에 대해 과도 상태 및 정상 상태에서 바람직한 성능을 제공할 수 있다. 본 발명의 제어 방법은 RFTAC 및 DOBPTATC보다 우수한 성능을 달성했지만, 세 컨트롤러에 필요한 초기 제어 노력은 거의 동일했다. 이 사실은 도 16에서 명확하게 관찰할 수 있다. 따라서, 새로운 제어 프레임워크는 추가적인 제어 에너지를 소비하지 않고도 훨씬 더 나은 성능을 보장한다.
도 17은 추정 오차를 나타내는 도면이고, 도 18은 총 교란과 그 추정값을 나타내는 도면이다.
도 17을 참조하면, 예상했던 대로, 추정 오차는 사전 정의된 시간 내에 0으로 빠르게 수렴되었다. t = 2와 t = 6에서 총 교란에 두 가지 불연속성이 있지만, 제안된 관찰자는 여전히 추정 오차를 안정화할 수 있었고 0으로 빠르게 수렴하도록 만들었다. 또한 전체 교란은 추정치와 함께 도 18에 설명되어 있다. 도 18은 제안된 관찰자에 의해 전체 교란이 정확하게 재구성되었음을 분명히 보여준다.
100: 임의 교란을 가진 강성 우주선을 위한 자세 추적 제어 장치
Claims (13)
- 임의 교란을 가진 강성 우주선을 위한 자세 추적 제어 방법으로서,
(a) 우주선에 대하여 측정된 자세값(이하 '측정 자세값'이라 한다)과, 목표 자세값과 차이값(이하 '자세 오차'라 한다)를 산출하는 단계;
(b) 상기 자세 오차 및, 성능 함수로부터 가상 자세 제어값을 산출하는 단계;
(c) 상기 측정 자세값 및, 우주선의 자세를 제어하기 위한 제어 토크의 피드백 값으로부터 우주선에 대하여 작용하는 외란(disturbance)의 추정값을 산출하는 단계; 및,
(d) 상기 측정 자세값, 상기 가상 자세 제어값 및 상기 외란 추정값으로부터 새로운 제어 토크 값을 산출하는 단계
를 포함하고,
상기 성능함수는,
시간 t=0으로부터 감소하여, 유한 시간 내에 특정 값으로 수렴하고, 이후 일정한 값으로 유지되는,
임의 교란을 가진 강성 우주선을 위한 자세 추적 제어 방법. - 청구항 1에 있어서,
상기 단계(d) 이후,
(e) 상기 단계(d)에서 산출된 제어 토크 값을, 우주선의 자세를 제어하는 제어부로 전달하는 단계
를 더 포함하는 것을 특징으로 하는, 임의 교란을 가진 강성 우주선을 위한 자세 추적 제어 방법. - 삭제
- 청구항 1에 있어서,
상기 가상 자세 제어값을 구하는 파라미터에는,
상기 제어 토크값 산출시, 해당 제어 토크값이 특정 수렴시간(이하, '제어 수렴시간'이라 한다) 내에 특정값으로 수렴하여 이후 일정한 값으로 유지되는 경우, 그 제어 수렴시간 값을 더 포함하는 것
을 특징으로 하는, 임의 교란을 가진 강성 우주선을 위한 자세 추적 제어 방법. - 청구항 4에 있어서,
상기 가상 자세 제어값은,
로 산출되고,
상기 ρi는 성능함수 값, , , β는 제어 파라미터이고,
이며,
는 제어 파라미터, Tc는 상기 제어 수렴시간인 것
을 특징으로 하는, 임의 교란을 가진 강성 우주선을 위한 자세 추적 제어 방법. - 청구항 1에 있어서,
상기 외란 추정값의 산출은,
(c1) 미지 입력 선형 시스템(unknown input linear system, UILS)을 정의하는 단계;
(c2) 상기 UILS의 상태변수의 추정값의 미분 식과, 상기 UILS 식으로부터 상기 UILS에 포함된 미지 입력값(unknown input)을 구하는 단계; 및,
(c3) 상기 미지 입력값으로부터 외란 추정값을 산출하는 단계
를 포함하는 것을 특징으로 하는, 임의 교란을 가진 강성 우주선을 위한 자세 추적 제어 방법. - 청구항 6에 있어서,
상기 UILS 식은,
이고, 여기서 는 상기 UILS의 상태변수, k1, k2는 상수 파라미터, 은 미지 입력값이고,
상기 상기 UILS의 상태변수의 추정값의 미분 식은,
이며, 여기서 로서 추정 오차이며,
상기 외란 추정값 은,
로 산출되며, 여기서 이고,
이며, 이때 로서, 목표 자세값()과 측정 자세값()의 오차를 나타내며,
이고,
상기 추정 오차는 수렴시간(이하 '외란 관찰 수렴시간'이라 한다) 내에 특정값으로 수렴함으로써 상기 외란 관찰 수렴시간 내에 상기 외란 추정값이 산출되는 것
을 특징으로 하는, 임의 교란을 가진 강성 우주선을 위한 자세 추적 제어 방법. - 청구항 5에 있어서,
상기 제어 토크값은,
로 산출되고, 여기서 는 상기 단계(c)에서 산출된 외란 추정값인 것
을 특징으로 하는, 임의 교란을 가진 강성 우주선을 위한 자세 추적 제어 방법. - 청구항 8에 있어서,
상기 제어 토크는 상기 제어 수렴시간 내에 특정값으로 수렴하는 것
을 특징으로 하는, 임의 교란을 가진 강성 우주선을 위한 자세 추적 제어 방법. - 청구항 9에 있어서,
상기 제어 수렴시간 는,
상기 외란 추정값 산출시의 추정 오차의 수렴시간(이하 '외란 관찰 수렴시간'이라 한다) 보다 큰 값인 것
을 특징으로 하는, 임의 교란을 가진 강성 우주선을 위한 자세 추적 제어 방법. - 임의 교란을 가진 강성 우주선의 자세를 추적 제어하기 위한 장치(이하 '자세 추적 제어 장치'라 한다)로서,
적어도 하나의 프로세서; 및
컴퓨터로 실행가능한 명령을 저장하는 적어도 하나의 메모리를 포함하되,
상기 적어도 하나의 메모리에 저장된 상기 컴퓨터로 실행가능한 명령은, 상기 적어도 하나의 프로세서에 의하여,
(a) 우주선에 대하여 측정된 자세값(이하 '측정 자세값'이라 한다)과, 목표 자세값과 차이값(이하 '자세 오차'라 한다)를 산출하는 단계;
(b) 상기 자세 오차 및, 성능 함수로부터 가상 자세 제어값을 산출하는 단계;
(c) 상기 측정 자세값 및, 우주선의 자세를 제어하기 위한 제어 토크의 피드백 값으로부터 우주선에 대하여 작용하는 외란(disturbance)의 추정값을 산출하는 단계; 및,
(d) 상기 측정 자세값, 상기 가상 자세 제어값 및 상기 외란 추정값으로부터 새로운 제어 토크 값을 산출하는 단계
가 실행되도록 하고,
상기 성능함수는,
시간 t=0으로부터 감소하여, 유한 시간 내에 특정 값으로 수렴하고, 이후 일정한 값으로 유지되는,
자세 추적 제어 장치. - 임의 교란을 가진 강성 우주선의 자세를 추적 제어하기 위한, 컴퓨터로 판독 가능한 비일시적 저장 매체에 저장된 컴퓨터 프로그램으로서,
비일시적 저장 매체에 저장되며, 프로세서에 의하여,
(a) 우주선에 대하여 측정된 자세값(이하 '측정 자세값'이라 한다)과, 목표 자세값과 차이값(이하 '자세 오차'라 한다)를 산출하는 단계;
(b) 상기 자세 오차 및, 성능 함수로부터 가상 자세 제어값을 산출하는 단계;
(c) 상기 측정 자세값 및, 우주선의 자세를 제어하기 위한 제어 토크의 피드백 값으로부터 우주선에 대하여 작용하는 외란(disturbance)의 추정값을 산출하는 단계; 및,
(d) 상기 측정 자세값, 상기 가상 자세 제어값 및 상기 외란 추정값으로부터 새로운 제어 토크 값을 산출하는 단계
가 실행되도록 하는 명령을 포함하고,
상기 성능함수는,
시간 t=0으로부터 감소하여, 유한 시간 내에 특정 값으로 수렴하고, 이후 일정한 값으로 유지되는,
임의 교란을 가진 강성 우주선의 자세를 추적 제어하기 위한, 컴퓨터로 판독 가능한 비일시적 저장 매체에 저장된 컴퓨터 프로그램. - 임의 교란을 가진 강성 우주선의 자세를 추적 제어하기 위한 장치(이하 '자세 추적 제어 장치'라 한다)로서,
우주선에 대하여 측정된 자세값(이하 '측정 자세값'이라 한다)과, 목표 자세값과 차이값(이하 '자세 오차'라 한다)를 산출하는 자세 오차 산출부;
상기 자세 오차 및, 성능 함수로부터 가상 자세 제어값을 산출하는 가상 제어부;
상기 측정 자세값 및, 우주선의 자세를 제어하기 위한 제어 토크의 피드백 값으로부터 우주선에 대하여 작용하는 외란(disturbance)의 추정값을 산출하는 외란 관찰부; 및,
상기 측정 자세값, 상기 가상 자세 제어값 및 상기 외란 추정값으로부터 새로운 제어 토크 값을 산출하는 자세 추적 제어부
를 포함하고,
상기 성능함수는,
시간 t=0으로부터 감소하여, 유한 시간 내에 특정 값으로 수렴하고, 이후 일정한 값으로 유지되는,
자세 추적 제어 장치.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
KR1020230045514A KR102605904B1 (ko) | 2023-04-06 | 2023-04-06 | 임의 교란을 가진 강성 우주선을 위한 자세 추적 제어 방법 및 장치 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
KR1020230045514A KR102605904B1 (ko) | 2023-04-06 | 2023-04-06 | 임의 교란을 가진 강성 우주선을 위한 자세 추적 제어 방법 및 장치 |
Publications (1)
Publication Number | Publication Date |
---|---|
KR102605904B1 true KR102605904B1 (ko) | 2023-11-27 |
Family
ID=88967908
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
KR1020230045514A KR102605904B1 (ko) | 2023-04-06 | 2023-04-06 | 임의 교란을 가진 강성 우주선을 위한 자세 추적 제어 방법 및 장치 |
Country Status (1)
Country | Link |
---|---|
KR (1) | KR102605904B1 (ko) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN118244790A (zh) * | 2024-05-28 | 2024-06-25 | 山东科技大学 | 一种具有时变误差约束的飞艇固定时间轨迹跟踪控制方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR19990063535A (ko) * | 1997-12-25 | 1999-07-26 | 가네꼬 히사시 | 이동 물체의 상태 제어 장치 및 방법 |
JP3859454B2 (ja) * | 2001-03-01 | 2006-12-20 | Nec東芝スペースシステム株式会社 | 人工衛星のマヌーバ制御装置 |
US9776741B1 (en) | 2016-06-02 | 2017-10-03 | Beihang University | Method for refined attitude control based on output feedback for flexible spacecraft |
CN112214890A (zh) * | 2020-10-10 | 2021-01-12 | 北京理工大学 | 一种小天体探测器绕飞轨道跟踪控制方法 |
CN114987800A (zh) * | 2022-05-19 | 2022-09-02 | 上海航天控制技术研究所 | 一种基于滑模干扰观测的航天器高精度姿态控制方法 |
-
2023
- 2023-04-06 KR KR1020230045514A patent/KR102605904B1/ko active IP Right Grant
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR19990063535A (ko) * | 1997-12-25 | 1999-07-26 | 가네꼬 히사시 | 이동 물체의 상태 제어 장치 및 방법 |
JP3859454B2 (ja) * | 2001-03-01 | 2006-12-20 | Nec東芝スペースシステム株式会社 | 人工衛星のマヌーバ制御装置 |
US9776741B1 (en) | 2016-06-02 | 2017-10-03 | Beihang University | Method for refined attitude control based on output feedback for flexible spacecraft |
CN112214890A (zh) * | 2020-10-10 | 2021-01-12 | 北京理工大学 | 一种小天体探测器绕飞轨道跟踪控制方法 |
CN114987800A (zh) * | 2022-05-19 | 2022-09-02 | 上海航天控制技术研究所 | 一种基于滑模干扰观测的航天器高精度姿态控制方法 |
Non-Patent Citations (1)
Title |
---|
"Funnel control for nonlinear systems with known strict relative degree (Automatica Volume 87, January 2018, Pages 345-357)" |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN118244790A (zh) * | 2024-05-28 | 2024-06-25 | 山东科技大学 | 一种具有时变误差约束的飞艇固定时间轨迹跟踪控制方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Schoellig et al. | Optimization-based iterative learning for precise quadrocopter trajectory tracking | |
Bu | Guaranteeing prescribed output tracking performance for air-breathing hypersonic vehicles via non-affine back-stepping control design | |
Dong et al. | Adaptive disturbance observer‐based finite‐time continuous fault‐tolerant control for reentry RLV | |
Zong et al. | Adaptive high‐order dynamic sliding mode control for a flexible air‐breathing hypersonic vehicle | |
Bu et al. | Minimal-learning-parameter based simplified adaptive neural back-stepping control of flexible air-breathing hypersonic vehicles without virtual controllers | |
Lu et al. | Finite‐time attitude control for rigid spacecraft‐based on adaptive super‐twisting algorithm | |
Tripathi et al. | Finite‐time super twisting sliding mode controller based on higher‐order sliding mode observer for real‐time trajectory tracking of a quadrotor | |
Wang et al. | Robust adaptive backstepping control for reentry reusable launch vehicles | |
KR102605904B1 (ko) | 임의 교란을 가진 강성 우주선을 위한 자세 추적 제어 방법 및 장치 | |
Chai et al. | Integrated missile guidance and control using optimization-based predictive control | |
CN103312256B (zh) | 基于干扰观测器的伺服电机网络化控制方法 | |
Chen et al. | Adaptive fixed-time backstepping control for three-dimensional trajectory tracking of underactuated autonomous underwater vehicles | |
Liu et al. | Adaptive saturated tracking control for spacecraft proximity operations via integral terminal sliding mode technique | |
Benosman et al. | Nonlinear control allocation for non-minimum phase systems | |
Bounemeur et al. | Active adaptive fuzzy fault-tolerant control for a class of nonlinear systems with actuator faults | |
Sun et al. | Augmented fixed‐time observer‐based continuous robust control for hypersonic vehicles with measurement noises | |
Yan et al. | Event-triggered output feedback sliding mode control of mechanical systems | |
Ma et al. | Observer-based event-triggered type-2 fuzzy control for uncertain steer-by-wire systems | |
Liu et al. | Application of the improved incremental nonlinear dynamic inversion in fixed-wing UAV flight tests | |
Xue et al. | Compound tracking control based on MPC for quadrotors with disturbances | |
Ataei-Esfahani et al. | Nonlinear control design of a hypersonic aircraft using sum-of-squares method | |
Wang et al. | Fixed-time event-triggered sliding mode cooperative path-following control with prescribed performance for USVs based on lumped disturbance observer | |
Khatri et al. | Accessible regions for controlled aircraft maneuvering | |
Ebirim et al. | Constrained Model Predictive Control With Integral Action for Twin Rotor MIMO Systems | |
Bi et al. | Attitude stabilization of rigid spacecraft implemented in backstepping control with input delay |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
E701 | Decision to grant or registration of patent right | ||
GRNT | Written decision to grant |