# 안정성 분석 방법

### 개요

비행체의 안정성 분석은 비행체의 거동을 이해하고, 비행체가 안정적인 비행을 유지할 수 있는지를 판단하는 데 중요한 역할을 한다. 안정성 분석 방법에는 여러 가지가 있으며, 주로 선형 시스템 이론을 바탕으로 한 수학적 방법이 많이 사용된다. 아래에서는 주요 안정성 분석 방법들을 설명한다.

### 선형화 및 작은 진동 해석

비행체의 운동 방정식은 보통 비선형이다. 그러나 안정성 분석을 위해 주로 선형화를 사용하여 작은 진동 해석을 수행한다.

#### 운동 방정식의 선형화

비선형 운동 방정식을 선형화하기 위해, 평형점 주변에서 테일러 전개를 이용하여 1차 항만을 고려한다. 예를 들어, 상태 벡터 $\mathbf{x}$와 제어 입력 $\mathbf{u}$에 대해 다음과 같은 비선형 시스템을 고려한다:

$$
\dot{\mathbf{x}} = \mathbf{f}(\mathbf{x}, \mathbf{u})
$$

이 시스템을 평형점 $(\mathbf{x}\_e, \mathbf{u}\_e)$ 주변에서 선형화하면 다음과 같다:

$$
\dot{\mathbf{x}} \approx \mathbf{A} \mathbf{x} + \mathbf{B} \mathbf{u}
$$

여기서 $\mathbf{A}$와 $\mathbf{B}$는 다음과 같이 정의된다:

$$
\mathbf{A} = \left. \frac{\partial \mathbf{f}}{\partial \mathbf{x}} \right|\_{(\mathbf{x}\_e, \mathbf{u}\_e)}
$$

$$
\mathbf{B} = \left. \frac{\partial \mathbf{f}}{\partial \mathbf{u}} \right|\_{(\mathbf{x}\_e, \mathbf{u}\_e)}
$$

#### 고유치 해석

선형화된 시스템의 안정성을 분석하기 위해 고유치 해석을 수행한다. 시스템 행렬 $\mathbf{A}$의 고유치 $\lambda\_i$를 계산하고, 이 고유치의 실수 부분이 음수인지 확인한다. 만약 모든 고유치의 실수 부분이 음수이면, 시스템은 안정적이다.

$$
\text{det}(\mathbf{A} - \lambda \mathbf{I}) = 0
$$

여기서 $\mathbf{I}$는 항등 행렬이다.

### 루스-허위츠 기준

루스-허위츠 기준은 시스템의 특성 방정식의 계수를 이용하여 시스템의 안정성을 판단하는 방법이다. 특성 방정식이 다음과 같을 때:

$$
a\_n \lambda^n + a\_{n-1} \lambda^{n-1} + \cdots + a\_1 \lambda + a\_0 = 0
$$

루스 행렬을 구성하고, 이 행렬의 첫 번째 열의 모든 요소가 양수인지 확인한다. 만약 첫 번째 열의 모든 요소가 양수이면, 시스템은 안정적이다.

### 루트 궤적법

루트 궤적법은 제어 시스템의 폐루프 특성 방정식의 근이 제어 이득 $K$에 따라 어떻게 변하는지를 시각적으로 보여주는 방법이다. 시스템의 개루프 전달 함수가 다음과 같을 때:

$$
G(s)H(s) = \frac{N(s)}{D(s)}
$$

폐루프 전달 함수의 특성 방정식은 다음과 같다:

$$
D(s) + K N(s) = 0
$$

이 방정식의 근을 $K$에 따라 플롯하여 시스템의 안정성을 분석한다.

### 나이퀴스트 및 보드 다이어그램

#### 나이퀴스트 다이어그램

나이퀴스트 다이어그램은 주파수 응답을 이용하여 시스템의 안정성을 분석하는 방법이다. 시스템의 개루프 전달 함수 $L(s) = G(s)H(s)$의 주파수 응답을 복소 평면에 플롯한다. 나이퀴스트 기준에 따라, $-1$점을 포함하지 않으면 시스템은 안정적이다.

