Partiële differentiaalvergelijkingen: scheiding der variabelen

Moderators: dirkwb, Xilvo

Reageer
Gebruikersavatar
Berichten: 524

Parti

Hallo allemaal,

Zij
\(a > 0\)
en
\(b > 0\)
. Beschouw voor
\(u = u(x,t)\)
de vergelijking:
\(au_{tt} + u_t = u_{xx}\)
, voor
\(0 < x < 1\)
.

Met de randvoorwaarden
\(u(0,t) = 0 = u_x(1,t) + bu(1,t)\)
heeft de differentiaalvergelijkingen oplossingen van de vorm
\(u(x,t) = T(t)X(x)\)
.

Leid een randwaardeprobleem af voor
\(X(x)\)
waarin een paramater
\(\lambda\)
voorkomt en laat zien dat er een rij
\(\lambda = \lambda_1 < \lambda_2 < ...\)
is waarvoor dit randwaardeprobleem een niet-triviale oplossing
\(X_n(x)\)
heeft én laat voor
\(m \neq n\)
zien dat:
\(\int_0^1 X_n(x)X_m(x) dx = 0\)
[/b]

-----

Het leek mij de bedoeling de differentiaalvergelijking te herschrijven door de scheiding van variabelen:
\(aT''(t)X(x) + T'(t)X(x) = T(t)X''(x)\)
.

Alles delen door
\(T(t)X(x)\)
:
\(a\frac{T''(t)}{T(t)} + \frac{T'(t)}{T(t)} = \frac{X''(x)}{X(x)}\)
De linkerkant is niet afhankelijk van
\(x\)
en de rechterkant niet van
\(t\)
, dus dat betekent dat beide kanten constant zijn. Dat levert op:
\(\frac{X''(x)}{X(x)} = \lambda\)
, dus
\(X''(x) = \lambda X(x)\)
Ik dacht dat
\(X(x) = A\cos(\sqrt{\lambda} x) + B\sin(\sqrt{\lambda} x)\)
, maar dan:
\(X'(x) = -A\sqrt{\lambda}\sin(\sqrt{\lambda} x) + B\cos(\sqrt{\lambda} x)\)
.
\(X''(x) = -A \lambda \cos(\sqrt{\lambda} x) - B\lambda \sin(\sqrt{\lambda} x)\)
Dus dan
\(X''(x) = -\lambda X(x)\)
en dat was niet de bedoeling.

Hoe vind ik
\(X(x)\)
en die rij
\(\lambda_1 < \lambda_2 < ...\)
?

Alvast bedankt!

- Fruitschaal.

Berichten: 7.068

Re: Parti

Fruitschaal schreef: wo 10 apr 2013, 10:34
\(\frac{X''(x)}{X(x)} = \lambda\)
, dus
\(X''(x) = \lambda X(x)\)
Ik dacht dat
\(X(x) = A\cos(\sqrt{\lambda} x) + B\sin(\sqrt{\lambda} x)\)
,
Dat is niet juist. Probeer eens een e-macht.

Gebruikersavatar
Berichten: 524

Re: Parti

\(X(x) = Ce^{\sqrt{\lambda}x}\)
is inderdaad een mogelijkheid, maar dan klopt een randvoorwaarde niet, want
\(X(0) = 0\)
en in dit geval
\(X(0) = C\)
, dus dan zou
\(X(x) = 0\)
moeten zijn, maar dat is de triviale oplossing.

Berichten: 4.246

Re: Parti

X(0) ka nooit nul zijn in het geval van die e-macht.
Quitters never win and winners never quit.

Berichten: 7.068

Re: Parti

\(X(x) = Ce^{\sqrt{\lambda}x}\)
is inderdaad een mogelijkheid,
Dat is niet de enige (e-macht) mogelijkheid. Als je de andere mogelijkheid ook meeneemt dan krijg je geen problemen met de randvoorwaarde.

Gebruikersavatar
Berichten: 524

Re: Parti

dirkwb schreef: wo 10 apr 2013, 18:59
X(0) ka nooit nul zijn in het geval van die e-macht.
Ja, ik bedoelde
\(X(0) \neq 0\)
.
EvilBro schreef: wo 10 apr 2013, 19:48
Dat is niet de enige (e-macht) mogelijkheid. Als je de andere mogelijkheid ook meeneemt dan krijg je geen problemen met de randvoorwaarde.
Welke dan?

Berichten: 7.068

Re: Parti

\(D e^{-\sqrt{\lambda} x}\)
en nu mag je jezelf voor de kop slaan... :)

Berichten: 7.068

Re: Parti

Misschien is het handig dat ik er toch iets meer woorden aan vuil maak. Je hebt:
\(\frac{X''(x)}{X(x)} = K^2\)
ofwel
\(X''(x) - K^2 X(x) = 0\)
Er zijn nu 3 mogelijkheden:
\(K^2 > 0\)
,
\(K^2 = 0\)
en
\(K^2 < 0\)
Alle drie moet je onderzoeken. In het laatste geval zou ik de substitutie doen
\(L^2 = -K^2\)
. Je krijgt dan:
\(X''(x) + L^2 X(x) = 0 \mbox{ met } L^2 > 0\)

Reageer