Hallo,
ich habe mal wieder ein Problem:
Ich moechte numerisch die DGL fuer eine 1-dimensionale
Wellengleichung loesen.
Fuer t=0 ist das System in Ruhe:
y[x,0]==0,
Derivative[0,1][y][x, 0]==0,
Ein Ende sei fest:
y[1,t]==0,
Wenn das andere Ende mit
y[0,t] ==Sin[2 Pi t]
peridisch bewegt wird, gibt es keine Probleme:
In[1]:=
welle1=
NDSolve[
{D[y[x, t],t, t] == D[y[x, t], x, x],
y[x,0]==0,
y[1,t]==0,
Derivative[0,1][y][x, 0]==0,
y[0,t]==Sin[2 Pi t]
},
y, {x, 0, 1}, {t, 0, 5}
]
Out[1]=
{{y\[Rule]InterpolatingFunction[{{0,1.},{0.,5.}},"<>"]}}
Wenn ich aber eine periodische KRAFT
Derivative[0,2][y][0,t]==Sin[2 Pi t]
auf ein Ende wirken lasse, findet Mathematica keine Loesung
In[2]:=
welle2=
NDSolve[
{D[y[x, t],t, t] == D[y[x, t], x, x],
y[x,0]==0,
y[1,t]==0,
Derivative[0,1][y][x, 0]==0,
Derivative[0,2][y][0,t]==Sin[2 Pi t]
},
y, {x, 0,1}, {t, 0, 5}
]
sondern es gibt folgende Fehlermeldung:
NDSolve::bcnorm : Boundary conditions
{MethodOfLinesDump`dv[1, 1](1, 0)[t, 0] == Sin[2 Pi t]}
include non-normal derivatives
Eigentlich sollte das doch zu loesen sein. Sehe ich nur den Wald
vor lauter Baeumen nicht, oder ist das tatsaechlich so ein
schwieriges Problem?
Mit freundlichen Gruessen
Jochen Skupin