hongyashan123 发表于 2017-9-4 16:39 
你把数据的和模型发过来
data <- read.csv("shuju1.csv",header = T)
data <- data[-which(dataB2grWindSpeed=="****"),] ## 剔除缺失值
Speed <- as.numeric(as.character(data
B2grWindSpeed=="****"),] ## 剔除缺失值 Speed <- as.numeric(as.character(data
B2grWindSpeed)) ## 风速转换为数值型
Process <- as.numeric(as.character(dataB2grGenPowerForProcess)) ## 功率转换为数值型
Direction <- 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 ## 显示拟合参数