역문제

편집 역사 토론
수치해석 통계 마지막 수정: 2026-07-12 04:07:33

1. 개요[편집]

역문제
Inverse Problem
대비 개념순문제(forward problem)
핵심 난점부적절성(ill-posedness)
대표 처방티호노프 정칙화(Tikhonov regularization)
응용토모그래피, 파라미터 추정, 지구물리 탐사

순문제: 원인을 알 때 결과를 구한다. 역문제: 결과를 보고 원인을 되짚는다. 후자가 훨씬 어렵다.

역문제(Inverse Problem)는 관측된 결과(데이터)로부터 그 결과를 만들어낸 원인(모델 파라미터·입력·계수)을 역추정하는 문제다. 정방향으로 “원인 → 결과”를 계산하는 것을 순문제(forward problem)라 부르는데, 물리 법칙과 초기·경계조건이 주어지면 순문제는 대체로 잘 정의되어 있고 계산도 안정적이다. 문제는 그 반대 방향이다. CT 스캐너가 측정한 X선 감쇠량으로부터 몸속 밀도 분포를 재구성하거나, 지진파 도착 시간으로 지구 내부 구조를 알아내거나, 관측된 온도장으로 미지의 열전도 계수를 찾아내는 일 — 전부 역문제다.

역문제가 악명 높은 이유는 순문제와 달리 부적절성(ill-posedness)에 시달리기 때문이다. 데이터의 미세한 잡음이 답을 폭발적으로 뒤흔들고, 답이 여러 개거나 아예 없을 수도 있다. 그래서 역문제 풀이의 절반은 “물리”가 아니라 “이 불안정성을 어떻게 길들일 것인가”라는 수치해석 싸움이다.

2. 아다마르의 적절성 조건[편집]

프랑스 수학자 자크 아다마르(Jacques Hadamard)는 잘 정의된 문제(well-posed problem)가 갖춰야 할 세 조건을 제시했다.1

  1. 존재성(existence): 해가 적어도 하나 존재한다.
  2. 유일성(uniqueness): 해가 오직 하나뿐이다.
  3. 안정성(stability): 데이터가 조금 바뀌면 해도 조금만 바뀐다(연속 의존성).

이 셋 중 하나라도 깨지면 그 문제는 부적절(ill-posed)하다고 부른다. 그리고 대부분의 역문제는 이 셋을 아주 신나게 위반한다. 해가 여러 개일 수도 있고(유일성 위반), 잡음 섞인 데이터에는 정확히 맞는 해가 없을 수도 있으며(존재성 위반), 무엇보다 안정성이 거의 항상 깨진다. 아다마르는 원래 이런 문제를 “물리적으로 무의미하다”며 배척했지만, 20세기 후반 의료영상·지구물리가 발전하면서 “무의미하기는커녕 돈이 되는” 문제임이 드러났다.

3. 왜 불안정한가 — 조건수의 저주[편집]

역문제의 불안정성은 조건수라는 개념으로 정량화된다. 이산화된 역문제는 대개 선형 시스템 Ax=bA\mathbf{x} = \mathbf{b}로 귀결되는데, 여기서 b\mathbf{b}가 관측 데이터, x\mathbf{x}가 추정하려는 미지수다. 문제는 이 행렬 AA의 조건수가 어마어마하게 크다는 것.

κ(A)=σmaxσmin\kappa(A) = \frac{\sigma_{\max}}{\sigma_{\min}}

특이값 분해로 보면 그 정체가 명확해진다. AA의 작은 특이값 σmin\sigma_{\min}은 순문제에서 “정보가 뭉개져 사라지는 방향”에 대응한다. 역문제는 이걸 되돌려야 하므로 1/σmin1/\sigma_{\min}으로 나누게 되는데, 이 과정에서 데이터에 섞인 미세한 잡음이 1/σmin1/\sigma_{\min}배로 증폭된다.2 조건수가 10810^{8}이면, 소수점 8자리 정밀도로 측정해도 답의 유효숫자가 통째로 날아갈 수 있다는 뜻이다. 순진하게 최소자승법으로 minAxb2\min \|A\mathbf{x}-\mathbf{b}\|^2을 풀면 잡음을 충실히 재현하는 쓰레기 해가 나온다.

4. 정칙화 — 티호노프의 처방[편집]

