11.08.2013 Views

Untitled - Cdm.unimo.it

Untitled - Cdm.unimo.it

Untitled - Cdm.unimo.it

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

178 Polynomial Approximation of Differential Equations<br />

of references is provided in chatelin(1983), where hints for the treatment of problems<br />

similar to the one considered here are given in chapter 4B.<br />

We can say something more about the Legendre case. Combining the proof of theo-<br />

rem 2.2.1 w<strong>it</strong>h that of theorem 8.2.2, one checks that, for α = β = 0, the eigenfunctions<br />

pn,m, 1 ≤ m ≤ n − 1, as well as their derivatives, are orthogonal w<strong>it</strong>h respect to the<br />

inner product (·, ·)w,n (see section 3.8). This is true if we assume that the eigenvalues<br />

λn,m, 1 ≤ m ≤ n − 1, are distinct.<br />

Let Xm, 1 ≤ m ≤ n − 1, denote a subspace of P 0 n of dimension m. Then, by<br />

virtue of the orthogonal<strong>it</strong>y of the eigenfunctions, the following characterization, known<br />

as Min-Max principle, holds:<br />

(8.6.7) λn,m = min<br />

Xm⊂P 0 n<br />

<br />

max<br />

p∈Xm<br />

p≡0<br />

p ′ 2 w<br />

p 2 w,n<br />

<br />

, 1 ≤ m ≤ n − 1.<br />

Applications of formula (8.6.7), in the field of fin<strong>it</strong>e-element approximations, are pre-<br />

sented for instance in strang and fix(1973), chapter six.<br />

Other s<strong>it</strong>uations may be considered. The problem<br />

⎧<br />

⎨ −φ<br />

(8.6.8)<br />

⎩<br />

′′ + µφ = λ φ in I =] − 1,1[, µ > 0,<br />

φ ′ (−1) = φ ′ (1) = 0,<br />

has a countable set of eigenvalues λm := π2<br />

4 m2 + µ, m ∈ N. These are effectively<br />

approximated by the eigenvalues of problem (8.2.13). We only remark that two of the<br />

λn,m’s are proportional to the parameter γ and are out of interest.<br />

Different arguments apply to problem (8.2.12). Here, the λn,m’s do not approach<br />

the eigenvalues of problem (8.6.8). If we replace the second equation in (8.2.12) by<br />

p ′ n,m(η (n)<br />

0 ) = λn,mpn,m(η (n)<br />

0 ), we get a negative eigenvalue, say λn,0. Another eigen-<br />

value, say λn,1, is such that limn→+∞(λn,0 + λn,1) = −1. The remaining eigenvalues<br />

approximate those of problem (8.6.8), w<strong>it</strong>h the exception of λ0 = µ.<br />

These considerations demonstrate the importance of the treatment of boundary<br />

cond<strong>it</strong>ions in spectral methods. A non correct specification may give raise to unexpected<br />

results. Though this is true in all numerical methods, the s<strong>it</strong>uation seems more crucial

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!