진동자 실험의 독립 자극과 계수 식별

두 번 측정해도 같은 정보일 수 있다

감쇠와 강성을 모르는 진동자의 가속도를 두 번 측정한다고 합시다. 한 번보다 정보가 많아 보이지만, 두 측정에서 변위와 속도의 비율이 같다면 같은 방정식을 반복했을 수 있습니다. 데이터 표의 행 개수와 독립적인 계수 제약의 개수는 다릅니다.

이 단원은 하나의 고정된 선형 진동자 모형에서 초기 상태를 제어한 실험을 사용합니다. 각 실험은 위치와 속도를 새로 설정한 뒤 순간 가속도를 측정합니다. 서로 다른 실험의 점들을 실제 시간 궤적으로 이어 붙이지 않습니다. 가능한 계수의 전체 해를 구하고 정확한 유일성, 물리적 허용과 잡음 민감도를 나눕니다.

계수와 단위를 먼저 정한다

알려진 스칼라 질량으로 정규화하고 외부 가속도 입력이 0인 모형을 사용합니다.

$$a=-\gamma v-\kappa x.$$

위치 x의 단위는 m, 속도 v는 m/s, 가속도 a는 m/s²입니다. 감쇠 계수 γ는 1/s, 단위 질량당 강성 κ는 1/s²입니다. 여기서 물리적으로 허용하는 영역은 γ≥0, κ>0입니다.

이들은 정규화된 계수입니다. 질량까지 보정되지 않은 미지수라면 γ와 κ를 알아도 원래의 질량·감쇠·강성을 각각 식별한 것은 아닐 수 있습니다. 입력의 물리적 크기가 왜 필요한지는 앞의 힘·질량 단원과 연결됩니다.

x와 v를 정한 실험에서 관측한 가속도는 γv+κx=−a라는 한 식을 줍니다. 두 실험은 같은 계수를 공유해야 합니다. 초기 상태를 다시 설정할 때 계수도 바꾸면 다른 역문제입니다.

두 실험을 행렬로 쓴다

두 초기 상태 실험의 관측을 다음처럼 모읍니다.

$$A\theta=y,\qquad A=\begin{pmatrix}v_1&x_1\\v_2&x_2\end{pmatrix},\quad \theta=\begin{pmatrix}\gamma\\\kappa\end{pmatrix},\quad y=\begin{pmatrix}-a_1\\-a_2\end{pmatrix}.$$

행렬식은 D=v1x2−v2x1입니다. 지원하는 유리수 입력에서는 0 여부를 부동소수점 허용 오차 없이 검사할 수 있습니다. D가 0이 아니면 다음처럼 대수적 해가 유일합니다.

$$\gamma=\frac{y_1x_2-y_2x_1}{D},\qquad \kappa=\frac{v_1y_2-v_2y_1}{D}.$$

이것은 지금 정한 데이터 행렬과 선형 모형에 대한 결론입니다. 실제 센서 오차, 모형에 없는 힘이나 모든 시간 궤적까지 설명했다는 뜻은 아닙니다.

서로 독립적인 자극을 손으로 계산한다

실제 계수가 γ=1/s, κ=2/s²라고 합시다. 첫 실험은 v1=1m/s, x1=0으로 시작하여 a1=−1m/s²를 줍니다. 둘째 실험은 v2=0, x2=1m로 시작하여 a2=−2m/s²를 줍니다.

선택한 SI 기준에서 행렬은 숫자로 단위행렬입니다. 첫 행은 감쇠를, 둘째 행은 강성을 분리하여 γ=1, κ=2를 얻습니다. 두 식을 모두 유지하면서 어느 계수도 바꿀 수 없습니다.

첫 실험을 두 번 복제하는 것과 다릅니다. 잡음 모형 아래에서 반복 측정으로 평균의 정확도를 높일 수는 있지만, 같은 행이 정확한 계수 정보의 새로운 방향을 만드는 것은 아닙니다.

