프로니 방법

1. 개요

프로니 방법 Prony's Method
모형	$x_k = \sum_{i=1}^{p} c_i z_i^{\,k}$ — 감쇠 진동 $p$ 개의 합
핵심 발상	비선형 맞춤을 선형 예측 + 다항식 근 + 선형 최소자승으로 쪼갠다
최소 데이터	잡음이 없으면 표본 $2p$ 개로 정확 복원
약점	잡음에 매우 민감, 최소자승판도 감쇠를 편향 추정
후계자	쿠마레산-터프츠(과차수+SVD), 행렬 펜슬, ESPRIT
발표	가스파르 드 프로니, 1795

지수함수 맞춤은 비선형 최적화 문제다. 프로니는 이걸 연립일차방정식 두 번과 다항식 근 찾기 한 번으로 바꿔 버렸다. 1795년에.

프로니 방법(Prony’s method)은 등간격 표본 $x_0, x_1, \dots, x_{N-1}$ 을 복소 지수함수 $p$ 개의 합

$$x_k = \sum_{i=1}^{p} c_i\, z_i^{\,k}, \qquad z_i = e^{(\sigma_i + j\omega_i)\Delta t}$$

로 맞추는 방법이다. 극 $z_i$ 의 크기가 감쇠율, 각도가 주파수, 진폭 $c_i$ 의 크기·위상이 모드의 세기·위상이다. 프랑스 공학자 가스파르 드 프로니(1755–1839)가 1795년 에콜 폴리테크니크 저널에 발표했는데, 원래 목적은 기체·증기의 압력-온도 실험 데이터를 지수함수로 보간하는 것이었다¹.

모형은 $z_i$ 에 대해 비선형이라 그대로 맞추려면 레벤버그-마쿼트 방법 같은 반복 최적화와 좋은 초기값이 필요하다. 프로니의 통찰은 지수함수의 합은 선형 점화식을 만족한다는 것이다. 이 점화식 계수를 먼저 선형으로 풀고, 그 특성다항식의 근으로 $z_i$ 를 얻고, 마지막으로 진폭을 선형으로 풀면 비선형 문제가 사라진다. 이 세 단계 구조는 자기회귀 모형의 선형 예측, 행켈 행렬의 크로네커 랭크 정리, 반데르몽드 행렬 최소자승을 한 줄에 꿰는 것이라, 신호처리·구조동역학·전력계통에서 「모드 추출」이라 불리는 거의 모든 기법의 조상이다.

2. 세 단계

1단계 — 선형 예측. 극 $z_1, \dots, z_p$ 를 근으로 갖는 다항식

$$P(z) = \prod_{i=1}^{p}(z - z_i) = z^p + a_1 z^{p-1} + \cdots + a_p$$

를 생각하자. 각 $z_i^{\,k-p}P(z_i) = 0$ 에 $c_i$ 를 곱해 더하면

$$x_k + a_1 x_{k-1} + \cdots + a_p x_{k-p} = 0, \qquad k = p, p+1, \dots$$

라는 선형 점화식이 나온다. 미지수는 $a_1, \dots, a_p$ 뿐이고 데이터에 대해 선형이다. $k = p, \dots, 2p-1$ 의 $p$ 개 식을 쓰면 정사각 토플리츠 행렬(열 순서를 뒤집으면 행켈) 계가 되고, 이것이 원조 프로니다. 표본이 더 있으면 $k = p, \dots, N-1$ 을 모두 써서 최소자승으로 푼다(최소자승 프로니).

2단계 — 근 찾기. $P(z) = 0$ 의 근이 극 $z_i$ 다. 실무에서는 동반 행렬(companion matrix)의 고유값으로 구한다. numpy.roots 가 정확히 그렇게 한다.

3단계 — 진폭. $z_i$ 를 알면 $x_k = \sum c_i z_i^{\,k}$ 는 $c_i$ 에 대해 선형이고, 계수 행렬이 $z_i$ 를 절점으로 하는 반데르몽드 행렬이다. 최소자승으로 $c_i$ 를 푼다.

