2e orde dv

Moderators: dirkwb, Xilvo

Berichten: 4.246

2e orde dv

Beschouw het Poisson probleem:
\( -D(x) \frac{d^2 T}{dx^2}=1\)
voor
\( x \in [0,1] \)


met D(x) = 10-3 in [0,0.5] en D(x) = 1 in (0.5,1] en homogene dirichlet randvoorwaarden.

Hoe los ik dit exact op?
Quitters never win and winners never quit.

Gebruikersavatar
Berichten: 24.578

Re: 2e orde dv

Je kan de dv al exact in de twee verschillende intervallen oplossen...

Ik zag je edit pas later, ik ben niet bekend met die term voor de randvoorwaarden.
"Malgré moi, l'infini me tourmente." (Alfred de Musset)

Berichten: 4.246

Re: 2e orde dv

Je kan de dv al exact in de twee verschillende intervallen oplossen...
Hoe zit het dan met x=0,5?

edit: homogene dirichlet beginvoorwaarden => x(0)=x(1)=0.
Quitters never win and winners never quit.

Gebruikersavatar
Berichten: 24.578

Re: 2e orde dv

Is er continuïteit van de oplossing vereist? Dan kan je in x=0.5 uitdrukken dat beide oplossingen in dat punt gelijk moeten zijn.
"Malgré moi, l'infini me tourmente." (Alfred de Musset)

Berichten: 4.246

Re: 2e orde dv

Is er continuïteit van de oplossing vereist?
Dat lijkt me wel.
Dan kan je in x=0.5 uitdrukken dat beide oplossingen in dat punt gelijk moeten zijn.
Dat is niet voldoende je hebt nog iets nodig.
Quitters never win and winners never quit.

Gebruikersavatar
Berichten: 24.578

Re: 2e orde dv

De oplossing ligt daar inderdaad niet uniek mee vast, maar ik zie ook niets anders gevraagd... Continuïteit van de afgeleide misschien (interpretatie: snelheid gelijk)?
"Malgré moi, l'infini me tourmente." (Alfred de Musset)

Berichten: 2.746

Re: 2e orde dv

Dat is niet voldoende je hebt nog iets nodig.
continuiteit van de eerste afgeleide?

Berichten: 4.246

Re: 2e orde dv

continuiteit van de eerste afgeleide?
Dat is het vreemde eraan: dat leidt naar een tegenspraak :D
Quitters never win and winners never quit.

Gebruikersavatar
Berichten: 24.578

Re: 2e orde dv

Ik vind geen tegenspraak, maar ik heb het ook maar snel berekend...
"Malgré moi, l'infini me tourmente." (Alfred de Musset)

Berichten: 2.746

Re: 2e orde dv

kan je hier eens zetten wat je tot nu toe hebt?

en die tegenspraak is misschien wel mogelijk, misschien omdat D(x) discontinu is?

Gebruikersavatar
Berichten: 24.578

Re: 2e orde dv

Ik vind (blauw voor x<1/2, groen voor x>1/2):

<!--graphstart--><script type="text/javascript">graph(0,1,0,100,300,300,600,600, '3001*x/8 - 500*x^2', '- x^2/2 - 995*x/8 + 999/8')</script><!--graphend-->
"Malgré moi, l'infini me tourmente." (Alfred de Musset)

Berichten: 2.746

Re: 2e orde dv

en is het analytisch ook continu, of enkel grafisch?

(want 10^-3 is een smerig klein getal!)

Berichten: 4.246

Re: 2e orde dv

Ik zag dat ik een fout maakte, ik had x(1) = 0 niet goed gebruikt, ik probeer het opnieuw.
Quitters never win and winners never quit.

Gebruikersavatar
Berichten: 24.578

Re: 2e orde dv

stoker schreef:en is het analytisch ook continu, of enkel grafisch?

(want 10^-3 is een smerig klein getal!)
Het is "analytisch berekend", dus het zou vreemd zijn moest het dan niet continu zijn :P

Wiskundig is "smerig klein" natuurlijk relatief, 10-100 is veel kleiner :D
"Malgré moi, l'infini me tourmente." (Alfred de Musset)

Berichten: 4.246

Re: 2e orde dv

TD's oplossing heb ik gevonden. Nu is er alleen nog een fout in de discretisatie. :D
Quitters never win and winners never quit.

Reageer