問題は、別の問題を解決するための初期条件として1つのPDEの解を使用したいので発生します。 NDSolveでは、解はInterpolatingFunctionによって与えられるので、2番目のPDEでInterpolatingFunctionを使う必要があります。これは可能ですか?なぜ私のコンピュータはそれを永遠に実行するのですか? マイコード:Mathematica; NDSolve; InterpolatingFunctionを初期条件として使用できますか?
(*********************Parameters********************)
\[Gamma] = 10^(-5);
T = 500;
tm = -250;
\[Sigma] = 100;
L = 1;
L0 = -(L/2);
L1 = L0 + L;
c = 29.979;
\[Omega] = 1.32949/10^8;
\[Eta] = 1539.1;
\[Mu] = 6.27;
\[Beta] = 0.1334;
k = 40895.3;
(*********************1st PDE********************)
solS = NDSolve[{D[sS[z, t], t] == (-(\[Gamma] + I*\[Beta]*z))*sS[z, t]-I*\[Omega]*aS[z, t],
D[aS[z, t], z] + (1/c)*D[aS[z, t], t] == (-I)*\[Eta]*k*sS[z, t], sS[z, -T] == 0,
aS[z, -T] == E^(-((-T - tm)^2/(2*\[Sigma]^2))),
aS[L0, t] == E^(-((t - tm)^2/(2*\[Sigma]^2)))}, {sS, aS}, {z, L0, L1}, {t, -T, 0},
MaxSteps -> Infinity]
(*************take the result and plot**************)
iniR = solS[[1]][[1]][[2]]
Plot3D[Re[iniR[z, t]], {z, L0, L1}, {t, -T, 0}, PlotRange -> All]
(***********2nd PDE, with the 1st result as ini. con.***********)
solR = NDSolve[{D[sR[z, t], t] == (-(\[Gamma] - I*\[Beta]*z))*sR[z, t] - I*\[Omega]*aR[z, t],
D[aR[z, t], z] + (1/c)*D[aR[z, t], t] == (-I)*\[Eta]*k*sR[z, t],
sR[z, 0] == iniR[z, 0], aR[z, 0] == 0, aR[L0, t] == 0}, {sR, aR}, {z, L0, L1},
{t, 0, T}, MaxSteps -> Infinity]
原則として、あるPDEの結果を使用して別のPDEを解決することができます。私はMathematicaでどのように処理されているのか分かりませんが、あなたの意図がこれを行うことであれば、補間を避けるために、両方のPDEに対して同じグリッド(またはより一般的には同じ表現)を持つことができます。 –
アレクサンドル:ありがとう。補間を避けるために、両方のPDEに同じグリッド(またはより一般的には同じ表現)を持たせることについて、より具体的に記述できますか?あなたは私の方程式を変えるのですか?どのようにMathematicaでそれを実現するのですか? – user1019155
これは[this physicsforums thread](http://www.physicsforums.com/showthread.php?t=544396)に関連しています。 – Simon