11.08.2013 Views

Untitled - Cdm.unimo.it

Untitled - Cdm.unimo.it

Untitled - Cdm.unimo.it

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

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!