ggplot2 で接線場や方向場を描く - 廿TT
接線場 微分方程式の解の振る舞いをみるために接線場を描いてみる。 以下で与えられるロジスティック方程式を考える。 接線場は t-x 平面の各点に傾き の小さな線分を描いたものである。 logis <- function(t,x,a=1){ a*x*(1-x) } #点(t1,x1)を通り, 傾きが slope の直線 TheLine <- function(t1,x1,slope,t){ slope(t1,x1)*(t-t1)+x1 } tv <- xv <- seq(-1.5,1.5,by=0.1) slopefield <- function(tv,xv,f){ tx <-expand.grid(x=tv,y=xv) data.frame(t=tx[,1], x=tx[,2], t2=tx[,1]+0.05, x2=TheLine(tx[,1],tx[,2],f,tx[,1]+0.05
1
リリース、障害情報などのサービスのお知らせ
最新の人気エントリーの配信
j次のブックマーク
k前のブックマーク
lあとで読む
eコメント一覧を開く
oページを開く
ggplot2 まとめ: 初歩から程よいレベルまで