Чувствительность параметрического волнового уравнения 

In[2]:=

X Panel@Grid[ Join[{{"", "a", "c"}}, Table[{StringJoin["t = ", ToString[\[Tau]]], Plot[Evaluate[(pfun[a, c][\[Tau], x] + .5 {0, 1, -1} D[pfun[a, c][\[Tau], x], a]) /. {a -> 1, c -> 1}], {x, -10, 10}, Filling -> {2 -> {3}}, PlotRange -> All], Plot[Evaluate[(pfun[a, c][\[Tau], x] + .5 {0, 1, -1} D[pfun[a, c][\[Tau], x], c]) /. {a -> 1, c -> 1}], {x, -10, 10}, Filling -> {2 -> {3}}, PlotRange -> All]}, {\[Tau], {1, 3, 5}}]]]

In[3]:=

X ColorSensitivityPlot3D[f_, fsense_, opts___] := Block[{fs = fsense}, Plot3D[f, {t, 0, 5}, {x, -10, 10}, Mesh -> None, ColorFunction -> Function[{t, x}, ColorData["LakeColors"][fs]], PlotPoints -> 60, ColorFunctionScaling -> False, opts]]; {ColorSensitivityPlot3D[pfun[1, 1][t, x], Abs[ifda[t, x]], PlotLabel -> a, PlotRange -> All], ColorSensitivityPlot3D[pfun[1, 1][t, x], .25 Abs[ifdc[t, x]], PlotLabel -> c, PlotRange -> All]}

Найдём и построим график параметрической чувствительности для параметрического волнового уравнения.

In[1]:=

X pfun = ParametricNDSolveValue[{D[u[t, x], t, t] == c^2 D[u[t, x], x, x], u[0, x] == Exp[-(a x)^2], Derivative[1, 0][u][0, x] == 0, u[t, -10] == u[t, 10]}, u, {t, 0, 5}, {x, -10, 10}, {a, c}]; ifun = pfun[1, 1]; ifda = D[pfun[a, 1], a] /. {a -> 1}; ifdc = D[pfun[1, c], c] /. {c -> 1};

Графики чувствительности к изменениям в высоте начального волнового профиля a и в скорости волны c.

Out[2]=

Ещё один способ визуализации чувствительности заключается в изменении цвета поверхности решения.

Out[3]=