Rank 1의 전체 계열과 물리 영역

이번에는 (v1,x1)=(1,1), (v2,x2)=(2,2)를 사용합니다. 같은 실제 계수에서 y=(3,6)이며 둘째 행은 첫째 행의 두 배입니다. 남는 식은 γ+κ=3 하나입니다.

전체 실수 해는 다음과 같습니다.

$$(\gamma,\kappa)=(3,0)+t(-1,1),\qquad t\in\mathbb R.$$

프로그램이 실제 계수 (1,2)만 반환하면 다른 가능한 해들을 숨긴 것입니다. 최소 노름이나 정규화 항이 한 해를 선택해도 관측 자체에는 직선 전체가 남습니다.

물리 영역에서는 γ=3−t≥0, κ=t>0이므로 0<t≤3입니다. 강성이 0인 아래 끝점은 제외하고 감쇠가 0인 위 끝점은 포함합니다. 몇 개 표본 점으로 대체하지 않고 열린 끝점과 닫힌 끝점을 보존해야 합니다.

Rank, 일관성과 물리적 허용은 따로 확인한다

서로 비례하는 행에 응답이 같은 비율로 비례하지 않으면 식은 모순입니다. 앞의 두 행에 y=(3,6.01)을 정확하게 맞추는 계수는 없습니다. Rank가 1이라는 정보만으로 해 직선이 존재한다고 결론낼 수 없으며 우변도 확인해야 합니다.

두 행이 모두 0이고 응답도 0이면 모든 실수 계수 쌍이 해입니다. 물리적 부분은 선언한 감쇠·강성 영역 전체입니다. 0인 행의 응답이 0이 아니면 그 식을 만족하는 계수는 없습니다.

Rank가 2이면 임의의 응답에 유일한 대수적 해가 있지만, 음의 감쇠나 0 이하 강성이 나올 수 있습니다. 실제 감쇠가 0인 예에서도 작은 가속도 오차가 음의 추정 감쇠를 만들 수 있습니다. 대수적 유일성과 물리적 허용은 다른 표시가 필요합니다.

물리 영역이 비었다는 결과는 장치의 실제 계수가 불가능하다는 뜻으로만 해석하지 않습니다. 잡음이나 모형 불일치를 정확한 등식으로 맞추는 것이 적절하지 않을 수 있습니다. 여기의 오차 조절은 결정론적 섭동이며 통계적인 신뢰 구간을 자동으로 주지 않습니다.

관측 기하를 무차원으로 만든다

A의 두 열은 서로 단위가 다릅니다. 원래 행렬식이나 열의 Euclidean 노름 숫자를 단위 선택과 무관한 실험 품질로 해석하면 안 됩니다.

양의 기준 길이 L과 기준 시간 T를 명시하여 다음을 정의합니다.

$$\bar A_i=\left(\frac{v_iT}{L},\frac{x_i}{L}\right),\qquad \bar\theta=(\gamma T,\kappa T^2),\qquad \bar y_i=\frac{y_iT^2}{L}.$$

이때 무차원 식은 Abar·thetabar=ybar입니다. 양의 기준 척도를 바꾸어도 정확한 rank와 복원한 물리 계수는 같습니다. 반면 특잇값과 조건수는 선택한 상대 계수 척도와 기하를 반영하므로 달라질 수 있습니다.

기본값 L=1m, T=1s는 익숙한 SI 숫자를 그대로 보여 주는 편리한 기준입니다. 그 편리함 때문에 정규화 가정을 생략하지 않아야 합니다.

정확한 유일성과 수치 조건은 다르다

선택한 무차원 행렬과 Euclidean 노름에서 2-노름 조건수는 rank가 2일 때 최대 특잇값/최소 특잇값입니다. 이 선형계의 특정 방향에서 오차가 얼마나 증폭될 수 있는지 보여 줍니다. Rank가 부족하면 유한한 역문제 조건수가 없으며 보통의 큰 유한 숫자로 숨기지 않습니다.

