두 변수 의 다항식 보간 에서, Padua 지점 은 그들 의 Lebesg 상수 의 최소 성장으로 증명된, O (로그 2 n ){\displaystyle O (\log ^{2}n)} 인, 비등결점 점 집합 의 첫 번째 예(그리고 지금까지 유일한 예)이다. [1] 이들의 이름은 원래 발견된 파두아 대학 덕분이다.[2]
포인트는 도메인 [ - 1 , 1 ] × [ - 1 ] ⊂ R 2 {\ displaystyle [-1,1]\times [- 1,1]\time [-1,1]\subset \mathb {R} ^{2 }}. 이후 90도 회전을 통해 얻은 포인트는 4개의 방향성으로 사용할 수 있다. 이렇게 하면 4개의 다른 패두아 포인트를 얻을 수 있다.
네 가족 첫 번째 패밀리와 5등급의 파두아 점들이 생성 곡선과 함께 표시된다. 첫 번째 패밀리의 Padua 포인트와 6도 포인트는 생성 곡선과 함께 표시된다. 우리는 파두아 점을 생성 곡선 이라 불리는 파라메트릭 곡선 의 "샘플링 "으로 볼 수 있는데, 이는 네 패밀리의 각 패밀리에 대해 약간씩 다르기 때문에 보간 정도 n[\displaystyle n} 과 패밀리 s {\displaysty s} 에 대한 점을 다음과 같이 정의할 수 있다.
패드 n s = { ξ = ( ξ 1 , ξ 2 ) } = { γ s ( k π n ( n + 1 ) ) , k = 0 , … , n ( n + 1 ) } . {\displaystyle {\text{Pad}}_{n}^{s}=\lbrace \mathbf {\xi } =(\xi _{1},\xi _{2})\rbrace =\left\lbrace \gamma _{s}\left({\frac {k\pi }{n(n+1)}}\right),k=0,\ldots ,n(n+1)\right\rbrace .} 실제로 파두아 지점은 곡선의 자기 교차점에 정확히 놓여 있고, 사각형의 경계와 함께 곡선의 교차점에 놓여 있다 [ - 1, 1 ] 2 {\ displaystyle [-1,1]^{2 }}. The cardinality of the set Pad n s {\displaystyle \scriptstyle {\text{Pad}}_{n}^{s}} is Pad n s = N = ( n + 1 ) ( n + 2 ) 2 {\displaystyle \scriptstyle {\text{Pad}}_{n}^{s} =N={\frac {(n+1)(n+2)}{2}}} . Moreover, for each family of Padua points, two points lie on consecutive vertices of the square [ - 1 , 1 ] 2 {\ displaystyle [-1,1]^{2 }}, 2n - 1 {\displaystyle 2n-1} 포인트는 정사각형의 가장자리에 놓여 있고, 나머지 포인트는 정사각형 내부 생성곡선의 자기 절편에 놓여 있다.[3] [4]
4개의 생성 곡선은 [ 0 , 2 π ]{\displaystyle [0,2\pi ] 의 구간 에서 닫힌 파라메트릭 곡선으로 리사주 곡선의 특별 한 경우다.
첫 번째 가족 첫 번째 패밀리의 Padua 포인트 생성 곡선은
γ 1 ( t ) = [ − cas ( ( n + 1 ) t ) , − cas ( n t ) ] , t ∈ [ 0 , π ] . {\displaystyle \pi _{1}(t)=[-\cos(nt)],-\cos(nt)],\cos(nt)],\caps t\in [0,\pi ]. } 위에서 설명한 대로 샘플을 채취하면 다음과 같은 결과를 얻을 수 있다.
패드 n 1 = { ξ = ( μ j , η k ) , 0 ≤ j ≤ n ; 1 ≤ k ≤ ⌊ n 2 ⌋ + 1 + δ j } , {\displaystyle {\text{Pad}}_{n}^{1}=\lbrace \mathbf {\xi } =(\mu _{j},\eta _{k}),0\leq j\leq n;1\leq k\leq \lfloor {\frac {n}{2}}\rfloor +1+\delta _{j}\rbrace ,} 여기 서 Δ j = 0 {\displaystyle \property_{ j}=0} 은(는) n {\ displaystyle n} 이 (가 ) 짝수이거나 홀수일 때, j {\displaystyle n}과 k {\displaysty k} 이(가 ) 모두 홀수일 경우 Δj = 1 {\daystytype _{j}
와 함께
μ j = cas ( j π n ) , η k = { cas ( ( 2 k − 2 ) π n + 1 ) j 기묘한 cas ( ( 2 k − 1 ) π n + 1 ) j 짝수의 {\displaystyle \mu _{j}=\cos \left({\frac {j\pi }{n}}\right),\eta _{k}={\begin{cases}\cos \left({\frac {(2k-2)\pi }{n+1}}\right)&j{\mbox{ odd}}\\\cos \left({\frac {(2k-1)\pi }{n+1}}\right)&j{\mbox{ even. }}}\end{case}} 여기서부터 첫 번째 패밀리의 Padua 지점은 n {\displaystyle n} 이 (가) 짝수일 경우 하단에 두 개의 꼭지점이 있고, n {\displaystyn} 이 (가) 홀수일 경우 왼쪽에 두 개의 정점이 있다.
두 번째 가족 두 번째 패밀리의 Padua 포인트 생성 곡선은
γ 2 ( t ) = [ − cas ( n t ) , − cas ( ( n + 1 ) t ) ] , t ∈ [ 0 , π ] , {\displaystyle \pi _{2}(t)=[-\cos(nt),-\cos(n+1)],\cos t\in [0,\pi ],} 따라서 n {\displaystyle n} 이 (가) 짝수일 경우 왼쪽에 정점이 있고, n {\displaystyle n} 이 (가) 홀수일 경우 아래쪽에 정점이 있다.
세 번째 가족 세 번째 패밀리의 Padua 포인트 생성 곡선은
γ 3 ( t ) = [ cas ( ( n + 1 ) t ) , cas ( n t ) ] , t ∈ [ 0 , π ] , {\displaystyle \pi \{3}(t)=[\cos(n+1)t],\cos(nt)],\cos(nt)],\caps t\in [0,\pi ],} 따라서 n {\displaystyle n} 이 (가) 짝수일 경우 상단에 정점이 있고, n {\displaystyle n} 이 (가) 홀수일 경우 오른쪽에 정점이 있다.
네번째 가족 네 번째 패밀리의 파두아 포인트 생성 곡선은
γ 4 ( t ) = [ cas ( n t ) , cas ( ( n + 1 ) t ) ] , t ∈ [ 0 , π ] , {\displaystyle \pi_{4}(t)=[\cos(nt),\cos((n+1)t)],\cos([0,\pi ],} n {\displaystyle n} 이 (가) 짝수일 경우 오른쪽에 정점이 있고, n {\displaystyle n} 이 (가) 홀수일 경우 상단에 정점이 있다.
보간식 The explicit representation of their fundamental Lagrange polynomial is based on the reproducing kernel K n ( x , y ) {\displaystyle \scriptstyle K_{n}(\mathbf {x} ,\mathbf {y} )} , x = ( x 1 , x 2 ) {\displaystyle \scriptstyle \mathbf {x} =(x_{1},x_{2})} and y = ( y 1 , y 2 ) {\displaystyle \s criptstyle \mathbf {y} = (y_{1},y_{2}) 공간 spacen 2 ([ - 1, 1 ] 2 ) {\ displaystyle \script style \Pi \{n}^{2}([-1,1]^{2})}
⟨ f , g ⟩ = 1 π 2 ∫ [ − 1 , 1 ] 2 f ( x 1 , x 2 ) g ( x 1 , x 2 ) d x 1 1 − x 1 2 d x 2 1 − x 2 2 {\displaystyle \langle f,g\rangle ={\frac {1}{\pi ^{2}}}\int _{[-1,1]^{2}}f(x_{1},x_{2})g(x_{1},x_{2}){\frac {dx_{1}}{\sqrt {1-x_{1}^{2}}}}{\frac {dx_{2}}{\sqrt {1-x_{2}^{2}}}}} 에 의해 정의된.
K n ( x , y ) = ∑ k = 0 n ∑ j = 0 k T ^ j ( x 1 ) T ^ k − j ( x 2 ) T ^ j ( y 1 ) T ^ k − j ( y 2 ) {\displaystyle K_{n}(\mathbf {x} ,\mathbf {y} )=\sum _{k=0}^{n}\sum _{j=0}^{k}{\hat {T}}_{j}(x_{1}){\hat {T}}_{k-j}(x_{2}){\hat {T}}_{j}(y_{1}){\hat {T}}_{k-j}(y_{2})} T ^ j {\ displaystyle \scriptstyle {\T}_{j} 이(가) 정규화된 Chebyshev polyomial of degree j {\displaystyle j}( T ^ 0 = T ^ 0 = T {\hat {T}_{0}=) 를 나타냄T_{0 }}, T ^ p = 2 T p {\ displaystyle \scriptstyle {\hat{T}_{p}={\sqrt{2}} T_ {p}} 여기서 T p ( ⋅ ) = cos (p arccos ( ⋅ ){\ displaystyle T_{p}(\cdot )=\cos(p\arccos(\cdot )}} 는 제1종 p {\displaystyp} 의 고전적인 체비셰프 다항식이다 .[3] For the four families of Padua points, which we may denote by Pad n s = { ξ = ( ξ 1 , ξ 2 ) } {\displaystyle \scriptstyle {\text{Pad}}_{n}^{s}=\lbrace \mathbf {\xi } =(\xi _{1},\xi _{2})\rbrace } , s = { 1 , 2 , 3 , 4 } {\displaystyle s=\lbrace 1,2,3,4\rbrace } , the interpolation formula of order n {\displa 함수 f 의 ystyle n} : [ - 1 , 1 ] 2 → R 2 {\ displaystyle \scriptstyle f\colon [-1]^{2}\to \mathb {R}{2}}: 일반 표적 x [ - 1 , 1 ] 2 {\ displaystyle \mathbf{x}{x}}}}}}}}}}}{ 2}}: 그 다음이다 .
L n s f ( x ) = ∑ ξ ∈ 패드 n s f ( ξ ) L ξ s ( x ) {\displaystyle {\mathcal {L}}_{n}^{s}f(\mathbf {x} )=\sum _{\mathbf {\xi } \in {\text{Pad}}_{n}^{s}}f(\mathbf {\xi } )L_{\mathbf {\xi } }^{s}(\mathbf {x} )} 여기서 L ξ s ( x ) {\ displaystyle \scriptstyle L_{\mathbf {\xi }^{s}(\mathbf {x} ) 는 기본 Lagrange 다항식이다.
L ξ s ( x ) = w ξ ( K n ( ξ , x ) − T n ( ξ i ) T n ( x i ) ) , s = 1 , 2 , 3 , 4 , i = 2 − ( s 모드의 2 ) . {\displaystyle L_{\mathbf {\xi}^{s}(\mathbf {x} )=w_{\\mathbf {n}(\mathbf {\xi},\mathbf {x})-T(\xi _{i})=w_{si} T_{n}(x_{i}),\quad s=1,2,3,4,\quad i=2-(s\mod 2). } w ξ {\ displaystyle \scriptstyle w_{\mathb{\xi}}}}}} 의 가중치는 다음과 같이 정의된다 .
w ξ = 1 n ( n + 1 ) ⋅ { 1 2 만일 ξ 꼭지점이다 1 만일 ξ 가장자리 포인트 2 만일 ξ 내부적인 지점이야. {\displaystyle w_{\mathbf {\xi } }={\frac {1}{n(n+1)}}\cdot {\begin{cases}{\frac {1}{2}}{\text{ if }}\mathbf {\xi } {\text{ is a vertex point}}\\1{\text{ if }}\mathbf {\xi } {\text{ is an edge point}}\\2{\text{ if }}\mathbf {\xi } {\text{ is an interior point. }}}\end{case}}
참조 ^ Caliari, Marco; Bos, Len; de Marchi, Stefano; Vianello, Marco; Xu, Yuan (2006), "Bivariate Lagrange interpolation at the Padua points: the generating curve approach", J. Approx. Theory , 143 (1): 15–25, arXiv :math/0604604 , doi :10.1016/j.jat.2006.03.008 ^ de Marchi, Stefano; Caliari, Marco; Vianello, Marco (2005), "Bivariate polynomial interpolation at new nodal sets", Appl. Math. Comput. , 165 (2): 261–274, doi :10.1016/j.amc.2004.07.001 ^ a b Caliari, Marco; de Marchi, Stefano; Vianello, Marco (2008), "Algorithm 886: Padua2D—Lagrange Interpolation at Padua Points on Bivariate Domains", ACM Transactions on Mathematical Software , 35 (3): 1–11, doi :10.1145/1391989.1391994 ^ Bos, Len; de Marchi, Stefano; Vianello, Marco; Xu, Yuan (2007), "Bivariate Lagrange interpolation at the Padua points: the ideal theory approach", Numerische Mathematik , 108 (1): 43–57, arXiv :math/0604604 , doi :10.1007/s00211-007-0112-z 외부 링크