定刻観測から差分を取る

  • こちらから
  • トリチェリの定理
    • 側面に穴の開いた容器の、液面変化の時間変化を観測して、液面の高さと流出量との関係について考える
    • 測定は変化する液面高の単位時間当たりの値
    • 単位時間当たりの値からは、差を取るだけで、時刻別の速度の近似値が出る(差分)

torricelli <- rbind(c(180, 147, 121, 101, 84, 68, 54, 42, 29, 17, 7),
c(181, 147, 121, 101, 83, 68, 53, 42, 28, 17, 6),
c(179, 147, 121, 102, 84, 67, 54, 41, 30, 18, 7),
c(179, 146, 120, 100, 83, 68, 53, 41, 28, 17, 6),
c(180, 146, 121, 101, 84, 68, 54, 41, 30, 18, 7),
c(179, 146, 120, 101, 83, 67, 54, 41, 29, 17, 6))

par(mfcol=c(1,2))
matplot(t(torricelli))

dd<-apply(torricelli,1,diff)
matplot(dd,type="l")
par(mfcol=c(1,1))