В качестве предисловия, я не пытаюсь делать что-то nls или подгонять данные с помощью функции, мы просто изучаем, как влияет модель. Ниже приведена функция, с которой я работаю: V0: скорость при t = 0 + (начальное значение = 100), r: время в секундах (0,1,5,10,15,20,30,40,50,60, 75,90,105,120), k: постоянная затухания (начальное значение = 0,2)
V=V0*exp(-k*r)
V=data.frame(V)
Цель состоит в том, чтобы увидеть влияние V0 и k на модель, но индивидуально. Я должен изменить V0 в диапазоне, оставляя k один, а затем изменить k, оставляя V0 в покое. После этого мне нужно сделать 2 участка из этих диапазонов, по одному для каждой переменной. У меня есть код для изменения k:
V0=100
r = c(0,1,5,10,15,20,30,40,50,60,75,90,105,120)
k=c(0,0.05,0.1,0.15,0.2,0.25,0.3,0.35,0.4,0.45,0.5,0.55,0.60,0.65)
V=NULL
V=data.frame(V)
K=cbind.data.frame("0"=V0*exp(-k[1]*r), "0.05"=V0*exp(-k[2]*r),"0.1"= V0*exp(-k[3]*r), "0.15" = V0*exp(-k[4]*r),".2"= V0*exp(-k[5]*r), ".25" = V0*exp(-k[6]*r),".3"= V0*exp(-k[7]*r),".35"= V0*exp(-k[8]*r),".4"= V0*exp(-k[9]*r),".45" =V0*exp(-k[10]*r),".5" = V0*exp(-k[11]*r),".55" =V0*exp(-k[12]*r), ".6"=V0*exp(-k[13]*r),".65"= V0*exp(-k[14]*r))
# x axis we want period
# y axis we want velocity
# want to change k in a range
# leave V0 constant
# graph velocity lines with each k as subplot
Сгенерированный df имеет скорости с min (r) / min (k) в верхнем левом углу и max (r) / max (k) скоростей справа внизу, с увеличением r вниз и увеличением k вправо. Таким образом, это означает, что каждый столбец имеет одинаковое k и разные r. Я застрял на том, как построить график для этого. Любые мысли будут с благодарностью.