잡음이 없으면 이 절차는 정확하다. $x_k = 2(0.98)^k\cos(0.6k) + 2(0.95)^k\cos(1.1k)$ (켤레쌍 2개, $p = 4$) 의 표본 8개만으로 극 크기 $0.98$, $0.95$, 각도 $0.6$, $1.1$ 과 진폭 크기 $1, 1, 1, 1$ 을 반올림 수준까지 복원한다. 같은 길이의 FFT 로는 주파수 해상도가 $2\pi/8 \approx 0.79$ rad 라 두 성분을 가르지도 못한다. 이 「초해상도」가 프로니 계열의 매력이다².

3. 잡음 — 원조 프로니가 무너지는 곳

같은 신호에 표준편차 0.01(진폭의 0.5%)의 백색 잡음을 섞고 $N = 64$ 표본을 준 뒤, 잡음 실현 300개로 통계를 냈다.

방법	쓴 표본	주파수 오차 중앙값	극 크기 평균 (참값 0.98 / 0.95)	극 크기 표준편차
원조 프로니 ($p = 4$, 정사각)	8	0.196 rad	— (88% 시행에서 주파수 오차 > 0.05)	—
최소자승 프로니 ($p = 4$)	64	$4.0\times10^{-3}$ rad	0.9746 / 0.9456	0.0021 / 0.0031
과차수 ($p = 20$) + 진폭 상위 4개	64	—	0.97999 / 0.94998	0.00014 / 0.00041

세 줄이 프로니 방법 200년사의 요약이다.

1. 원조 프로니는 잡음에서 쓸 수 없다. 진폭의 0.5% 잡음에 주파수가 0.2 rad 씩 틀린다. 점화식 계수 $a_i$ 를 정확히 $p$ 개 식으로 맞추면 잡음까지 그대로 「지수함수」로 해석해 버린다.
2. 최소자승으로 바꾸면 분산은 줄지만 편향이 남는다. 극 크기 평균이 0.9746 으로 참값 0.98 보다 표준편차의 2.5배만큼 작다 — 감쇠를 체계적으로 과대 추정한다. 원인은 1단계 회귀 행렬 자체에 잡음이 들어 있다는 것(설명변수 오차, errors-in-variables)이다. 보통 최소자승은 우변의 오차만 가정하므로 일관 추정량이 아니다. 총최소제곱이 이 상황을 위한 도구이고, 오스본·스미스(1995)의 수정 프로니는 반복으로 이 편향을 없앤 최대우도 추정이다.
3. 차수를 일부러 크게 잡으면 오히려 정확해진다. $p = 20$ 으로 점화식을 맞추면 잡음을 흡수할 여분의 극 16개가 생기고, 진짜 극 4개는 진폭이 커서 골라낼 수 있다. 결과는 편향이 사라지고 표준편차가 7~15배 줄었다. 이것을 체계화한 것이 쿠마레산·터프츠(1982)다 — 과차수 선형 예측 행렬을 SVD 로 랭크 $p$ 로 잘라 잡음 부분공간을 버리고, 역방향 선형 예측을 쓰면 여분의 극이 단위원 안쪽에, 신호 극(의 역수)은 바깥쪽에 떨어져 둘을 가를 수 있다는 것을 보였다.

이 마지막 아이디어 — 「행켈/토플리츠 구조의 데이터 행렬을 SVD 로 신호 부분공간과 잡음 부분공간으로 가른다」 — 를 끝까지 밀면 다항식 근 찾기 단계 자체가 사라지고 이동 불변성으로 극을 바로 고유값으로 읽는 행렬 펜슬 방법(후아·사르카르 1990)과 ESPRIT 이 된다. 펜슬의 구체적 절차와 같은 신호에서의 실측은 행켈 행렬 문서에 있다.

4. 어디에 쓰이는가