행이 (1,1), (1,1.1)인 설계를 봅시다. 행렬식은 0.1이라 해는 정확히 유일하지만 두 방향이 가까워 일부 계수 조합을 분리하기 어렵습니다.

실제 계수 (1,2)의 잡음 없는 응답은 (3,3.2)입니다. 첫 가속도가 0.01만큼 커져 응답이 2.99가 되고 둘째 응답이 3.2로 유지되면, 추정값은 γ=0.89, κ=2.1입니다. 독립적인 단위행렬 설계에서 같은 응답 섭동은 γ만 1에서 0.99로 바꿉니다. 두 경우 모두 해가 유일하지만 민감도는 다릅니다.

정확한 판정과 근사 그림을 구분한다

전체 아핀 해 계열과 물리 영역의 교집합은 유리수 소거와 정확한 부등식으로 구합니다. 특잇값과 조건수는 수치적인 기하 출력입니다. 그 작은 오차를 정확한 rank나 해 존재 증명으로 대신하지 않습니다.

계산 공식도 검산해야 합니다. 특잇값이 같을 때 거의 같은 숫자를 빼는 공식은 정확도를 잃을 수 있습니다. 독립적인 고정밀 SVD 기준값으로 구현을 확인하면 이런 오류를 찾을 수 있습니다. 코드에서도 수치 조건과 유리수 rank 판정을 나누어 유지합니다.

그림은 무한한 제약 직선의 보이는 일부를 표시합니다. 전체 해는 샘플 곡선이 교차해 보인다는 사실이 아니라 아핀 계열과 정의역 경계로 설명합니다. 행 벡터 그림은 명시한 정규화 아래의 무차원 좌표입니다.

물리 계수를 학습할 때의 의미

정보가 부족한 자극에서 방정식 오차가 작아도 계수가 크게 달라질 수 있습니다. 같은 관측을 더 복제해도 영공간 방향이 없어지지는 않습니다. 사전 가정이나 정규화 항은 한 답을 안정시킬 수 있지만 데이터의 유일한 식별 정보와 구분해야 합니다.

평가에서는 특성 설계, rank, 물리 계수 영역, 정규화와 잡음 가정을 기록합니다. 전체 운동 예측이 목표라면 순간 방정식의 적합과 장시간 시뮬레이션 정확성도 구분합니다. 두 순간 실험이 모든 시간의 동역학을 검증하는 것은 아닙니다.

연습문제와 해설

1. 문제: 행 (1,0), (0,1)과 응답 (1,2)의 계수는? 해설: γ=1, κ=2이며 각 행이 하나의 계수를 분리합니다.
2. 문제: 행 (1,1), (2,2)과 응답 (3,6)의 전체 해는? 해설: 모든 실수 t에 대해 (γ,κ)=(3,0)+t(−1,1)입니다.
3. 문제: 그중 γ≥0, κ>0을 만족하는 부분은? 해설: 0<t≤3이며 아래 끝점은 열리고 위 끝점은 닫힙니다.
4. 문제: Rank 1이면 응답 (3,6.01)에도 해가 있나요? 해설: 아닙니다. 비례하는 행에는 비례하는 응답이 필요합니다.
5. 문제: Rank 2의 유일한 추정값이 물리적으로 허용되지 않을 수 있나요? 해설: 가능합니다. 관측 오차로 음의 감쇠나 0 이하 강성이 나올 수 있습니다.
6. 문제: 기준 L이나 T를 바꾸면 새로운 식별 정보를 얻나요? 해설: 아닙니다. 수치 표현과 조건 기하는 달라져도 정확한 rank와 물리 계수는 같습니다.

참고 자료

MIT Underactuated Robotics의 system identification 장은 구조가 알려진 기계 모형의 계수 추정과 데이터 행렬의 배경을 제공합니다. 초기 상태 실험, 전체 계열, 물리 구간과 무차원 조건수 예제는 위 모형에서 직접 유도했습니다.