WOLFRAM

Wolfram言語 ＆ システムドキュメントセンター

パラメータ感度を計算する

    
パラメータを微分方程式の中で摂動すると，解の族が得られる：
Wolfram Language code: pfun = ParametricNDSolveValue[{y''[t] + y[t] == a Sin[2t ] + Cos[2y[t]], y[0] == 1, y'[0] == 0}, y, {t, 0, 10}, {a}]; Show[Plot[Evaluate@Table[pfun[a][t], {a, -.6, .6, .1}], {t, 0, 5}, PlotStyle -> Lighter[Gray, 0.3], PlotRange -> All], Plot[pfun[0][t], {t, 0, 5}, PlotStyle -> {Green, Thick}], ImageSize -> Medium]
感度の解は，パラメータの小さな変化の影響をはっきりと示す：
Wolfram Language code: Plot[Evaluate[pfun[0][t] + {-.5, 0, .5}pfun'[0][t]], {t, 0, 5}, PlotStyle -> {{Dashed, Lighter[Blue, .4]}, {Thick, Green}, {Dashed, Lighter[Red, .4]}}, Filling -> {1 -> {2}, 2 -> {3}}, ImageSize -> Medium]