# 加载 ggplot2 库
library(ggplot2)
# 固定参数赋值
# 定义表达式 y 的函数
y_function <- function(t,n) {
Ka <- 0.1
F <- 0.8
Dose <- 3
ktr <- 0.5
Input_Rate<- Ka*F*Dose*((ktr * t)^n / factorial(n)) * exp(-ktr*t)
}
# 定义时间范围
t <- seq(0, 24, by = 0.5) # 时间范围从 0 到 5，步长为 0.01
# 定义 n 的取值范围
n_values <- c(0, 1, 2, 3, 4, 5,6,7,8) # n 的取值
# 创建数据框以便 ggplot 使用
data <- data.frame()
for (n in n_values) {
y <- sapply(t, y_function, n = n) # 计算不同 t 下的 y 值
temp_data <- data.frame(t = t, y = y, n = factor(n)) # 将 n 转换为因子类型
data <- rbind(data, temp_data) # 合并数据
}
# 使用 ggplot 绘制图像
ggplot(data, aes(x = t, y = y, color = n)) +
geom_line(size = 1) + # 绘制曲线
labs(title = "y vs t for different n", x = "time", y = "Input Rate", color = "n") + # 添加标题和标签
theme_minimal() + # 使用简洁主题
scale_color_brewer(palette = "Set1") # 设置颜色样式
