|
data <- read.csv("shuju1.csv",header = T)
data <- data[-which(data$B2grWindSpeed=="****"),] ## 剔除缺失值
Speed <- as.numeric(as.character(data$B2grWindSpeed)) ## 风速转换为数值型
Process <- as.numeric(as.character(data$B2grGenPowerForProcess)) ## 功率转换为数值型
Direction <- as.numeric(as.character(data$B2grWindDirection)) ## 风向转换为数值型
Speed = Speed[-which(is.na (Process))]
Direction = Direction[-which(is.na (Process))]
Process = Process[-which(is.na (Process))]
p_max = max(p)
start = list(a = 1, b=4)
Mode2 <- nls(p~p_max*(1-exp(-(v/b)^a)), start = start)
M2 <- summary(Mode2)
M2$coefficients ## 显示拟合参数
|