Savitzky‑Golay 필터 차수 선택을 위한 N‑fold 교차검증 기법

초록

본 논문은 Savitzky‑Golay(SG) 스무딩 필터의 다항식 차수를 N‑fold(Leave‑One‑Out) 교차검증으로 자동 선택하는 방법을 제안한다. SG 필터가 최소 노름(minimum‑norm) 문제로 표현될 수 있다는 점을 이용해, 전체 예측 오차와 SG 투영 공간 사이의 새로운 수학적 관계를 도출함으로써 기존 BIC 기반 방법보다 낮은 SNR·샘플 수 환경에서도 정확한 차수 선택이 가능함을 보였다. 또한 행렬 연산을 QR 기반 직교화로 최적화해 계산 복잡도를 O(N³)→O(N²) 수준으로 크게 낮추었다. MATLAB 구현 코드와 실험 결과도 제공한다.

상세 분석

Savitzky‑Golay(SG) 필터는 입력 시퀀스를 일정 길이 N 의 윈도우 안에서 P 차 다항식 부분공간에 최소제곱 투영(projection)한 뒤, 중앙 샘플을 추출하는 FIR 구조로 이해된다. 이때 필터 계수 h 는 Aᵀ(AAᵀ)⁻¹c 형태의 최소 노름 해이며, A 는 0차부터 P 차까지의 다항식 값을 행렬로 나열한 것이다. 논문은 이러한 구조를 이용해 “예측 필터”(prediction filter)와 “스무딩 필터”(smoothing filter)의 차이를 정확히 수식화한다. 예측 필터는 중앙 샘플을 제외한 N‑1 개의 관측값만을 사용해 중앙 값을 복원하도록 설계되며, 이는 A 의 k‑번째 열을 0으로 만든 A′ 와 동일한 선형 시스템으로 표현된다.

핵심적인 수학적 변환은 hₚₖ = hₛₖ – γₖeₖ (식 7) 형태로, 여기서 γₖ = 1/(1–