확률 공명

1. 개요

확률 공명 Stochastic Resonance
재료	문턱(쌍안정) 비선형계 + 약한 주기 신호 + 잡음
현상	출력 신호 대 잡음비가 잡음 세기에 대해 최대를 가짐
표준 모형	과감쇠 이중우물 랑주뱅 방정식
핵심 척도	크라머스 탈출률 rK 와 구동 주기의 비교
제안	벤치·수테라·불피아니 (1981), 니콜리스 (1981)
관측	슈미트 트리거, 링 레이저, 가재 기계수용기 등

잡음을 줄이면 신호가 잘 보인다. 보통은. 그런데 문턱이 있는 시스템에서는 잡음을 조금 더 넣었더니 신호가 더 잘 보이는 구간이 있다.

확률 공명(stochastic resonance)은 문턱이나 쌍안정 구조를 가진 비선형계에 혼자서는 문턱을 못 넘는 약한 주기 신호를 넣었을 때, 적당한 세기의 잡음이 더해지면 출력에서 그 신호 성분의 신호 대 잡음비(SNR)가 오히려 커지고, 잡음 세기에 대해 최댓값을 갖는 현상이다. 잡음이 너무 약하면 문턱을 못 넘고, 너무 세면 신호와 상관없이 마구 넘으니, 그 사이 어딘가에 최적점이 있다.

1981년 벤치·수테라·불피아니와 니콜리스가 각각 빙하기 주기 문제에서 제안했다. 지구 공전 궤도 이심률이 약 10만 년 주기로 바뀌어 일사량이 아주 조금 흔들리는데, 그 약한 강제만으로는 빙하기-간빙기 전환이 설명되지 않으니 기후의 단기 변동(잡음)이 전환을 「도와준다」는 가설이었다¹. 기후 가설로서의 지위는 여전히 논쟁 중이지만, 메커니즘 자체는 전자 회로·레이저·감각 뉴런에서 실험으로 확인됐다.

크라머스 탈출률 문서가 잡음에 의한 장벽 넘기의 속도를 다뤘다. 이 문서는 그 속도가 바깥 리듬과 박자를 맞추는 이야기다.

2. 표준 모형

과감쇠 이중우물 랑주뱅 방정식이 교과서 모형이다.

$$\dot x = -U'(x) + A\cos\Omega t + \xi(t),\qquad U(x) = -\frac{x^2}{2} + \frac{x^4}{4},\qquad \langle \xi(t)\xi(s)\rangle = 2D\,\delta(t-s).$$

우물 바닥은 $x = \pm1$, 장벽 높이 $\Delta V = 1/4$, 곡률은 바닥 $U''(\pm1) = 2$, 장벽 꼭대기 $U''(0) = -1$ 이다. 구동 $A\cos\Omega t$ 는 퍼텐셜을 주기적으로 기울여 두 우물의 장벽 높이를 번갈아 $\Delta V \mp A$ 쯤으로 바꾼다. 잡음 없이 우물을 넘으려면 $A > 2/(3\sqrt3) \approx 0.385$ 가 필요한데, 아래 검산은 $A = 0.1$ 로 그 4분의 1 수준이다. 혼자서는 절대 못 넘는다.

잡음이 있으면 크라머스 탈출률(과감쇠 극한)

$$r_K = \frac{\sqrt{U''(\pm1)\,\lvert U''(0)\rvert}}{2\pi}\,e^{-\Delta V/D} = \frac{\sqrt2}{2\pi}\,e^{-1/(4D)}$$

로 우물 사이를 뛴다. 구동이 장벽을 번갈아 낮추니 탈출률도 주기적으로 오르내리고, 그 결과 점프가 구동 위상에 몰리기 시작한다.

3. 두 상태 이론 — SNR 의 봉우리

맥나마라-위젠펠드(1989)는 연속 변수를 두 상태 $x = \pm1$ 로 줄이고, 구동이 느리고 약하다는 가정(단열·선형 응답)에서 출력 파워 스펙트럼의 구동 주파수 성분과 배경을 계산했다.

$$\mathrm{SNR} \;\approx\; \pi \Bigl(\frac{A\,x_m}{D}\Bigr)^{2} r_K \;\propto\; \frac{1}{D^2}\, e^{-\Delta V/D}.$$

($x_m = 1$ 은 우물 위치.) $D$ 가 작으면 지수 $e^{-\Delta V/D}$ 가 죽이고, 크면 $1/D^2$ 가 죽인다. 미분해서 0으로 놓으면 최적점은

$$D_{\mathrm{opt}} = \frac{\Delta V}{2}$$

이다. 이 모형에서는 $D_{\mathrm{opt}} = 0.125$. 단열 근사라 이 위치는 구동 주파수 $\Omega$ 에 의존하지 않는다는 점도 기억해 두자.

3.1. 시간척도 정합

또 하나 자주 인용되는 기준이 있다. 평균 체류 시간 $1/r_K$ 가 구동 반주기 $\pi/\Omega$ 와 같아지면, 즉

$$r_K = \frac{\Omega}{\pi}$$

이면 매 반주기마다 평균 한 번씩 뛰어 「동기화」가 가장 잘 된다는 직관이다. 이건 SNR 최댓값과는 다른 점이다. 아래 검산에서 정합 조건은 $D \approx 0.080$ 인데 SNR 봉우리는 $D \approx 0.125$ 근처에 있다. 체류 시간 분포의 $T/2$ 봉우리처럼 동기화 자체를 재는 척도는 정합 근처에서 최대가 된다는 결과가 있어, 「무엇을 재느냐에 따라 최적 잡음이 다르다」가 정확한 표현이다².