이 폭주를 막는 표준 처방이 정칙화(regularization)다. 핵심 아이디어는 “데이터에 잘 맞으면서 동시에 얌전한(smooth·small) 해”를 고르도록 벌점(penalty)을 추가하는 것. 가장 유명한 티호노프 정칙화(Tikhonov regularization)는 최소자승 목적함수에 해의 크기 항을 더한다.

minx  {Axb2+λ2Lx2}\min_{\mathbf{x}} \;\left\{\, \|A\mathbf{x} - \mathbf{b}\|^2 + \lambda^2\, \|L\mathbf{x}\|^2 \,\right\}

여기서 λ\lambda는 정칙화 파라미터, LL은 항등행렬이거나 미분 연산자(해를 매끄럽게)를 쓴다. λ\lambda가 크면 잡음에 둔감하지만 답이 뭉개지고(과도한 평활화), 작으면 데이터엔 잘 맞지만 잡음이 되살아난다. 이 저울질이 역문제의 영원한 딜레마다. 적절한 λ\lambda를 고르는 방법으로 L-곡선(L-curve), 일반화 교차검증(GCV), 오차 원리(Morozov’s discrepancy principle) 등이 쓰인다.3

통계적으로 보면 티호노프 정칙화는 베이즈 추정의 사전분포(prior)와 정확히 대응한다. L2L^2 벌점은 가우시안 사전분포, L1L^1 벌점(LASSO·전변동 정칙화)은 라플라스 사전분포에 해당한다. 그래서 역문제는 결정론적 관점(정칙화)과 통계적 관점(베이즈 추론) 양쪽에서 접근할 수 있으며, 후자는 답뿐 아니라 그 불확실성까지 정량화해준다는 장점이 있다.

5. 응용 분야[편집]

역문제는 이름은 낯설어도 응용은 우리 삶 깊숙이 들어와 있다.

  • 토모그래피: CT, MRI, PET 등 의료영상은 전부 역문제다. 몸을 관통한 신호로 내부 단면을 재구성한다. 지구물리의 지진 토모그래피, 산업용 비파괴검사도 같은 계보.
  • 파라미터 추정: 시뮬레이션 모델의 미지 물성치(열전도율, 투과율, 재료 상수)를 실험 데이터에 맞춰 역산한다. CFD·구조해석 모델 보정의 핵심 절차이며, 데이터 동화(data assimilation)로 기상예보에도 쓰인다.
  • 비파괴 탐상·원격 탐사: 표면 측정만으로 내부 결함이나 지하 자원 분포를 추정.
  • 디컨볼루션: 흐릿한 영상, 번진 신호를 원래대로 복원하는 문제 역시 대표적 부적절 역문제.

시뮬레이션 관점에서 특히 중요한 것은 파라미터 추정이다. 검증 및 확인(V&V) 절차에서 모델을 현실에 맞추는 보정(calibration) 단계가 사실상 역문제 풀이이기 때문이다. 다만 여기엔 함정이 있다. 부적절한 역문제로 파라미터를 억지로 맞추면, 물리적으로 틀린 파라미터 조합이 우연히 데이터를 재현하는 과적합이 생긴다. “실험과 잘 맞으니 검증 완료”라는 착각이 가장 위험한 지점이며, 불확실성 정량화와 교차검증으로 방어해야 한다.4

6. 관련 문서[편집]

7. Footnotes[편집]

  1. Hadamard, J. (1902). 아다마르는 물리 문제라면 당연히 이 세 조건을 만족해야 한다고 믿었다. 반세기 뒤 CT 스캐너가 부적절 문제로 노벨상(1979년 코맥·하운스필드)을 받을 줄은 몰랐을 것이다.

  2. 그래서 역문제 연구자들 사이엔 “순문제는 정보를 지우고, 역문제는 지워진 정보를 상상해내는 일”이라는 말이 있다. 없는 정보를 만들어내야 하니 잡음에 취약할 수밖에.

  3. L-곡선은 잔차 노름과 해 노름을 로그-로그로 그리면 L자 모양이 나오고, 그 꺾이는 모서리가 최적 λ\lambda라는 휴리스틱이다. 예쁘게 꺾이면 다행이고, 안 꺾이면 그냥 감으로 찍는다.

  4. “파라미터 3개만 튜닝하면 코끼리도 그리고, 4개면 코를 흔들게 할 수 있다”는 폰 노이만의 농담이 정확히 이 위험을 지적한다. 자유도가 많은 역문제는 무엇이든 맞출 수 있다는 게 오히려 저주다.