Error Analysis for Iterative Methods
Order of Convergence
Definition 2.7
Suppose {pn}n=0∞ is a sequence that converges to p, with pn=p for all n. If positive constants λ and α exist with
{pn}n=0∞가 모든 n에 대해 pn=p인 p로 수렴되는 시퀀스라고 가정한다. 양의 상수 λ와 α가 존재하는 경우
n→∞lim∣pn−p∣α∣pn+1−p∣=λ
then {pn}n=0∞ converges to p of order α, with asymptotic error constant λ.
그런 다음 {pn}n=0∞는 점근 오류 상수 α의 p로 수렴한다.
An iterative technique of the form pn=g(pn−1) is said to be of order α if the sequence {pn}n=0∞ converges to the solution p=g(p) of order α.
{pn}n=0∞ 시퀀스가 α 차수의 솔루션 p=g(p)로 수렴하는 경우 pn=g(pn−1) 형태의 반복 기술은 α 차수라고 한다.
Order α
- In general, a sequence with a high order of convergence converges more rapidly than a sequence with a lower order.
- The asymptotic constant affects the speed of convergence but not to the extent of the order. Two cases of order are given special attention.
- Two cases of order are given special attention.
(i) If α=1 (and λ<1 ), the sequence is linearly convergent.
(ii) If α=2, the sequence is quadratically convergent.
일반적으로 수렴 순서가 높은 수열은 낮은 수열보다 빠르게 수렴한다.
점근 상수는 수렴 속도에 영향을 미치지만 차수의 범위까지는 영향을 미치지 않는다.
두 가지 order에 특별한 주의를 기울인다.
(i) α=1(및 α<1)이면 시퀀스는 선형 수렴한다.
(ii) α=2이면 시퀀스는 2차적으로 수렴한다.
Illustration
Suppose that {pn}n=0∞ is linearly convergent to 0 with
n→∞lim∣pn∣∣pn+1∣=0.5
and that {pn}n=0∞ is quadratically convergent to 0 with the same asymptotic error constant,
n→∞lim∣p~n∣2∣p~n+1∣=0.5
For simplicity we assume that for each n we have
∣pn∣∣pn+1∣≈0.5 and ∣p~n∣2∣p~n+1∣≈0.5
- For the linearly convergent scheme, this means that
∣pn−0∣=∣pn∣≈0.5∣pn−1∣≈(0.5)2∣pn−2∣≈⋯≈(0.5)n∣p0∣,
- whereas the quadratically convergent procedure has
∣p~n−0∣=∣p~n∣≈0.5∣p~n−1∣2≈0.5[0.5∣p~n−2∣2]2=(0.5)3∣p~n−2∣4≈(0.5)3[0.5∣p~n−3∣2]4=(0.5)7∣p~n−3∣8≈⋯≈(0.5)2n−1∣p~0∣2
- The quadratically convergent sequence is within 10−38 of 0 by the seventh term. At least 126 terms are needed to ensure this accuracy for the linearly convergent sequence
이차 수렴 시퀀스는 7번째 항까지 0의 10−38 내에 있다. 선형 수렴 시퀀스에 대한 이 정확도를 보장하기 위해 최소 126개의 항이 필요하다.
Theorem 2.8
Let g∈C[a,b] be such that g(x)∈[a,b], for all x∈[a,b]. Suppose, in addition, that g′ is continuous on (a,b) and a positive constant k<1 exists with
g∈C[a,b]가 모든 x∈[a,b]에 대해 g(x)∈[a,b]가 되도록 하자. 또한 g′가 (a,b)에서 연속되고 양의 상수 k<1가 존재한다고 가정하자.
∣g′(x)∣≤k, for all x∈(a,b).
If g′(p)=0, then for any number p0=p in [a,b], the sequence
pn=g(pn−1), for n≥1,
converges only linearly to the unique fixed point p in [a,b].
[a,b]에서 고유한 고정점 p에만 선형적으로 수렴한다.
Theorem 2.9
Let p be a solution of the equation x=g(x). Suppose that g′(p)=0 and g′′ is continuous with ∣g′′(x)∣<M on an open interval I containing p. Then there exists a δ>0 such that, for p0∈[p−δ,p+δ], the sequence defined by pn=g(pn−1), when n≥1, converges at least quadratically to p. Moreover, for sufficiently large values of n,
p를 x=g(x) 방정식의 해라고 하자. g′(p)=0 및 g′′(x)가 p를 포함하는 열린 간격 I에서 ∣g′′(x)∣<M와 연속된다고 가정한다.
그런 다음 n≥1일 때 p0∈[p−δ,p+δ]에 대해 pn=g(pn−1)에 의해 정의된 시퀀스가 적어도 p로 사분면적으로 수렴하는 δ>0가 존재한다. 또한 n이 충분히 큰 값의 경우,
∣pn+1−p∣<2M∣pn−p∣2