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.

160 Polynomial Approximation of Differential Equations<br />

Theorem 8.2.3 - Let α = β = 0. Then, for any n ≥ 2, the eigenvalues of problem<br />

(8.2.1) are real and strictly pos<strong>it</strong>ive.<br />

Proof - Given that w ≡ 1, <strong>it</strong> is sufficient to recall the equal<strong>it</strong>y<br />

(8.2.11)<br />

1<br />

p<br />

−1<br />

′ n,m(¯pn,mw) ′ dx =<br />

1<br />

|p<br />

−1<br />

′ n,m| 2 dx > 0, 1 ≤ m ≤ n − 1.<br />

For different values of the parameters α and β, problem (8.2.1) still furnishes real<br />

eigenvalues, though the proof of this statement is more involved than for the Legendre<br />

case. By studying the roots of the characteristic polynomial, gottlieb and lust-<br />

man(1983) give the proof in the Chebyshev case. To this end, they analyze the eigen-<br />

values of two auxiliary first-order problems. Other cases are still an open question. For<br />

−1 < α < 1 and −1 < β < 1 , the numerical experiments generate real eigenvalues,<br />

and some complex eigenvalues w<strong>it</strong>h relatively small imaginary part. In general, the<br />

analysis of the eigenvalue problem related to (7.4.12) does not lead to real eigenvalues,<br />

due to the presence of the first-order derivative (see previous section). However, when<br />

A and B are sufficiently close to zero, the second order operator dominates, and the<br />

eigenvalues are a small perturbation of those in (8.2.1). Thus, their real part is pos<strong>it</strong>ive.<br />

We examine the eigenvalue problems associated to (7.4.13) and (7.4.16). Namely,<br />

for 0 ≤ m ≤ n, we have<br />

⎧<br />

(8.2.12)<br />

(8.2.13)<br />

⎪⎨<br />

⎪⎩<br />

⎧<br />

⎪⎨<br />

⎪⎩<br />

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

i ) + µpn,m(η (n)<br />

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

i ), 1 ≤ i ≤ n − 1,<br />

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

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

0 ),<br />

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

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

n ),<br />

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

i ) + µpn,m(η (n)<br />

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

i ), 1 ≤ i ≤ n − 1,<br />

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

0 ) + µpn,m(η (n)<br />

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

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

0 ),<br />

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

n ) + µpn,m(η (n)<br />

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

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

n ),<br />

where µ > 0 and γ ∈ R w<strong>it</strong>h γ = 0.

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

Saved successfully!

Ooh no, something went wrong!