Gauss-Legendre quadrature와 Gauss-Lobatto quadrature 비교하기

두 방법은 모두 수치적분을 위한 Gaussian 방법의 한 종류이지만, 서로 다른 특성을 갖습니다.

• Gauss-Legendre quadrature:
  • 적분 구간 내에서 적절한 위치(대개 직교 다항식의 근)에서 노드를 선택하여, 가능한 한 높은 차수의 다항식에 대해 적분을 정확하게 계산할 수 있습니다.
  • 일반적으로 구간의 양 끝점은 노드에 포함되지 않습니다.
• Gauss-Lobatto quadrature:
  • 노드 선택 시 구간의 양 끝점을 반드시 포함합니다.
  • 이로 인해 선택할 수 있는 노드의 위치가 제한되므로, 동일한 노드 수를 사용했을 때 정확도가 Gauss quadrature보다 낮을 수 있지만, 경계 값을 반드시 고려해야 하는 문제에 유리합니다.

따라서, 두 quadrature는 기본 아이디어는 유사하지만 노드 선택 조건이 달라 서로 다른 개념으로 이해할 수 있습니다.

1. Gauss-Legendre Quadrature

기본 아이디어

• 목적: 구간 $[a, b]$에서의 적분
  $$
  I[f] = \int_a^b f(x),dx
  $$
  를 $n$개의 샘플 점과 가중치를 사용하여 근사합니다.
• 표현식:
  $$
  I[f] \approx \sum_{i=1}^{n} w_i, f(x_i),
  $$
  여기서 $x_i$는 노드(node), $w_i$는 가중치(weight)입니다.

노드 선정

• 원리:
  적분 구간과 관련된 직교 다항식 (예를 들어, 구간 $[-1, 1]$에서는 Legendre 다항식)의 근(zeros)을 노드로 선택합니다.
• 특징:
  • 노드들은 구간의 내부에 위치하며, 양 끝점은 일반적으로 포함하지 않습니다.
  • 예를 들어, Gauss–Legendre quadrature에서는 $n$개의 노드 $x_1, x_2, \dots, x_n$가 Legendre 다항식 $P_n(x)$의 근으로 결정됩니다.

가중치 결정

• 계산 방법:
  각 노드에 대응하는 가중치 $w_i$는 정해진 정밀도 조건(정확하게 적분 가능한 다항식의 최대 차수를 만족하도록)과 직교성 조건을 이용해 결정됩니다.
• 예시:
  $$
  w_i = \frac{2}{\left(1-x_i^2\right)[P_n'(x_i)]^2},
  $$
  여기서 $P_n'(x_i)$는 Legendre 다항식 $P_n(x)$의 도함수입니다.

정확도

• 정확하게 통합되는 다항식 차수:
  $n$개의 노드를 사용할 때, Gauss quadrature는 최대 $2n-1$차 다항식에 대해 정확합니다.
• 오차항:
  오차항은 함수 $f$의 $2n$차 미분과 관련되어 있으며, 적분 구간의 성질과 사용된 직교 다항식의 특성에 따라 결정됩니다.

장점 및 적용 분야

• 장점:
  동일한 노드 수로 최대의 다항식 차수를 정확하게 적분할 수 있어 매우 효율적입니다.
• 적용 분야:
  고차원 적분, 정밀 수치 해석, 스펙트럴 방법 등에서 널리 사용됩니다.

---

2. Gauss–Lobatto Quadrature

기본 아이디어

• 목적:
  구간의 양 끝점의 함수값을 반드시 포함하는 적분 근사법입니다. 경계 조건이 중요한 문제나, 경계에서의 값이 이미 알려진 경우에 유용합니다.
• 표현식:
  $$
  I[f] = \int_a^b f(x),dx \approx w_0, f(a) + \sum_{i=1}^{n-2} w_i, f(x_i) + w_{n-1}, f(b),
  $$
  여기서 $x_0 = a$와 $x_{n-1} = b$는 고정된 노드입니다.

노드 선정

• 원리:
  • 구간의 양 끝점 $a$와 $b$는 반드시 포함됩니다.
  • 나머지 $n-2$개의 내부 노드는 특정 조건(예를 들어, 잔차 항을 최소화하도록) 하에 선택되며, 보통 관련 직교 다항식의 도함수의 근으로 결정됩니다.
• 예시 (구간 $[-1, 1]$의 경우):
  • 양 끝점은 $x_0 = -1$과 $x_{n-1} = 1$로 고정됩니다.
  • 내부 노드 $x_1, x_2, \dots, x_{n-2}$는 Legendre 다항식 $P_{n-1}(x)$의 도함수 $P'_{n-1}(x)$의 근으로 결정됩니다.

가중치 결정

• 계산 방법:
  내부 노드와 양 끝점에 대한 가중치는, 주어진 노드들이 다항식의 최대 차수에 대해 정확히 적분되도록 결정됩니다.
• 특징:
  내부 노드에 대한 가중치는 Gauss quadrature의 경우와 유사하게 계산되지만, 양 끝점에 대한 가중치는 고정된 형태를 갖습니다.
  예를 들어, Gauss–Lobatto quadrature에서는 양 끝점의 가중치가 특정한 상수 형태(예: $\frac{2}{n(n-1)}$)를 가질 수 있습니다.

정확도

• 정확하게 통합되는 다항식 차수:
  $n$개의 노드를 사용할 때, Gauss–Lobatto quadrature는 최대 $2n-3$차 다항식에 대해 정확합니다.
• 비교:
  동일한 노드 수를 사용했을 때 Gauss quadrature가 $2n-1$차까지 정확한 것에 비해, Gauss–Lobatto quadrature는 한두 차수가 낮습니다.

장점 및 적용 분야

• 장점:
  경계에서의 함수값을 반드시 포함하므로, 경계 조건이 중요한 문제에서 유리합니다.
  함수가 경계에서 특이한 거동을 보이거나 경계 값이 이미 알려진 경우, 추가적인 정보를 활용할 수 있습니다.
• 적용 분야:
  유한 요소법, 스펙트럴 요소법, 경계 조건을 엄격하게 반영해야 하는 문제 등에서 활용됩니다.

---

요약 비교

항목	Gauss Quadrature	Gauss–Lobatto Quadrature
노드	구간 내부의 직교 다항식 근	구간 양 끝점 포함, 내부 노드는 $P'_{n-1}(x)$의 근
정확도	최대 $2n-1$차 다항식까지 정확	최대 $2n-3$차 다항식까지 정확
가중치	직교 다항식의 도함수를 이용해 계산	양 끝점은 고정, 내부는 최적 조건에 따라 결정
적용	고차 정확성이 필요한 경우	경계 조건 반영이 중요한 문제

---

이와 같이, 두 quadrature 기법은 기본 아이디어는 유사하지만 노드 선정 방식과 정확도 측면에서 차이가 있습니다. 문제의 성격(예를 들어 경계 조건의 중요성)에 따라 적절한 방법을 선택할 수 있습니다.