Slope2

測定した信号が対象ポイントx*で示す勾配を計算します。

この勾配の計算には有限差分法を使用します。次の図に示すように、計算対象とする信号の左右にある2つの信号の値を使用します。

1. 有限差分による曲線の勾配の評価


時間T=5秒における並進非線形スプリングの剛性を計算します。このスプリングの特性は、I MARKER=66、J MARKER=55、自由長L0=100mmです。

この計算を実行するコードを以下に示します。
>>> # Force at T=5.1
>>> y2 = ValueAtTime (
      label         = "Force @ T=5.1",
      measuredValue = "FM(66,55,55)",
      atTime        = 5.1,
      delta         = 1e-1,
    )

>>> # Displacement at T=5.1
>>> x2 = ValueAtTime (
      label         = "Displacement @ T=5.1",
      measuredValue = "DM(66,55,55)",
      atTime        = 5.1,
      delta         = 1e-1,
    )

>>> # Force at T=4.9
>>> y1 = ValueAtTime (
      label         = "Force @ T=4.9",
      measuredValue = "FM(66,55,55)",
      atTime        = 4.9,
      delta         = 1e-1,
    )

>>> # Displacement at T=4.9
>>> x1 = ValueAtTime (
      label         = "Displacement @ T=4.9",
      measuredValue = "DM(66,55,55)",
      atTime        = 4.9,
      delta         = 1e-1,
    )

>>> # The Slope at T=5.0
>>> k5 = Slope2 (
      label = "Slope@5", 
      points = (y2,x2,y1,x1),
    )

この例ではSlope2が次のように計算されます。

k =   y 2  -  y 1 x 2  -  x 1 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbdfwBIj xAHbstHrhAaerbuLwBLnhiov2DGi1BTfMBaeXafv3ySLgzGmvETj2B Sbqefm0B1jxALjhiov2Daebbfv3ySLgzGueE0jxyaibaieYdf9irVe eu0dXdh9vqqj=hEeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVk FHe9pgea0dXdar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaiaaciWadm aadaGaciaaeeGabuGaaOqaaabbaaaaaG+acXwDLbWdbiaadUgacaqG 1aGaaeiiaiabg2da9iaabccadaWcaaqaaiaadMhadaWgaaWcbaGaaG OmaaqabaGccaqGGaGaaeylaiaabccacaWG5bWaaSbaaSqaaiaaigda aeqaaaGcbaGaamiEamaaBaaaleaacaaIYaaabeaakiaabccacaqGTa GaaeiiaiaadIhadaWgaaWcbaGaaGymaaqabaaaaaaa@50BA@