3 차 스플라인 보간이 보간 다항식보다 나은 경우는 언제입니까?


9

다음 그림은 교과서의 예를 약간 변형 한 것입니다. 저자는 동일한 간격으로 배치 된 샘플에 대한 보간 다항식이 보간 간격의 끝 근처에서 큰 진동을 가짐을 설명하기 위해이 예제를 사용했습니다. 물론 입방 스플라인 보간은 전체 간격에 걸쳐 좋은 근사치를 제공합니다. 몇 년 동안, 여기에 설명 된 이유로 같은 간격의 샘플에 대한 고차 다항식 보간을 피해야한다고 생각했습니다.

여기에 이미지 설명을 입력하십시오

그러나 최근에 고차 보간 다항식이 입방 스플라인 보간보다 근사 오차가 적은 대역 제한 신호의 많은 예를 발견했습니다. 일반적으로 보간 다항식은 샘플 속도가 충분히 높을 때 전체 보간 간격에 대해 더 정확합니다. 이것은 샘플이 신호의 나이키 스트 주파수보다 적어도 3 배 더 큰 샘플 속도로 동일한 간격으로있을 때 유지되는 것으로 보입니다. 또한, 큐빅 스플라인 보간에 대한 이점은 (샘플 속도) / (Nyquist frequency)가 증가함에 따라 개선됩니다.

예를 들어, 나이 퀴 스트 주파수가 2Hz이고 샘플 속도가 6.5Hz 인 사인파에 대해 3 차 스플라인 보간과 보간 다항식을 비교합니다. 샘플 포인트 사이에서 보간 다항식은 실제 신호와 정확히 동일하게 보입니다. 여기에 이미지 설명을 입력하십시오


아래에서는 두 근사치의 오류를 비교합니다. 첫 번째 예에서와 같이 다항식 보간은 샘플 간격의 시작과 끝 근처에서 최악입니다. 그러나 보간 다항식은 전체 샘플 간격에서 3 차 스플라인보다 오류가 적습니다. 보간 다항식은 작은 간격으로 외삽 할 때도 오류가 적습니다. 잘 알려진 사실을 발견 했습니까? 그렇다면 어디서 읽을 수 있습니까?

여기에 이미지 설명을 입력하십시오


수식이나 데이터를 근사합니까? 공식과 마찬가지로 주어진 고차 파생 상품도 고려되는 고급 스플라인을 항상 사용할 수 있습니다. 또한 큐빅 스플라인이 특정 "에너지"기능을 최소화한다는 사실을 확인해야합니다. wikipedia en.wikipedia.org/wiki/Spline_interpolation을보십시오 . 어떤 의미에서는 곡률을 최소화하면 더 잘 할 수 없습니다. 다른 해석은 입방 스플라인이 피팅에 사용되었다는 것입니다. 근사하지 않습니다. "피팅"은 특정 메트릭이 최적화됨을 의미합니다.
rrogers

@ rrogers, 보간 다항식은 측정 된 샘플에서 함수를 추정하고 신호의 대역폭이 샘플 속도의 1/6 미만인 것으로 알려진 경우 더 나은 접근 방법이 될 것이라고 생각했습니다. It
테드 Ersek

@TedErsek : 하나의 정 성적 고려 : 다항식 함수는 본질적 으로 가로 좌표 변수 로 로 분기됩니다 . 이 효과는 다항식 차수가 증가함에 따라 악화됩니다. 첫 번째 예에서 근사 될 신호는 보간 간격의 끝 근처에서 0으로 감소합니다. 이것은 보간 기의 점근 적 행동과 호환되지 않습니다. 두 번째 그림에는 구간의 가장자리 근처에 가파른 기울기와 0이 아닌 값이 있으므로 더 나은 근사값을 얻을 수 있습니다. 여기서 이론적 인 것은 아니고 단지 관측 일뿐입니다. ±
Jason R

@TedErsek Ted Ersek의 의견을 다루는 실질적인 측면으로서; 합리적인 다항식 근사를 시도 했습니까? BTW : 저는 1 년 전부터 무료로 제공되는 곡선 공식 추정 프로그램의 무료 사본을 가지고 있습니다. 프로그램이 베타에서 결제로 전환되었으므로 현재 버전이 없습니다.
rrogers

@JasonR 나는 당신에게 나의 마지막 코멘트를 다룰 예정이었습니다. 주제를 되돌아 가면, 어떤 경우 에도 함수를 알고 있다면 다항식에서 균일 한 오차 (최소 / 최대)를 제공 하는 en.wikipedia.org/wiki/Chebyshev_polynomials 가 있습니다. 그러나 함수를 알고 있다면 항상 "일치 필터"를 합성 할 수 있습니다.
rrogers

답변:


4

논의되는 현상은 Runge 's 현상 이다.

의 최대 절대 값 의 파생 상품 (ω) 이다 ω. 대한 룽게의 기능 1252+1 의 최대 절대 값 th (짝수) 도함수는 5!, 어디 !계승을 나타냅니다. 이것은 훨씬 빠른 성장입니다. 파생 상품이 증가하여 너무 빠르게 성장하는 경우에만보간 순서가 증가함에 따라 보간 오차가 발산 될 수있다. 지수아직 너무 빠르지 않습니다. James F. Epperson , The Runge 예제 , The American Mathematical Monthly , vol. 94, 1987, 329-341 쪽.

함수에 연속 도함수 만있는 경우, 경쟁 방식 인 부분 별 다항식 스플라인 보간은 소수의 초기 도함수가 관심 구간에 묶여 있으면 항상 수렴됩니다 ( 예 : 선형 보간에 대한 Wikipedia 기사 참조) .

두 방법이 수렴하면 (비 조각 별) 다항식 보간은 많은 샘플이 사용되는 경우 다항식도가 더 높으며 사인 예제에서 보듯이 더 나은 근사치를 제공 할 수 있습니다. 또한 LN Trefethen, 균등 한 점에서의 다항식 보간에 대한 두 가지 결과 , Journal of Approximation Theory 65, Issue 3, 1991 년 6 월, 페이지 247-260에 관심이있을 수 있습니다. 인용문:

복잡한 지수 함수의 대역 제한 보간에서 [...] 이자형나는α엑스(α아르 자형), 오류는 0 같이 만약에 α 파장 당 최소 6 개의 포인트를 제공 할만큼 작습니다.

파장 당 6.5 개의 샘플이 있습니다.

당사 사이트를 사용함과 동시에 당사의 쿠키 정책개인정보 보호정책을 읽고 이해하였음을 인정하는 것으로 간주합니다.
Licensed under cc by-sa 3.0 with attribution required.