• 전력계통 저주파 진동 해석. 송전망의 발전기 군 사이 진동(0.1–2 Hz 대역, 감쇠비 몇 %)을 PMU 측정이나 시뮬레이션 과도 응답에서 뽑아 감쇠가 충분한지 판정한다. 하우어 등(1990)이 이 분야에 프로니 분석을 정착시켰고, 지금도 계통 운영 도구의 표준 기능이다.
• 구조물 모드 식별. 충격 응답이나 자유 진동 기록에서 고유진동수·감쇠비를 뽑는다. 구조동역학에서는 다채널로 확장한 복소 지수법(CE)·최소자승 복소 지수법(LSCE)이라는 이름으로 쓰이고, 상태공간판이 ERA 다(시스템 식별 · 모드 해석 · 비행 플러터 시험).
• CFD·시뮬레이션 후처리. 와류 방출·한계 사이클의 주파수와 성장률을 짧은 시계열에서 추정하거나, 안정성 해석에서 과도 응답의 지배 고유값을 뽑는다. 스냅숏 전체 장(field)으로 일반화한 것이 동적 모드 분해이고, DMD 의 고유값-진폭 구조는 정확히 프로니의 2·3단계를 고차원으로 한 것이다.
• 전자기 해석의 시간영역 외삽. FDTD 같은 시간영역 솔버가 공진기 응답을 다 감쇠할 때까지 돌리는 대신, 앞부분에서 극을 뽑아 나머지를 외삽해 계산 시간을 줄인다. 고-Q 공진기일수록 이득이 크다.

5. 실무 체크리스트

• 차수 선택이 전부다. 데이터 행렬의 특이값을 그려 낭떠러지를 보고 $p$ 를 정하되, 1단계 회귀에는 그보다 넉넉한 차수를 써라. 위 실험에서 보듯 차수가 모자라면 편향, 넉넉하면 여분의 극이 생기니 진폭·위치로 골라낸다.
• 표본 간격. 최고 주파수가 나이퀴스트 밑이어야 하는 건 당연하고, 반대로 너무 촘촘히 표본을 찍으면 $z_i$ 가 모두 1 근처에 몰려 근 찾기와 반데르몽드 진폭 단계가 악조건이 된다. 관심 모드가 한 주기에 수십 표본 정도가 되도록 데시메이션하는 것이 국룰이다.
• 고차 다항식 근은 민감하다. 차수 20 이상 다항식의 근은 계수 섭동에 민감하다(윌킨슨 다항식의 교훈). 과차수 프로니를 쓸 땐 근 찾기 대신 펜슬·ESPRIT 처럼 고유값으로 직접 가는 쪽이 안정하다.
• 결과는 잔차로 검증하라. 뽑은 극·진폭으로 신호를 재구성해 원 신호와 비교하고, 잔차가 백색 잡음처럼 보이는지 확인한다. 극 하나를 뺐을 때 잔차가 얼마나 늘어나는지 보면 그 모드가 진짜인지 잡음 흡수용인지 갈린다³.

6. 관련 문서

• 행켈 행렬 · 토플리츠 행렬 · 반데르몽드 행렬 · 특이값 분해 · 저랭크 근사
• 자기회귀 모형 · 레빈슨-더빈 알고리즘 · 총최소제곱 · 최소자승법 · 레벤버그-마쿼트 방법
• 시스템 식별 · 모드 해석 · 동적 모드 분해 · 파데 근사 · 고속 푸리에 변환
• 행렬 펜슬 방법 · ESPRIT · 선형 예측 부호화

7. Footnotes

1. 프로니는 같은 시기에 프랑스 혁명 정부의 미터법 개편용 거대 로그·삼각함수표 편찬을 지휘했다. 「계산원」 수십 명에게 덧셈·뺄셈만 시키는 분업 체계였는데, 이게 배비지의 차분기관 구상에 영감을 줬다는 이야기가 유명하다. 지수함수 맞춤도, 표 계산도 결국 「비싼 연산을 싼 연산으로 쪼개기」다. ↩

2. 「표본 8개로 FFT 가 못 가르는 두 주파수를 가른다」는 잡음 0 이라는 조건 아래의 이야기다. 모형이 정확히 맞을 때 얻는 해상도는 모형이 틀리는 순간 같이 증발한다. 초해상도는 공짜가 아니라 모형에 대한 외상이다. ↩

3. 프로니가 뽑아 준 극 20개를 전부 「물리 모드」로 보고서에 쓰면 리뷰어에게 털린다. 잡음 흡수용 극은 진폭이 작고 감쇠가 크며, 데이터 구간을 바꾸면 위치가 멋대로 춤춘다. 진짜 모드는 구간을 바꿔도 제자리에 있다. ↩