연주 시 연동되는 흰색 잡음이 포함된 리드 악기의 지연 호프 분기 예측

초록

본 논문은 입구 압력(블로잉 압력)이 시간에 따라 선형적으로 증가하면서 정적 호프 분기점을 통과할 때, 흰색 잡음이 가해진 단일 리드 악기 모델의 동적 거동을 분석한다. 확률 평균법을 이용해 느린 동역학을 도출하고, 잡음이 무시되는 경우와 무시될 수 없는 경우를 각각 해석하여 동적 분기점(소리가 실제로 시작되는 압력값)을 구한다. 이론적 결과는 수치 시뮬레이션과 비교해 좋은 일치를 보이며, 잡음 수준과 압력 상승 속도가 분기 지연에 미치는 영향을 명확히 제시한다.

상세 분석

이 연구는 전통적인 정적 호프 분기점(정상 상태에서 진동이 시작되는 임계 압력)과는 달리, 실제 연주 상황에서 블로잉 압력이 서서히 증가할 때 나타나는 ‘동적 분기점’(소리가 실제로 발생하는 순간)을 정량화한다. 모델은 하나의 공명 모드만을 고려한 최소화된 연속시간 ODE 시스템으로, 리드의 비선형 흐름‑압력 관계를 3차 다항식으로 전개하고, 공기 컬럼의 임피던스를 단일 고유진동수와 감쇠비로 표현한다. 블로잉 압력 γ는 시간에 따라 γ(t)=γ₀+ε̂t 로 선형 증가하며, ε̂는 매우 작은 비율(ε̂≪1)이다. 또한, 실제 연주에서 발생하는 난류와 수치적 오차를 흰색 잡음 νξ(t) 형태로 추가한다.

확률 평균법을 적용해 빠른 진동 성분을 평균화하고, 진폭‑위상 변수를 도입해 1차 비자코프 방정식 형태의 확률 미분 방정식으로 축소한다. 결과적으로 얻어지는 비자코프 방정식은 비자율적(시간 의존) 이토형 SDE이며, 진폭 A(t)²의 평균 동역학을 기술한다. 여기서 핵심은 두 가지 근사 시나리오이다. 첫 번째는 잡음 항이 무시 가능한 경우(ν→0)로, 시스템은 결정론적 지연을 보이며, 초기 압력 y₀(=γ₀−γ̂_st)의 크기에 따라 동적 분기점 ŷ_dyn_det이 y₀보다 크게 이동한다. 두 번째는 잡음이 유의미한 경우로, 잡음 강도 ν와 압력 상승 속도 ε가 동적 분기점 ŷ_dyn_stoch에 직접적인 영향을 미친다. 이 경우, 평균 진폭이 초기값으로 복귀하는 순간을 확률적 첫 통과 시간(first‑passage time) 문제로 풀어, ŷ_dyn_stoch≈(ν²/2ε)·ln(·) 형태의 근사식을 얻는다. 논문은 이 근사식을 두 단계(첫 번째, 두 번째 근사)로 구분해 ŷ_dyn_stoch,a와 ŷ_dyn_stoch,b를 제시하고, 각각 잡음이 약한 경우와 강한 경우에 적합함을 보인다.

수치 검증에서는 Euler‑Maruyama 방법으로 SDE를 직접 통합하고, 결정론적 경우에는 고전적인 Runge‑Kutta를 사용한다. 파라미터(α₁, ω₁, ζ 등)를 실제 클라리넷에 맞춰 설정한 뒤, ε와 ν를 다양하게 바꾸어 동적 분기점의 변화를 관찰한다. 결과는 이론식과 거의 일치하며, 특히 잡음이 없는 경우에는 ŷ_dyn_det이 초기 조건에 민감하게 변하고, 잡음이 존재하면 ŷ_dyn_stoch는 초기 조건에 거의 독립적이면서 ν와 ε에 비례·반비례 관계를 보인다. 이는 기존 이산시간 클라리넷 모델에서 보고된 현상과 연속시간 모델에서도 동일하게 적용됨을 확인한다. 마지막으로, 논문은 이러한 분석이 실제 연주자에게 블로잉 압력 제어 전략을 제공하고, 디지털 악기 시뮬레이션에서 잡음 모델링의 중요성을 강조한다는 점을 강조한다.