2-4. Error Analysis for Iterative Methods

공부하자·2022년 9월 19일
0

수치해석

목록 보기
4/9
post-thumbnail

Error Analysis for Iterative Methods

Order of Convergence

Definition 2.7

Suppose {pn}n=0\left\{p_n\right\}_{n=0}^{\infty} is a sequence that converges to pp, with pnpp_n \neq p for all nn. If positive constants λ\lambda and α\alpha exist with

{pn}n=0\left\{p_n\right\}_{n=0}^{\infty}가 모든 nn에 대해 pnpp_n \neq ppp로 수렴되는 시퀀스라고 가정한다. 양의 상수 λ\lambdaα\alpha가 존재하는 경우

limnpn+1ppnpα=λ\lim _{n \rightarrow \infty} \frac{\left|p_{n+1}-p\right|}{\left|p_n-p\right|^\alpha}=\lambda

then {pn}n=0\left\{p_n\right\}_{n=0}^{\infty} converges to pp of order α\alpha, with asymptotic error constant λ\lambda.

그런 다음 {pn}n=0\left\{p_n\right\}_{n=0}^{\infty}는 점근 오류 상수 α\alphapp로 수렴한다.

An iterative technique of the form pn=g(pn1)p_n=g\left(p_{n-1}\right) is said to be of order α\alpha if the sequence {pn}n=0\left\{p_n\right\}_{n=0}^{\infty} converges to the solution p=g(p)p=g(p) of order α\alpha.

{pn}n=0\left\{p_n\right\}_{n=0}^{\infty} 시퀀스가 α\alpha 차수의 솔루션 p=g(p)p=g(p)로 수렴하는 경우 pn=g(pn1)p_n=g\left(p_{n-1}\right) 형태의 반복 기술은 α\alpha 차수라고 한다.

Order α\alpha

  • 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\alpha=1 (and λ<1\lambda<1 ), the sequence is linearly convergent.
    (ii) If α=2\alpha=2, the sequence is quadratically convergent.

일반적으로 수렴 순서가 높은 수열은 낮은 수열보다 빠르게 수렴한다.
점근 상수는 수렴 속도에 영향을 미치지만 차수의 범위까지는 영향을 미치지 않는다.
두 가지 order에 특별한 주의를 기울인다.
(i) α=1\alpha=1(및 α<1\alpha<1)이면 시퀀스는 선형 수렴한다.
(ii) α=2\alpha=2이면 시퀀스는 2차적으로 수렴한다.

Illustration

Suppose that {pn}n=0\left\{p_n\right\}_{n=0}^{\infty} is linearly convergent to 0 with

limnpn+1pn=0.5\lim _{n \rightarrow \infty} \frac{\left|p_{n+1}\right|}{\left|p_n\right|}=0.5

and that {pn}n=0\left\{p_n\right\}_{n=0}^{\infty} is quadratically convergent to 0 with the same asymptotic error constant,

limnp~n+1p~n2=0.5\lim _{n \rightarrow \infty} \frac{\left|\tilde{p}_{n+1}\right|}{\left|\tilde{p}_n\right|^2}=0.5

For simplicity we assume that for each nn we have

pn+1pn0.5 and p~n+1p~n20.5\frac{\left|p_{n+1}\right|}{\left|p_n\right|} \approx 0.5 \text { and } \frac{\left|\tilde{p}_{n+1}\right|}{\left|\tilde{p}_n\right|^2} \approx 0.5
  • For the linearly convergent scheme, this means that
    pn0=pn0.5pn1(0.5)2pn2(0.5)np0,\left|p_n-0\right|=\left|p_n\right| \approx 0.5\left|p_{n-1}\right| \approx(0.5)^2\left|p_{n-2}\right| \approx \cdots \approx(0.5)^n\left|p_0\right|,
  • whereas the quadratically convergent procedure has
    p~n0=p~n0.5p~n120.5[0.5p~n22]2=(0.5)3p~n24(0.5)3[0.5p~n32]4=(0.5)7p~n38(0.5)2n1p~02\begin{aligned} \left|\tilde{p}_n-0\right|=\left|\tilde{p}_n\right| & \approx 0.5\left|\tilde{p}_{n-1}\right|^2 \approx 0.5\left[0.5\left|\tilde{p}_{n-2}\right|^2\right]^2=(0.5)^3\left|\tilde{p}_{n-2}\right|^4 \\ & \approx(0.5)^3\left[0.5\left|\tilde{p}_{n-3}\right|^2\right]^4=(0.5)^7\left|\tilde{p}_{n-3}\right|^8 \\ & \approx \cdots \approx(0.5)^{2^n-1}\left|\tilde{p}_0\right|^2 \end{aligned}

  • The quadratically convergent sequence is within 103810^{-38} of 0 by the seventh term. At least 126 terms are needed to ensure this accuracy for the linearly convergent sequence

    이차 수렴 시퀀스는 7번째 항까지 0의 103810^{-38} 내에 있다. 선형 수렴 시퀀스에 대한 이 정확도를 보장하기 위해 최소 126개의 항이 필요하다.

Theorem 2.8

Let gC[a,b]g \in C[a, b] be such that g(x)[a,b]g(x) \in[a, b], for all x[a,b]x \in[a, b]. Suppose, in addition, that gg^{\prime} is continuous on (a,b)(a, b) and a positive constant k<1k<1 exists with

gC[a,b]g \in C[a, b]가 모든 x[a,b]x \in[a, b]에 대해 g(x)[a,b]g(x) \in[a, b]가 되도록 하자. 또한 gg^{\prime}(a,b)(a, b)에서 연속되고 양의 상수 k<1k<1가 존재한다고 가정하자.

g(x)k,\left|g^{\prime}(x)\right| \leq k, \quad for all x(a,b)x \in(a, b).
If g(p)0g^{\prime}(p) \neq 0, then for any number p0pp_0 \neq p in [a,b][a, b], the sequence

pn=g(pn1), for n1,p_n=g\left(p_{n-1}\right), \quad \text { for } n \geq 1,

converges only linearly to the unique fixed point pp in [a,b][a, b].

[a,b][a, b]에서 고유한 고정점 pp에만 선형적으로 수렴한다.

Theorem 2.9

Let pp be a solution of the equation x=g(x)x=g(x). Suppose that g(p)=0g^{\prime}(p)=0 and gg^{\prime \prime} is continuous with g(x)<M\left|g^{\prime \prime}(x)\right|<M on an open interval II containing pp. Then there exists a δ>0\delta>0 such that, for p0[pδ,p+δ]p_0 \in[p-\delta, p+\delta], the sequence defined by pn=g(pn1)p_n=g\left(p_{n-1}\right), when n1n \geq 1, converges at least quadratically to pp. Moreover, for sufficiently large values of nn,

ppx=g(x)x=g(x) 방정식의 해라고 하자. g(p)=0g^{\prime}(p)=0g(x)g^{\prime \prime}(x)pp를 포함하는 열린 간격 II에서 g(x)<M|g^{\prime \prime}(x)|<M와 연속된다고 가정한다.
그런 다음 n1n \geq 1일 때 p0[pδ,p+δ]p_0 \in[p-\delta, p+\delta]에 대해 pn=g(pn1)p_n=g\left(p_{n-1}\right)에 의해 정의된 시퀀스가 적어도 pp로 사분면적으로 수렴하는 δ>0\delta>0가 존재한다. 또한 nn이 충분히 큰 값의 경우,

pn+1p<M2pnp2\left|p_{n+1}-p\right|<\frac{M}{2}\left|p_n-p\right|^2
profile
아주대학교 수업 기록

0개의 댓글