スコットランドの高校生(16〜18歳)対象の数学問題
1匹のワニが20メートル上流の川の対岸にいる獲物に忍び寄ろうとしている。
![](https://blogimg.goo.ne.jp/user_image/09/3e/f4b3822b49743b6ca86cec1a2ee97cb5.jpg)
ワニの移動速度は地上と水中では異なる。
図に示されているxメートル上流にある対岸のP地点へ向けてワニが泳いだ場合,獲物に到達するまでの移動時間が最小となる。
この所要時間は,以下の等式によって求められる。(引用者注:時間の単位は特に定める必要がないので省略し,無名数とする)
T(x) = 5*sqrt(36+x^2)+4*(20-x)
(a) (i) ワニが地上を移動しなかった場合の所要時間を求めよ。
(ii) ワニの泳いだ距離が最小だった場合の所要時間を求めよ。
(b) これらの両極端の2つの選択肢の中間に,所要時間を最小化するxの値が1つある。そのxの値を求めよ。
この問題を R を用いて解く。
まず,T(x) の図を描いておこう(0 ≦ x ≦ 20)
> T = function(x) 5*sqrt(36+x^2)+4*(20-x)
> x = seq(0, 20, by=0.01)
> y = T(x)
> plot(x, y, type="l")
![](https://blogimg.goo.ne.jp/user_image/54/89/a9a7d824b558498403fcf2e9bc2ffb18.png)
(a) (i) 「ワニが地上を移動しなかった場合」というのは,川を斜めに渡ったということ。すなわち,x = 20 のときの T(x) が答え。
> T(0)
[1] 110
答え 110
(a) (ii) 「ワニの泳いだ距離が最小だった場合」というのは,川を垂直に渡ったということ。すなわち,x = 0 のときの T(x) が答え。
> T(20)
[1] 104.4031
答え 104.4031
(b) T(x) が最小になるときの x を求めるということは,T(x) の微分係数が 0 になるときの x を求めるということ。
> g = deriv(~5*sqrt(36+x^2)+4*(20-x), "x", function.arg="x")
> g # 微分係数の関数
function (x)
{
.expr2 <- 36 + x^2
.value <- 5 * sqrt(.expr2) + 4 * (20 - x)
.grad <- array(0, c(length(.value), 1L), list(NULL, c("x")))
.grad[, "x"] <- 5 * (0.5 * (2 * x * .expr2^-0.5)) - 4
attr(.value, "gradient") <- .grad
.value
}
微分係数のグラフを描いておこう。
> plot(x, attr(g(x), "gradient"))
> plot(x, attr(g(x), "gradient"), type="l")
> abline(h=0, v=8, col="red", lty=2) # 結果論ではあるが x = 8 のとき,微分係数が 0 になることを図示しておく
![](https://blogimg.goo.ne.jp/user_image/63/b7/30a59a849e6642b3a5a04b391046e54a.png)
関数の最小値を求めるのは optimize 関数が使える
> optimize(g, lower=0, upper=20, tol = .Machine$double.eps) # deriv 関数の戻り値 g を与える場合
$minimum
[1] 8
$objective
[1] 98
attr(,"gradient")
x
[1,] -5.039269e-08
または
> optimize(T, lower=0, upper=20, tol=.Machine$double.eps) # もとの関数を与える場合
$minimum
[1] 8
$objective
[1] 98
いずれにせよ $minimum = 8 が解である(最短時間は $objective = 98)。
なお,deriv 関数が返す導関数 5 * (0.5 * (2 * x * (36 + x^2)^-0.5)) - 4 が 0 になる値は,uniroot でも求められる。
> h = function(x) 5 * (0.5 * (2 * x * (36 + x^2)^-0.5)) - 4
> uniroot(h, lower=0, upper=20, tol=.Machine$double.eps)
$root
[1] 8
$f.root
[1] 0
$iter
[1] 8
$init.it
[1] NA
$estim.prec
[1] 2.993504e-05
$root = 8 が解である(そのとき,$f.root = 0 となる)
> h(8)
[1] 0