⚡ electrical engineering

Fractional Order Filters with √f Characteristic

이 논문은 스킨 효과를 모델링하기 위해 맥스웰 방정식으로부터 유도된 편미분 방정식을 공간적으로 이산화함으로써 $\sqrt{f}$ 및 $1/\sqrt{f}$ 전달 함수를 갖는 분수 차수 필터를 유도하고, Oustaloup 재귀 필터와 벤치마크를 수행하는 쌍선형 변환을 통한 그에 상응하는 디지털 구현을 제시한다.

원저자: Hans Georg Brachtendorf

게시일 2026-09-14

📖 4 분 읽기☕ 가벼운 읽기

CC BY 4.0

원저자: Hans Georg Brachtendorf

원본 논문은 CC BY 4.0 (https://creativecommons.org/licenses/by/4.0/) 라이선스로 제공됩니다. ✨ 이것은 아래 논문에 대한 AI 생성 설명입니다. 저자가 작성하거나 승인한 것이 아닙니다. 기술적 정확성을 위해서는 원본 논문을 참조하세요. 전체 면책 조항 읽기

금속을 통해 흐르는 전기의 숨겨진 세계에는 신호의 속도에 따라 신호의 거동을 변화시키는 독특한 현상이 존재합니다. 전류가 도체를 통과할 때, 전류는 항상 전선 전체를 균일하게 채우지는 않습니다. 매우 높은 속도에서 전류는 외표면 쪽으로 밀려나는데, 이러한 현상을 표피 효과(skin effect)라고 합니다. 이 효과는 단순한 호기심의 대상이 아니라, 전기 신호가 어떻게 필터링되고 처리되는지를 근본적으로 변화시킵니다. 엔지니어들은 신호의 강도가 주파수의 제곱근에 비례하여 변화하는 특정 수학적 관계를 모방하는 전자 부품을 만들기 위해 오랫동안 노력해 왔습니다. 이러한 부품들은 민감한 측정에서 노이즈를 제거하고, 복잡한 기계 시스템을 제어하며, 생체 조직이 전자기파와 상호작용하는 방식을 모델링하는 데 필수적입니다. 수십 년 동안 이를 달성하기 위한 최선의 도구는 제한된 범위 내에서만 잘 작동하는 복잡한 근사법들이었으며, 이는 위상 응답을 평탄하게 만들기 위해 정교한 극점(poles)과 영점(zeros)의 배치를 요구하곤 했습니다.

오버 오스트리아 응용과학대학교(University of Applied Sciences Upper Austria)의 한 연구자는 이제 회로 근사법에서 시작하는 것이 아니라, 물리학의 근본 법칙 자체에서 시작하는 다른 경로를 제안했습니다. 저자는 금속 막대 내부에서 전기장과 자기장이 어떻게 상호작용하는지를 설명하는 방정식들을 분석하고 이를 작고 다룰 수 있는 단계들로 분해함으로써, 이러한 특수한 필터를 구축하는 새로운 방법을 유도해 냈습니다. 표준 회로가 제곱근 필터처럼 작동하도록 강제하는 대신, 이 방법은 물리적인 표피 효과를 직접적인 디지털 표현으로 삼아 자연스럽게 그렇게 작동하는 필터를 구축합니다. 이 연구는 계산 단계들을 세심하게 배치함으로써, 즉 현상이 발생하는 막대의 표면에 더 가깝게 배치함으로써, 이전 방식보다 훨씬 넓은 주파수 범위에서 매우 정확하게 작동하는 모델을 만들 수 있음을 보여줍니다.

여정은 단순한 물리적 형상, 즉 긴 직사각형 금속 막대로부터 시작됩니다. 연구자는 재료가 균일하며 내부의 전기장과 자기장이 전자기학의 표준 규칙, 구체적으로 맥스웰 방정식(Maxwell's equations)을 따른다고 가정합니다. 표피 효과가 강한 영역에서는 전류가 금속 내부로 침투하는 깊이가 막대의 두께에 비해 매우 작습니다. 이러한 물리적 실체는 주파수와 신호 손실 사이의 특정한 관계를 규정합니다. 이 연속적인 물리 법칙을 컴퓨터가 사용할 수 있는 형태로 바꾸기 위해, 연구자는 막대를 일련의 작은 조각들로 나눕니다. 공간 이산화(spatial discretization)라고 알려진 이 과정은 물리학의 매끄럽게 흐르는 방정식들을 상미분 방정식의 사슬로 변환합니다. 이 방정식들은 저항과 커패시터로 구성된 전기 회로에 직접 매핑될 수 있으며, 이를 통해 들어오는 신호에 대한 막대의 응답을 모방하는 행동 모델을 생성합니다.

이 연구의 핵심적인 통찰은 이 조각들이 어떻게 배치되는가에 있습니다. 만약 조각들의 크기가 모두 같다면, 금속 표면 바로 근처에서 발생하는 급격한 변화를 포착하기 위해 비현실적으로 많은 수의 조각이 필요할 것입니다. 대신, 연구자는 표면에 가까워질수록 조각들이 점진적으로 작아지는 비균일 격자(non-uniform grid)를 사용합니다. 이러한 계산 지점의 밀집은 가장 중요한 변화가 도체의 가장 얇은 층에서 일어난다는 물리적 실체를 반영합니다. 연구는 인접한 조각들 사이의 크기 비율이 결정적이라는 것을 발견했습니다. 만약 표면에서 멀어짐에 따라 조각의 크기가 너무 빠르게 커지면, 모델은 불안정해지고 원치 않는 링잉(ringing) 현상을 발생시킵니다. 시뮬레이션 결과에 따르면, 인접한 간격의 비율은 특정 수학적 상수 미만으로 유지되어야 하며, 약 2.5라는 값이 정확도와 커버 가능한 주파수 범위 사이의 최적의 균형을 제공하는 것으로 나타났습니다.

아날로그 모델이 설정되면, 다음 단계는 이를 현대적인 프로세서에서 실행할 수 있는 디지털 형식으로 변환하는 것입니다. 연구자는 연속 시간 시스템을 이산 시간 시스템으로 변환하는 몇 가지 표준적인 방법들을 테스트했습니다. 쌍선형 변환(bilinear transformation)으로 알려진 한 대중적인 방법이 명확한 승자로 떠올랐습니다. 인위적인 댐핑을 도입하거나 하드웨어 구현을 복잡하게 만드는 고차 계산을 요구하는 다른 기술들과 달리, 쌍선형 변환은 원래의 물리 모델의 안정성과 정확성을 보존합니다. 결과물인 디지털 필터는 원하는 제곱근 관계를 넓은 스펙트럼에 걸쳐 유지하면서, 물리적인 표피 모델과 정확히 동일하게 작동합니다. 시뮬레이션에서 이 새로운 접근 방식은 분수 차수 시스템(fractional-order systems)을 근사하는 데 널리 쓰이는 표준 기술인 우스탈루프 재귀 필터(Oustaloup recursive filter)보다 우수함을 입증했습니다. 우스탈루프 방식은 위상 응답을 평탄하게 하기 위해 로그 스케일 상에 극점과 영점을 배치하지만, 여기서 유도된 표피 효과 모델은 자연스럽게 더 높은 위상 정확도를 달상하며 우스탈루프 모델보다 더 큰 주파수 범위에서 유효합니다.

시뮬레이션 결과는 놀랍습니다. 단 8개의 격점(grid points)을 사용하는 모델의 경우, 이 새로운 필터는 시작 주파수의 약 6,500배에 달하는 범위 동안 그 특성을 유지합니다. 이는 정확도를 유지하는 데 어려움을 겪는 기존 방식들에 비해 상당한 개선입니다. 이 연구는 표피 효과 뒤에 숨겨진 물리적 직관이 이러한 필터를 설계하는 데 견고한 토대를 제공하며, 복잡한 최적화 알고리즘의 필요성을 우회한다는 것을 확인시켜 줍니다. 저자는 현재의 연구가 주파수의 제곱근이라는 특정 사례에 집중하고 있지만, 편미분 방정식을 이산화하는 근본적인 방법론은 향후 다른 분수 차수로 확장될 가능성이 있다고 언급했습니다. 현재로서는, 이 연구는 정밀한 신호 처리가 필요한 엔지니어들에게 물리적으로 근거가 확실한 대안을 제공하며, 때로는 디지털 문제를 해결하는 가장 좋은 방법이 실제 세계의 물리학을 면밀히 살펴보는 것임을 증명하고 있습니다.

기술 요약: $\sqrt{f}$ 특성을 갖는 분수 차수 필터

문제 정의

본 논문은 전달 함수가 $s^{\alpha}$ (여기서 $\alpha = \pm 1/2$ 에 해당하여 $\sqrt{f}$ 및 $1/\sqrt{f}$ 주파수 거동을 나타냄)에 비례하는 특성을 갖는 분수 차수 필터의 설계 및 분석을 다룬다. 이러한 필터는 프리엠퍼시스(preemphasis)/디엠퍼시스(deemphasis), $1/f$ (플리커) 노이즈 전력 스펙트럼 모델링, 분수 차수 제어, 그리고 생체 조직과 같은 매질의 분산 특성 시뮬레이션 등의 응용 분야에서 필수적이다.

Oustaloup 재귀 필터나 Charef 방법과 같은 기존의 근사 기법들은 위상 또는 크기 응답을 평탄하게 만들기 위해 로그 스케일 상에 극점(pole)과 영점(zero)을 배치하는 방식에 의존한다. 본 논문은 순수하게 수학적인 근사에 기반한 방식이 아닌, 물리 법칙에 뿌리를 둔 대안적인 유도 방식을 모색하며, 맥스웰 방정식의 이산화를 통해 원하는 분수 차수 특성을 자연스럽게 나타내는 행동 모델을 도출하는 것을 목표로 한다.

방법론

1. 맥스웰 방정식으로부터의 물리적 유도

핵심 방법론은 선형, 등방성 및 균질한 도전성 막대(conducting rod)에 대한 맥스웰 방정식에서 시작된다. 변위 전류가 무시 가능하다는 가정( $J \gg \partial D/\partial t$ ) 하에, 그리고 길이 $l$ 이 폭 $b$ 및 두께 $d$ 보다 훨씬 큰( $l \gg b, d$ ) 직사각형 막대에 대한 준정적(quasistationary) 접근법을 활용하면, 방정식은 전기장( $E_z$ )과 자기장( $H_x$ )을 기술하는 편미분 방정식(PDE) 시스템으로 축소된다.

본 논문은 스킨 효과가 강한 영역(스킨 깊이가 기하학적 구조에 비해 작은 경우)에서 도전성 물질의 스킨 깊이( $\Delta x$ )가 $1/\sqrt{f}$ 법칙을 따른다는 알려진 물리적 현상을 활용한다. 목표는 이러한 거동을 모사하는 필터를 만들기 위해 상기 PDE를 이산화하는 것이다.

2. 공간 이산화 및 회로 모델링

PDE는 유한 차분(Finite Difference, FD) 기법을 사용하여 이산화된다.

균등 격자(Uniform Grid): 초기에는 등간격 격자가 고려되며, 이는 상미분-대수 방정식(DAE) 시스템으로 이어진다. 이 시스템은 인덕터와 저항(LR) 또는 커패시터와 도체(GC)로 구성된 수동 행동 회로 모델(래더 네트워크)로 매핑된다.
비균등 격자(Non-Uniform Grid): 균등 격자는 고주파수에서 스킨 효과를 포착하기 위해 실용적이지 않을 정도로 많은 점을 필요로 한다는 점을 인식하여, 본 논문은 비균등 격자를 제안한다. 격자 간격은 스킨 효과가 가장 두드러지는 표면( $y = d/2$ )에서 가장 작다. 인접한 간격 사이의 비율은 파라미터 $c$ 에 의해 정의된다.
회로 구현: 이산화된 방정식은 등가 수동 회로로 구현된다. 본 논문은 LR 모델이 자연스럽기는 하지만, 수정 노드 해석법(MNA)에서 미지수의 수를 줄여주는 GC 모델(커패시터와 도체 사용)이 SPICE 유사 시뮬레이터에 더 적합하다고 언급한다.

3. 디지털 변환

디지털 등가물을 얻기 위해, 아날로그 필터는 쌍선형 변환(bilinear transformation, 사다리꼴 공식)을 사용하여 변환되며, 후방 차분 공식(Backward Differentiation Formulas, BDF)과 비교된다.

쌍선형 변환: 허수축을 단위 원으로 매핑하면서도 BDF 방법에서 관찰되는 수치적 댐핑(numerical damping)을 도입하지 않기 때문에 선호된다. 또한, 이는 $p \ge 2$ 인 고차 BDF 방법보다 낮은 필터 차수를 결과로 가져온다.
구현: 특정 재귀 단계에서 곱셈기를 피하기 위해 보조 벡터를 도입하는 것을 포함하여, 결과적으로 도출된 차분 방정식의 효율적인 구현 전략을 논의한다.

주요 기여

분수 차수 필터의 물리적 유도: 본 논문은 극점-영점 배치 휴리스틱에만 의존하는 대신, 도전성 막대에 대한 맥스벨 방정식의 이산화로부터 $\sqrt{f}$ 및 $1/\sqrt{f}$ 필터를 직접 유도하여 분수 차수 거동에 대한 물리적 근거를 제공한다.
비균등 격자 최적화: 특정 간격 비율( $c$ )을 갖는 비균등 격자를 사용하는 것이 정확도를 유지하면서도 시스템 차수를 크게 줄일 수 있음을 입증한다. 본 논문은 링잉(ringing) 현상을 피하기 위해 인접 격자 간격의 비율이 오일러 상수보다 작아야 한다고 명시하며, $c=2$ 에서 $c=2.5$ 가 최적임을 찾아냈다.
행동 회로 모델: 유도된 모델은 수동 래더 회로(LR 및 GC 모두)이며, 오믹 손실(Ohmic losses) 덕분에 극점이 좌반평면에 머물게 되어 본질적으로 안정적이다.
디지털 구현: 쌍선형 변환을 통해 디지털 대응물을 구현하는 명확한 경로를 제공하며, 이는 안정적이고 효율적인 디지털 필터 구조를 제공한다.

결과

Oustaloup와의 비교: 필터 차수 $K=11$ 에 대한 시뮬레이션 결과, 본 논문에서 유도된 스킨 깊이 모델이 Oustaloup 재귀 필터에 비해 더 넓은 근사 유효 주파수 범위와 높은 위상 정확도를 제공함을 보여준다.
주파수 범위: $K=8$ 개의 점과 $c=2$ 를 갖는 비균등 격자에 대해, 모델은 약 4개 자릿수(decades)의 주파수 범위에서 예상된 $\sqrt{f}$ 거동(임피던스의 이중 로그 플롯에서 0.5의 선형 기울기)을 보인다. 이론적 추정치는 더 넓은 범위를 시사했으나, 시뮬레이션 결과 실질적인 유효 범위는 약 $f_1/10$ 인 것으로 나타났다.
파라미터 민감도: 테스트를 통해 격자 간격 비율 $c$ 의 값이 오일러 상수( $\approx 2.718$ )와 $2 $사이일 때 시스템 차수, 유효 범위 및 정확도 사이의 최적의 절충안을 제공함을 확인했다.$ c$가 오일러 상수를 초과하는 값은 바람직하지 않은 링잉 현상을 초래한다.
디지털 성능: 쌍선형 변환을 통해 구현된 디지털 대응물은 정규화된 각주파수에 대해 약 4개 자릿수의 제곱근 거동을 유지함으로써 변환의 효과를 입증했다.

의의 및 주장

본 논문은 제안된 방법이 유효 범위 및 위상 정확도 측면에서 Oustaloup와 같은 기존의 재귀적 근사법보다 우수한 대안을 제공한다고 주장한다. 도전체의 물리적 스킨 효과에 기반을 둠으로써, 결과 모델은 수학적으로 타당할 뿐만 아니라 수동 네트워크로서 물리적으로 해석 가능하다.

저자는 비균등 이산화 전략을 통해 (균등 이산화보다) 훨씬 적은 격자점으로 스킨 효과를 정확하게 모델링할 수 있음을 강조한다. 쌍선형 변환의 선호도는 안정성을 보존하고 수치적 댐핑을 피할 수 있는 능력에 의해 정당화되며, 이는 하드웨어 구현에 적합하다.

마지막으로, 본 논문은 이 연구를 분수 차수 필터링( $\alpha = \pm 1/2$ )의 특수한 사례로 위치시키며, 편미분 방정식을 이산화하는 기본 방법론이 향후 $\pm 1/2$ 을 넘어 임의의 차수 $\alpha$ 를 갖는 분수 차수 필터로 일반화될 수 있음을 시사한다.

기술 요약: f\sqrt{f}f​ 특성을 갖는 분수 차수 필터