4. 검산 — 오일러-마루야마 앙상블

위 모형을 $A = 0.1$, 주기 $T = 200$ ($\Omega = 0.0314$), $dt = 0.05$ 오일러-마루야마 방법(확률미분방정식 참고)으로 궤적 256개 × 40주기 적분했다(앞 1/8 은 과도구간으로 버림). 각 궤적의 주기도표(periodogram)를 평균한 뒤 구동 주파수 빈의 파워를 이웃 빈 배경과 비교해 SNR 을 쟀다.

$D$	실측 SNR (빈 단위)	이론 $\pi A^2 r_K/D^2$	$1/r_K$
0.040	5.0	0.0085	2302
0.060	13.4	0.0305	287
0.080	23.9	0.0485	101
0.100	28.3	0.0580	54
0.125	30.4	0.0612	33
0.150	29.7	0.0594	24
0.200	26.3	0.0506	16
0.300	20.2	0.0341	10

두 열의 절대값은 스펙트럼 규약(단측·양측, 주파수 빈 폭, 관측 길이) 때문에 상수배(여기선 약 440~590배) 차이 나므로 모양만 비교한다. 실측 SNR 은 $D = 0.125$ 에서 최대이고, 이론의 $D_{\mathrm{opt}} = \Delta V/2 = 0.125$ 와 정확히 같은 자리다. 잡음을 0.04 → 0.125 로 세 배 키웠더니 SNR 이 여섯 배 좋아졌다. 반대로 $D$ 를 더 올리면 천천히 떨어진다 — 봉우리가 비대칭인 것도 $D^{-2}e^{-\Delta V/D}$ 의 모양 그대로다.

$D = 0.04$ 에서 $1/r_K \approx 2300$ 은 관측 시간(약 7000) 동안 우물당 몇 번 못 뛴다는 뜻이다. 이 영역은 「잡음이 모자라」 신호가 안 보이는 쪽이다.

5. 오해와 한계

• “잡음이 정보를 만든다”가 아니다. 입력에 이미 있던 약한 신호를, 문턱 비선형성이 잃어버리지 않도록 잡음이 돕는 것이다. 신호가 작은 선형 응답 영역에서는 출력 SNR 이 입력 SNR 을 넘을 수 없다는 결과(드위즈-비알렉 1995)도 있다. 확률 공명은 「비선형 검출기가 잡아먹는 손실을 잡음이 줄인다」로 읽는 게 정확하다.
• 비선형성이 필수다. 선형 시스템에서는 잡음을 더하면 SNR 이 단조 감소한다. 쌍안정·문턱·흥분성(뉴런의 발화 문턱) 중 하나가 있어야 한다.
• 잡음 없이 신호를 키울 수 있으면 그게 낫다. 확률 공명은 신호를 못 키우는 상황(감각 뉴런, 이미 주어진 하드웨어)에서의 차선책이다.

6. 어디서 보이나

• 전자 회로. 포브-에슬로(1983)가 슈미트 트리거로 처음 실험 측정했다. 이중 문턱 회로 그 자체가 두 상태 모형이다.
• 레이저. 맥나마라·위젠펠드·로이(1988)가 양방향 링 레이저의 두 회전 모드 사이 전환에서 관측했다.
• 감각 생물학. 더글러스·모스 등(1993)이 가재 꼬리의 기계수용 뉴런에 약한 주기 자극과 잡음을 함께 주어, 잡음 세기에 따라 발화의 SNR 이 봉우리를 갖는 것을 보였다. 이온 채널의 확률적 개폐와 호지킨-헉슬리 모형의 발화 문턱이 바로 그 비선형성이다.
• 신호 처리. 양자화기에 일부러 잡음을 넣어 문턱 아래 신호를 살리는 디더링은 같은 원리를 공학적으로 쓰는 경우다³.
• 수치 실험. 확률미분방정식 적분기를 시험할 때 좋은 벤치마크다. 결과가 지수 $e^{-\Delta V/D}$ 에 민감해 시간 간격이 거칠면 탈출률이 틀어지고 봉우리 위치가 밀린다. 분포 쪽은 포커-플랑크 방정식을 주기 계수로 풀어 대조할 수 있다.

7. 관련 문서

• 크라머스 탈출률 · 랑주뱅 동역학 · 확률미분방정식 · 포커-플랑크 방정식
• 이중우물 퍼텐셜 · 아레니우스 방정식 · 평균 첫 통과 시간 · 파워 스펙트럼
• 이온 채널 · 호지킨-헉슬리 모형 · 디더링 · 위너 필터

8. Footnotes

1. 10만 년 주기의 빙하기를 설명하려고 나온 이론이 지금은 가재 꼬리와 보청기 연구에서 더 많이 인용된다. 논문의 운명은 저자도 모른다. ↩

2. 이 혼동 때문에 1990년대 중반에 “그게 진짜 공명이냐”는 논쟁이 붙었고, 체류 시간 분포로 잰 쪽을 ‘진짜(bona fide) 확률 공명’이라 부르는 논문까지 나왔다. 이름에 공명이 붙으면 물리학자들은 일단 공명 조건부터 따지고 본다. ↩

3. 오디오 엔지니어는 16비트로 줄일 때 디더를 꼭 넣으라고 한다. 잡음을 넣으면 소리가 좋아진다는 말이 처음엔 사기처럼 들리지만, 이 문서를 읽고 나면 꽤 그럴듯하게 들린다. ↩