#### 보드 다이어그램

보드 다이어그램은 주파수 응답을 크기와 위상으로 나누어 로그 스케일로 플롯한 것이다. 주파수 응답을 이용하여 시스템의 이득 여유(gain margin)와 위상 여유(phase margin)를 분석한다. 충분한 이득 여유와 위상 여유가 있으면 시스템은 안정적이다.

### 리아푸노프 안정성

리아푸노프 안정성 이론은 비선형 시스템의 안정성을 분석하는 데 매우 유용한 도구이다. 리아푸노프 함수는 에너지 함수와 유사한 개념으로, 시스템의 상태가 평형점에서 멀어질수록 값이 증가하고, 평형점에 가까워질수록 값이 감소하는 함수이다.

#### 리아푸노프 함수

리아푸노프 함수 $V(\mathbf{x})$는 다음 조건을 만족해야 한다:

1. $V(\mathbf{x})$는 평형점에서 0이 된다 ($V(\mathbf{x}\_e) = 0$).
2. $V(\mathbf{x})$는 평형점 주변에서 양의 정의이다 ($V(\mathbf{x}) > 0$ for $\mathbf{x} \neq \mathbf{x}\_e$).
3. $V(\mathbf{x})$의 시간에 대한 도함수 $\dot{V}(\mathbf{x})$는 평형점 주변에서 음의 정의이다 ($\dot{V}(\mathbf{x}) < 0$ for $\mathbf{x} \neq \mathbf{x}\_e$).

#### 리아푸노프의 직접 방법

리아푸노프의 직접 방법을 사용하면 시스템의 운동 방정식 $\dot{\mathbf{x}} = \mathbf{f}(\mathbf{x})$에 대해 다음과 같은 리아푸노프 함수를 선택하고, 시간에 따른 변화율 $\dot{V}(\mathbf{x})$를 계산하여 시스템의 안정성을 분석할 수 있다:

$$
\dot{V}(\mathbf{x}) = \frac{\partial V}{\partial \mathbf{x}} \mathbf{f}(\mathbf{x})
$$

만약 $\dot{V}(\mathbf{x})$가 평형점 주변에서 음의 정의이면, 시스템은 안정적이다.

### 상태 공간 분석

상태 공간 분석 방법은 비행체의 운동 방정식을 상태 변수 형태로 표현하고, 이 상태 변수의 시간에 따른 변화를 분석하여 안정성을 평가하는 방법이다.

#### 상태 공간 표현

비행체의 운동 방정식을 상태 변수 $\mathbf{x}$와 제어 입력 $\mathbf{u}$로 나타내면 다음과 같다:

$$
\dot{\mathbf{x}} = \mathbf{A}\mathbf{x} + \mathbf{B}\mathbf{u}
$$

$$
\mathbf{y} = \mathbf{C}\mathbf{x} + \mathbf{D}\mathbf{u}
$$

여기서 $\mathbf{A}$, $\mathbf{B}$, $\mathbf{C}$, $\mathbf{D}$는 시스템 매트릭스이다.

#### 안정성 평가

상태 공간 표현을 이용하여 시스템의 고유치, 제어 가능성, 관측 가능성을 분석함으로써 시스템의 안정성을 평가할 수 있다. 예를 들어, 시스템 매트릭스 $\mathbf{A}$의 고유치가 모두 실수 부분이 음수이면, 시스템은 안정적이다.

***

비행체의 안정성 분석은 여러 가지 방법을 통해 수행될 수 있으며, 각 방법은 특정 상황과 조건에서 유용하다. 선형화 및 작은 진동 해석, 루스-허위츠 기준, 루트 궤적법, 나이퀴스트 및 보드 다이어그램, 리아푸노프 안정성 이론, 상태 공간 분석 등 다양한 방법을 종합적으로 사용하여 비행체의 안정성을 평가하는 것이 중요하다. 각 방법의 장단점을 이해하고, 특정 비행체와 조건에 적합한 방법을 선택하여 정확한 안정성 분석을 수행하는 것이 필요하다.